Changeset 2519
- Timestamp:
- Mar 9, 2006, 5:32:48 PM (19 years ago)
- Location:
- production/onslow_2006
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
production/onslow_2006/project.py
r2470 r2519 9 9 10 10 from pmesh.create_mesh import convert_points_from_latlon_to_utm 11 12 from coordinate_transforms.redfearn import degminsec2decimal_degrees 11 13 12 14 #Making assumptions about the location of scenario data … … 16 18 # 250m data to be provided 17 19 coarsename = 'onsl_bathydem250' # get from Neil/Ingo (DEM or topo data) 18 """ 20 19 21 # 30m data to be provided 20 finename = 'onshore_30' # get from Neil/Ingo (DEM or topo data) 21 """ 22 onshore_name = 'onslow_onshore_30m_dted' # get from Neil/Ingo (DEM or topo data) 23 24 25 boundary_basename = 'SU-AU' 22 26 23 27 # clipping region for fine elevation data 24 west = 240000 25 east = 340000 26 south = 7570000 27 north = 7645000 28 eastingmin = 250000 29 eastingmax = 330000 30 northingmin = 7580000 31 northingmax = 7635000 32 33 south = degminsec2decimal_degrees(-22,00,0) 34 north = degminsec2decimal_degrees(-21,10,0) 35 west = degminsec2decimal_degrees(114,30,0) 36 east = degminsec2decimal_degrees(115,30,0) 37 28 38 29 39 #swollen/ all data output … … 45 55 46 56 coarsedemname = datadir + coarsename 47 """ 48 finedemname = datadir + finename 57 58 onshore_dem_name = datadir + onshore_name 59 49 60 combineddemname = datadir + 'onslow_combined_elevation' 50 """51 61 52 62 outputname = outputdir + basename #Used by post processing -
production/onslow_2006/run_onslow.py
r2470 r2519 41 41 # filenames 42 42 coarsedemname = project.coarsedemname 43 ''' 44 #finedemname = project.finedemname45 ''' 43 44 onshore_dem_name = project.onshore_dem_name 45 46 46 meshname = project.meshname+'.msh' 47 48 source_dir = project.boundarydir 47 49 48 50 # coarse data 49 51 convert_dem_from_ascii2netcdf(coarsedemname, use_cache=True, verbose=True) 50 52 dem2pts(coarsedemname, use_cache=True, verbose=True) 51 ''' 53 52 54 # fine data (clipping the points file to smaller area) 53 convert_dem_from_ascii2netcdf( finedemname, use_cache=True, verbose=True)54 dem2pts( finedemname,55 easting_min=project.eastingmin,56 easting_max=project.eastingmax,57 northing_min=project.northingmin,58 northing_max= project.northingmax,55 convert_dem_from_ascii2netcdf(onshore_dem_name, use_cache=True, verbose=True) 56 dem2pts(onshore_dem_name, 57 # easting_min=project.eastingmin, 58 # easting_max=project.eastingmax, 59 # northing_min=project.northingmin, 60 # northing_max= project.northingmax, 59 61 use_cache=True, 60 62 verbose=True) … … 62 64 63 65 # combining the coarse and fine data 64 combine_rectangular_points_files(project. finedemname + '.pts',66 combine_rectangular_points_files(project.onshore_dem_name + '.pts', 65 67 project.coarsedemname + '.pts', 66 68 project.combineddemname + '.pts') 67 ''' 69 68 70 69 71 #------------------------------------------------------------------------------- … … 109 111 domain.set_name(project.basename) 110 112 domain.set_datadir(project.outputdir) 111 domain.set_quantities_to_be_stored(['stage' , 'xmomentum', 'ymomentum'])113 domain.set_quantities_to_be_stored(['stage']) 112 114 113 115 … … 136 138 domain.set_quantity('stage', tide) 137 139 domain.set_quantity('friction', 0.0) 140 138 141 domain.set_quantity('elevation', 139 # filename = project.combineddemname + '.pts', 140 filename = project.coarsedemname + '.pts', 142 # 0. 143 filename = project.combineddemname + '.pts', 144 # filename = project.coarsedemname + '.pts', 141 145 use_cache = True, 142 146 verbose = True … … 147 151 # Setup boundary conditions (all reflective) 148 152 #------------------------------------------------------------------------------- 149 153 ''' 150 154 from pyvolution.data_manager import ferret2sww 151 155 … … 163 167 'minlon': west-1, 164 168 'maxlon': east+1, 165 'origin': project.mesh_origin, 169 # 'origin': project.mesh_origin, 170 'origin': domain.geo_reference.get_origin(), 166 171 'mean_stage': tide, 167 172 'zscale': 1, #Enhance tsunami … … 171 176 verbose = True) 172 177 173 178 ''' 174 179 print 'Available boundary tags', domain.get_boundary_tags() 175 180 … … 194 199 t0 = time.time() 195 200 196 for t in domain.evolve(yieldstep = 50, finaltime = 100):201 for t in domain.evolve(yieldstep = 100, finaltime = 100): 197 202 domain.write_time() 198 203 domain.write_boundary_statistics(tags = 'top')
Note: See TracChangeset
for help on using the changeset viewer.