1 | from anuga.utilities.numerical_tools import ensure_numeric |
---|
2 | from Scientific.IO.NetCDF import NetCDFFile |
---|
3 | from Numeric import asarray,transpose,sqrt,argmax,argmin,arange,Float,\ |
---|
4 | compress,zeros,fabs,allclose,ones |
---|
5 | from anuga.utilities.polygon import inside_polygon,read_polygon |
---|
6 | from os import sep |
---|
7 | from time import localtime, strftime, gmtime |
---|
8 | from anuga.shallow_water.data_manager import urs2sts,create_sts_boundary |
---|
9 | import os |
---|
10 | import project |
---|
11 | |
---|
12 | |
---|
13 | ###################################### |
---|
14 | # Create urs boundary from mux2files # |
---|
15 | ###################################### |
---|
16 | dir= os.path.join(project.muxhome,'mux') |
---|
17 | print dir |
---|
18 | prefix=os.path.join(dir,'Java-00') |
---|
19 | suffix='-z.grd' |
---|
20 | |
---|
21 | urs_filenames=[] |
---|
22 | for i in range(49): |
---|
23 | if i<10: |
---|
24 | urs_filenames.append(prefix+'0'+str(i)+suffix) |
---|
25 | else: |
---|
26 | urs_filenames.append(prefix+str(i)+suffix) |
---|
27 | print urs_filenames |
---|
28 | |
---|
29 | weights=10.4544*ones(len(urs_filenames),Float) |
---|
30 | |
---|
31 | base_name='tide_polyline' |
---|
32 | order_filename='boxingday_boundary_order.txt' |
---|
33 | |
---|
34 | # Create ordered sts file |
---|
35 | print 'creating sts file' |
---|
36 | urs2sts(urs_filenames,basename_out=base_name, |
---|
37 | ordering_filename=order_filename, |
---|
38 | weights=weights, |
---|
39 | mean_stage=project.tide, |
---|
40 | verbose=False) |
---|
41 | |
---|
42 | # Read in boundary from ordered sts file |
---|
43 | urs_boundary=create_sts_boundary(base_name) |
---|
44 | |
---|