source: anuga_work/production/onslow_2006/project_urs.py @ 4856

"""Common filenames and locations for topographic data, meshes and outputs.
"""

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

# file and system info
#---------------------------------
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)


# INUNDATIONHOME is the inundation directory, not the data directory.
#home += sep +'data'

#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 = -1.5
#tide = 0.0
tide = 1.5

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

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=1
    yieldstep=60
if setup =='store':
    print'store'
    res_factor=10
    time_thinning=48
    yieldstep=240    

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

#Making assumptions about the location of scenario data
state = 'western_australia'
scenario_name = 'onslow'
scenario = 'onslow_tsunami_scenario_2006'

# onshore data provided by WA DLI
onshore_name = 'onslow_onshore_30m_dted2' # original
onshore_name1 = 'onslow_onshore_20m_dli' # original

# AHO + DPI data + colin French coastline
coast_name = 'onslow_coast'
islands_name = 'onslow_islands_dted2'
offshore_name = 'onslow_offshore_points'
offshore1_name = '1'
offshore2_name = '2'

#final topo name
combined_name ='onslow_combined_elevation'
combined_smaller_name = 'onslow_combined_elevation_small'
combined_smallest_name = 'onslow_combined_elevation_smallest'

anuga_dir = home+state+sep+scenario+sep+'anuga'+sep
#anuga_dir_name = temp_dir + temp_time+sep

topographies_in_dir = home+state+sep+scenario+sep+'elevation_final'+sep+'070112'+sep+'points'+sep
topographies_dir = anuga_dir+'topographies'+sep
topographies_time_dir = topographies_dir+build_time+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_smaller_dir_name = topographies_dir + combined_smaller_name
#combined_smallest_dir_name = topographies_dir + combined_smallest_name
#combined_time_dir_final_name = topographies_time_dir + combined_final_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


if source =='dampier':
    boundaries_name = 'onslow_3854_17042007' #Dampier gun
    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'dampier'+sep+'1_10000'+sep

if source=='onslow':
    boundaries_name = 'onslow_3859_16052007' #onslow_hedland_broome gun
    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'onslow_hedland_broome'+sep+'1_10000'+sep
    
if source=='exmouth':
    boundaries_name = 'onslow_3103_18052007' #exmouth gun
    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep

#boundaries locations
boundaries_source = ''
#boundaries_in_dir = anuga_dir+'boundaries'+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+'_'+source+sep
#output_run_time_dir = output_dir +run_time+dir_comment+sep
output_run_time_dir = output_dir +'20070619_042140_run_final_0_onslow_nbartzis'+sep
output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing

#gauges
gauges_dir = anuga_dir+'gauges'+sep

gauges_dir_name = gauges_dir + 'gauge_location_onslow.csv'
gauges_dir_name_bindi = gauges_dir + 'gauge_location_bindi.csv'
gauges_comparison_dir_name_50 = gauges_dir + '50_gauges.xya'
gauges_comparison_dir_name = gauges_dir + 'MOST_comparison_gauges.xya'
gauges_comparison_dir_name_3d = gauges_dir + 'MOST_comparison_gauges_3d.xya'

community_filename = gauges_dir + 'CHINS_v2.csv'
community_scenario = gauges_dir + 'community_onslow.csv'
buildings_filename = gauges_dir + 'Onslow_res_Project.csv'
buildings_filename_out =output_run_time_dir + 'Onslow_res_Project_modified.csv'
###############################
# Domain definitions
###############################

south_boundary = degminsec2decimal_degrees(-22,15,0)
north_boundary = degminsec2decimal_degrees(-20,30,0)
west_boundary = degminsec2decimal_degrees(114,15,0)
east_boundary = degminsec2decimal_degrees(115,50,0)

#Updated Main Domain of Onslow: first run NB 6/4/06
d0 = [310000, 7690000]
d1 = [280000, 7690000]
d2 = [270000, 7645000]
d3 = [240000, 7625000]
d4 = [270000, 7580000]
d5 = [300000, 7590000]
d6 = [340000, 7610000]

poly_all = [d0, d1, d2, d3, d4, d5, d6]
res_poly_all = 100000*res_factor
print 'bounding polygon area', polygon_area(poly_all)/1000000.0
#Interior region - Onslow town

#for visualisation
#i0 = [304000, 7609000]
#i1 = [301500, 7605000]
#i2 = [304000, 7603000]
#i3 = [307000, 7602000]
#i4 = [309000, 7603000]
#i5 = [311000, 7607000]

i0 = [304000, 7607000]
i1 = [302000, 7605000]
i2 = [304000, 7603000]
i3 = [307000, 7602000]
i4 = [309000, 7603000]
i5 = [307000, 7606000]
poly_onslow = [i0, i1, i2, i3, i4, i5]
res_onslow = 500*res_factor
print 'onslow polygon area', polygon_area(poly_onslow)/1000000.0

#Thevenard Island
j0 = [294000, 7629000]
j1 = [285000, 7625000]
j2 = [294000, 7621000]
j3 = [299000, 7625000]
poly_thevenard = [j0, j1, j2, j3]
res_thevenard = 2500*res_factor
print 'thevenard polygon area', polygon_area(poly_thevenard)/1000000.0

#med res around onslow
l0 = [300000, 7610000]
l1 = [285000, 7600000]
l2 = [300000, 7597500]
l3 = [310000, 7600000]
l4 = [315000, 7610000]
poly_coast = [l0, l1, l2, l3, l4]
res_coast = 2500*res_factor
print 'coast polygon area', polygon_area(poly_coast)/1000000.0

#general coast and local area to onslow region
m0 = [270000, 7581000]
m1 = [300000, 7591000]
m2 = [339000, 7610000]
m3 = [330000, 7630000]
m4 = [290000, 7640000]
m5 = [260000, 7600000]
poly_region = [m0, m1, m2, m3, m4, m5]
res_region = 20000*res_factor
print 'region polygon area', polygon_area(poly_region)/1000000.0

interior_regions = [[poly_onslow,res_onslow],[poly_coast,res_coast],
                    [poly_region,res_region]
#                   ,[poly_thevenard,res_thevenard]
                    ]

trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)

print 'min number triangles', trigs_min


###################################################################
# Clipping regions for export to asc and regions for clipping data
###################################################################

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

# region to export
eastingmin = 300000
eastingmax = 310000
northingmin = 7600000
northingmax = 7610000