Changeset 5551


Ignore:
Timestamp:
Jul 22, 2008, 1:39:46 PM (11 years ago)
Author:
ole
Message:

Work towards using the STS boundary with the Perth production script

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/perth/build_boundary.py

    r5468 r5551  
    2121suffix='-z.grd'
    2222
    23 urs_filenames=[]
    24 for i in range(49):
    25     if i<10:
    26         urs_filenames.append(prefix+'0'+str(i)+suffix)
    27     else:
    28         urs_filenames.append(prefix+str(i)+suffix)
    29 print urs_filenames
     23
     24
     25
     26urs_filenames={}
     27for filename in os.listdir(dir):
     28    basename = filename[:-7]
     29
     30    urs_filenames[basename] = True
     31
     32
     33# Make list of basenames 
     34urs_filenames = [os.path.join(dir, basename) for basename in urs_filenames.keys()]
     35
     36# Reduce this list slightly
     37urs_filenames = urs_filenames[:30]
     38for name in urs_filenames:
     39    print name
    3040
    3141weights=10.4544*ones(len(urs_filenames),Float)
    3242
    3343base_name=project.scenario_name
    34 order_filename=os.path.join(project.boundaries_dir,'boundary_ordering.txt')
     44order_filename=os.path.join(project.boundaries_dir, 'boundary_ordering.txt')
    3545
     46print 'reading', order_filename
    3647# Create ordered sts file
    3748print 'creating sts file'
     
    4051        weights=weights,
    4152        mean_stage=project.tide,
    42         verbose=False)
     53        verbose=True)
    4354
    4455# Read in boundary from ordered sts file
Note: See TracChangeset for help on using the changeset viewer.