Changeset 4533

Timestamp:

Jun 5, 2007, 1:39:30 PM (17 years ago)

Author:

nick

Message:

update to WA model scripts

Location:

anuga_work/production

Files:

• broome_2006/build_broome.py (modified) (7 diffs)
• broome_2006/project_urs.py (modified) (4 diffs)
• broome_2006/run_broome_urs.py (modified) (4 diffs)
• exmouth_2006/project.py (modified) (3 diffs)
• exmouth_2006/run_exmouth.py (modified) (2 diffs)
• onslow_2006/project_urs.py (modified) (4 diffs)
• onslow_2006/run_onslow_urs.py (modified) (2 diffs)
• pt_hedland_2006/build_pt_hedland.py (modified) (3 diffs)
• pt_hedland_2006/project_urs.py (modified) (4 diffs)

anuga_work/production/broome_2006/build_broome.py

@@ -2 +2 @@
 
 Source data such as elevation and boundary data is assumed to be available in
-directories specified by project.py
-The output sww file is stored in project.output_time_dir 
+directories specified by project_urs.py
+The output sww file is stored in project_urs.output_time_dir 
 
-The scenario is defined by a triangular mesh created from project.polygon,
+The scenario is defined by a triangular mesh created from project_urs.polygon,
 the elevation data and a simulated submarine landslide.
 
@@ -27 +27 @@
 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
 from anuga.geospatial_data.geospatial_data import *
-from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files
+from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files
+from anuga_parallel.parallel_abstraction import get_processor_name
 
 # Application specific imports
-import project   # Definition of file names and polygons
+import project_urs   # Definition of file names and polygons
 
 #------------------------------------------------------------------------------
@@ -37 +38 @@
 #------------------------------------------------------------------------------
 
-copy_code_files(project.output_build_time_dir,__file__, 
-               dirname(project.__file__)+sep+ project.__name__+'.py' )
+copy_code_files(project_urs.output_build_time_dir,__file__, 
+               dirname(project_urs.__file__)+sep+ project_urs.__name__+'.py' )
 
-start_screen_catcher(project.output_build_time_dir)
+start_screen_catcher(project_urs.output_build_time_dir)
 
-print 'USER:    ', project.user
+print 'USER:    ', project_urs.user
 
 #-------------------------------------------------------------------------------
@@ -51 +52 @@
 # Fine pts file to be clipped to area of interest
 #-------------------------------------------------------------------------------
-print"project.combined_dir_name",project.combined_dir_name
+print"project_urs.combined_dir_name",project_urs.combined_dir_name
 
 # topography directory filenames
-onshore_in_dir_name = project.onshore_in_dir_name
-coast_in_dir_name = project.coast_in_dir_name
-#island_in_dir_name = project.island_in_dir_name
-offshore_in_dir_name = project.offshore_in_dir_name
-offshore_in_dir_name1 = project.offshore_in_dir_name1
-offshore_in_dir_name2 = project.offshore_in_dir_name2
+onshore_in_dir_name = project_urs.onshore_in_dir_name
+coast_in_dir_name = project_urs.coast_in_dir_name
+#island_in_dir_name = project_urs.island_in_dir_name
+offshore_in_dir_name = project_urs.offshore_in_dir_name
+offshore_in_dir_name1 = project_urs.offshore_in_dir_name1
+offshore_in_dir_name2 = project_urs.offshore_in_dir_name2
 
