source: anuga_work/production/Broome_2008/project.py @ 6874

"""Common filenames and locations for elevation, meshes and outputs.
This script is the heart of all scripts in folder
"""
#------------------------------------------------------------------------------
# Import necessary modules
#------------------------------------------------------------------------------

from os import sep, environ, getenv, getcwd
from os.path import expanduser
import sys
from time import localtime, strftime, gmtime
from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, number_mesh_triangles
from anuga.utilities.system_tools import get_user_name, get_host_name
from anuga.shallow_water.data_manager import urs2sts,create_sts_boundary
from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon

#------------------------------------------------------------------------------
# Directory setup
#------------------------------------------------------------------------------
##Note: INUNDATIONHOME is the inundation directory, not the data directory.

home = getenv('INUNDATIONHOME') + sep +'data'+sep #Sandpit's parent diruser = get_user_name()
muxhome = getenv('MUXHOME')
user = get_user_name()
host = get_host_name()

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

#------------------------------------------------------------------------------
# Initial Conditions
#------------------------------------------------------------------------------

##Changed to reflect the modeled community (also sets up the directory system)
state = 'western_australia'
scenario_name = 'broome'
scenario = 'broome_tsunami_scenario_2008'
##One or all can be changed each time the run_scenario script is excuted
tide = 0.6
event_number = 27255
#event_number = 27283
alpha = 0.1
friction=0.01
starttime=0
finaltime=80000 

setup='final'  ## Can replace with trial or basic, this action will thin the mesh
               ## so that the run script will take less time (hence less acurate) 

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

#------------------------------------------------------------------------------
# Output Filename
#------------------------------------------------------------------------------
##Important to distiquish each run
dir_comment='_'+setup+'_'+str(tide)+'_'+str(event_number)+'_'+ 'alpha' +str(alpha)+'_'+str(user)

#------------------------------------------------------------------------------
# INPUT DATA
#------------------------------------------------------------------------------

##ELEVATION DATA## - used in build_broome.py
## onshore data: format ascii grid with accompaning projection file
onshore_name = 'town_topo_10m' 
## island: format ascii grid with accompaning projection file
#island_name = 'rott_dli_ext' 
#island_name1 = 'gard_dli_ext'
#island_name2 = 'carnac_island_dli_ext'
#island_name3 = 'penguin_dli_ext'
## coastline: format x,y,elevation (with title)
coast_name = 'Broome_coastline.txt'
## bathymetry: format x,y,elevation (with title)
offshore_name = 'Broome_Bathymetry.txt'
offshore_name1 = 'cable_250m'
offshore_name2 = 'inferred_north'
offshore_name3 = 'inferred_south'
offshore_name4 = 'north_250m'
offshore_name5 = 'other_topo_250m'
offshore_name6 = 'south_250m'

##GAUGES## - used in get_timeseries.py
#gauge_name = 'Broome.csv'
#gauge_name2 = 'thinned_MGA50.csv'

##BOUNDING POLYGON## -used in build_boundary.py and run_broome.py respectively
#NOTE: when files are put together must be in sequence - for ease go clockwise!
#Check the run_broome.py for boundary_tags
##thinned ordering file from Hazard Map: format index,latitude,longitude (with title)
order_filename = 'thinned_boundary_ordering.csv'
##landward bounding points
landward = 'landward_bounding_polygon.csv'

#------------------------------------------------------------------------------
# OUTPUT ELEVATION DATA
#------------------------------------------------------------------------------
##Output filename for elevation
## its a combination of all the data put together (utilisied in build_boundary)
combined_name ='broome_combined_elevation'
combined_smaller_name = 'broome_combined_elevation_smaller'

#------------------------------------------------------------------------------
# Directory Structure
#------------------------------------------------------------------------------
anuga_dir = home+state+sep+scenario+sep+'anuga'+sep
topographies_in_dir = home+state+sep+scenario+sep+'elevation_final'+sep+'points'+sep
topographies_dir = anuga_dir+'topographies'+sep
polygons_dir = anuga_dir+'polygons'+sep
tide_dir = anuga_dir+'tide_data'+sep
boundaries_dir = anuga_dir+'boundaries'+ sep 
output_dir = anuga_dir+'outputs'+sep
gauges_dir = anuga_dir+'gauges'+sep
meshes_dir = anuga_dir+'meshes'+sep

