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

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 #
######################################

#print project.muxhome
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=project.scenario_name
order_filename=os.path.join(project.boundaries_dir,'boundary_ordering.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)