Changeset 5415 for anuga_work/production/busselton

Timestamp:

Jun 24, 2008, 10:24:21 AM (17 years ago)

Author:

kristy

Message:

 

Location:

anuga_work/production/busselton

Files:

• build_busselton.py (modified) (4 diffs)
• export_results.py (modified) (1 diff)
• project.py (modified) (10 diffs)
• run_busselton.py (modified) (4 diffs)

anuga_work/production/busselton/build_busselton.py

@@ -37 +37 @@
 #------------------------------------------------------------------------------
 
-start_screen_catcher(project.output_build_time_dir)
-
-print 'time stamp: ',project.time
-print 'USER: ', project.user
-
-
 copy_code_files(project.output_build_time_dir,__file__, 
                dirname(project.__file__)+sep+ project.__name__+'.py' )
 
+start_screen_catcher(project.output_build_time_dir)
 
+print 'USER: ', project.user
 
 #-------------------------------------------------------------------------------
@@ -55 +51 @@
 # Fine pts file to be clipped to area of interest
 #-------------------------------------------------------------------------------
+print"project.poly_all",project.poly_all
+print"project.combined_dir_name",project.combined_dir_name
 
 # topography directory filenames
@@ -60 +58 @@
 coast_in_dir_name = project.coast_in_dir_name
 coast_in_dir_name1 = project.coast_in_dir_name1
-#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_dir_name = project.onshore_dir_name
 coast_dir_name = project.coast_dir_name
 coast_dir_name1 = project.coast_dir_name1
-#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
 
 # creates DEM from asc data
@@ -91 +92 @@
 G3 = Geospatial_data(file_name = coast_in_dir_name1 + '.txt',verbose=True)
 G_off = Geospatial_data(file_name = offshore_in_dir_name + '.txt',verbose=True)
+G_off1 = Geospatial_data(file_name = offshore_in_dir_name1 + '.txt',verbose=True)
+G_off2 = Geospatial_data(file_name = offshore_in_dir_name2 + '.txt',verbose=True)
 
 print'add all geospatial objects'
-G = G1 + G2 + G3 + G_off 
+G = G1 + G2 + G3 + G_off + G_off1 + G_off2 
 
 print'clip combined geospatial object by bounding polygon'
-#G_clipped = G.clip(project.bounding_polygon)
-#FIXME: add a clip function to pts
-#print'shape of clipped data', G_clipped.get_data_points().shape
+G_clipped = G.clip(project.poly_all)
 
 print'export combined DEM file'
 if access(project.topographies_dir,F_OK) == 0:
     mkdir (project.topographies_dir)
-G.export_points_file(project.combined_dir_name + '.txt')
-#G_clipped.export_points_file(project.combined_dir_name + '.xya')
+G_clipped.export_points_file(project.combined_dir_name + '.txt')
 
 print'project.combined_dir_name + .txt',project.combined_dir_name + '.txt'

anuga_work/production/busselton/export_results.py

@@ -6 +6 @@
 
 time_dir = '20080530_170705_run_final_0.6_exmouth_kvanputt'
+time_dir = '20080619_115523_run_trial_0.6_exmouth_current_kvanputt'
+
 cellsize = 25
 #cellsize = 150
 #timestep = 0
 directory = project.output_dir
-name = directory+time_dir+sep+project.scenario_name
+name = directory+sep+'test'+sep+time_dir+sep+project.scenario_name
 
 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts

anuga_work/production/busselton/project.py

@@ -13 +13 @@
 # file and system info
 #---------------------------------
+#codename = 'project.py'
 
 home = getenv('INUNDATIONHOME') + sep +'data'+sep #Sandpit's parent dir   
 user = get_user_name()
 host = get_host_name()
+
 # INUNDATIONHOME is the inundation directory, not the data directory.
-
-#needed when running using mpirun, mpirun doesn't inherit umask from .bashrc
-umask(002)
 
 #time stuff
 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'
 run_time = time+'_run'
-
-tide = 0.6
+print 'gtime: ', gtime
 
 #Making assumptions about the location of scenario data
@@ -34 +33 @@
 scenario = 'busselton_tsunami_scenario'
 
