source: anuga_core/documentation/user_manual/examples/project.py @ 3944

Last change on this file since 3944 was 3944, checked in by sexton, 16 years ago

(i) update sww2timeseries so can handle gauges which don't fall within sww domain (ii) script for generating report for dampier based on sww files from parallel setup

File size: 2.3 KB
Line 
1"""Common filenames and locations for topographic data, meshes and outputs.
2Also includes origin for slump scenario.
3"""
4
5from os import sep, environ
6from os.path import expanduser
7import sys
8
9# We'll do this for the final thing
10home = '.'
11
12if sys.platform == 'win32':
13    home = environ['INUNDATIONHOME']     #Sandpit's parent dir
14else:   
15    home = expanduser('~')
16   
17basename = 'Sydney'
18
19#Setting the location of scenario data
20scenario_dir_name = 'sydney_tsunami_example_2006'
21# revised 100m data
22#coarsename = 'bathyland100'
23# revised 25m data
24#finename = 'bathy_dem25'
25# revised 25m data
26dem_name = 'sydney_MGA56'
27
28home += sep + 'data' + sep + 'new_south_wales'
29
30#Derive subdirectories and filenames
31meshdir = home+sep+scenario_dir_name+sep+'meshes'+sep
32datadir = home+sep+scenario_dir_name+sep+'topographies'+sep
33outputdir = home+sep+scenario_dir_name+sep+'output'+sep
34polygondir = home+sep+scenario_dir_name+sep+'polygons'+sep
35
36meshname = meshdir + basename
37#coarsedemname = datadir + coarsename
38#finedemname = datadir + finename
39demname = datadir + dem_name
40combineddemname = datadir + 'sydneytopo'
41
42refzone = 56
43print 'refzone', refzone
44
45# demo poly
46j0 = [385000, 6280000]
47j1 = [360000, 6273000]
48j2 = [335000, 6273000]
49j3 = [330000, 6265000]
50j31 = [325000, 6260000]
51j4 = [316000, 6260000]
52j5 = [316000, 6246750]
53j6 = [350000, 6246750]
54j7 = [385000, 6238000]
55
56demopoly = [j0, j1, j2, j3, j31, j4, j5, j6, j7]
57
58from anuga.utilities.polygon import read_polygon, plot_polygons
59#polygonptsfile4 = polygondir + 'poly1'
60#polygonptsfile0 = polygondir + 'poly2'
61#polygonptsfile1 = polygondir + 'poly3'
62#polygonptsfile2 = polygondir + 'poly4'
63#polygonptsfile3 = polygondir + 'poly5'
64#northern_polygon = read_polygon(polygonptsfile0 + '.csv')
65#manly_polygon = read_polygon(polygonptsfile1 + '.csv')
66#harbour_polygon = read_polygon(polygonptsfile2 + '.csv')
67#southern_polygon = read_polygon(polygonptsfile3 + '.csv')
68#top_polygon = read_polygon(polygonptsfile4 + '.csv')
69coastal_polygon = read_polygon(polygondir+'coastal'+'.csv')
70shallow_polygon = read_polygon(polygondir+'shallow'+'.csv')
71
72#plot_polygons([demopoly,northern_polygon,manly_polygon,harbour_polygon,southern_polygon,top_polygon],'model_setup',verbose=False)
73plot_polygons([demopoly,coastal_polygon,shallow_polygon],'new_model_setup',verbose=False)
74
75slump_origin = [372500.0, 6255000.0] #Absolute UTM
Note: See TracBrowser for help on using the repository browser.