source: anuga_work/production/busselton/build_boundary.py @ 5756

Last change on this file since 5756 was 5756, checked in by kristy, 16 years ago
File size: 3.4 KB
Line 
1from anuga.utilities.numerical_tools import ensure_numeric
2from Scientific.IO.NetCDF import NetCDFFile
3from Numeric import asarray,transpose,sqrt,argmax,argmin,arange,Float,\
4    compress,zeros,fabs,allclose,ones
5from anuga.utilities.polygon import inside_polygon,read_polygon
6from os import sep
7from time import localtime, strftime, gmtime
8from anuga.shallow_water.data_manager import urs2sts,create_sts_boundary
9import os
10import project
11
12######################################
13# Create urs boundary from mux2files #
14######################################
15
16dir=os.path.join(project.muxhome,'mux')
17
18# taken from David's event list
19urs_filenames = [
20    os.path.join(dir,'Java-0000-z.grd'),
21    os.path.join(dir,'Java-0001-z.grd'),
22    os.path.join(dir,'Java-0002-z.grd'),
23    os.path.join(dir,'Java-0003-z.grd'),
24    os.path.join(dir,'Java-0005-z.grd'), 
25    os.path.join(dir,'Java-0006-z.grd'), 
26    os.path.join(dir,'Java-0007-z.grd'), 
27    os.path.join(dir,'Java-0008-z.grd'), 
28    os.path.join(dir,'Java-0010-z.grd'), 
29    os.path.join(dir,'Java-0011-z.grd'), 
30    os.path.join(dir,'Java-0012-z.grd'), 
31    os.path.join(dir,'Java-0013-z.grd'), 
32    os.path.join(dir,'Java-0015-z.grd'), 
33    os.path.join(dir,'Java-0016-z.grd'), 
34    os.path.join(dir,'Java-0017-z.grd'), 
35    os.path.join(dir,'Java-0018-z.grd'), 
36    os.path.join(dir,'Java-0020-z.grd'), 
37    os.path.join(dir,'Java-0021-z.grd'), 
38    os.path.join(dir,'Java-0022-z.grd'), 
39    os.path.join(dir,'Java-0023-z.grd'), 
40    os.path.join(dir,'Java-0025-z.grd'), 
41    os.path.join(dir,'Java-0026-z.grd'), 
42    os.path.join(dir,'Java-0027-z.grd'), 
43    os.path.join(dir,'Java-0028-z.grd'), 
44    os.path.join(dir,'Java-0029-z.grd'), 
45    os.path.join(dir,'Java-0030-z.grd'), 
46    os.path.join(dir,'Java-0031-z.grd'), 
47    os.path.join(dir,'Java-0032-z.grd'), 
48    os.path.join(dir,'Java-0033-z.grd'), 
49    os.path.join(dir,'Java-0034-z.grd'), 
50    os.path.join(dir,'Java-0035-z.grd'), 
51    os.path.join(dir,'Java-0036-z.grd'), 
52    os.path.join(dir,'Java-0037-z.grd'), 
53    os.path.join(dir,'Java-0038-z.grd'), 
54    os.path.join(dir,'Java-0039-z.grd'), 
55    os.path.join(dir,'Java-0040-z.grd'), 
56    os.path.join(dir,'Java-0041-z.grd'), 
57    os.path.join(dir,'Java-0042-z.grd'), 
58    os.path.join(dir,'Java-0043-z.grd'), 
59    os.path.join(dir,'Java-0044-z.grd'), 
60    os.path.join(dir,'Java-0045-z.grd'), 
61    os.path.join(dir,'Java-0046-z.grd'), 
62    os.path.join(dir,'Java-0047-z.grd'),
63    os.path.join(dir,'Java-0048-z.grd')] 
64 
65print 'number of sources', len(urs_filenames)
66
67# AS per David Burbidge email on Friday 4th July the mag 9.3 event
68# has 1m worth of slip on each sub fault therefore mutliple each unit
69# source by the slip (10.4544) and sum the 44 time series together to
70# get the time series for this event at the points on your boundary.
71
72weight_factor = 10.4544
73weights=weight_factor*ones(len(urs_filenames),Float)
74
75scenario_name=project.scenario_name+'_all'
76order_filename=os.path.join(project.boundaries_dir, 'thinned_bounding_order.txt')
77
78print 'reading', order_filename
79# Create ordered sts file
80print 'creating sts file'
81
82urs2sts(urs_filenames,
83        basename_out=os.path.join(project.boundaries_dir,scenario_name),
84        ordering_filename=order_filename,
85        weights=weights,
86        mean_stage=project.tide,
87        verbose=True)
Note: See TracBrowser for help on using the repository browser.