Changeset 6398


Ignore:
Timestamp:
Feb 24, 2009, 12:28:04 PM (15 years ago)
Author:
myall
Message:

ceduna working using older scripts - going to upgrade to the latest ones anyway, so this isn't very important

Location:
anuga_work/production/australia_ph2/ceduna
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/ceduna/build_boundary_64468.py

    r6342 r6398  
    6363            ordering_filename = urs_order,
    6464            weights = weights,
     65            zone=53,
    6566            mean_stage = project.tide,
    6667            verbose = True)
     
    9394
    9495    maxname = 'max_sts_stage.csv'
    95     fid_max = open(project.boundaries_dir_event+sep+maxname,'w')
     96    fid_max = open(project.event_folder+sep+maxname,'w')
    9697    s = 'index, x, y, max_stage \n'
    9798    fid_max.write(s)   
     
    110111
    111112    minname = 'min_sts_stage.csv'
    112     fid_min = open(project.boundaries_dir_event+sep+minname,'w')
     113    fid_min = open(project.event_folder+sep+minname,'w')
    113114    s = 'index, x, y, max_stage \n'
    114115    fid_min.write(s)   
     
    124125
    125126
    126         fid_sts = open(project.boundaries_dir_event+sep+basename+'_'+ str(index)+'.csv', 'w')
     127        fid_sts = open(project.event_folder+sep+basename+'_'+ str(index)+'.csv', 'w')
    127128        s = 'time, stage, xmomentum, ymomentum \n'
    128129        fid_sts.write(s)
     
    142143    return quantities,elevation,time
    143144
    144 quantities,elevation,time=get_sts_gauge_data(os.path.join(project.boundaries_dir_event,project.scenario_name),verbose=False)
     145quantities,elevation,time=get_sts_gauge_data(project.event_sts,verbose=False)
    145146
    146147print len(elevation), len(quantities['stage'][0,:])
  • anuga_work/production/australia_ph2/ceduna/project.py

    r6359 r6398  
    4040# run_scenario script is executed
    4141tide = 0
    42 event_number = 27255 # Java 9.3 worst case for Perth
     42#event_number = 27255 # Java 9.3 worst case for Perth
     43event_number = 64468 # event for SA coast
    4344alpha = 0.1             # smoothing parameter for mesh
    4445friction = 0.01           # manning's friction coefficient
     
    104105# Thinned ordering file from Hazard Map (geographic)
    105106# Format is index,latitude,longitude (with header)
    106 urs_order_filename = 'urs_order.csv'
     107urs_order_filename = 'urs_order_austtg.csv'
    107108
    108109# Landward bounding points
  • anuga_work/production/australia_ph2/ceduna/run_model.py

    r6370 r6398  
    6060# Read in boundary from ordered sts file
    6161event_sts = create_sts_boundary(project.event_sts)
     62print 'HELLO event_sts', event_sts
    6263
    6364# Reading the landward defined points, this incorporates the original clipping
    6465# polygon minus the 100m contour
    6566landward_boundary = read_polygon(project.landward_boundary)
     67print 'landward boundary', landward_boundary
    6668
    6769# Combine sts polyline with landward points
    6870bounding_polygon_sts = event_sts + landward_boundary
     71print 'bounding polygon', bounding_polygon_sts
    6972
    7073# Number of boundary segments
Note: See TracChangeset for help on using the changeset viewer.