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

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

updates (i) export to xya for Broome (ii) gauge investigation for Hobart (compare MOST and ANUGA) (iii) script to compare onslow outputs with change in parameters

File size: 9.0 KB
Line 
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
9from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area
10from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_from_latlon_to_utm
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
24#Making assumptions about the location of scenario data
25state = 'tasmania'
26scenario_dir_name = 'hobart_tsunami_scenario_2006'
27
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'
37
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'
55
56# developed by NM&I
57coast_name = 'coastline_points'
58
59boundary_basename = 'puysegur' # Mw 8.7
60#boundary_basename = 'puysegur_clip' # Mw 8.5
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
74boundarydir = home+sep+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep
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
79#gauge_filename = gaugedir + 'hobart_gauges_final.csv'
80#gauge_filename = gaugedir + 'hobartGauges20061019.csv'
81gauge_filename = gaugedir + 'gauges_MRT.csv'
82gauge_filename_bom = gaugedir + 'gauges_bruny_fine.csv'
83#gauge_filename = gaugedir + 'check_boundary.csv'
84
85codedir = getcwd()+sep                           
86codedirname = codedir + 'project.py'
87meshname = outputtimedir + 'mesh_' + basename
88
89# Necessary if using point datasets, rather than grid
90onshore_dem_name = datadir + onshore_name
91onshore_dem_name_25 = datadir + onshore_name_25
92all_onshore_dem_name = datadir + 'combined_onshore'
93offshore_dem_name_local1 = datadir + offshore_name_tas1
94offshore_dem_name_local2 = datadir + offshore_name_tas2
95offshore_dem_name_local3 = datadir + offshore_name_tas3
96offshore_dem_name_local4 = datadir + offshore_name_tas4
97offshore_dem_name_aho1 = datadir + offshore_name1
98offshore_dem_name_aho2 = datadir + offshore_name2
99offshore_dem_name_aho3 = datadir + offshore_name3
100offshore_dem_name_aho4 = datadir + offshore_name4
101offshore_dem_name_aho5 = datadir + offshore_name5
102offshore_dem_name_aho6 = datadir + offshore_name6
103offshore_dem_name_aho7 = datadir + offshore_name7
104offshore_dem_name_aho8 = datadir + offshore_name8
105offshore_dem_name_aho9 = datadir + offshore_name9
106offshore_dem_name_aho10 = datadir + offshore_name10
107offshore_dem_name_aho11 = datadir + offshore_name11
108offshore_dem_name_aho12 = datadir + offshore_name12
109offshore_dem_name_aho13 = datadir + offshore_name13
110offshore_dem_name_aho14 = datadir + offshore_name14
111offshore_dem_name_aho15 = datadir + offshore_name15
112offshore_dem_name_aho16 = datadir + offshore_name16
113coast_dem_name = datadir + coast_name
114
115# addition once total grid delivered
116onshore_offshore_dem_name_25 = datadir + '25m_se_tas' #25m grid
117onshore_offshore_dem_name    = datadir + '50m_se_tas' #50m grid
118
119# output names
120bruny_dem_name_25   = datadir + 'bruny_25_dem'
121hobart_dem_name_25  = datadir + 'hobart_25_dem'
122combined_dem_name   = datadir + 'hobart_combined_elevation'
123combined_dem_name_2 = datadir + 'hobart_combined_elevation_2'
124combined_dem_name_3 = datadir + 'hobart_combined_elevation_jane'
125combined_dem_name_4 = datadir + 'hobart_combined_elevation_test'
126
127#outputname = outputtimedir + basename  #Used by post processing
128
129###############################
130# Domain definitions
131###############################
132
133# bounding box for clipping MOST output (much bigger than study area)
134south = degminsec2decimal_degrees(-44,45,0)
135north = degminsec2decimal_degrees(-42,0,0)
136west  = degminsec2decimal_degrees(146,45,0)
137east  = degminsec2decimal_degrees(148,25,0) 
138
139###Main Domain of Hobart:
140d0 = [south, west]
141d1 = [south, east]
142d2 = [north, east]
143d3 = [north, west]
144polyAll2, zone = convert_from_latlon_to_utm([d0, d1, d2, d3])
145refzone = zone
146
147# Mark run - morning Fri 29 Sep; surrounds -100m and 20mish elevation
148polyAll = read_polygon(polygondir+'new_extent_2.csv')
149plot_polygons([polyAll, polyAll2],'boundingpoly',verbose=False)
150#print 'Area of bounding polygon', polygon_area(polyAll)
151polyAll_refine = read_polygon(polygondir+'extent_refine.csv')
152#print 'Area of refined bounding polygon', polygon_area(polyAll_refine)
153#print 'change', polygon_area(polyAll_refine)/polygon_area(polyAll)
154
155polyAll_refine2 = read_polygon(polygondir+'export_refine_2.csv')
156#print 'Area of refined bounding polygon 2', polygon_area(polyAll_refine2)
157#print 'change 2', polygon_area(polyAll_refine2)/polygon_area(polyAll)
158
159###################################################################
160# Clipping regions for export to asc and regions for clipping data
161###################################################################
162
163# clipping 12.5m onshore data set
164eastingmin = 520000
165eastingmax = 536000
166northingmin = 5245000
167northingmax = 5260000
168# for original model setup
169# clipping 25m data set - Hobart
170eastingmin25 = 524208.387
171eastingmax25 = 554867.24
172northingmin25 = 5229154.555
173northingmax25 = 5258511.857
174hob_clip = [[eastingmin25,northingmin25],[eastingmin25,northingmax25],[eastingmax25,northingmax25],[eastingmax25,northingmin25]]
175# clipping 25m data set - Bruny
176eastingmin25_2 = 523672.502
177eastingmax25_2 = 536020.057
178northingmin25_2 = 5203120.035
179northingmax25_2 = 5212052.309
180bruny_clip = [[eastingmin25_2,northingmin25_2],[eastingmin25_2,northingmax25_2],[eastingmax25_2,northingmax25_2],[eastingmax25_2,northingmin25_2]]
181# overall 25m clipping region - covering Hobart and Bruny
182eastingmin25_3 = 523672.502
183eastingmax25_3 = 554867.24
184northingmin25_3 = 5203120.035
185northingmax25_3 = 5258511.857
186
187# export for refined grids at paleo sites
188eastingmin_site13 = 532835.94
189eastingmax_site13 = 533848.23
190northingmin_site13 = 5237384.18
191northingmax_site13 = 5238195.03
192
193eastingmin_kingston = 525500 #525905.49
194eastingmax_kingston = 527773.42
195northingmin_kingston = 5240256.71
196northingmax_kingston = 5242300 #5242148.58
197
198eastingmin_bruny = 523638.9
199eastingmax_bruny = 533000.88
200northingmin_bruny = 5204967.41
201northingmax_bruny = 5212177.69
202
203###############################
204# Interior region definitions
205###############################
206
207# Hobart digitized polygons
208poly_hobart1 = read_polygon(polygondir+'Hob_poly1.csv')
209poly_hobart2 = read_polygon(polygondir+'Hob_poly2.csv')
210poly_hobart3 = read_polygon(polygondir+'Hob_poly3.csv')
211poly_hobart4 = read_polygon(polygondir+'Hob_poly4.csv')
212poly_hobart5 = read_polygon(polygondir+'new_interior_region.csv')
213
214plot_polygons([polyAll, poly_hobart1,poly_hobart2,poly_hobart3,poly_hobart5],'orig_model',verbose=False)
215
216# to refine around paleo sites
217poly_kingston = read_polygon(polygondir+'Kingston.csv')
218poly_bruny = read_polygon(polygondir+'bruny_refine.csv')
219poly_site13 = read_polygon(polygondir+'Site13.csv')
220#print 'areas of refined polys', polygon_area(poly_kingston), polygon_area(poly_bruny), polygon_area(poly_site13)
221
222plot_polygons([polyAll, poly_kingston,poly_bruny,poly_site13,poly_hobart5],'orig_model',verbose=False)
223#print 'Area of interior regions', (polygon_area(poly_hobart1)+\
224#      polygon_area(poly_hobart2)+polygon_area(poly_hobart3)+\
225#      polygon_area(poly_hobart4)+polygon_area(poly_hobart5))/1000000.
226
227#print 'Area of original 25m data clipping', (polygon_area(hob_clip)+polygon_area(bruny_clip))/1000000.
228#print 'Area of new interior region', polygon_area(poly_hobart5)/1000000.
Note: See TracBrowser for help on using the repository browser.