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

updated sts boundary

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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
Note: See TracChangeset for help on using the changeset viewer.