Changeset 4246
- Timestamp:
- Feb 8, 2007, 10:01:26 AM (18 years ago)
- Location:
- anuga_work/production/dampier_2006
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/dampier_2006/build_dampier.py
r4212 r4246 127 127 G_small.export_points_file(project.combined_dir_name + '_small' + '.pts') 128 128 G_other.export_points_file(project.combined_dir_name + '_other' + '.pts') 129 ''' 129 130 130 print 'start reading:',project.combined_dir_name + '.pts' 131 131 G = Geospatial_data(file_name = project.combined_dir_name + '.pts') 132 132 print 'start export',project.combined_dir_name + '.txt' 133 133 G.export_points_file(project.combined_dir_name + '.txt') 134 ''' 135 #------------------------------------------------------------------------- 136 # Convert URS to SWW file for boundary conditions 137 #------------------------------------------------------------------------- 138 print 'starting to create boundary conditions' 139 boundaries_in_dir_name = project.boundaries_in_dir_name 140 141 from anuga.shallow_water.data_manager import urs2sww 142 143 print 'minlat=project.north_boundary, maxlat=project.south_boundary',project.north_boundary, project.south_boundary 144 print 'minlon= project.west_boundary, maxlon=project.east_boundary',project.west_boundary, project.east_boundary 145 146 #import sys; sys.exit() 147 148 #if access(project.boundaries_dir,F_OK) == 0: 149 # mkdir (project.boundaries_dir) 150 151 from caching import cache 152 cache(urs2sww, 153 (boundaries_in_dir_name, 154 project.boundaries_dir_name1), 155 {'verbose': True, 156 'minlat': project.south_boundary, 157 'maxlat': project.north_boundary, 158 'minlon': project.west_boundary, 159 'maxlon': project.east_boundary, 160 'mint': 0, 'maxt': 40000, 161 # 'origin': domain.geo_reference.get_origin(), 162 'mean_stage': project.tide, 163 'zscale': 1, #Enhance tsunami 164 'fail_on_NaN': False}, 165 verbose = True, 166 ) 167 # dependencies = source_dir + project.boundary_basename + '.sww') 134 168 135 169 -
anuga_work/production/dampier_2006/project.py
r4228 r4246 78 78 boundaries_source = '' 79 79 boundaries_name = 'o' 80 boundaries_name1 = 'o_new' 81 80 82 #boundaries locations 81 83 boundaries_in_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries' … … 84 86 boundaries_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries'+sep 85 87 boundaries_dir_name = boundaries_dir + boundaries_name 88 boundaries_dir_name1 = boundaries_dir + boundaries_name1 86 89 #boundaries_time_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries'+sep+build_time+sep 87 90 #boundaries_time_dir_name = boundaries_time_dir + boundaries_name #Used by post processing … … 125 128 print 'Area of bounding polygon', polygon_area(poly_all)/1000000.0 126 129 127 #res_poly_all = 100000128 res_poly_all = 500000130 res_poly_all = 100000 131 #res_poly_all = 500000 129 132 130 133 ############################### … … 136 139 137 140 poly_dampier = read_polygon(polygons_dir+'dampier_town.csv') 138 #res_dampier = 500139 res_dampier = 5000141 res_dampier = 500 142 #res_dampier = 5000 140 143 141 144 poly_karratha = read_polygon(polygons_dir+'karrathav2.csv') … … 143 146 144 147 poly_karratha_town = read_polygon(polygons_dir+'karratha_townv2.csv') 145 #res_karratha_town = 500146 res_karratha_town = 5000148 res_karratha_town = 500 149 #res_karratha_town = 5000 147 150 148 151 poly_facility = read_polygon(polygons_dir+'facility.csv') … … 153 156 154 157 poly_coast = read_polygon(polygons_dir+'coastpoly.csv') 155 #res_coast = 1000156 res_coast = 10000158 res_coast = 1000 159 #res_coast = 10000 157 160 158 161 poly_NWislands = read_polygon(polygons_dir+'nw_islands_area.csv') … … 171 174 res_island0 = res_poly_all 172 175 173 #res_islands = 5000174 res_islands = 15000176 res_islands = 5000 177 #res_islands = 15000 175 178 176 179 poly_ref_nw4 = read_polygon(polygons_dir+'ref_nw4.csv') -
anuga_work/production/dampier_2006/run_dampier.py
r4212 r4246 33 33 from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files 34 34 from anuga_parallel.parallel_api import distribute, numprocs, myid, barrier 35 35 from anuga_parallel.parallel_abstraction import get_processor_name 36 36 # Application specific imports 37 37 import project # Definition of file names and polygons … … 43 43 44 44 start_screen_catcher(project.output_run_time_dir, myid, numprocs) 45 print "Processor Name:",get_processor_name() 45 46 46 47 # filenames … … 148 149 # MUST USE TXT FILES FOR CACHING TO WORK! 149 150 filename = project.combined_dir_name + '.txt', 150 use_cache = True,151 use_cache = False, 151 152 verbose = True, 152 153 alpha = 0.1) … … 181 182 print 'Reading Boundary file' 182 183 Bf = File_boundary(boundaries_dir_name + '.sww', 183 domain, time_thinning= 1, use_cache=True, verbose=True)184 domain, time_thinning=4, use_cache=True, verbose=True) 184 185 185 186 print 'finished reading boundary file'
Note: See TracChangeset
for help on using the changeset viewer.