[4297] | 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 |
---|
| 6 | import sys |
---|
| 7 | from time import localtime, strftime, gmtime |
---|
| 8 | from anuga.utilities.polygon import read_polygon, plot_polygons, \ |
---|
| 9 | polygon_area, is_inside_polygon |
---|
| 10 | from anuga.coordinate_transforms.redfearn import convert_from_latlon_to_utm |
---|
| 11 | from time import localtime, strftime |
---|
| 12 | |
---|
| 13 | ############################### |
---|
| 14 | # Domain definitions |
---|
| 15 | ############################### |
---|
[5241] | 16 | ''' |
---|
[4297] | 17 | # Bathymetry and topography filenames |
---|
| 18 | bathymetry = 'phuket_utm47' |
---|
| 19 | |
---|
| 20 | # bounding polygon for study area |
---|
| 21 | bounding_polygon_latlon = read_polygon('phuket_boundary.csv') |
---|
| 22 | # only works for zone 47 |
---|
| 23 | south = 1000 |
---|
| 24 | east = -1000 |
---|
| 25 | west = 1000 |
---|
| 26 | north = -1000 |
---|
| 27 | for points in bounding_polygon_latlon: |
---|
| 28 | south = min(points[1],south) |
---|
| 29 | north = max(points[1],north) |
---|
| 30 | east = max(points[0],east) |
---|
| 31 | west = min(points[0],west) |
---|
| 32 | |
---|
| 33 | print "south, north, east, west", south, north, east, west |
---|
| 34 | |
---|
| 35 | time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir |
---|
| 36 | #combined_dir_name = 'phuket_bathymetry'+time |
---|
| 37 | combined_dir_name = 'phuket_bathymetry' |
---|
| 38 | meshname = 'phuket_mesh'+time |
---|
| 39 | mesh_elevname = 'phuket_mesh'+time+'_elev' |
---|
| 40 | bathymetry_points = 'bathymetry_points'+time |
---|
| 41 | boundary_most_in = 'out' |
---|
| 42 | boundary_most_out = 'most_boundary_condition' |
---|
| 43 | ############################### |
---|
| 44 | # Interior region definitions |
---|
| 45 | ############################### |
---|
[5241] | 46 | ''' |
---|
| 47 | # location of polygons developed by John Jakeman April 2008 |
---|
| 48 | home = getenv('INUNDATIONHOME') + sep +'data'+sep |
---|
| 49 | anuga_dir = home+'anuga_validation'+sep+'boxing_day_sceanario'+sep+'anuga'+sep |
---|
| 50 | poly_dir = anuga_dir+'outputs'+sep+'john'+sep |
---|
| 51 | gauges_dir = anuga_dir+'gauges'+sep |
---|
[4297] | 52 | # interior polygons |
---|
[5241] | 53 | #poly_island1 = read_polygon('island1.csv') |
---|
| 54 | #poly_north_island = read_polygon('north_island.csv') |
---|
| 55 | #poly_50m_contour = read_polygon('phuket_50m_contour_small.csv') |
---|
[4297] | 56 | |
---|
[5241] | 57 | # bounding polygon |
---|
| 58 | poly_all = read_polygon(poly_dir + 'extent.csv')#'domain.csv') |
---|
| 59 | |
---|
[4297] | 60 | gauge_filename = 'gauges.csv' |
---|