source: anuga_work/production/geraldton/run_boundary_points.py @ 5498

Last change on this file since 5498 was 5382, checked in by kristy, 16 years ago

urs boundary points

File size: 2.2 KB
Line 
1"""Script for running tsunami inundation scenario for Dampier, WA, Australia.
2
3Source data such as elevation and boundary data is assumed to be available in
4directories specified by project.py
5The output sww file is stored in project.output_run_time_dir
6
7The scenario is defined by a triangular mesh created from project.polygon,
8the elevation data and a simulated tsunami generated with URS code.
9
10Ole Nielsen and Duncan Gray, GA - 2005 and Jane Sexton, Nick Bartzis, GA - 2006
11"""
12
13#------------------------------------------------------------------------------
14# Import necessary modules
15#------------------------------------------------------------------------------
16
17# Standard modules
18
19# # Related major packages
20# from anuga.shallow_water import Domain
21from anuga.shallow_water.data_manager import URS_points_needed_to_file, \
22     urs_ungridded2sww
23import project                 # Definition of file names and polygons
24
25#-----------------------------------------------------------------------------
26# run_boundary_points
27#-----------------------------------------------------------------------------
28 
29LL_LAT = -50.0          #ll_lat - lower left latitude, in decimal degrees
30LL_LONG = 80.0          #ll-long - lower left longitude, in decimal degrees
31GRID_SPACING = 2.0/60.0 #grid_spacing - in deciamal degrees
32LAT_AMOUNT = 4800       #lat_amount - number of latitudes
33LONG_AMOUNT = 3600      #long_amount- number of longs
34ZONE = 49
35
36#------------------------------------------------------------------------------
37# Copy scripts to time stamped output directory and capture screen
38# output to file
39#------------------------------------------------------------------------------
40poly = project.poly_all
41print "poly", poly
42print "project.boundaries_in_dir_name",project.boundaries_in_dir_name
43URS_points_needed_to_file(project.boundaries_in_dir_name,poly, ZONE,
44                          LL_LAT, LL_LONG,GRID_SPACING,LAT_AMOUNT, LONG_AMOUNT,export_csv=True)
45#URS_points_needed_to_file(project.boundaries_dir_name5,poly, export_csv=True, zone=50)
46#URS_points_needed_to_file(project.boundaries_dir_name5,poly, export_csv=True)
47#urs_ungridded2sww(project.boundaries_in_dir_name2, project.boundaries_dir_name5, mean_stage= project.tide,  verbose=True)
Note: See TracBrowser for help on using the repository browser.