Changeset 6217

Timestamp:

Jan 20, 2009, 2:04:09 PM (15 years ago)

Author:

kristy

Message:

updated scripts, new poly_all (no land)

Location:

anuga_work/production/australia_ph2/dampier

Files:

• export_results_max.py (modified) (4 diffs)
• get_timeseries.py (modified) (3 diffs)
• project_250m.py (modified) (7 diffs)
• run_dampier_ph2_250m.py (modified) (3 diffs)

anuga_work/production/australia_ph2/dampier/export_results_max.py

@@ -20 +20 @@
 directory = project_250m.output_dir
 
-##time_dir1 = '20090115_151520_run_final_0_250m_kvanputt'
-time_dir2 = '20090115_164756_run_final_0_250m_small_kvanputt'
-time_dirs = [time_dir2]   
+time_dir1 = '20090115_151520_run_final_0_250m_kvanputt'
+##time_dir2 = '20090115_164756_run_final_0_250m_small_kvanputt'
+time_dirs = [time_dir1]   
 
 #cellsize = 5
@@ -34 +34 @@
 ######
 
-#area = ['Sorrento', 'Geordie' ]  # strings must match keys in var_equations below
 area = ['All']      # 'All' means no special areas - the whole thing
 
@@ -78 +77 @@
 
             name1 = directory+time_dir+sep+project_250m.scenario_name
-##            name2 = directory+time_dir+sep+project_250m.scenario_name+'_time_39900_0'
-##            name3 = directory+time_dir+sep+project_250m.scenario_name+'_time_79800_0'
+            name2 = directory+time_dir+sep+project_250m.scenario_name+'_time_16200_0'
+            name3 = directory+time_dir+sep+project_250m.scenario_name+'_time_32400_0'
+            name4 = directory+time_dir+sep+project_250m.scenario_name+'_time_48600_0'
+            name5 = directory+time_dir+sep+project_250m.scenario_name+'_time_64800_0'
 
-            names = [name1] #, name2, name3]
+            names = [name1, name2, name3, name4, name5]
       
-##            asc_name = []   
+            asc_name = []   
                
             for name in names:
@@ -104 +105 @@
                             format = 'asc')
 
-##                asc_name.append(outname + '.asc')
-##    
-##            maxasc_outname = directory+time_dir+sep+project_250m_250m.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
-##
-##            print 'max asc outname ', maxasc_outname
-##            print 'asc_name ', str(asc_name)
-##            
-##            MaxAsc(maxasc_outname, asc_name)
+                asc_name.append(outname + '.asc')
+    
+            maxasc_outname = directory+time_dir+sep+project_250m.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
+            print 'max asc outname ', maxasc_outname
+            print 'asc_name ', str(asc_name)
+            
+            MaxAsc(maxasc_outname, asc_name)

anuga_work/production/australia_ph2/dampier/get_timeseries.py

@@ -2 +2 @@
 Generate time series of nominated "gauges".
 Input: sww file from run_perth.py
-       gauges project.gauge_filename
+       gauges project_250m.gauge_filename
 Output: csv files stage,speed,depth,elevation over time
 Stored in the 'outputs_dir' folder for respective .sww file
@@ -14 +14 @@
 from os import sep, rename, listdir, system
 from anuga.abstract_2d_finite_volumes.util import sww2csv_gauges,csv2timeseries_graphs
-import project
+import project_250m
 
-directory = project.output_dir
+directory = project_250m.output_dir
 
-time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'+sep+'gauge'
-time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'+sep+'gauge'
-time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'+sep+'gauge'
-time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'+sep+'gauge'
-time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'+sep+'gauge'
-time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'+sep+'gauge'
+time_dir1 = '20090115_151520_run_final_0_250m_kvanputt'
+time_dir2 = '20090115_164756_run_final_0_250m_small_kvanputt'
 
-time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1, time_dir2]
 
 for time_dir in time_dirs:
