source: anuga_work/production/pt_hedland_2006/project_urs.py @ 4467

"""Common filenames and locations for topographic data, meshes and outputs.
Also includes origin for slump scenario.
"""

from os import sep, environ, getenv, getcwd,umask
from os.path import expanduser, basename
from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles
import sys
from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
from time import localtime, strftime, gmtime
from anuga.utilities.system_tools import get_user_name, get_host_name

codename = 'project.py'

home = getenv('INUNDATIONHOME') + sep +'data'+sep #Sandpit's parent dir   
user = get_user_name()
host = get_host_name()
#needed when running using mpirun, mpirun doesn't inherit umask from .bashrc
umask(002)

#Making assumptions about the location of scenario data
state = 'western_australia'
scenario_name = 'pt_hedland'
scenario = 'pt_hedland_tsunami_scenario_2006'

#time stuff
time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
build_time = time+'_build'
run_time = time+'_run'

#tide = -3.9
tide = 0.0
#tide = 3.6

#Maybe will try to make project a class to allow these parameters to be passed in.
alpha = 0.1
friction=0.01
finaltime=25000
starttime=3600
setup='final'

#dir_comment='_trial_'+str(tide)+'_'+str(user)
#dir_comment='_final_'+str(tide)+'_'+str(user)
#time_thinning=4
#yieldstep=60
#res_factor = 1

if setup =='trial':
    print'trial'
    res_factor=10
    time_thinning=48
    yieldstep=240
if setup =='basic': 
    print'basic'
    res_factor=4
    time_thinning=12
    yieldstep=120
if setup =='final': 
    print'final'
    res_factor=1
    time_thinning=4
    yieldstep=60


dir_comment='_'+setup+'_'+str(tide)+'_'+str(user)

# onshore data from 30m DTED level 2
onshore_name = 'dli_dem' # get from Neil/Ingo (DEM or topo data)
offshore_name = 'pt_hedland_bathy_edit'
offshore1_name = '1'
offshore2_name = '2'
coast_name = 'Coastline'

#final topo name
combined_name ='pt_hedland_combined_elevation'
combined_small_name = 'pt_hedland_combined_elevation_small'
combined_smallest_name = 'pt_hedland_export'

anuga_dir = home+state+sep+scenario+sep+'anuga'+sep

topographies_in_dir = home+state+sep+scenario+sep+'elevation_final'+sep+'20070416'+sep+'points'+sep
topographies_dir = anuga_dir+'topographies'+sep

# input topo file location
onshore_in_dir_name = topographies_in_dir + onshore_name
coast_in_dir_name = topographies_in_dir + coast_name
offshore_in_dir_name = topographies_in_dir + offshore_name
offshore1_in_dir_name = topographies_in_dir + offshore1_name
offshore2_in_dir_name = topographies_in_dir + offshore2_name

onshore_dir_name = topographies_dir + onshore_name
coast_dir_name = topographies_dir + coast_name
offshore_dir_name = topographies_dir + offshore_name
offshore1_dir_name = topographies_dir + offshore1_name
offshore2_dir_name = topographies_dir + offshore2_name

#final topo files
combined_dir_name = topographies_dir + combined_name
#combined_time_dir_name = topographies_time_dir + combined_name
combined_small_dir_name = topographies_dir + combined_small_name
combined_smallest_dir_name = topographies_dir + combined_smallest_name

meshes_dir = anuga_dir+'meshes'+sep
meshes_dir_name = meshes_dir + scenario_name

polygons_dir = anuga_dir+'polygons'+sep
tide_dir = anuga_dir+'tide_data'+sep

#boundaries locations
boundaries_name = 'pt_hedland_3854_17042007'

boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'1_10000'+sep
boundaries_in_dir_name = boundaries_in_dir + boundaries_name
boundaries_dir = anuga_dir+'boundaries'+sep
boundaries_dir_name = boundaries_dir + boundaries_name

#output locations
output_dir = anuga_dir+'outputs'+sep
output_build_time_dir = output_dir+build_time+sep
output_run_time_dir = output_dir +run_time+dir_comment+sep
output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing

#gauges
gauge_name = 'gauge_location_port_hedland.csv' 
gauge_name_test = 'gauge_checking_test.csv' 
gauges_dir = anuga_dir+'gauges'+sep
gauges_dir_name = gauges_dir + gauge_name
gauges_dir_name_test = gauges_dir + gauge_name_test

tide_dir = anuga_dir+'tide_data'+sep

