Changeset 5381
- Timestamp:
- May 30, 2008, 10:52:58 AM (16 years ago)
- Location:
- anuga_work/production
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/busselton/build_busselton.py
r5149 r5381 72 72 print "creates DEMs from ascii data" 73 73 convert_dem_from_ascii2netcdf(onshore_in_dir_name, basename_out=onshore_dir_name, use_cache=True, verbose=True) 74 #convert_dem_from_ascii2netcdf(island_in_dir_name, basename_out=island_dir_name, use_cache=True, verbose=True)75 74 76 75 #creates pts file for onshore DEM … … 111 110 #G_all=Geospatial_data(file_name = project.combined_dir_name + '.txt') 112 111 print'split' 113 G_all_1 ,= G.split(.10)112 G_all_1 = G.split(.10) 114 113 print'export 1' 115 114 G_all_1.export_points_file(project.combined_dir_name+'_small' + '.txt') -
anuga_work/production/busselton/project.py
r5160 r5381 72 72 # AHO + DPI data 73 73 coast_name = 'coastline_excluding_beach_survey' 74 coast_name1 = 'b each_survey_final'75 offshore_name = ' busselton'74 coast_name1 = 'bussleton_beach_data' 75 offshore_name = 'Busselton_bathymetry' 76 76 77 77 #final topo name … … 84 84 topographies_dir = home+state+sep+scenario+sep+'anuga'+sep+'topographies'+sep 85 85 86 #input topo file location 87 onshore_in_dir_name = topographies_in_dir + onshore_name 88 #island_in_dir_name = topographies_in_dir + island_name 86 # input topo file location 87 onshore_in_dir_name = topographies_in_dir + onshore_name #topo 89 88 90 coast_in_dir_name = topographies_in_dir + coast_name 91 coast_in_dir_name1 = topographies_in_dir + coast_name1 92 offshore_in_dir_name = topographies_in_dir + offshore_name 89 coast_in_dir_name = topographies_in_dir + coast_name #coastline 90 coast_in_dir_name1 = topographies_in_dir + coast_name1 #beach survey 91 offshore_in_dir_name = topographies_in_dir + offshore_name #bathymetry 93 92 93 #output to anuga from build file 94 94 onshore_dir_name = topographies_dir + onshore_name 95 #island_dir_name = topographies_dir + island_name 95 96 96 coast_dir_name = topographies_dir + coast_name 97 97 coast_dir_name1 = topographies_dir + coast_name1 … … 177 177 res_busselton1 = 10000*res_factor 178 178 179 poly_busselton2 = read_polygon(polygons_dir+'neg5_pos5_poly _.csv')179 poly_busselton2 = read_polygon(polygons_dir+'neg5_pos5_poly.csv') 180 180 res_busselton2 = 500*res_factor 181 181 … … 183 183 184 184 interior_regions = [[poly_busselton1,res_busselton1],[poly_busselton2,res_busselton2]] 185 print 'Hello',interior_regions 185 186 186 boundary_tags={'back': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1 5],187 'side': [10 , 14], 'ocean': [11, 12, 13]}187 boundary_tags={'back': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 14], 188 'side': [10], 'ocean': [11, 12, 13]} 188 189 189 190 trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all) 190 191 191 poly_mainland=read_polygon( polygons_dir+'initial_condition_south.csv')192 poly_mainland=read_polygon(topographies_in_dir+'initial_condition_south.csv') 192 193 193 194 print 'min number triangles', trigs_min -
anuga_work/production/busselton/run_busselton.py
r5157 r5381 87 87 88 88 print 'start create mesh from regions' 89 89 90 90 create_mesh_from_regions(project.poly_all, 91 91 boundary_tags=project.boundary_tags, … … 94 94 filename=project.meshes_dir_name+'.msh', 95 95 use_cache=False, 96 verbose= True)96 verbose=False) 97 97 barrier() 98 98 … … 175 175 Bo = Dirichlet_boundary([kwargs['tide']+10.0,0,0]) 176 176 177 if project.source != 'other': 177 if project.source != 'test': 178 print 'start reading boundary file' 179 178 180 Bf = Field_boundary(kwargs['boundary_file'], 179 181 domain, time_thinning=kwargs['time_thinning'], mean_stage=kwargs['tide'], … … 189 191 'ocean': Bd}) 190 192 191 #kwargs['input_start_time']=domain.starttime193 kwargs['input_start_time']=domain.starttime 192 194 193 195 print'finish set boundary' … … 198 200 t0 = time.time() 199 201 200 for t in domain.evolve(yieldstep = 240, finaltime = kwargs[' finaltime']):202 for t in domain.evolve(yieldstep = 240, finaltime = kwargs['starttime']): 201 203 domain.write_time() 202 204 domain.write_boundary_statistics(tags = 'ocean') … … 298 300 export_model(**kwargs) 299 301 barrier 302 -
anuga_work/production/carnarvon/project.py
r5005 r5381 31 31 #Making assumptions about the location of scenario data 32 32 state = 'western_australia' 33 scenario_name = ' busselton'34 scenario = ' busselton_tsunami_scenario_2006'33 scenario_name = 'carnarvon' 34 scenario = 'carnarvon_tsunami_scenario' 35 35 36 36 #Maybe will try to make project a class to allow these parameters to be passed in. … … 73 73 74 74 #final topo name 75 combined_name =' busselton_combined_elevation.pts'76 combined_smaller_name = ' busselton_combined_elevation_smaller'75 combined_name ='carnarvon_combined_elevation.pts' 76 combined_smaller_name = 'carnarvon_combined_elevation_smaller' 77 77 78 78 anuga_dir = home+state+sep+scenario+sep+'anuga'+sep … … 137 137 gauges_dir_name = gauges_dir + gauge_name 138 138 139 buildings_filename = gauges_dir + ' Busselton_res_Project.csv'140 buildings_filename_out = ' Busselton_res_Project_modified.csv'139 buildings_filename = gauges_dir + 'carnarvon_res_Project.csv' 140 buildings_filename_out = 'carnarvon_res_Project_modified.csv' 141 141 142 142 community_filename = gauges_dir +'' … … 169 169 170 170 #digitized polygons 171 poly_ busselton1 = read_polygon(polygons_dir+'neg20_pos10_polygon.csv')172 res_ busselton1 = 10000*res_factor173 poly_ busselton2 = read_polygon(polygons_dir+'neg5_pos5_poly_.csv')174 res_ busselton2 = 500*res_factor171 poly_carnarvon1 = read_polygon(polygons_dir+'neg20_pos10_polygon.csv') 172 res_carnarvon1 = 10000*res_factor 173 poly_carnarvon2 = read_polygon(polygons_dir+'neg5_pos5_poly_.csv') 174 res_carnarvon2 = 500*res_factor 175 175 176 176 #plot_polygons([polyAll,poly_broome1,poly_broome2,poly_broome3],figname='boundingpoly2',verbose=False) 177 177 178 interior_regions = [[poly_ busselton1,res_busselton1],[poly_busselton2,res_busselton2]]178 interior_regions = [[poly_carnarvon1,res_carnarvon1],[poly_carnarvon2,res_carnarvon2]] 179 179 180 180 boundary_tags={'back': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 14], -
anuga_work/production/geraldton/project.py
r5150 r5381 104 104 tide_dir = anuga_dir+'tide_data'+sep 105 105 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 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' 125 126 boundaries_dir = anuga_dir+'boundaries'+sep 126 127 boundaries_dir_name = boundaries_dir + scenario_name … … 168 169 #digitized polygons 169 170 # bounding polygon for study area 170 poly_all = read_polygon(polygons_dir+' extent_new_points.csv')171 poly_all = read_polygon(polygons_dir+'poly_all_z50.csv') 171 172 res_poly_all = 100000*res_factor 172 173 … … 177 178 res_neg10_pos10 = 2000*res_factor 178 179 179 poly_cbd = read_polygon(polygons_dir+'cbd_ points.csv')180 poly_cbd = read_polygon(polygons_dir+'cbd_500m.csv') 180 181 res_cbd = 500*res_factor 181 182 -
anuga_work/production/onslow_2006/export_results.py
r4887 r5381 19 19 #time_dir = '20070904_235118_run_final_1.5_jsexton_exmouth_original' # revised model original data 20 20 #time_dir = '20070905_052311_run_final_1.5_jsexton_exmouth_revised' # revised model original + survey data 21 time_dir = '20070906_061418_run_final_1.5_jsexton_exmouth_original' # refined again21 #time_dir = '20070906_061418_run_final_1.5_jsexton_exmouth_original' # refined again 22 22 #time_dir = '20070907_055936_run_final_1.5_jsexton_exmouth_revised' # refined again 23 cellsize = 50 23 time_dir = '20080526_043334_run_final_1.5_kvanputt_exmouth_revised' # New 24 cellsize = 25 24 25 #cellsize = 150 25 26 timestep = None … … 43 44 print 'output dir:', name 44 45 45 #var = [2,3,4] # depth and speed46 var = [2,3,4] # depth and speed 46 47 #var = [2] # depth 47 48 #var = [4] 48 var = [0]49 #var = [0] 49 50 50 51 for which_var in var: -
anuga_work/production/onslow_2006/run_onslow_grad.py
r4698 r5381 106 106 interior_regions=project_grad.interior_regions, 107 107 filename=project_grad.mesh_name+'.msh', 108 use_cache= False,108 use_cache=True, 109 109 verbose=True) 110 110 … … 117 117 #above don't work 118 118 domain = Domain(project_grad.mesh_name+'.msh', 119 use_cache= False, verbose=True)119 use_cache=True, verbose=True) 120 120 121 121 print domain.statistics() … … 150 150 use_cache = True, 151 151 verbose = True, 152 alpha = alpha)152 alpha = project_grad.alpha) 153 153 154 154 … … 173 173 print 'domain id', id(domain) 174 174 domain.beta_h = 0 175 #domain.tight_slope_limiters = 1175 domain.tight_slope_limiters = 0 176 176 177 177 … … 187 187 time_thinning=time_thinning, 188 188 mean_stage=tide, 189 use_cache= False,189 use_cache=True, 190 190 verbose=True) 191 191 -
anuga_work/production/perth/project.py
r5360 r5381 155 155 gauges_dir_name = gauges_dir + gauge_name 156 156 157 buildings_filename = gauges_dir + 'Perth_res _Project.csv'157 buildings_filename = gauges_dir + 'Perth_resA.csv' 158 158 buildings_filename_out = 'Perth_res_Project_modified.csv' 159 159
Note: See TracChangeset
for help on using the changeset viewer.