source: anuga_work/production/australia_ph2/eucla_motel/project.py @ 6316

Last change on this file since 6316 was 6316, checked in by myall, 15 years ago

correcting name in project.py

File size: 7.0 KB
Line 
1"""Common filenames and locations for elevation, meshes and outputs.
2This script is the heart of all scripts in the folder
3"""
4#------------------------------------------------------------------------------
5# Import necessary modules
6#------------------------------------------------------------------------------
7
8import os
9from os.path import join
10from os import sep, getenv
11from time import localtime, strftime, gmtime
12from anuga.utilities.polygon import read_polygon, number_mesh_triangles
13from anuga.utilities.system_tools import get_user_name, get_host_name
14
15#------------------------------------------------------------------------------
16# Directory setup
17#------------------------------------------------------------------------------
18# Note: INUNDATIONHOME is the inundation directory, not the data directory.
19
20home = getenv('INUNDATIONHOME')+sep+'data'+sep # Absolute path for data folder
21muxhome = getenv('MUXHOME')
22user = get_user_name()
23host = get_host_name()
24
25# determines time for setting up output directories
26time = strftime('%Y%m%d_%H%M%S',localtime()) 
27gtime = strftime('%Y%m%d_%H%M%S',gmtime()) 
28build_time = time+'_build'
29run_time = time+'_run'
30
31# this section needs to be updated to reflect the modelled community.
32# Note, the user needs to set up the directory system accordingly
33state = 'australia_ph2'
34scenario_name = 'eucla_motel'
35
36#------------------------------------------------------------------------------
37# Initial Conditions
38#------------------------------------------------------------------------------
39# Model specific parameters. One or all can be changed each time the
40# run_scenario script is executed
41tide = 0 
42event_number = 27255 # Java 9.3 worst case for Perth
43alpha = 0.1             # smoothing parameter for mesh
44friction = 0.01           # manning's friction coefficient
45starttime = 0             
46finaltime = 1000 #80000         # final time for simulation
47
48setup = 'final'  # Final can be replaced with trial or basic.
49               # Either will result in a coarser mesh that will allow a
50               # faster, but less accurate, simulation.
51
52if setup =='trial':
53    print'trial'
54    scale_factor=100
55    time_thinning=96
56    yieldstep=240
57if setup =='basic': 
58    print'basic'
59    scale_factor=4
60    time_thinning=12
61    yieldstep=120
62if setup =='final': 
63    print'final'
64    scale_factor=1
65    time_thinning=4
66    yieldstep=60
67
68
69#------------------------------------------------------------------------------
70# Output Filename
71#------------------------------------------------------------------------------
72# Important to distinguish each run - ensure str(user) is included!
73# Note, the user is free to include as many parameters as desired
74output_comment= ('_' + setup + '_' + str(tide)+ '_' + str(event_number) +
75                 '_' + str(user))
76
77#------------------------------------------------------------------------------
78# Input Data
79#------------------------------------------------------------------------------
80# ELEVATION DATA
81# Used in build_elevation.py
82# Format for ascii grids, as produced in ArcGIS + a projection file
83##ascii_grid_filenames = ['grid250m'] # 250m grid 2005
84
85# Format for point is x,y,elevation (with header)
86point_filenames = ['grid250m_pts.txt'] # To extend boundary
87
88# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
89# Used in build_elevation.py
90# Format for points easting,northing (no header)
91bounding_polygon_filename = 'bounding_polygon.csv'
92
93# GAUGES - for creating timeseries at a specific point
94# Used in get_timeseries.py
95# Format easting,northing,name,elevation (with header)
96##gauges_filename = 'gauges.csv'
97
98# BOUNDING POLYGON
99# used in build_boundary.py and run_model.py respectively
100# NOTE: when files are put together the points must be in sequence
101# For ease go clockwise!
102# Check the run_model.py for boundary_tags
103
104# Thinned ordering file from Hazard Map (geographic)
105# Format is index,latitude,longitude (with header)
106urs_order_filename = 'urs_order.csv'
107
108# Landward bounding points
109# Format easting,northing (no header)
110landward_boundary_filename = 'landward_boundary.csv'
111
112#------------------------------------------------------------------------------
113# Output Elevation Data
114#------------------------------------------------------------------------------
115# Output filename for elevation
116# this is a combination of all the data generated in build_elevation.py
117combined_elevation_basename = scenario_name + '_combined_elevation'
118
119#------------------------------------------------------------------------------
120# Directory Structure
121#------------------------------------------------------------------------------
122anuga_folder = join(home, state, scenario_name, 'anuga')
123topographies_folder = join(anuga_folder, 'topographies')
124polygons_folder = join(anuga_folder, 'polygons')
125boundaries_folder = join(anuga_folder, 'boundaries')
126output_folder = join(anuga_folder, 'outputs')
127gauges_folder = join(anuga_folder,'gauges')
128meshes_folder = join(anuga_folder, 'meshes')
129
130#------------------------------------------------------------------------------
131# Location of input and output data
132#------------------------------------------------------------------------------
133
134# The absolute pathname of the all elevation, generated in build_elevation.py
135combined_elevation = join(topographies_folder, combined_elevation_basename)
136
137# The absolute pathname of the mesh, generated in run_model.py
138meshes = join(meshes_folder, scenario_name) + '.msh'
139
140# The absolute pathname for the urs order points, used within build_boundary.py
141urs_order = join(boundaries_folder, urs_order_filename)
142
143# The absolute pathname for the landward points of the bounding polygon,
144# Used within run_model.py)
145landward_boundary = join(boundaries_folder, landward_boundary_filename)
146
147# The absolute pathname for the .sts file, generated in build_boundary.py
148event_sts = join(boundaries_folder, str(event_number), scenario_name)
149
150# The absolute pathname of the event folder
151event_folder = join(boundaries_folder, str(event_number))
152
153# The absolute pathname for the output folder names
154# Used for build_elevation.py
155output_build = join(output_folder, build_time) + '_' + str(user) 
156# Used for run_model.py
157output_run = join(output_folder, run_time) + output_comment
158# Used by post processing
159output_run_time = join(output_run, scenario_name) 
160
161# The absolute pathname for the gauges file
162# Used for get_timeseries.py
163##gauges = join(gauges_folder, gauges_filename)       
164
165#------------------------------------------------------------------------------
166# Reading polygons and creating interior regions
167#------------------------------------------------------------------------------
168
169# Initial bounding polygon for data clipping
170bounding_polygon = read_polygon(join(polygons_folder,
171                                     bounding_polygon_filename))
172bounding_maxarea = 100000*scale_factor
173
174interior_regions = []
175
176# Estimate the number of triangles                     
177trigs_min = number_mesh_triangles(interior_regions,
178                                  bounding_polygon, bounding_maxarea)
179print 'min estimated number of triangles', trigs_min
180   
181
Note: See TracBrowser for help on using the repository browser.