[3627] | 1 | """Common filenames and locations for topographic data, meshes and outputs. |
---|
| 2 | """ |
---|
| 3 | |
---|
| 4 | from os import sep, environ, getenv, getcwd |
---|
| 5 | from os.path import expanduser, basename |
---|
| 6 | #from anuga.utilities.polygon import read_polygon |
---|
| 7 | import sys |
---|
| 8 | from anuga.coordinate_transforms.redfearn import\ |
---|
| 9 | degminsec2decimal_degrees,\ |
---|
| 10 | convert_points_from_latlon_to_utm |
---|
| 11 | |
---|
| 12 | from time import localtime, strftime |
---|
| 13 | from anuga.geospatial_data.geospatial_data import * |
---|
| 14 | |
---|
| 15 | # Location and naming of scenario data |
---|
| 16 | state = 'western_australia' |
---|
| 17 | scenario_dir_name = 'karratha_tsunami_scenario_2005' |
---|
| 18 | |
---|
| 19 | basename = 'elevation50m' |
---|
| 20 | boundary_basename = 'SU-AU' |
---|
| 21 | |
---|
| 22 | |
---|
| 23 | |
---|
| 24 | # Data output |
---|
| 25 | codename = 'project.py' |
---|
| 26 | |
---|
| 27 | home = getenv('INUNDATIONHOME') #Sandpit's parent dir |
---|
| 28 | |
---|
| 29 | if sys.platform == 'win32': |
---|
| 30 | user = getenv('USERPROFILE') |
---|
| 31 | else: |
---|
| 32 | user = getenv('LOGNAME') |
---|
| 33 | |
---|
| 34 | |
---|
| 35 | #Derive subdirectories and filenames |
---|
| 36 | time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir |
---|
| 37 | outputtimedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'output'+sep+time+sep |
---|
| 38 | |
---|
| 39 | meshdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'meshes'+sep |
---|
| 40 | datadir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'topographies'+sep |
---|
| 41 | gaugedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'gauges'+sep |
---|
| 42 | polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep |
---|
| 43 | boundarydir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep |
---|
| 44 | outputdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'output'+sep |
---|
| 45 | tidedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'tide_data'+sep |
---|
| 46 | |
---|
| 47 | #gauge_filename = gaugedir + 'gauge_location_broome.csv' |
---|
| 48 | |
---|
| 49 | meshname = meshdir + basename |
---|
| 50 | outputname = outputtimedir + basename #Used by post processing |
---|
| 51 | #boundaryname = outputtimedir + basename #Used by post processing |
---|
| 52 | |
---|
| 53 | |
---|
| 54 | |
---|
| 55 | # Regions |
---|
| 56 | |
---|
| 57 | refzone = 50 |
---|
| 58 | south = degminsec2decimal_degrees(-20,55,0) |
---|
| 59 | north = degminsec2decimal_degrees(-20,15,0) |
---|
| 60 | west = degminsec2decimal_degrees(116,17,0) |
---|
| 61 | east = degminsec2decimal_degrees(117,10,0) |
---|
| 62 | |
---|
| 63 | p0 = [south, degminsec2decimal_degrees(116,32,0)] |
---|
| 64 | p1 = [south, west] |
---|
| 65 | p2 = [degminsec2decimal_degrees(-20,23,0), west] |
---|
| 66 | p3 = [north, degminsec2decimal_degrees(116,45,0)] |
---|
| 67 | p4 = [north, degminsec2decimal_degrees(117,0,0)] |
---|
| 68 | p5 = [p2[0], degminsec2decimal_degrees(117,8,0)] |
---|
| 69 | p6 = [degminsec2decimal_degrees(-20,30,0), east] |
---|
| 70 | p7 = [degminsec2decimal_degrees(-20,38,0), east] |
---|
| 71 | p8 = [south, east] |
---|
| 72 | |
---|
| 73 | |
---|
| 74 | bounding_polygon, zone = convert_points_from_latlon_to_utm([p0, p1, p2, p3, p4, p5, p6, p7, p8]) |
---|
| 75 | refzone = zone |
---|
| 76 | |
---|
| 77 | |
---|
| 78 | #Interior regions - NOT CONVERTED YET TO UTM |
---|
| 79 | karratha_south = degminsec2decimal_degrees(-20,44,0) |
---|
| 80 | karratha_north = degminsec2decimal_degrees(-20,42,0) |
---|
| 81 | karratha_west = degminsec2decimal_degrees(116,48,0) |
---|
| 82 | karratha_east = degminsec2decimal_degrees(116,53,30) |
---|
| 83 | |
---|
| 84 | k0 = [karratha_south, karratha_west] |
---|
| 85 | k1 = [karratha_south, karratha_east] |
---|
| 86 | k2 = [karratha_north, karratha_east] |
---|
| 87 | k3 = [karratha_north, karratha_west] |
---|
| 88 | |
---|
| 89 | karratha_polygon = [k0, k1, k2, k3] |
---|
| 90 | |
---|
| 91 | |
---|
| 92 | #Interior regions |
---|
| 93 | dampier_south = degminsec2decimal_degrees(-20,40,0) |
---|
| 94 | dampier_north = degminsec2decimal_degrees(-20,38,10) |
---|
| 95 | dampier_west = degminsec2decimal_degrees(116,43,0) |
---|
| 96 | dampier_east = degminsec2decimal_degrees(116,45,0) |
---|
| 97 | |
---|
| 98 | d0 = [dampier_south, dampier_west] |
---|
| 99 | d1 = [dampier_south, dampier_east] |
---|
| 100 | d2 = [dampier_north, dampier_east] |
---|
| 101 | d3 = [dampier_north, dampier_west] |
---|
| 102 | |
---|
| 103 | dampier_polygon = [d0, d1, d2, d3] |
---|
| 104 | |
---|
| 105 | |
---|
| 106 | #Interior regions |
---|
| 107 | refinery_south = degminsec2decimal_degrees(-20,37,50) |
---|
| 108 | refinery_north = degminsec2decimal_degrees(-20,36,0) |
---|
| 109 | refinery_west = degminsec2decimal_degrees(116,44,0) |
---|
| 110 | refinery_east = degminsec2decimal_degrees(116,46,10) |
---|
| 111 | |
---|
| 112 | d0 = [refinery_south, refinery_west] |
---|
| 113 | d1 = [refinery_south, refinery_east] |
---|
| 114 | d2 = [refinery_north, refinery_east] |
---|
| 115 | d3 = [refinery_north, refinery_west] |
---|
| 116 | |
---|
| 117 | refinery_polygon = [d0, d1, d2, d3] |
---|
| 118 | |
---|
| 119 | |
---|
| 120 | #Interior region around 468899, 7715177: |
---|
| 121 | #lat (-20, 39, 44.93753), lon (116, 42, 5.09106) |
---|
| 122 | |
---|
| 123 | point_south = degminsec2decimal_degrees(-20,39,46) |
---|
| 124 | point_north = degminsec2decimal_degrees(-20,39,42) |
---|
| 125 | point_west = degminsec2decimal_degrees(116,42,0) |
---|
| 126 | point_east = degminsec2decimal_degrees(116,42,10) |
---|
| 127 | |
---|
| 128 | d0 = [point_south, point_west] |
---|
| 129 | d1 = [point_south, point_east] |
---|
| 130 | d2 = [point_north, point_east] |
---|
| 131 | d3 = [point_north, point_west] |
---|
| 132 | |
---|
| 133 | point_polygon = [d0, d1, d2, d3] |
---|
| 134 | |
---|
| 135 | |
---|
| 136 | #Neils areas around interesting points |
---|
| 137 | neil1_point1 = [degminsec2decimal_degrees(-20,35,34), |
---|
| 138 | degminsec2decimal_degrees(116,45,18)] |
---|
| 139 | neil1_point2 = [degminsec2decimal_degrees(-20,36,15), |
---|
| 140 | degminsec2decimal_degrees(116,46,18)] |
---|
| 141 | neil1_point3 = [degminsec2decimal_degrees(-20,35,9), |
---|
| 142 | degminsec2decimal_degrees(116,47,17)] |
---|
| 143 | neil1_point4 = [degminsec2decimal_degrees(-20,34,26), |
---|
| 144 | degminsec2decimal_degrees(116,46,17)] |
---|
| 145 | |
---|
| 146 | neil1_polygon = [neil1_point1, neil1_point2, neil1_point3, neil1_point4] |
---|
| 147 | |
---|
| 148 | |
---|
| 149 | |
---|
| 150 | neil2_point1 = [degminsec2decimal_degrees(-20,39,36), |
---|
| 151 | degminsec2decimal_degrees(116,41,33)] |
---|
| 152 | neil2_point2 = [degminsec2decimal_degrees(-20,40,10), |
---|
| 153 | degminsec2decimal_degrees(116,42,13)] |
---|
| 154 | neil2_point3 = [degminsec2decimal_degrees(-20,38,39), |
---|
| 155 | degminsec2decimal_degrees(116,43,49)] |
---|
| 156 | neil2_point4 = [degminsec2decimal_degrees(-20,38,5), |
---|
| 157 | degminsec2decimal_degrees(116,43,9)] |
---|
| 158 | |
---|
| 159 | neil2_polygon = [neil2_point1, neil2_point2, neil2_point3, neil2_point4] |
---|
| 160 | |
---|
| 161 | |
---|
| 162 | |
---|
| 163 | |
---|
| 164 | #Withnell bay |
---|
| 165 | wb_point1 = [degminsec2decimal_degrees(-20,35,34), |
---|
| 166 | degminsec2decimal_degrees(116,45,18)] |
---|
| 167 | wb_point2 = [degminsec2decimal_degrees(-20,36,15), |
---|
| 168 | degminsec2decimal_degrees(116,46,18)] |
---|
| 169 | wb_point3 = [degminsec2decimal_degrees(-20,35,9), |
---|
| 170 | degminsec2decimal_degrees(116,47,17)] |
---|
| 171 | wb_point4 = [degminsec2decimal_degrees(-20,34,26), |
---|
| 172 | degminsec2decimal_degrees(116,46,17)] |
---|
| 173 | |
---|
| 174 | wb_polygon = [wb_point1, wb_point2, wb_point3, wb_point4] |
---|
| 175 | |
---|
| 176 | |
---|
| 177 | |
---|
| 178 | |
---|
| 179 | #Larger Withnell bay |
---|
| 180 | lwb_point1 = [degminsec2decimal_degrees(-20,35,59), |
---|
| 181 | degminsec2decimal_degrees(116,42,00)] |
---|
| 182 | lwb_point2 = [degminsec2decimal_degrees(-20,36,50), |
---|
| 183 | degminsec2decimal_degrees(116,46,50)] |
---|
| 184 | lwb_point3 = [degminsec2decimal_degrees(-20,34,00), |
---|
| 185 | degminsec2decimal_degrees(116,47,39)] |
---|
| 186 | lwb_point4 = [degminsec2decimal_degrees(-20,33,00), |
---|
| 187 | degminsec2decimal_degrees(116,42,50)] |
---|
| 188 | |
---|
| 189 | lwb_polygon = [lwb_point1, lwb_point2, lwb_point3, lwb_point4] |
---|
| 190 | |
---|
| 191 | |
---|
| 192 | |
---|