Changeset 4045 for anuga_work/production/pt_hedland_2006
- Timestamp:
- Nov 30, 2006, 7:03:43 PM (18 years ago)
- Location:
- anuga_work/production/pt_hedland_2006
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/pt_hedland_2006/project.py
r3788 r4045 7 7 #from anuga.utilities.polygon import read_polygon 8 8 import sys 9 from anuga.coordinate_transforms.redfearn import convert_points_from_latlon_to_utm10 9 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees 11 10 from time import localtime, strftime … … 25 24 # offshore data from AHO fairsheets 26 25 offshore_name2 = 'pt_hedland_offshore_points_fairsheet_orig' 26 27 # reissued offshore data from PMD and NM&I Nov 2006 28 offshore_name0 = 'hydrographic_charts' 29 offshore_name1 = 'xy100003817' 30 offshore_name2 = 'xy100003819' 31 offshore_name3 = 'xy100003822' 32 offshore_name4 = 'xy100003824' 33 offshore_name5 = 'xy100003825' 34 offshore_name6 = 'xy100003826' 35 offshore_name7 = 'xy100003827' 36 offshore_name8 = 'xy100003828' 37 offshore_name9 = 'xy100003829' 38 offshore_name10 = 'xy100003830' 39 offshore_name11 = 'xy100003831' 40 offshore_name12 = 'xy100003832' 41 offshore_name13 = 'xy100003837' 42 offshore_name14 = 'xy100003838' 43 offshore_name15 = 'xy100003840' 44 offshore_name16 = 'xy100003841' 45 offshore_name17 = 'xy100003843' 46 offshore_name18 = 'xy100003845' 47 offshore_name19 = 'xy100003846' 48 offshore_name20 = 'xy100003847' 49 offshore_name21 = 'xy100003945' 50 offshore_name22 = 'xy100003946' 51 offshore_name23 = 'xy100003947' 52 offshore_name24 = 'xy100003948' 53 offshore_name25 = 'xy100004071' 54 offshore_name26 = 'xy100004072' 55 offshore_name27 = 'xy100004074' 56 offshore_name28 = 'xy100004078' 57 offshore_name29 = 'xy100004079' 58 offshore_name30 = 'xy100004080' 59 offshore_name31 = 'xy100004081' 60 offshore_name32 = 'xy100004082' 61 offshore_name33 = 'xy100004083' 62 offshore_name34 = 'xy100004084' 63 offshore_name35 = 'xy100004085' 64 offshore_name36 = 'xy100004086' 65 offshore_name37 = 'xy100004087' 66 offshore_name38 = 'xy100004088' 67 offshore_name39 = 'xy100004089' 68 offshore_name40 = 'xy100004090' 69 offshore_name41 = 'xy100011753' 70 offshore_name42 = 'xy100014967' 71 offshore_name43 = 'xy100014968' 72 offshore_name44 = 'xy100014969' 73 offshore_name45 = 'xy100021030' 74 offshore_name46 = 'xy100021031' 75 offshore_name47 = 'xy100021032' 76 offshore_name48 = 'xy100021033' 77 offshore_name49 = 'xy100021035' 78 offshore_name50 = 'xy100022767' 27 79 28 80 # coastline developed from aerial photography and 1.5m DLI contour … … 86 138 #onshore_dem_name = datadir + onshore_name_dted 87 139 onshore_dem_name = datadir + onshore_name_dli 140 offshore_dem_name0 = datadir + offshore_name0 88 141 offshore_dem_name1 = datadir + offshore_name1 89 142 offshore_dem_name2 = datadir + offshore_name2 143 offshore_dem_name3 = datadir + offshore_name3 144 offshore_dem_name4 = datadir + offshore_name4 145 offshore_dem_name5 = datadir + offshore_name5 146 offshore_dem_name6 = datadir + offshore_name6 147 offshore_dem_name7 = datadir + offshore_name7 148 offshore_dem_name8 = datadir + offshore_name8 149 offshore_dem_name9 = datadir + offshore_name9 150 offshore_dem_name10 = datadir + offshore_name10 151 offshore_dem_name11 = datadir + offshore_name11 152 offshore_dem_name12 = datadir + offshore_name12 153 offshore_dem_name13 = datadir + offshore_name13 154 offshore_dem_name14 = datadir + offshore_name14 155 offshore_dem_name15 = datadir + offshore_name15 156 offshore_dem_name16 = datadir + offshore_name16 157 offshore_dem_name17 = datadir + offshore_name17 158 offshore_dem_name18 = datadir + offshore_name18 159 offshore_dem_name19 = datadir + offshore_name19 160 offshore_dem_name20 = datadir + offshore_name20 161 offshore_dem_name21 = datadir + offshore_name21 162 offshore_dem_name22 = datadir + offshore_name22 163 offshore_dem_name23 = datadir + offshore_name23 164 offshore_dem_name24 = datadir + offshore_name24 165 offshore_dem_name25 = datadir + offshore_name25 166 offshore_dem_name26 = datadir + offshore_name26 167 offshore_dem_name27 = datadir + offshore_name27 168 offshore_dem_name28 = datadir + offshore_name28 169 offshore_dem_name29 = datadir + offshore_name29 170 offshore_dem_name30 = datadir + offshore_name30 171 offshore_dem_name31 = datadir + offshore_name31 172 offshore_dem_name32 = datadir + offshore_name32 173 offshore_dem_name33 = datadir + offshore_name33 174 offshore_dem_name34 = datadir + offshore_name34 175 offshore_dem_name35 = datadir + offshore_name35 176 offshore_dem_name36 = datadir + offshore_name36 177 offshore_dem_name37 = datadir + offshore_name37 178 offshore_dem_name38 = datadir + offshore_name38 179 offshore_dem_name39 = datadir + offshore_name39 180 offshore_dem_name40 = datadir + offshore_name40 181 offshore_dem_name41 = datadir + offshore_name41 182 offshore_dem_name42 = datadir + offshore_name42 183 offshore_dem_name43 = datadir + offshore_name43 184 offshore_dem_name44 = datadir + offshore_name44 185 offshore_dem_name45 = datadir + offshore_name45 186 offshore_dem_name46 = datadir + offshore_name46 187 offshore_dem_name47 = datadir + offshore_name47 188 offshore_dem_name48 = datadir + offshore_name48 189 offshore_dem_name49 = datadir + offshore_name49 190 offshore_dem_name50 = datadir + offshore_name50 191 90 192 coast_dem_name = datadir + coast_name 91 193 combined_dem_name = datadir + 'pt_hedland_combined_elevation' -
anuga_work/production/pt_hedland_2006/run_pt_hedland.py
r3788 r4045 6 6 7 7 The scenario is defined by a triangular mesh created from project.polygon, 8 the elevation data and a simulated submarine landslide.8 the elevation data and a tsunami generated by a subduction zone earthquake. 9 9 10 10 Ole Nielsen and Duncan Gray, GA - 2005 and Nick Bartzis, GA - 2006 … … 28 28 from anuga.geospatial_data.geospatial_data import * 29 29 import sys 30 from anuga.abstract_2d_finite_volumes.util import Screen_Catcher30 from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files 31 31 32 32 # Application specific imports … … 38 38 #------------------------------------------------------------------------------- 39 39 40 # creates copy of code in output dir if dir doesn't exist 41 if access(project.outputtimedir,F_OK) == 0 : 42 mkdir (project.outputtimedir) 43 copy (dirname(project.__file__) +sep+ project.__name__+'.py', project.outputtimedir + project.__name__+'.py') 44 copy (__file__, project.outputtimedir + basename(__file__)) 45 print 'project.outputtimedir',project.outputtimedir 46 47 # normal screen output is stored in 48 screen_output_name = project.outputtimedir + "screen_output.txt" 49 screen_error_name = project.outputtimedir + "screen_error.txt" 50 51 # used to catch screen output to file 52 sys.stdout = Screen_Catcher(screen_output_name) 53 sys.stderr = Screen_Catcher(screen_error_name) 40 copy_code_files(project.outputtimedir,__file__,dirname(project.__file__)+sep+ project.__name__+'.py' ) 41 myid = 0 42 numprocs = 1 43 start_screen_catcher(project.outputtimedir, myid, numprocs) 44 54 45 print 'USER: ', project.user 55 46 … … 79 70 verbose=True) 80 71 81 print 'create G1' 82 G1 = Geospatial_data(file_name = project.offshore_dem_name1 + '.xya') 83 print 'create G2' 84 G2 = Geospatial_data(file_name = project.offshore_dem_name2 + '.xya') 85 print 'create G3' 86 G3 = Geospatial_data(file_name = project.onshore_dem_name + '.pts') 87 print 'create G4' 88 G4 = Geospatial_data(file_name = project.coast_dem_name + '.xya') 89 print 'add G1+G2+G3+G4' 90 G = G1 + G2 + G3 + G4 72 print 'create offshore' 73 G11= Geospatial_data(file_name = project.offshore_dem_name0 + '.xya')+\ 74 Geospatial_data(file_name = project.offshore_dem_name1 + '.xya')+\ 75 Geospatial_data(file_name = project.offshore_dem_name2 + '.xya')+\ 76 Geospatial_data(file_name = project.offshore_dem_name3 + '.xya')+\ 77 Geospatial_data(file_name = project.offshore_dem_name4 + '.xya')+\ 78 Geospatial_data(file_name = project.offshore_dem_name5 + '.xya')+\ 79 Geospatial_data(file_name = project.offshore_dem_name6 + '.xya')+\ 80 Geospatial_data(file_name = project.offshore_dem_name7 + '.xya')+\ 81 Geospatial_data(file_name = project.offshore_dem_name8 + '.xya')+\ 82 Geospatial_data(file_name = project.offshore_dem_name9 + '.xya')+\ 83 Geospatial_data(file_name = project.offshore_dem_name10 + '.xya') 84 G12= Geospatial_data(file_name = project.offshore_dem_name11 + '.xya')+\ 85 Geospatial_data(file_name = project.offshore_dem_name12 + '.xya')+\ 86 Geospatial_data(file_name = project.offshore_dem_name13 + '.xya')+\ 87 Geospatial_data(file_name = project.offshore_dem_name14 + '.xya')+\ 88 Geospatial_data(file_name = project.offshore_dem_name15 + '.xya')+\ 89 Geospatial_data(file_name = project.offshore_dem_name16 + '.xya')+\ 90 Geospatial_data(file_name = project.offshore_dem_name17 + '.xya')+\ 91 Geospatial_data(file_name = project.offshore_dem_name18 + '.xya')+\ 92 Geospatial_data(file_name = project.offshore_dem_name19 + '.xya')+\ 93 Geospatial_data(file_name = project.offshore_dem_name20 + '.xya') 94 G13= Geospatial_data(file_name = project.offshore_dem_name21 + '.xya')+\ 95 Geospatial_data(file_name = project.offshore_dem_name22 + '.xya')+\ 96 Geospatial_data(file_name = project.offshore_dem_name23 + '.xya')+\ 97 Geospatial_data(file_name = project.offshore_dem_name24 + '.xya')+\ 98 Geospatial_data(file_name = project.offshore_dem_name25 + '.xya')+\ 99 Geospatial_data(file_name = project.offshore_dem_name26 + '.xya')+\ 100 Geospatial_data(file_name = project.offshore_dem_name27 + '.xya')+\ 101 Geospatial_data(file_name = project.offshore_dem_name28 + '.xya')+\ 102 Geospatial_data(file_name = project.offshore_dem_name29 + '.xya') 103 G14= Geospatial_data(file_name = project.offshore_dem_name30 + '.xya')+\ 104 Geospatial_data(file_name = project.offshore_dem_name31 + '.xya')+\ 105 Geospatial_data(file_name = project.offshore_dem_name32 + '.xya')+\ 106 Geospatial_data(file_name = project.offshore_dem_name33 + '.xya')+\ 107 Geospatial_data(file_name = project.offshore_dem_name34 + '.xya')+\ 108 Geospatial_data(file_name = project.offshore_dem_name35 + '.xya')+\ 109 Geospatial_data(file_name = project.offshore_dem_name36 + '.xya')+\ 110 Geospatial_data(file_name = project.offshore_dem_name37 + '.xya')+\ 111 Geospatial_data(file_name = project.offshore_dem_name38 + '.xya')+\ 112 Geospatial_data(file_name = project.offshore_dem_name39 + '.xya') 113 G15= Geospatial_data(file_name = project.offshore_dem_name40 + '.xya')+\ 114 Geospatial_data(file_name = project.offshore_dem_name41 + '.xya')+\ 115 Geospatial_data(file_name = project.offshore_dem_name42 + '.xya')+\ 116 Geospatial_data(file_name = project.offshore_dem_name43 + '.xya')+\ 117 Geospatial_data(file_name = project.offshore_dem_name44 + '.xya')+\ 118 Geospatial_data(file_name = project.offshore_dem_name45 + '.xya')+\ 119 Geospatial_data(file_name = project.offshore_dem_name46 + '.xya')+\ 120 Geospatial_data(file_name = project.offshore_dem_name47 + '.xya')+\ 121 Geospatial_data(file_name = project.offshore_dem_name48 + '.xya')+\ 122 Geospatial_data(file_name = project.offshore_dem_name49 + '.xya')+\ 123 Geospatial_data(file_name = project.offshore_dem_name50 + '.xya') 124 125 print 'create onshore' 126 G2 = Geospatial_data(file_name = project.onshore_dem_name + '.pts') 127 print 'create coast' 128 G3 = Geospatial_data(file_name = project.coast_dem_name + '.xya') 129 print 'add' 130 G = G11 + G12 + G13 + G14 + G15 + G2 + G3 91 131 print 'export G' 92 132 G.export_points_file(project.combined_dem_name + '.pts')
Note: See TracChangeset
for help on using the changeset viewer.