-#Maybe will try to make project a class to allow these parameters to be passed in.
+tide = 0.6
+
 alpha = 0.1
 friction=0.01
 starttime=0
-midtime=21600
-#finaltime=25000
-finaltime=86000
+finaltime=80000
 export_cellsize=25
 setup='final'
@@ -62 +60 @@
     yieldstep=60
 
-dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+str(user)
-
-
-# onshore data provided by WA DLI
-#onshore_name = 'DLI_orthophoto_DEM' # original
-onshore_name = 'dli_dem_clipped' # original
-#islands
+rev_num = 'current'
+##rev_num = '4695' # 2nd Sept 2007
+##rev_num = '4743' # 3nd Oct 2007
+##rev_num = '4777' # 1st Nov 2007
+##rev_num = '4874' # 3rd Dec 2007
+##rev_num = '4901' # 3rd Jan 2007
+##rev_num = '4990' # 5th Feb 2007
+##rev_num = '5103' # 3rd March 2007
+##rev_num = '5185' # 1st April 2007
+##rev_num = '5273' # 2nd May 2007
+
+
+dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+str(rev_num)+'_'+str(user)
+
+
+# onshore data provided by WA DLI - provided by Hamish on the 17th June 2008
+
+onshore_name = 'busselton_v2_gda94_mga50' # original
 
 # AHO + DPI data
-coast_name = 'coastline_excluding_beach_survey'
-coast_name1 = 'bussleton_beach_data'
-offshore_name = 'Busselton_bathymetry'
+coast_name = 'Contour_0' # provided by hamish, represent better coastline than the 100km as compared to charts
+coast_name1 = 'beach_survey_final'
+offshore_name = 'Bathymetry_final'
+offshore_name1 = 'Busselton_Chart'
+offshore_name2 = 'Busselton_Digitised'
 
 #final topo name
@@ -82 +93 @@
 
 topographies_in_dir = home+state+sep+scenario+sep+'elevation_final'+sep+'points'+sep
-topographies_dir = home+state+sep+scenario+sep+'anuga'+sep+'topographies'+sep
+topographies_dir = anuga_dir+'topographies'+sep
 
 # input topo file location
@@ -89 +100 @@
 coast_in_dir_name = topographies_in_dir + coast_name #coastline
 coast_in_dir_name1 = topographies_in_dir + coast_name1 #beach survey
+
 offshore_in_dir_name = topographies_in_dir + offshore_name #bathymetry
+offshore_in_dir_name1 = topographies_in_dir + offshore_name1 #bathymetry Charts
+offshore_in_dir_name2 = topographies_in_dir + offshore_name2 #Digitised Fairsheet
 
 #output to anuga from build file
@@ -96 +110 @@
 coast_dir_name = topographies_dir + coast_name
 coast_dir_name1 = topographies_dir + coast_name1
+
 offshore_dir_name = topographies_dir + offshore_name
+offshore_dir_name1 = topographies_dir + offshore_name1
+offshore_dir_name2 = topographies_dir + offshore_name2
 
 #final topo files
@@ -108 +125 @@
 tide_dir = anuga_dir+'tide_data'+sep
 
+#boundaries_source = '1'
+
 if source=='exmouth':
     boundaries_name = 'busselton_3103_30052008' #exmouth gun
@@ -114 +133 @@
 if source=='test':
     boundaries_name = 'other' #exmouth gun
-    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
+    boundaries_in_dir = anuga_dir+'boundaries'+sep
 
 
@@ -124 +143 @@
 #output locations
 output_dir = anuga_dir+'outputs'+sep
-output_build_time_dir = output_dir+build_time+dir_comment+sep
+output_build_time_dir = output_dir +'test'+sep + build_time + dir_comment + sep
 #output_run_time_dir = output_dir +run_time+dir_comment+sep
-output_run_time_dir = anuga_dir+sep+'outputs'+sep+run_time+dir_comment+sep
+output_run_time_dir = output_dir +'test'+sep + run_time + dir_comment +sep # take test out
 output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
 
 #gauges
-gauge_name = '???.csv'
+gauge_name = 'busselton.csv'
 gauges_dir = home+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
 gauges_dir_name = gauges_dir + gauge_name