#------------------------------------------------------------------------------
# Location of input and output Data
#------------------------------------------------------------------------------
##where the input data sits
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
offshore_in_dir_name1 = topographies_in_dir + offshore_name1
offshore_in_dir_name2 = topographies_in_dir + offshore_name2
offshore_in_dir_name3 = topographies_in_dir + offshore_name3
offshore_in_dir_name4 = topographies_in_dir + offshore_name4
offshore_in_dir_name5 = topographies_in_dir + offshore_name5
offshore_in_dir_name6 = topographies_in_dir + offshore_name6

##where the output data sits
onshore_dir_name = topographies_dir + onshore_name
coast_dir_name = topographies_dir + coast_name
offshore_dir_name = topographies_dir + offshore_name
offshore_dir_name1 = topographies_dir + offshore_name1
offshore_dir_name2 = topographies_dir + offshore_name2
offshore_dir_name3 = topographies_dir + offshore_name3
offshore_dir_name4 = topographies_dir + offshore_name4
offshore_dir_name5 = topographies_dir + offshore_name5
offshore_dir_name6 = topographies_dir + offshore_name6

##where the combined file sits
combined_dir_name = topographies_dir + combined_name
combined_smaller_name_dir = topographies_dir + combined_smaller_name

##where the mesh sits (this is created during the run_broome.py)
meshes_dir_name = meshes_dir + scenario_name+'.msh'

#where the boundary order files sit (this is used within build_boundary.py)
order_filename_dir = boundaries_dir + order_filename

#where the landward points of boundary extent sit (this is used within run_broome.py)
landward_dir = boundaries_dir + landward

#where the event sts files sits (this is created during the build_boundary.py)
boundaries_dir_event = boundaries_dir + str(event_number) + sep
boundaries_dir_mux = muxhome

#where the directory of the output filename sits
output_build_time_dir = output_dir+build_time+dir_comment+sep #used for build_broome.py
output_run_time_dir = output_dir+run_time+dir_comment+sep #used for run_broome.py
output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing

#where the directory of the gauges sit
#gauges_dir_name = gauges_dir + gauge_name #used for get_timeseries.py
#gauges_dir_name2 = gauges_dir + gauge_name2 #used for get_timeseries.py

#------------------------------------------------------------------------------
# Interior region definitions
#------------------------------------------------------------------------------

#Land, to set the stage/water to be offcoast only
poly_mainland = read_polygon(polygons_dir+'initial_conditions_mainland.csv')

#Initial bounding polygon for data clipping 
poly_all = read_polygon(polygons_dir+'poly_all.csv')
res_poly_all = 100000*res_factor

#Area of Interest 1 (Broome)
poly_aoi1 = read_polygon(polygons_dir+'Area_Of_Interest.csv')
res_aoi1 = 500*res_factor

#Area of Significance 1 (Broome)
poly_aos1 = read_polygon(polygons_dir+'Area_Of_Significance.csv')
res_aos1 = 1000*res_factor

#Shallow water 1
poly_sw1 = read_polygon(polygons_dir+'Shallow_Water.csv')
res_sw1 = 25000*res_factor

##Deep water (land of Rottnest, ANUGA does not do donuts!)
#poly_dw1 = read_polygon(polygons_dir+'rottnest_internal.csv')
#res_dw1 = 100000*res_factor

# Combined all regions, must check that all are included!
interior_regions = [[poly_aoi1,res_aoi1],
                    [poly_aos1,res_aos1],
                    [poly_sw1,res_sw1]]

    
trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
print 'min estimated number of triangles', trigs_min
   

###------------------------------------------------------------------------------
### Clipping regions for export to asc and regions for clipping data
###------------------------------------------------------------------------------
##
###Geordie Bay extract ascii grid
##xminGeordie = 358000
##xmaxGeordie = 362000
##yminGeordie = 6458500
##ymaxGeordie = 6461000
##
###Sorrento extract ascii grid
##xminSorrento = 379000
##xmaxSorrento = 382500
##yminSorrento = 6477000
##ymaxSorrento = 6480000
##
###Fremantle extract ascii grid
##xminFremantle = 376000
##xmaxFremantle = 388000
##yminFremantle = 6449000
##ymaxFremantle = 6461000
##
###Rockingham extract ascii grid
##xminRockingham = 373500
##xmaxRockingham = 385500
##yminRockingham = 6424000
##ymaxRockingham = 6433000
##
##