buildings_filename = gauges_dir + 'pt_hedland_res.csv'
buildings_filename_out = 'pt_hedland_res_modified.csv'
buildings_filename_damage_out = 'pt_hedland_res_modified_damage.csv'
tidal_filename = tide_dir + 'pt_hedland_tide.txt'
community_filename = gauges_dir + 'CHINS_v2.csv'
community_scenario = gauges_dir + 'community_pt_hedland.csv'

# clipping region to make DEM (pts file) from onshore data
eastingmin = 594000
eastingmax = 715000
northingmin = 7720000
northingmax = 7880000

# for ferret2sww
south = degminsec2decimal_degrees(-20,30,0)
north = degminsec2decimal_degrees(-17,10,0)
west = degminsec2decimal_degrees(117,00,0)
east = degminsec2decimal_degrees(120,00,0)

# region to export (used from export_results.py)
e_min_area = 648000
e_max_area = 675000
n_min_area = 7745000
n_max_area = 7761000

export_region = [[e_min_area,n_min_area],[e_min_area,n_max_area],
                 [e_max_area,n_max_area],[e_max_area,n_min_area]]
                 
refzone = 50 

from anuga.coordinate_transforms.redfearn import redfearn
# boundary up to 50 m contour
lat1_50 = degminsec2decimal_degrees(-19,20,0)
lat2_50 = degminsec2decimal_degrees(-19,30,0)
lat3_50 = degminsec2decimal_degrees(-19,45,0)
lon1_50 = degminsec2decimal_degrees(119,05,0)
lon2_50 = degminsec2decimal_degrees(118,20,0)
lon3_50 = degminsec2decimal_degrees(117,45,0)
z, easting, northing = redfearn(lat1_50, lon1_50)
d0_50 = [easting, northing]
z, easting, northing = redfearn(lat2_50, lon2_50)
d1_50 = [easting, northing]
z, easting, northing= redfearn(lat3_50, lon3_50)
d2_50 = [easting, northing]

d4_50 = [285000, 7585000]
d6_50 = [330000, 7605000]
#bounding_poly50 = [p0_50, p1_50, p2_50, d6_50, d5, d4_50]

d0 = [763852.0, 7934358.0]
d1 = [710987.0, 7925797.0]
d2 = [658264.0, 7926314.0]
d3 = [552686.0, 7871580.0]
#d4 = [604415.81, 7733013.56]
d4 = [638000.0, 7733013.56]
#d5 = [656561.15, 7732615.11]
d5 = [662000.0, 7732615.11]
#d6 = [708940.32, 7750510.33]
d6 = [690000.0, 7740510.33]
#polyAll = [d0, d1, d2, d3, d4, d5, d6]
#polyAll = [d0_50, d1_50, d2_50, d4, d5, d6]
# from Hamish
h0=[629262.17, 7747205.47]
h1=[552686.00, 7871579.99] #d3
h2=[658264.00, 7926314.00] #d2
h3=[710986.99, 7925796.99] #d1
h4=[763851.99, 7934357.99] #d0
h5=[701485.21, 7770656.86]
h6=[698273.75, 7762227.38]
h7=[698194.23, 7762018.65]
h8=[691627.41, 7744781.98]
h9=[679220.75, 7743604.59]
h10=[653512.59, 7740528.56]
h11=[634777.71, 7738247.17]
h12=[629443.86, 7746910.37]
h13=[629396.84, 7746986.75]
h14=[629352.32, 7747059.06]
h15=[629276.24, 7747182.63]
h16=[629262.17, 7747205.47] #repeat of h0
# using Hamish's new bounding polygon
#polyAll = [d0_50, d1_50, d2_50, h16,h15,h14,h13,h12,h11,h10,h9,h8,h7,h6,h5]
poly_all = [d0_50, d1_50, d2_50, h16,h11,h8,h6, h5]
res_poly_all = 250000*res_factor

#Interior region - Pt Hedland town
i0 = [668000, 7757000]
i1 = [659000, 7755000]
i2 = [660000, 7749000]
i3 = [667000, 7746000]
i4 = [678000, 7751000]

poly_pt_hedland = [i0, i1, i2, i3, i4]
res_pt_hedland=500*res_factor

#Are there other significant features?
j0 = [670000, 7760000]
j1 = [633000, 7745000]
j2 = [665000, 7743000]
j3 = [690000, 7755000]

poly_region = [j0, j1, j2, j3]
res_region = 60000*res_factor

poly_mainland = read_polygon(polygons_dir+'Initial_Condition.csv')

poly_topo_clip = [[664000,7752000],[664000,7754000],
                  [666000,7754000],[666000,7752000]]

interior_regions = [[poly_pt_hedland,res_pt_hedland],[poly_region,res_region]
#                    ,[poly_region,res_region]
                    ]

trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)

print 'min number triangles', trigs_min