-onshore_dir_name = project.onshore_dir_name
-coast_dir_name = project.coast_dir_name
-#island_dir_name = project.island_dir_name
-offshore_dir_name = project.offshore_dir_name
-offshore_dir_name1 = project.offshore_dir_name1
-offshore_dir_name2 = project.offshore_dir_name2
+onshore_dir_name = project_urs.onshore_dir_name
+coast_dir_name = project_urs.coast_dir_name
+#island_dir_name = project_urs.island_dir_name
+offshore_dir_name = project_urs.offshore_dir_name
+offshore_dir_name1 = project_urs.offshore_dir_name1
+offshore_dir_name2 = project_urs.offshore_dir_name2
 '''
 # creates DEM from asc data
@@ -101 +102 @@
 
 print'clip combined geospatial object by bounding polygon'
-G_clipped = G.clip(project.poly_all)
+G_clipped = G.clip(project_urs.poly_all)
 #FIXME: add a clip function to pts
 #print'shape of clipped data', G_clipped.get_data_points().shape
 
 print'export combined DEM file'
-if access(project.topographies_dir,F_OK) == 0:
-    mkdir (project.topographies_dir)
-print'export',project.combined_dir_name+'_unclipped' + '.txt'
-G.export_points_file(project.combined_dir_name+'_unclipped' + '.txt')
+if access(project_urs.topographies_dir,F_OK) == 0:
+    mkdir (project_urs.topographies_dir)
+print'export',project_urs.combined_dir_name+'_unclipped' + '.txt'
+G.export_points_file(project_urs.combined_dir_name+'_unclipped' + '.txt')
 
 print'split'
@@ -115 +116 @@
 
 
-print 'export',project.combined_dir_name + '.txt'
-G_clipped.export_points_file(project.combined_dir_name + '.txt')
-print 'export', project.combined_small_dir_name + '.txt'
-G_small.export_points_file(project.combined_small_dir_name + '.txt')
-#G_clipped.export_points_file(project.combined_dir_name + '.xya')
+print 'export',project_urs.combined_dir_name + '.txt'
+G_clipped.export_points_file(project_urs.combined_dir_name + '.txt')
+print 'export', project_urs.combined_small_dir_name + '.txt'
+G_small.export_points_file(project_urs.combined_small_dir_name + '.txt')
+#G_clipped.export_points_file(project_urs.combined_dir_name + '.xya')
 
 
 
-print'project.combined_dir_name + 1.xya',project.combined_dir_name + '1.xya'
-G_all=Geospatial_data(file_name = project.combined_dir_name + '1.xya')
+print'project_urs.combined_dir_name + 1.xya',project_urs.combined_dir_name + '1.xya'
+G_all=Geospatial_data(file_name = project_urs.combined_dir_name + '1.xya')
 print'split'
 G_all_1, G_all_2 = G_all.split(.10)
 print'export 1'
-G_all_1.export_points_file(project.combined_dir_name+'_small1' + '.xya')
+G_all_1.export_points_file(project_urs.combined_dir_name+'_small1' + '.xya')
 print'export 2'
-G_all_2.export_points_file(project.combined_dir_name+'_other1' + '.xya')
+G_all_2.export_points_file(project_urs.combined_dir_name+'_other1' + '.xya')
 
 '''
@@ -137 +138 @@
 #-------------------------------------------------------------------------
 print 'starting to create boundary conditions'
-boundaries_in_dir_name = project.boundaries_in_dir_name
-
 from anuga.shallow_water.data_manager import urs2sww, urs_ungridded2sww
 
-print 'boundaries_in_dir_name',boundaries_in_dir_name
-print 'project.boundaries_dir_name',project.boundaries_dir_name
+print 'boundaries_in_dir_name',project_urs.boundaries_in_dir_name
 
