source: anuga_work/production/perth/build_boundary.py @ 5467

from anuga.utilities.numerical_tools import ensure_numeric
from Scientific.IO.NetCDF import NetCDFFile
from Numeric import asarray,transpose,sqrt,argmax,argmin,arange,Float,\
    compress,zeros,fabs,allclose,ones
from anuga.utilities.polygon import inside_polygon,read_polygon
from os import sep
from time import localtime, strftime, gmtime
from anuga.shallow_water.data_manager import urs2sts,create_sts_boundary
import os
import project


######################################
# Create urs boundary from mux2files #
######################################
dir= os.path.join(project.muxhome,'mux')
print dir 
prefix=os.path.join(dir,'Java-00')
suffix='-z.grd'

urs_filenames=[]
for i in range(49):
    if i<10:
        urs_filenames.append(prefix+'0'+str(i)+suffix)
    else:
        urs_filenames.append(prefix+str(i)+suffix)
print urs_filenames

weights=10.4544*ones(len(urs_filenames),Float)

base_name='tide_polyline'
order_filename='boxingday_boundary_order.txt'

# Create ordered sts file
print 'creating sts file'
urs2sts(urs_filenames,basename_out=base_name,
ordering_filename=order_filename,
weights=weights,
mean_stage=project.tide,
verbose=False)

# Read in boundary from ordered sts file
urs_boundary=create_sts_boundary(base_name)