Changeset 4587 for anuga_work/production/dampier_2006/run_dampier_urs.py

Timestamp:

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

Author:

nick

Message:

update model scripts

File:

• anuga_work/production/dampier_2006/run_dampier_urs.py (modified) (4 diffs)

anuga_work/production/dampier_2006/run_dampier_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
@@ -69 +70 @@
     kwargs['boundary_file']=project_urs.boundaries_in_dir_name + '.sww'
 #    kwargs['Completed']=''
-
+    
     start_screen_catcher(kwargs['output_dir'], myid, numprocs)
 
@@ -223 +224 @@
     
     kwargs['completed']=str(time.time()-t0)
-    
-    if myid == 0:
-        store_parameters(**kwargs)
-
-    swwfile = kwargs['output_dir']+kwargs['aa_scenario_name']
-    export_grid(swwfile, extra_name_out = 'town',
+     
+    if myid == 0:
+ #       store_parameters(**kwargs)
+        
+        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']
+        export_grid(swwfile, extra_name_out = 'town',
                 quantities = ['elevation','depth','stage','speed'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation
                 timestep = None,
@@ -243 +248 @@
                 format = 'asc')
     
-    buildings_filename = project_urs.buildings_filename
-    buildings_filename_out = project_urs.buildings_filename_out
+        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
+        inundation_damage(swwfile+'.sww', buildings_filename, buildings_filename_out)
+        print '\n Augmented building file written to %s \n' %buildings_filename_out
     
     barrier()