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

# -*- coding: cp1252 -*-
"""Common filenames and locations for topographic data, meshes and outputs.
"""

from os import sep, environ, getenv, getcwd
from os.path import expanduser
import sys
from time import localtime, strftime, gmtime
                
#Making assumptions about the location of scenario data
state = 'tasmania'
scenario_dir_name = 'hobart_tsunami_scenario_2006'

# data provided by Tas SES and checked by NM&I
#onshore_name = 'hob3_topo' # original
onshore_name = 'hob5_topo' # 12.5m grid and clipped to 100m elevation or 3000m from coast
onshore_name_25 = 'hob5_topo_25' # 25m grid and clipped to 100m elevation or 3000m from coast
offshore_name_tas1 = 'derwent_2m'
offshore_name_tas2 = 'derwent_5m'
offshore_name_tas3 = 'south_east_tas' #actually this is AHO
offshore_name_tas4 = 'hobart_1m'

# AHO data and checked by NM&I
offshore_name1 = 'xy100003760'
offshore_name2 = 'xy100003761'
offshore_name3 = 'xy100003762'
offshore_name4 = 'xy100003907'
offshore_name5 = 'xy100003908'
offshore_name6 = 'xy100003909'
offshore_name7 = 'xy100003910'
offshore_name8 = 'xy100003932'
offshore_name9 = 'xy100003933'
offshore_name10 = 'xy100003934'
offshore_name11 = 'xy100003935'
offshore_name12 = 'xy100003936'
offshore_name13 = 'xy100003964'
offshore_name14 = 'xy100014250'
offshore_name15 = 'xy100014253'
offshore_name16 = 'xy100016142'

# developed by NM&I
coast_name = 'coastline_points'

boundary_basename = 'SU-AU_clip'

#swollen/ all data output
basename = 'source'

codename = 'project.py'

if sys.platform == 'win32':
    home = getenv('INUNDATIONHOME')
    user = getenv('USERPROFILE')    
else:    
    home = getenv('INUNDATIONHOME', sep+'d'+sep+'xrd'+sep+'gem'+sep+'2'+sep+'ramp'+'risk_assessment_methods_project'+sep+'inundation'+sep+'data')     
    user = getenv('LOGNAME')
    print 'USER:', user

#Derive subdirectories and filenames
#time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
local_time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir

meshdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'meshes'+sep
datadir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'topographies'+sep
gaugedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'gauges'+sep
polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
boundarydir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep
#output dir without time
outputdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'outputs'+sep
outputtimedir = outputdir + local_time + sep
polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep

gauge_filename = gaugedir + 'hobart_gauges.xya'
buildings_filename = gaugedir + 'hobart_res.csv'
buildings_filename_damage_out = 'hobart_res_modified.csv'

gaugetimeseries = gaugedir + 'hobart'

# boundary source data
#MOST_dir = 'f:'+sep+'3'+sep+'ehn'+sep+'users'+sep+'davidb'+sep+'tsunami'+sep+'WA_project'+sep+'SU-AU_90'+sep+'most_2'+sep+'detailed'+sep

codedir = getcwd()+sep
                                
codedirname = codedir + 'project.py'

meshname = outputtimedir + 'mesh_' + basename

