Changeset 4422

Timestamp:

Apr 24, 2007, 12:03:30 PM (18 years ago)

Author:

nick

Message:

good version of dampier and update the export_results.py

Location:

anuga_work/production/dampier_2006

Files:

• export_results.py (modified) (4 diffs)
• get_timeseries.py (modified) (2 diffs)
• project.py (modified) (1 diff)
• run_dampier.py (modified) (5 diffs)

anuga_work/production/dampier_2006/export_results.py

@@ -5 +5 @@
 from os import sep
 
-time_dir = '20061127_055714_run' #
+time_dir = '20070418_035124_run' #
+cellsize = 25
+timestep = None
 directory = project.output_dir
 #name = directory + time_dir + sep + project.scenario_name
-name = directory+sep+project.scenario_name
+name = directory+time_dir+sep+project.scenario_name
 
-is_parallel = True
+is_parallel = False
 if is_parallel == True: nodes = 4
 print 'output dir:', name
 
-which_var = 1
+which_var = 2
 if which_var == 0:  # Stage
     outname = name + '_stage'
@@ -35 +37 @@
     quantityname = 'elevation'  #Elevation
 
+
+
 if is_parallel == True:
+#    print 'is_parallel',is_parallel
     for i in range(0,nodes):
         namei = name + '_P%d_%d' %(i,nodes)
@@ -42 +47 @@
         sww2dem(namei, basename_out = outnamei,
                     quantity = quantityname,
-                    timestep = 1,
-                    cellsize = 100,      
-                    #easting_min = project.e_min_area,
-                    #easting_max = project.e_max_area,
-                    #northing_min = project.n_min_area,
-                    #northing_max = project.n_max_area,        
+                    timestep = timestep,
+                    cellsize = cellsize,      
+                    easting_min = project.e_min_area,
+                    easting_max = project.e_max_area,
+                    northing_min = project.n_min_area,
+                    northing_max = project.n_max_area,        
                     reduction = max, 
                     verbose = True,
@@ -55 +60 @@
     sww2dem(name, basename_out = outname,
                 quantity = quantityname,
-                cellsize = 20,      
+                timestep = timestep,
+                cellsize = cellsize,      
                 easting_min = project.e_min_area,
                 easting_max = project.e_max_area,

anuga_work/production/dampier_2006/get_timeseries.py

@@ -24 +24 @@
 # sww file created from URS boundary
 production_dirs = {#'boundaries': '10000 yr wave height'
-                   '20070413_033455_run': 'URS boundary normal'
+                   '20070418_035124_run': 'URS boundary normal'
                    }
 
     
-is_parallel = True
+is_parallel = False
 if is_parallel == True:
     nodes = 8
@@ -89 +89 @@
         
     texname, elev_output = sww2timeseries(swwfiles,
-                                          project.gauges_dir_name_simple,
+                                          project.gauges_dir_name,
+#                                          project.gauges_dir_name_simple,
                                           production_dirs,
                                           report = False,

anuga_work/production/dampier_2006/project.py

@@ -21 +21 @@
 
 #time stuff
-time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
+time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
+#time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
 gtime = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
 build_time = time+'_build'

anuga_work/production/dampier_2006/run_dampier.py

@@ -100 +100 @@
 #domain = cache(Domain, (meshes_dir_name), {'use_cache':True, 'verbose':True}, verbose=True)
 #above don't work
-domain = Domain(meshes_dir_name, use_cache=True, verbose=True)
+domain = Domain(meshes_dir_name, use_cache=False, verbose=True)
 print 'domain id', id(domain) 
 print 'myhash', myhash(domain)      
@@ -106 +106 @@
 print domain.statistics()
 
-boundaries_dir_name=project.boundaries_dir_name
-
-print 'starting to create boundary conditions'
-
-from anuga.shallow_water.data_manager import urs2sww
+#boundaries_dir_name=project.boundaries_dir_name
+
+#print 'starting to create boundary conditions'
+
+#from anuga.shallow_water.data_manager import urs2sww
 
 
@@ -122 +122 @@
     from polygon import *
     #following sets the stage/water to be offcoast only
-    IC = Polygon_function( [(project.poly_mainland, 0.)], default = tide)
-#    IC = Polygon_function( [(project.poly_mainland, 0.)], default = 200)
+#    IC = Polygon_function( [(project.poly_mainland, 0.)], default = tide)
+    IC = Polygon_function( [(project.poly_mainland, 0.)], default = 200)
     domain.set_quantity('stage', IC)
     domain.set_quantity('friction', 0.01) 
@@ -159 +159 @@
 domain.set_maximum_allowed_speed(0.1) # Allow a little runoff (0.1 is OK)
 print 'domain id', id(domain)
-#domain.beta_h = 0
+domain.beta_h = 0
 #domain.limit2007 = 1
 
@@ -170 +170 @@
 
 Bf = Field_boundary(project.boundaries_dir_namea + '.sww',
-                    domain, time_thinning=12, mean_stage=tide, 
+                    domain, time_thinning=1, mean_stage=tide, 
                     use_cache=True, verbose=True)