-urs_ungridded2sww(boundaries_in_dir_name, project.boundaries_dir_name,
+urs_ungridded2sww(project_urs.boundaries_in_dir_name, project_urs.boundaries_in_dir_name,
                   verbose=True, mint=4000, maxt=35000, zscale=1)
 

anuga_work/production/broome_2006/project_urs.py

@@ -18 +18 @@
 host = get_host_name()
 #needed when running using mpirun, mpirun doesn't inherit umask from .bashrc
-umask(002)
+#umask(002)
 
 
@@ -41 +41 @@
 finaltime=25000
 setup='final'
+extra='onslow'
 
 if setup =='trial':
@@ -58 +59 @@
     yieldstep=60
 
-dir_comment='_'+setup+'_'+str(tide)+'_'+str(user)
+dir_comment='_'+setup+'_'+str(tide)+'_'+str(extra)+'_'+str(user)
 
 # onshore data from 30m DTED level 2
@@ -105 +106 @@
 
 #boundaries locations
-boundaries_name = 'broome_3854_17042007'
+#boundaries_name = 'broome_3854_17042007' #dampier
+#boundaries_name = 'broome_3103_18052007' #exmouth
+boundaries_name = 'broome_3859_16052007' #onslow
 
-boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'1_10000'+sep
+#boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'1_10000'+sep
+#boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'dampier'+sep+'1_10000'+sep
+#boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
+boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'onslow_hedland_broome'+sep+'1_10000'+sep
 boundaries_in_dir_name = boundaries_in_dir + boundaries_name
 boundaries_dir = anuga_dir+'boundaries'+sep

anuga_work/production/broome_2006/run_broome_urs.py

@@ -30 +30 @@
 from anuga.shallow_water import Field_boundary
 from Numeric import allclose
+from anuga.shallow_water.data_manager import export_grid
 
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
@@ -36 +37 @@
 from anuga_parallel.parallel_abstraction import get_processor_name
 from anuga.caching import myhash
+from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage
+
 # Application specific imports
 import project_urs                 # Definition of file names and polygons
@@ -67 +70 @@
     kwargs['bathy_file']=project_urs.combined_dir_name_unclipped1+'.txt'
 #    kwargs['bathy_file']=project_urs.combined_small_dir_name + '.pts'
-    kwargs['boundary_file']=project_urs.boundaries_dir_name + '.sww'
+    kwargs['boundary_file']=project_urs.boundaries_in_dir_name + '.sww'
 
     print 'output_dir',kwargs['output_dir']
@@ -232 +235 @@
         store_parameters(**kwargs)
     barrier()
+    
+    swwfile = kwargs['output_dir']+kwargs['aa_scenario_name']
+    export_grid(swwfile, extra_name_out = 'town',
+                quantities = ['elevation','stage','speed'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation
+                timestep = None,
+                reduction = max,
+                cellsize = 25,
+                NODATA_value = -9999,
+                easting_min = project.eastingmin,
+                easting_max = project.eastingmax,
+                northing_min = project.northingmin,
+                northing_max = project.northingmax,
+                verbose = False,
+                origin = None,
+                datum = 'WGS84',
+                format = 'asc')
+    
+    buildings_filename = project.buildings_filename
+    buildings_filename_out = project.buildings_filename_out
+                
+    inundation_damage(swwfile+'.sww', buildings_filename, buildings_filename_out)
+    print '\n Augmented building file written to %s \n' %buildings_filename_out
 
 

anuga_work/production/exmouth_2006/project.py

@@ -17 +17 @@
 host = get_host_name()
 #needed when running using mpirun, mpirun doesn't inherit umask from .bashrc
-#umask(002)
+umask(002)
 
 #Making assumptions about the location of scenario data
@@ -39 +39 @@
 finaltime=25000
 setup='final'
-extra='exmouth'
+extra='onslow'
 
 if setup =='trial':
@@ -126 +126 @@
 output_dir = anuga_dir+'outputs'+sep
 output_build_time_dir = output_dir+build_time+sep
-#output_run_time_dir = output_dir +run_time+dir_comment+sep
-output_run_time_dir = output_dir +'20070528_061822_run_final_1.4_nbartzis'+sep
+output_run_time_dir = output_dir +run_time+dir_comment+sep
+#output_run_time_dir = output_dir +'20070530_063246_run_final_0_dampier_nbartzis'+sep
 output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
 

anuga_work/production/exmouth_2006/run_exmouth.py

@@ -66 +66 @@
     kwargs['boundary_file']=project.boundaries_in_dir_name + '.sww'
     
-    
-    '''
     print 'output_dir',kwargs['output_dir']
     if myid == 0:
@@ -231 +229 @@
 
 
-    '''
-    print'hello'
     swwfile = kwargs['output_dir']+kwargs['aa_scenario_name']
     export_grid(swwfile, extra_name_out = 'town',

anuga_work/production/onslow_2006/project_urs.py

@@ -18 +18 @@
 host = get_host_name()
 #needed when running using mpirun, mpirun doesn't inherit umask from .bashrc
-#umask(002)
+umask(002)
 
 
@@ -39 +39 @@
 finaltime=25000
 setup='final'
-extra='onslow'
+extra='exmouth'
 
 if setup =='trial':
@@ -121 +121 @@
 boundaries_source = ''
 #boundaries_name = 'onslow_3854_17042007' #Dampier gun
-boundaries_name = 'onslow_3859_16052007' #onslow_hedland_broome gun
-#boundaries_name = 'onslow_3103_18052007' #exmouth gun
+#boundaries_name = 'onslow_3859_16052007' #onslow_hedland_broome gun
+boundaries_name = 'onslow_3103_18052007' #exmouth gun
 
 #boundaries locations
 #boundaries_in_dir = anuga_dir+'boundaries'+sep
 #boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'1_10000'+sep
-boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'onslow_hedland_broome'+sep+'1_10000'+sep
-#boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
+#boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'dampier'+sep+'1_10000'+sep
+#boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'onslow_hedland_broome'+sep+'1_10000'+sep
+boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
 boundaries_in_dir_name = boundaries_in_dir + boundaries_name
 
@@ -137 +138 @@
 output_dir = anuga_dir+'outputs'+sep
 output_build_time_dir = output_dir+build_time+sep
-#output_run_time_dir = output_dir +run_time+dir_comment+sep
-output_run_time_dir = output_dir +'20070530_082920_run_final_1.5_nbartzis'+sep
+output_run_time_dir = output_dir +run_time+dir_comment+sep
+#output_run_time_dir = output_dir +'20070531_002753_run_final_0.0_onslow_nbartzis'+sep
 output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
 

anuga_work/production/onslow_2006/run_onslow_urs.py

@@ -66 +66 @@
 #    kwargs['bathy_file']=project_urs.combined_small_dir_name + '.txt'
     kwargs['boundary_file']=project_urs.boundaries_in_dir_name + '.sww'
-    '''
+    
     print 'output_dir',kwargs['output_dir']
     if myid == 0:
@@ -228 +228 @@
         store_parameters(**kwargs)
     barrier()
-    '''
+    
     swwfile = kwargs['output_dir']+kwargs['aa_scenario_name']
     print 'swwfile',swwfile

anuga_work/production/pt_hedland_2006/build_pt_hedland.py

@@ -25 +25 @@
 # Related major packages
 from anuga.shallow_water import Domain
-from anuga.shallow_water import Dirichlet_boundary
-from anuga.shallow_water import File_boundary
-from anuga.shallow_water import Reflective_boundary
-from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
-from anuga.pmesh.mesh_interface import create_mesh_from_regions
-from anuga.geospatial_data.geospatial_data import Geospatial_data
-from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files
+from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts,urs_ungridded2sww
+from anuga.geospatial_data.geospatial_data import *
+from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files
 from anuga_parallel.parallel_abstraction import get_processor_name
 
@@ -105 +101 @@
 #G_other.export_points_file(project_urs.combined_dir_name + '_other' + '.pts')
 
-'''
+
 print 'start reading:',project_urs.combined_dir_name + '.txt'
 G = Geospatial_data(file_name = (project_urs.combined_dir_name + '.txt'))
@@ -116 +112 @@
 #G.export_points_file(project_urs.combined_smaller_dir_name + '.txt')
 G_clipped.export_points_file(project_urs.combined_smallest_dir_name + '.txt')
-
+'''
 #-------------------------------------------------------------------------
 # Convert URS to SWW file for boundary conditions 
 #-------------------------------------------------------------------------
 print 'starting to create boundary conditions'
-#boundaries_in_dir_name = project_urs.boundaries_in_dir_name
+     
+print 'boundaries_in_dir_name',project_urs.boundaries_in_dir_name
 
-#print 'minlat=project_urs.north_boundary, maxlat=project_urs.south_boundary',project_urs.north_boundary, project_urs.south_boundary
-#print 'minlon= project_urs.west_boundary, maxlon=project_urs.east_boundary',project_urs.west_boundary, project_urs.east_boundary
-'''       
-from anuga.shallow_water.data_manager import urs_ungridded2sww
-
-print 'boundaries_in_dir_name',project_urs.boundaries_in_dir_name
-print 'project.boundaries_dir_name',project_urs.boundaries_dir_name
-
-urs_ungridded2sww(project_urs.boundaries_in_dir_name, project_urs.boundaries_dir_name,
+urs_ungridded2sww(project_urs.boundaries_in_dir_name, project_urs.boundaries_in_dir_name,
                   verbose=True, mint=5000, maxt=35000, zscale=1)
-'''
 
 
 
 
+

anuga_work/production/pt_hedland_2006/project_urs.py

@@ -36 +36 @@
 alpha = 0.1
 friction=0.01
+starttime=3600
 finaltime=25000
-starttime=3600
 setup='final'
-
-#dir_comment='_trial_'+str(tide)+'_'+str(user)
-#dir_comment='_final_'+str(tide)+'_'+str(user)
-#time_thinning=4
-#yieldstep=60
-#res_factor = 1
+extra='dampier'
 
 if setup =='trial':
@@ -62 +57 @@
     yieldstep=60
 
-
-dir_comment='_'+setup+'_'+str(tide)+'_'+str(user)
+dir_comment='_'+setup+'_'+str(tide)+'_'+str(extra)+'_'+str(user)
 
 # onshore data from 30m DTED level 2
@@ -108 +102 @@
 
 #boundaries locations
-boundaries_name = 'pt_hedland_3854_17042007'
-
-boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'1_10000'+sep
+boundaries_name = 'pt_hedland_3854_17042007' #dampier
+#boundaries_name = 'pt_hedland_3103_18052007' #exmouth
+#boundaries_name = 'pt_hedland_3859_16052007' #onslow
+
+#boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'1_10000'+sep
+boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'dampier'+sep+'1_10000'+sep
+#boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
+#boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'onslow_hedland_broome'+sep+'1_10000'+sep
 boundaries_in_dir_name = boundaries_in_dir + boundaries_name
 boundaries_dir = anuga_dir+'boundaries'+sep
@@ -117 +116 @@
 #output locations
 output_dir = anuga_dir+'outputs'+sep
-output_build_time_dir = output_dir+build_time+sep
+output_build_time_dir = output_dir+build_time+'_'+extra+sep
 output_run_time_dir = output_dir +run_time+dir_comment+sep
 output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing