Changeset 3808
- Timestamp:
- Oct 17, 2006, 12:05:03 PM (18 years ago)
- Location:
- anuga_work/production/dampier_2006
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/dampier_2006/project.py
r3802 r3808 15 15 # Location and naming of scenario data 16 16 state = 'western_australia' 17 #scenario_dir_name = 'dampier_tsunami_scenario_2006'18 scenario_dir_name = 'karratha_tsunami_scenario_2005' # Tmp location17 scenario_dir_name = 'dampier_tsunami_scenario_2006' 18 #scenario_dir_name = 'karratha_tsunami_scenario_2005' # Tmp location 19 19 20 20 basename = 'elevation50m' 21 boundary_basename = ' SU-AU'21 boundary_basename = 'dampier' 22 22 23 23 … … 43 43 gaugedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'gauges'+sep 44 44 polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep 45 boundarydir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep 45 boundarydir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep+'urs'+sep+'test'+sep 46 46 outputdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'output'+sep 47 47 tidedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'tide_data'+sep … … 51 51 meshname = meshdir + basename 52 52 outputname = outputtimedir + basename #Used by post processing 53 #boundaryname = outputtimedir + basename #Used by post processing 54 55 53 boundaryname = boundarydir + boundary_basename #Used by post processing 56 54 57 55 … … 60 58 refzone = 50 61 59 south = degminsec2decimal_degrees(-20,55,0) 62 north = degminsec2decimal_degrees(- 20,15,0)60 north = degminsec2decimal_degrees(-19,15,0) 63 61 west = degminsec2decimal_degrees(116,17,0) 64 east = degminsec2decimal_degrees(11 7,10,0)62 east = degminsec2decimal_degrees(118,10,0) 65 63 66 64 p0 = [south, degminsec2decimal_degrees(116,32,0)] -
anuga_work/production/dampier_2006/run_dampier.py
r3802 r3808 47 47 meshname = project.meshname+'.msh' 48 48 source_dir = project.boundarydir 49 boundary_file = project.boundaryname 50 51 from anuga.shallow_water.data_manager import urs2sww 52 53 urs2sww(boundary_file, verbose='true', 54 minlat=project.south, 55 maxlat=project.north, 56 minlon=project.west, 57 maxlon=project.east) 58 59 import sys; sys.exit() 60 Bf = File_boundary(source_dir + project.boundary_basename + '.sww', 61 domain, verbose = True) 62 63 Br = Reflective_boundary(domain) 64 Bd = Dirichlet_boundary([tide,0,0]) 65 domain.set_boundary({'back': Br, 66 'side': Bd, 67 'ocean': Bf}) 49 68 50 69 … … 75 94 76 95 96 #--------------------------- 97 # this is check that no interior polygon is outside the bounding poly 98 #------------------------------ 99 count = 0 100 for i in range(len(interior_regions)): 101 region = interior_regions[i] 102 interior_polygon = region[0] 103 if len(inside_polygon(interior_polygon, bounding_polygon, 104 closed = True, verbose = False)) <> len(interior_polygon): 105 print 'WARNING: interior polygon %d is outside bounding polygon' %(i) 106 count += 1 107 108 if count == 0: 109 print 'interior regions OK' 110 else: 111 print 'check out your interior polygons' 112 print 'check %s in production directory' %figname 113 import sys; sys.exit() 114 115 #------------------------------------- 77 116 78 117 print 'start create mesh from regions' … … 119 158 print 'Available boundary tags', domain.get_boundary_tags() 120 159 160 161 162 urs2sww(boundary_file, 163 minlat=project.south, 164 maxlat=project.north, 165 minlon=project.west, 166 maxlon=project.east) 167 121 168 Bf = File_boundary(source_dir + project.boundary_basename + '.sww', 122 169 domain, verbose = True)
Note: See TracChangeset
for help on using the changeset viewer.