Changeset 4587 for anuga_work/production/pt_hedland_2006/run_pt_hedland_urs.py

Timestamp:

Jul 4, 2007, 5:22:18 PM (17 years ago)

Author:

nick

Message:

update model scripts

File:

• anuga_work/production/pt_hedland_2006/run_pt_hedland_urs.py (modified) (6 diffs)

anuga_work/production/pt_hedland_2006/run_pt_hedland_urs.py

@@ -38 +38 @@
 from anuga.caching import myhash
 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage
+from anuga.fit_interpolate.benchmark_least_squares import mem_usage 
 # Application specific imports
 import project_urs                 # Definition of file names and polygons
@@ -53 +54 @@
     
     kwargs['output_dir']=project_urs.output_run_time_dir
-    kwargs['bathy_file']=project_urs.combined_dir_name + '.txt'
+    kwargs['bathy_file']=project_urs.combined_dir_name + '.pts'
     kwargs['boundary_file']=project_urs.boundaries_in_dir_name + '.sww'
 
@@ -61 +62 @@
     #------------------------------------------------------------------------------
 
-    start_screen_catcher(kwargs['output_dir'], myid, numprocs)
-
     print 'output_dir',kwargs['output_dir']
     if myid == 0:
@@ -71 +70 @@
 
     barrier()
+    
+    start_screen_catcher(kwargs['output_dir'], myid, numprocs)
+    
 
     print "Processor Name:",get_processor_name()
@@ -124 +126 @@
         from polygon import Polygon_function
         #following sets the stage/water to be offcoast only
-        IC = Polygon_function( [(project_urs.poly_mainland, -1.0)], default = kwargs['tide'],
-                                 geo_reference = domain.geo_reference)
-        domain.set_quantity('stage', IC)
+        #IC = Polygon_function( [(project_urs.poly_mainland, -50.0)], default = kwargs['tide'],
+        #                         geo_reference = domain.geo_reference)
+        #domain.set_quantity('stage', IC)
+        domain.set_quantity('stage', kwargs['tide'])
         domain.set_quantity('friction', kwargs['friction']) 
         
-        print 'Start Set quantity'
-
+        print 'Start Set quantity', kwargs['bathy_file'],kwargs['alpha']
+        
         domain.set_quantity('elevation', 
                             filename = kwargs['bathy_file'],
@@ -217 +220 @@
     if myid == 0:
         store_parameters(**kwargs)
-
-    swwfile = kwargs['output_dir']+kwargs['aa_scenario_name']
-    export_grid(swwfile, extra_name_out = 'town',
-                quantities = ['elevation','depth','stage','speed'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation
-                timestep = None,
-                reduction = max,
-                cellsize = 25,
-                NODATA_value = -9999,
-                easting_min = project_urs.eastingmin,
-                easting_max = project_urs.eastingmax,
-                northing_min = project_urs.northingmin,
-                northing_max = project_urs.northingmax,
-                verbose = True,
-                origin = None,
-                datum = 'WGS84',
-                format = 'asc')
-    
-    buildings_filename = project_urs.buildings_filename
-    buildings_filename_out = project_urs.buildings_filename_out
-                
-    inundation_damage(swwfile+'.sww', buildings_filename, buildings_filename_out)
-    print '\n Augmented building file written to %s \n' %buildings_filename_out
+        
+        print 'memory usage before del domain',mem_usage()
+        del domain
+        print 'memory usage after del domain',mem_usage()
+    
+        swwfile = kwargs['output_dir']+kwargs['aa_scenario_name']
+        print'swwfile',swwfile
+
+        export_grid(swwfile, extra_name_out = 'town',
+                    quantities = ['speed','depth','stage','elevation'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation
+                    timestep = None,
+                    reduction = max,
+                    cellsize = 25,
+                    NODATA_value = -9999,
+                    easting_min = project_urs.eastingmin,
+                    easting_max = project_urs.eastingmax,
+                    northing_min = project_urs.northingmin,
+                    northing_max = project_urs.northingmax,
+                    verbose = True,
+                    origin = None,
+                    datum = 'WGS84',
+                    format = 'asc')
+        
+        buildings_filename = project_urs.buildings_filename
+        buildings_filename_out = project_urs.buildings_filename_out
+                    
+        inundation_damage(swwfile+'.sww', buildings_filename, buildings_filename_out)
+        print '\n Augmented building file written to %s \n' %buildings_filename_out
     
     barrier()