@@ -151 +170 @@
 res_poly_all = 100000*res_factor
 
+
+###############################
+# Interior region definitions
+###############################
+
+#digitized polygons
+poly_large = read_polygon(polygons_dir+'coast_5km_d20m.csv')
+res_large = 1000*res_factor
+
+poly_busselton = read_polygon(polygons_dir+'busselton_1km.csv')
+res_busselton = 500*res_factor
+
+poly_bunbury = read_polygon(polygons_dir+'bunbury_1km.csv')
+res_bunbury = 500*res_factor
+
+
+interior_regions = [[poly_large,res_large],[poly_busselton,res_busselton],[poly_bunbury,res_bunbury]]
+
+
+boundary_tags={'back': [0, 1, 2, 3, 4, 9],
+               'side': [5], 'ocean': [6, 7, 8]}
+
+trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
+
+poly_mainland=read_polygon(polygons_dir+'initial_condition.csv')
+
+print 'min number triangles', trigs_min
+
+
 ###################################################################
 # Clipping regions for export to asc and regions for clipping data
 ###################################################################
 
-# exporting asc grid
-eastingmin = 340000
-eastingmax = 350000
-northingmin = 6273400
-northingmax = 6277700
-
-###############################
-# Interior region definitions
-###############################
-
-#digitized polygons
-poly_busselton1 = read_polygon(polygons_dir+'neg20_pos10_polygon.csv')
-res_busselton1 = 10000*res_factor
-
-poly_busselton2 = read_polygon(polygons_dir+'neg5_pos5_poly.csv')
-res_busselton2 = 500*res_factor
-
-#plot_polygons([polyAll,poly_broome1,poly_broome2,poly_broome3],figname='boundingpoly2',verbose=False)
-
-interior_regions = [[poly_busselton1,res_busselton1],[poly_busselton2,res_busselton2]]
-print 'Hello',interior_regions
-
-boundary_tags={'back': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 14],
-               'side': [10], 'ocean': [11, 12, 13]}
-
-trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
-
-poly_mainland=read_polygon(polygons_dir+'initial_condition.csv')
-
-print 'min number triangles', trigs_min
-
-
+# exporting asc grid for Busselton
+xminBusselton = 325000
+xmaxBusselton = 360000
+yminBusselton = 6273000
+ymaxBusselton = 6290000
+
+# exporting asc grid for Busselton
+xminBunbury = 370000
+xmaxBunbury = 386000
+yminBunbury = 6305000
+ymaxBunbury = 6334000

anuga_work/production/busselton/run_busselton.py

@@ -173 +173 @@
     Br = Reflective_boundary(domain)
     Bd = Dirichlet_boundary([kwargs['tide'],0,0])
-   
+    Bo = Dirichlet_boundary([kwargs['tide']+10.0,0,0]) # To be deleted for FESA runs
 
     Bf = Field_boundary(kwargs['boundary_file'],
@@ -181 +181 @@
     domain.set_boundary({'back': Bd,
                              'side': Bd,
-                             'ocean': Bf})
+                             'ocean': Bd}) #changed from Bf to Bd for large wave
 
     kwargs['input_start_time']=domain.starttime
@@ -192 +192 @@
     t0 = time.time()
 
-    for t in domain.evolve(yieldstep = 240, finaltime = kwargs['finaltime']): 
+    for t in domain.evolve(yieldstep = project.yieldstep, finaltime = kwargs['finaltime']
+                            ,skip_initial_step = False ): 
         domain.write_time()
-        domain.write_boundary_statistics(tags = 'ocean')     
+        domain.write_boundary_statistics(tags = 'ocean')
+
+        # To be deleted for FESA runs
+        if allclose(t, 240):
+            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bo})
+
+        if allclose(t, 1440):
+            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bd})
+
+
 
     x, y = domain.get_maximum_inundation_location()
@@ -222 +232 @@
     kwargs['starttime']=project.starttime
     kwargs['yieldstep']=project.yieldstep
-    kwargs['midtime']=project.midtime
     kwargs['finaltime']=project.finaltime