-    name = directory+time_dir+sep+project.scenario_name
-    gauge = project.gauges_dir_name
+    name = directory+time_dir+sep+project_250m.scenario_name
+    gauge = project_250m.gauges_dir_name
     sww2csv_gauges(name+'.sww',gauge,
                    quantities = ['stage','speed','depth','elevation'],
@@ -35 +31 @@
 
     
-##    name1 = directory+time_dir+sep+project.scenario_name 
-####    name2 = directory+time_dir+sep+project.scenario_name+'_time_43380_0'
-####    name3 = directory+time_dir+sep+project.scenario_name+'_time_78360_0'
+##    name1 = directory+time_dir+sep+project_250m.scenario_name 
+####    name2 = directory+time_dir+sep+project_250m.scenario_name+'_time_43380_0'
+####    name3 = directory+time_dir+sep+project_250m.scenario_name+'_time_78360_0'
 ##
 ##    names = [name1]#, name2, name3]
 ##
 ##    for i, name in enumerate(names):
-##        gauge = project.gauges_dir_name
+##        gauge = project_250m.gauges_dir_name
 ####        out_name = 'gauge_%d_' % i
 ##        print name

anuga_work/production/australia_ph2/dampier/project_250m.py

@@ -40 +40 @@
 scenario_name = 'dampier_ph2'
 scenario = 'dampier'
+area = 'large_b10' #'small_b10'
 
 # Model specific parameters. One or all can be changed each time the
@@ -75 +76 @@
 # Important to distinguish each run - ensure str(user) is included!
 # Note, the user is free to include as many parameters as desired
-dir_comment='_'+setup+'_'+str(tide)+'_250m_small_'+str(user)
+dir_comment='_'+setup+'_'+str(tide)+'_250m_' + str(area)+'_'+str(user)
 
 #------------------------------------------------------------------------------
@@ -86 +87 @@
 
 # gauges - used in get_timeseries.py
-#gauge_name = 'perth.csv'
+gauge_name = 'gauges.csv'
 
 # BOUNDING POLYGON - used in build_boundary.py and run_perth.py respectively
@@ -101 +102 @@
 # Output filename for elevation
 # this is a combination of all the data (utilisied in build_boundary)
-combined_name = scenario_name + '_combined_elevation_250m_small'
+combined_name = scenario_name + '_combined_elevation_250m_' + area
 
 #------------------------------------------------------------------------------
@@ -129 +130 @@
 
 # where the mesh sits (this is created during the run_perth.py)
-meshes_dir_name = meshes_dir + scenario_name+ '_small.msh'
+meshes_dir_name = meshes_dir + scenario_name+ '_' + area + '.msh'
 
 # where the boundary ordering files sit (this is used within build_boundary.py)
@@ -146 +147 @@
 output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
 
-#w here the directory of the gauges sit
-##gauges_dir_name = gauges_dir + gauge_name       #used for get_timeseries.py
+# where the directory of the gauges sit
+gauges_dir_name = gauges_dir + gauge_name       #used for get_timeseries.py
 ##gauges_dir_name2 = gauges_dir + gauge_name2     #used for get_timeseries.py
 
@@ -154 +155 @@
 #------------------------------------------------------------------------------
 
-#Land, to set the initial stage/water to be offcoast only
-poly_mainland = read_polygon(polygons_dir+'initial_condition_small.csv')
-
 # Initial bounding polygon for data clipping 
-poly_all = read_polygon(polygons_dir+'poly_all_small.csv')
+poly_all = read_polygon(polygons_dir+'poly_all_' +area + '.csv')
 res_poly_all = 125000*res_factor
 

anuga_work/production/australia_ph2/dampier/run_dampier_ph2_250m.py

@@ -96 +96 @@
 ##    boundary_tags={'back': [N+1,N+2,N+3], 'side': [N,N+4], 'ocean': range(N)}
 
-##    land = range (0,11)
-##    sea = range (12,122)
-##    bounding_polygon = project_250m.poly_all
-##    boundary_tags={'back': land, 'side': [11,122], 'ocean': sea}
-   
-    sea = range (4,22)
-    bounding_polygon = project_250m.poly_all
-    boundary_tags={'back': [0,1,2,23], 'side': [3,22], 'ocean': sea}
+    if project_250m.area =='large_b10':
+        land = range (0,516)
+        sea = range (516,626)
+        bounding_polygon = project_250m.poly_all
+        boundary_tags={'back': land, 'ocean': sea}
+    elif project_250m.area =='small_b10':
+        land = range (0,276)
+        sea = range (276,294)
+        bounding_polygon = project_250m.poly_all
+        boundary_tags={'back': land, 'ocean': sea}
+    else:
+        print 'area not defined therefore bondary_tags not defined'
 
     #--------------------------------------------------------------------------
@@ -144 +148 @@
 
     # sets the initial stage in the offcoast region only
-    IC = Polygon_function( [(project_250m.poly_mainland, 0)], default = kwargs['tide'],
-                             geo_reference = domain.geo_reference)
-    domain.set_quantity('stage', IC)
+##    IC = Polygon_function( [(project_250m.poly_mainland, 0)], default = kwargs['tide'],
+##                             geo_reference = domain.geo_reference)
+    domain.set_quantity('stage', 0)
     #domain.set_quantity('stage',kwargs['tide'] )
     domain.set_quantity('friction', kwargs['friction']) 
@@ -200 +204 @@
     Rb = Time_boundary (domain, f=lambda t: [(sin(t*pi/600) * exp(-(t/3600)**2)), 0.0, 0.0])   
 
-    domain.set_boundary({'back': Br,
-                         'side': Bt,
+    domain.set_boundary({'back': Bt,
                          'ocean': Rb})