Changeset 5578


Ignore:
Timestamp:
Jul 25, 2008, 5:21:30 PM (11 years ago)
Author:
kristy
Message:

updated sts boundary

Location:
anuga_work/production
Files:
4 edited

Legend:

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

    r5575 r5578  
    5252        verbose=True)
    5353
    54 # Read in boundary from ordered sts file
    55 urs_boundary_polygon=create_sts_boundary(scenario_name)
    5654
    5755
     
    6058
    6159
    62 
  • anuga_work/production/busselton/run_busselton.py

    r5575 r5578  
    1818from os import sep
    1919from os.path import dirname, basename
     20import os
    2021from os import mkdir, access, F_OK
    2122from shutil import copy
     
    3031from anuga.shallow_water import Field_boundary
    3132from Numeric import allclose
    32 from anuga.shallow_water.data_manager import export_grid
     33from anuga.shallow_water.data_manager import export_grid, create_sts_boundary
    3334
    3435from anuga.pmesh.mesh_interface import create_mesh_from_regions
     
    7778
    7879    # Read in boundary from ordered sts file
    79     urs_bounding_polygon=create_sts_boundary(project.scenario_name)
     80    urs_bounding_polygon=create_sts_boundary(os.path.join(project.boundaries_dir,project.scenario_name))
    8081
    8182    # Reading the landward defined points, this incorporates the original clipping
    8283    # polygon minus the 100m contour
    83     landward_bounding_polygon = read_polygon(project.polygons_dir+'landward_bounding_polygon.txt')
     84    landward_bounding_polygon = read_polygon(project.boundaries_dir+'landward_bounding_polygon.txt')
    8485
    8586    # Combine sts polyline with landward points
     
    104105        print 'start create mesh from regions'
    105106       
    106         create_mesh_from_regions(project.poly_all,
    107                              boundary_tags=project.boundary_tags,
     107        create_mesh_from_regions(bounding_polygon,
     108                             boundary_tags=boundary_tags,
    108109                             maximum_triangle_area=project.res_poly_all,
    109110                             interior_regions=project.interior_regions,
     
    255256    kwargs['bathy_file']=project.combined_dir_name+'.txt'
    256257#    kwargs['bathy_file']=project.combined_small_dir_name + '.pts'
    257     kwargs['boundary_file']=project.boundaries_in_dir_name + '.sww'
     258 
    258259    kwargs['file_name']=project.home+'detail.csv'
    259260    kwargs['aa_scenario_name']=project.scenario_name
  • anuga_work/production/perth/build_boundary.py

    r5576 r5578  
    2121suffix='-z.grd'
    2222
    23 
    24 
    25 
    2623urs_filenames={}
    2724for filename in os.listdir(dir):
     
    2926
    3027    urs_filenames[basename] = True
    31 
    3228
    3329# Make list of basenames 
     
    5551        verbose=True)
    5652
    57 # Read in boundary from ordered sts file
    58 urs_boundary_polygon=create_sts_boundary(scenario_name)
    5953
    6054
    6155
    6256
    63 
    64 
    65 
  • anuga_work/production/perth/run_perth.py

    r5569 r5578  
    1717# Standard modules
    1818from os import sep
     19import os
    1920from os.path import dirname, basename
    2021from os import mkdir, access, F_OK
     
    7778
    7879    # Read in boundary from ordered sts file
    79     urs_bounding_polygon=create_sts_boundary(project.scenario_name)
     80    urs_bounding_polygon=create_sts_boundary(os.path.join(project.boundaries_dir,project.scenario_name))
    8081
    8182    # Reading the landward defined points, this incorporates the original clipping
Note: See TracChangeset for help on using the changeset viewer.