source: anuga_work/production/hobart_2006/project.py @ 3897

Last change on this file since 3897 was 3897, checked in by sexton, 17 years ago

polygon updates to Broome

File size: 9.0 KB
RevLine 
[3559]1# -*- coding: cp1252 -*-
2"""Common filenames and locations for topographic data, meshes and outputs.
3"""
4
5from os import sep, environ, getenv, getcwd
6from os.path import expanduser
7import sys
8from time import localtime, strftime, gmtime
[3721]9from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area
[3779]10from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_from_latlon_to_utm
[3671]11
12if sys.platform == 'win32':
13    home = getenv('INUNDATIONHOME')
14    user = getenv('USERPROFILE')
15
16else:   
17    home = getenv('INUNDATIONHOME', sep+'d'+sep+'xrd'+sep+'gem'+sep+'2'+sep+'ramp'+sep+'risk_assessment_methods_project'+sep+'inundation')     
18    user = getenv('LOGNAME')
19    print 'USER:', user
20
21# INUNDATIONHOME is the inundation directory, not the data directory.
22home += sep +'data'
23
[3559]24#Making assumptions about the location of scenario data
25state = 'tasmania'
26scenario_dir_name = 'hobart_tsunami_scenario_2006'
27
[3721]28# data provided by Tas SES and checked by NM&I
29#onshore_name = 'hob3_topo' # original
30onshore_name_25 = 'hob5_topo_25m' # 25m grid and clipped to 100m elevation or 3000m from coast
31onshore_name = 'hob5_topo' # 12.5m grid and clipped to 100m elevation or 3000m from coast
32#onshore_name_25 = 'hob6_topo_25m' # 25m grid NOT clipped
33offshore_name_tas1 = 'derwent_2m'
34offshore_name_tas2 = 'derwent_5m'
35offshore_name_tas3 = 'south_east_tas' #actually this is AHO
36offshore_name_tas4 = 'hobart_1m'
[3559]37
[3721]38# AHO data and checked by NM&I
39offshore_name1 = 'xy100003760'
40offshore_name2 = 'xy100003761'
41offshore_name3 = 'xy100003762'
42offshore_name4 = 'xy100003907'
43offshore_name5 = 'xy100003908'
44offshore_name6 = 'xy100003909'
45offshore_name7 = 'xy100003910'
46offshore_name8 = 'xy100003932'
47offshore_name9 = 'xy100003933'
48offshore_name10 = 'xy100003934'
49offshore_name11 = 'xy100003935'
50offshore_name12 = 'xy100003936'
51offshore_name13 = 'xy100003964'
52offshore_name14 = 'xy100014250'
53offshore_name15 = 'xy100014253'
54offshore_name16 = 'xy100016142'
[3626]55
[3721]56# developed by NM&I
57coast_name = 'coastline_points'
[3559]58
[3779]59boundary_basename = 'puysegur' # Mw 8.7
60#boundary_basename = 'puysegur_clip' # Mw 8.5
[3559]61
62#swollen/ all data output
63basename = 'source'
64
65codename = 'project.py'
66
67#Derive subdirectories and filenames
68#time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
69local_time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
70meshdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'meshes'+sep
71datadir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'topographies'+sep
72gaugedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'gauges'+sep
73polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
[3669]74boundarydir = home+sep+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep
[3559]75outputdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'outputs'+sep
76outputtimedir = outputdir + local_time + sep
77polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
78
[3779]79#gauge_filename = gaugedir + 'hobart_gauges_final.csv'
[3867]80#gauge_filename = gaugedir + 'hobartGauges20061019.csv'
81gauge_filename = gaugedir + 'gauges_MRT.csv'
[3852]82#gauge_filename = gaugedir + 'check_boundary.csv'
[3559]83
[3671]84codedir = getcwd()+sep                           
[3559]85codedirname = codedir + 'project.py'
86meshname = outputtimedir + 'mesh_' + basename
87
[3701]88# Necessary if using point datasets, rather than grid
[3721]89onshore_dem_name = datadir + onshore_name
90onshore_dem_name_25 = datadir + onshore_name_25
91all_onshore_dem_name = datadir + 'combined_onshore'
92offshore_dem_name_local1 = datadir + offshore_name_tas1
93offshore_dem_name_local2 = datadir + offshore_name_tas2
94offshore_dem_name_local3 = datadir + offshore_name_tas3
95offshore_dem_name_local4 = datadir + offshore_name_tas4
96offshore_dem_name_aho1 = datadir + offshore_name1
97offshore_dem_name_aho2 = datadir + offshore_name2
98offshore_dem_name_aho3 = datadir + offshore_name3
99offshore_dem_name_aho4 = datadir + offshore_name4
100offshore_dem_name_aho5 = datadir + offshore_name5
101offshore_dem_name_aho6 = datadir + offshore_name6
102offshore_dem_name_aho7 = datadir + offshore_name7
103offshore_dem_name_aho8 = datadir + offshore_name8
104offshore_dem_name_aho9 = datadir + offshore_name9
105offshore_dem_name_aho10 = datadir + offshore_name10
106offshore_dem_name_aho11 = datadir + offshore_name11
107offshore_dem_name_aho12 = datadir + offshore_name12
108offshore_dem_name_aho13 = datadir + offshore_name13
109offshore_dem_name_aho14 = datadir + offshore_name14
110offshore_dem_name_aho15 = datadir + offshore_name15
111offshore_dem_name_aho16 = datadir + offshore_name16
112coast_dem_name = datadir + coast_name
[3701]113
[3671]114# addition once total grid delivered
[3683]115onshore_offshore_dem_name_25 = datadir + '25m_se_tas' #25m grid
[3701]116onshore_offshore_dem_name    = datadir + '50m_se_tas' #50m grid
117
[3683]118# output names
[3701]119bruny_dem_name_25   = datadir + 'bruny_25_dem'
120hobart_dem_name_25  = datadir + 'hobart_25_dem'
121combined_dem_name   = datadir + 'hobart_combined_elevation'
[3683]122combined_dem_name_2 = datadir + 'hobart_combined_elevation_2'
[3721]123combined_dem_name_3 = datadir + 'hobart_combined_elevation_jane'
124combined_dem_name_4 = datadir + 'hobart_combined_elevation_test'
[3559]125
[3701]126#outputname = outputtimedir + basename  #Used by post processing
[3559]127
[3671]128###############################
129# Domain definitions
130###############################
[3650]131
[3701]132# bounding box for clipping MOST output (much bigger than study area)
[3679]133south = degminsec2decimal_degrees(-44,45,0)
134north = degminsec2decimal_degrees(-42,0,0)
[3701]135west  = degminsec2decimal_degrees(146,45,0)
136east  = degminsec2decimal_degrees(148,25,0) 
[3559]137
[3701]138###Main Domain of Hobart:
[3615]139d0 = [south, west]
140d1 = [south, east]
141d2 = [north, east]
142d3 = [north, west]
[3769]143polyAll2, zone = convert_from_latlon_to_utm([d0, d1, d2, d3])
[3559]144refzone = zone
145
[3679]146# Mark run - morning Fri 29 Sep; surrounds -100m and 20mish elevation
[3721]147polyAll = read_polygon(polygondir+'new_extent_2.csv')
[3679]148plot_polygons([polyAll, polyAll2],'boundingpoly',verbose=False)
[3809]149#print 'Area of bounding polygon', polygon_area(polyAll)
[3779]150polyAll_refine = read_polygon(polygondir+'extent_refine.csv')
[3809]151#print 'Area of refined bounding polygon', polygon_area(polyAll_refine)
152#print 'change', polygon_area(polyAll_refine)/polygon_area(polyAll)
[3671]153
[3779]154polyAll_refine2 = read_polygon(polygondir+'export_refine_2.csv')
[3809]155#print 'Area of refined bounding polygon 2', polygon_area(polyAll_refine2)
156#print 'change 2', polygon_area(polyAll_refine2)/polygon_area(polyAll)
[3779]157
[3671]158###################################################################
159# Clipping regions for export to asc and regions for clipping data
160###################################################################
161
[3721]162# clipping 12.5m onshore data set
163eastingmin = 520000
164eastingmax = 536000
165northingmin = 5245000
166northingmax = 5260000
[3809]167# for original model setup
[3683]168# clipping 25m data set - Hobart
[3701]169eastingmin25 = 524208.387
170eastingmax25 = 554867.24
[3683]171northingmin25 = 5229154.555
172northingmax25 = 5258511.857
[3779]173hob_clip = [[eastingmin25,northingmin25],[eastingmin25,northingmax25],[eastingmax25,northingmax25],[eastingmax25,northingmin25]]
[3683]174# clipping 25m data set - Bruny
175eastingmin25_2 = 523672.502
176eastingmax25_2 = 536020.057
177northingmin25_2 = 5203120.035
178northingmax25_2 = 5212052.309
[3779]179bruny_clip = [[eastingmin25_2,northingmin25_2],[eastingmin25_2,northingmax25_2],[eastingmax25_2,northingmax25_2],[eastingmax25_2,northingmin25_2]]
[3732]180# overall 25m clipping region - covering Hobart and Bruny
181eastingmin25_3 = 523672.502
182eastingmax25_3 = 554867.24
183northingmin25_3 = 5203120.035
184northingmax25_3 = 5258511.857
185
[3809]186# export for refined grids at paleo sites
187eastingmin_site13 = 532835.94
188eastingmax_site13 = 533848.23
189northingmin_site13 = 5237384.18
190northingmax_site13 = 5238195.03
191
[3886]192eastingmin_kingston = 525500 #525905.49
[3809]193eastingmax_kingston = 527773.42
[3897]194northingmin_kingston = 5240256.71
[3886]195northingmax_kingston = 5242300 #5242148.58
[3809]196
197eastingmin_bruny = 523638.9
198eastingmax_bruny = 533000.88
199northingmin_bruny = 5204967.41
200northingmax_bruny = 5212177.69
201
[3671]202###############################
203# Interior region definitions
204###############################
205
[3679]206# Hobart digitized polygons
207poly_hobart1 = read_polygon(polygondir+'Hob_poly1.csv')
208poly_hobart2 = read_polygon(polygondir+'Hob_poly2.csv')
209poly_hobart3 = read_polygon(polygondir+'Hob_poly3.csv')
[3721]210poly_hobart4 = read_polygon(polygondir+'Hob_poly4.csv')
[3795]211poly_hobart5 = read_polygon(polygondir+'new_interior_region.csv')
[3559]212
[3809]213plot_polygons([polyAll, poly_hobart1,poly_hobart2,poly_hobart3,poly_hobart5],'orig_model',verbose=False)
[3779]214
[3809]215# to refine around paleo sites
216poly_kingston = read_polygon(polygondir+'Kingston.csv')
217poly_bruny = read_polygon(polygondir+'bruny_refine.csv')
218poly_site13 = read_polygon(polygondir+'Site13.csv')
[3852]219#print 'areas of refined polys', polygon_area(poly_kingston), polygon_area(poly_bruny), polygon_area(poly_site13)
[3809]220
221plot_polygons([polyAll, poly_kingston,poly_bruny,poly_site13,poly_hobart5],'orig_model',verbose=False)
222#print 'Area of interior regions', (polygon_area(poly_hobart1)+\
223#      polygon_area(poly_hobart2)+polygon_area(poly_hobart3)+\
224#      polygon_area(poly_hobart4)+polygon_area(poly_hobart5))/1000000.
225
226#print 'Area of original 25m data clipping', (polygon_area(hob_clip)+polygon_area(bruny_clip))/1000000.
227#print 'Area of new interior region', polygon_area(poly_hobart5)/1000000.
Note: See TracBrowser for help on using the repository browser.