Changeset 5652 for anuga_work/production/geraldton/project.py
- Timestamp:
- Aug 13, 2008, 4:40:45 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/geraldton/project.py
r5381 r5652 10 10 #from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm 11 11 from anuga.utilities.system_tools import get_user_name, get_host_name 12 from anuga.shallow_water.data_manager import urs2sts,create_sts_boundary 13 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon 12 14 13 15 # file and system info … … 15 17 16 18 home = getenv('INUNDATIONHOME') + sep +'data'+sep #Sandpit's parent dir 19 muxhome = getenv('MUXHOME') 17 20 user = get_user_name() 18 21 host = get_host_name() 22 19 23 # INUNDATIONHOME is the inundation directory, not the data directory. 20 21 #needed when running using mpirun, mpirun doesn't inherit umask from .bashrc22 umask(002)23 24 24 25 #time stuff 25 26 time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir 27 gtime = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir 26 28 build_time = time+'_build' 27 29 run_time = time+'_run' 28 29 tide = 0.75 #??? must check!!! 30 print 'gtime: ', gtime 30 31 31 32 #Making assumptions about the location of scenario data … … 34 35 scenario = 'geraldton_tsunami_scenario' 35 36 36 #Maybe will try to make project a class to allow these parameters to be passed in. 37 38 tide = 0.75 #??? must check!!! 39 37 40 alpha = 0.1 38 41 friction=0.01 39 starttime=10000 40 midtime=21600 41 #finaltime=25000 42 finaltime=10000 43 export_cellsize=50 42 starttime=0 43 finaltime=80000 44 export_cellsize=25 44 45 setup='final' 45 source=' other'46 source='polyline' 46 47 47 48 … … 62 63 yieldstep=60 63 64 64 dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+ str(user)65 dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+ 'alpha' +str(alpha)+'_'+str(user) 65 66 66 67 … … 70 71 71 72 # AHO + DPI data 72 coast_name = 'coastline' 73 offshore_name = 'geraldton_bathy' 73 coast_name = 'XYcoastline_KVP' 74 offshore_name = 'Geraldton_bathymetry' 75 offshore_name1 = 'DPI_Data' 76 offshore_name2 = 'grid250' 77 offshore_name3 = 'Top_bathymetry' 78 74 79 75 80 #final topo name … … 82 87 topographies_dir = home+state+sep+scenario+sep+'anuga'+sep+'topographies'+sep 83 88 84 # input topo file location89 ##input topo file location 85 90 onshore_in_dir_name = topographies_in_dir + onshore_name 86 91 island_in_dir_name = topographies_in_dir + island_name … … 88 93 coast_in_dir_name = topographies_in_dir + coast_name 89 94 offshore_in_dir_name = topographies_in_dir + offshore_name 95 offshore_in_dir_name1 = topographies_in_dir + offshore_name1 96 offshore_in_dir_name2 = topographies_in_dir + offshore_name2 97 offshore_in_dir_name3 = topographies_in_dir + offshore_name3 90 98 91 99 onshore_dir_name = topographies_dir + onshore_name … … 93 101 coast_dir_name = topographies_dir + coast_name 94 102 offshore_dir_name = topographies_dir + offshore_name 103 offshore_dir_name1 = topographies_dir + offshore_name1 104 offshore_dir_name2 = topographies_dir + offshore_name2 105 offshore_dir_name3 = topographies_dir + offshore_name3 95 106 96 107 #final topo files … … 104 115 tide_dir = anuga_dir+'tide_data'+sep 105 116 106 ##if source =='dampier': 107 ## boundaries_name = 'broome_3854_17042007' #Dampier gun 108 ## boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'dampier'+sep+'1_10000'+sep 109 ## 110 ##if source=='onslow': 111 ## boundaries_name = 'broome_3859_16052007' #onslow_hedland_broome gun 112 ## boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'onslow_hedland_broome'+sep+'1_10000'+sep 113 ## 114 ##if source=='exmouth': 115 ## boundaries_name = 'broome_3103_18052007' #exmouth gun 116 ## boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep 117 ## 118 ##if source=='other': 119 ## boundaries_name = 'other' #exmouth gun 120 ## boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep 121 ## 122 ## 123 ###boundaries locations 124 ##boundaries_in_dir_name = boundaries_in_dir + boundaries_name 125 boundaries_in_dir_name = anuga_dir + sep + 'boundaries' + sep + scenario_name + '49' 117 #boundaries_source = '1' 118 119 #boundaries locations 126 120 boundaries_dir = anuga_dir+'boundaries'+sep 127 121 boundaries_dir_name = boundaries_dir + scenario_name 122 boundaries_dir_mux = muxhome 128 123 129 124 #output locations 130 125 output_dir = anuga_dir+'outputs'+sep 131 output_build_time_dir = output_dir+build_time+dir_comment+sep 132 output_run_time_dir = output_dir +run_time+dir_comment+sep 133 #output_run_time_dir = output_dir+sep+run_time+sep 126 output_build_time_dir = anuga_dir+'outputs'+sep+build_time+dir_comment+sep 127 output_run_time_dir = anuga_dir+'outputs'+sep+run_time+dir_comment+sep 134 128 output_run_time_dir_name = output_run_time_dir + scenario_name #Used by post processing 129 130 vertex_filename = output_run_time_dir + 'mesh_vertex.csv' 135 131 136 132 #gauges … … 147 143 148 144 ############################### 149 # Domain definitions145 # Interior region definitions 150 146 ############################### 151 147 152 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon 148 # bounding polygon for study area 149 poly_all = read_polygon(polygons_dir+'poly_all.csv') 150 res_poly_all = 100000*res_factor 153 151 152 poly_neg20_pos20 = read_polygon(polygons_dir+'neg20_pos20.csv') 153 res_neg20_pos20 = 20000*res_factor 154 155 poly_CBD_1km= read_polygon(polygons_dir+'CBD_1km.csv') 156 res_CBD_1km = 800*res_factor 157 158 poly_cbd = read_polygon(polygons_dir+'CBD_500m.csv') 159 res_cbd = 500*res_factor 160 161 poly_wallabi = read_polygon(polygons_dir+'island_wallabi_poly.csv') 162 res_wallabi = 1000*res_factor 163 164 poly_dingiville = read_polygon(polygons_dir+'island_dingiville_poly.csv') 165 res_dingiville = 1000*res_factor 166 167 poly_pelsaert = read_polygon(polygons_dir+'island_pelsaert_poly.csv') 168 res_pelsaert = 1000*res_factor 169 170 171 #plot_polygons([polyAll,poly_broome1,poly_broome2,poly_broome3],figname='boundingpoly2',verbose=False) 172 173 interior_regions = [[poly_neg20_pos20,res_neg20_pos20],[poly_CBD_1km,res_CBD_1km], 174 [poly_cbd,res_cbd],[poly_wallabi,res_wallabi], 175 [poly_dingiville,res_dingiville],[poly_pelsaert, res_pelsaert]] 176 177 trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all) 178 print 'min number triangles', trigs_min 179 180 poly_mainland=read_polygon(topographies_in_dir +'initial_condition.csv') 154 181 155 182 ################################################################### … … 157 184 ################################################################### 158 185 159 # exporting asc grid160 eastingmin = 340000161 eastingmax = 350000162 northingmin = 6273400163 northingmax = 6277700164 186 165 ############################### 166 # Interior region definitions 167 ############################### 187 # Geraldton CBD extract ascii grid 188 xminCBD = 262700 189 xmaxCBD = 267600 190 yminCBD = 6811500 191 ymaxCBD = 6816400 168 192 169 #digitized polygons170 # bounding polygon for study area171 poly_all = read_polygon(polygons_dir+'poly_all_z50.csv')172 res_poly_all = 100000*res_factor173 174 poly_neg20_pos20 = read_polygon(polygons_dir+'neg20_pos20_points.csv')175 res_neg20_pos20 = 20000*res_factor176 177 poly_neg10_pos10= read_polygon(polygons_dir+'neg10_pos10_points.csv')178 res_neg10_pos10 = 2000*res_factor179 180 poly_cbd = read_polygon(polygons_dir+'cbd_500m.csv')181 res_cbd = 500*res_factor182 183 poly_wallabi = read_polygon(polygons_dir+'island_wallabi_poly_points.csv')184 res_wallabi = 20000*res_factor185 186 poly_dingiville = read_polygon(polygons_dir+'island_dingiville_poly_points.csv')187 res_dingiville = 20000*res_factor188 189 poly_pelsaert = read_polygon(polygons_dir+'island_pelsaert_poly_points.csv')190 res_pelsaert = 20000*res_factor191 192 193 #plot_polygons([polyAll,poly_broome1,poly_broome2,poly_broome3],figname='boundingpoly2',verbose=False)194 195 interior_regions = [[poly_neg20_pos20,res_neg20_pos20],[poly_neg10_pos10,res_neg10_pos10],196 [poly_cbd,res_cbd],[poly_wallabi,res_wallabi],197 [poly_dingiville,res_dingiville],[poly_pelsaert, res_pelsaert]]198 199 boundary_tags={'back': [2, 3, 4, 5],200 'side': [1, 6], 'ocean': [0, 7, 8, 9, 10]}201 202 trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)203 204 poly_mainland=read_polygon(polygons_dir+'initial_condition.csv')205 206 print 'min number triangles', trigs_min207 208
Note: See TracChangeset
for help on using the changeset viewer.