onshore_dem_name = datadir + onshore_name
onshore_dem_name_25 = datadir + onshore_name_25
all_onshore_dem_name = datadir + 'combined_onshore'
offshore_dem_name_local1 = datadir + offshore_name_tas1
offshore_dem_name_local2 = datadir + offshore_name_tas2
offshore_dem_name_local3 = datadir + offshore_name_tas3
offshore_dem_name_local4 = datadir + offshore_name_tas4
offshore_dem_name_aho1 = datadir + offshore_name1
offshore_dem_name_aho2 = datadir + offshore_name2
offshore_dem_name_aho3 = datadir + offshore_name3
offshore_dem_name_aho4 = datadir + offshore_name4
offshore_dem_name_aho5 = datadir + offshore_name5
offshore_dem_name_aho6 = datadir + offshore_name6
offshore_dem_name_aho7 = datadir + offshore_name7
offshore_dem_name_aho8 = datadir + offshore_name8
offshore_dem_name_aho9 = datadir + offshore_name9
offshore_dem_name_aho10 = datadir + offshore_name10
offshore_dem_name_aho11 = datadir + offshore_name11
offshore_dem_name_aho12 = datadir + offshore_name12
offshore_dem_name_aho13 = datadir + offshore_name13
offshore_dem_name_aho14 = datadir + offshore_name14
offshore_dem_name_aho15 = datadir + offshore_name15
offshore_dem_name_aho16 = datadir + offshore_name16
coast_dem_name = datadir + coast_name
combined_dem_name = datadir + 'hobart_combined_elevation'

outputname = outputtimedir + basename  #Used by post processing

# clipping 12.5m onshore data set
eastingmin = 520000
eastingmax = 536000
northingmin = 5245000
northingmax = 5260000

# region to export for inundation map
e_min_area = 500000
e_max_area = 580000
n_min_area = 5200000
n_max_area = 5275000

from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
                        
# bounding box
south = degminsec2decimal_degrees(-43,45,0)
north = degminsec2decimal_degrees(-42,30,0)
west = degminsec2decimal_degrees(146,45,0)
east = degminsec2decimal_degrees(148,15,0)

#Main Domain of Hobart: first run JS 18/9/06
d0 = [south, west]
d1 = [south, east]
d2 = [north, east]
d3 = [north, west]

polyAll, zone = convert_points_from_latlon_to_utm([d0, d1, d2, d3])
refzone = zone
NWcorner = polyAll[0]
SEcorner = polyAll[2]

#Main Domain of Hobart: updated region from Alex JS 22/9/06
j0 = [583819.54,5193158.37]
j1 = [585025.05,5172148.01]
j2 = [583991.76,5155443.05]
j3 = [486000.80,5155443.05]
j4 = [486173.02,5165603.80]
j5 = [478939.94,5170770.28]
j6 = [479628.80,5290787.67]
j7 = [570731.12,5281488.00]
j8 = [583991.76,5266160.77]
j9 = [590363.75,5268399.58]
j10 = [597252.40,5276838.17]
#j11 = [602418.88,5276838.17]
j11 = [600000,5276838.17]
#j12 = [600868.93,5201924.17]
j12 = [597600, 5201914]
#j13 = [593119.21,5196241.04]
j13 = [590000,5196241]
#j14 = [583819.54,5193158.37]
j14 = [580000, 5193160]ZZ
polyAll = [j0,j1,j2,j3,j4,j5,j6,j7,j8,j9,j10,j11,j12,j13,j14]

# region to export for Alex to make bathymetry map
e_min_area = 480000#NWcorner[0]
e_max_area = 600000#SEcorner[0]
n_min_area = 5160000#NWcorner[1]
n_max_area = 5280000#SEcorner[1]

#Interior region - Hobart city area + Glenorchy, Kingston
i0 = [517000, 5267000]
i1 = [517000, 5255000]
i2 = [520000, 5250000]
i3 = [522000, 5239000]
i4 = [524000, 5238000]
i5 = [526000, 5236000]
i6 = [530000, 5244000]
i7 = [530000, 5250000]
i8 = [534000, 5254000]
i9 = [520000, 5270000]

poly_hobart = [i0, i1, i2, i3, i4, i5, i6, i7, i8, i9]

# Tasman Peninsula
l0 = [550000, 5247000]
l1 = [550000, 5211000]
l2 = [583000, 5211000]
l3 = [583000, 5247000]

poly_tasman_peninsula = [l0, l1, l2, l3]

# Bruny Island
from anuga.utilities.polygon import read_polygon
poly_bruny = read_polygon(polygondir+'bruny.csv')