Changeset 4613


Ignore:
Timestamp:
Jul 9, 2007, 3:11:24 PM (17 years ago)
Author:
nick
Message:

update onslow

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/onslow_2006/boundary_plot.py

    r4594 r4613  
    1010time_dir2 = '20070608_062811_run_final_1.5_exmouth_nbartzis'# HAT exmouth
    1111directory = project.outputdir
     12boundary_dir=project.boundaries_in_dir
    1213
    1314def get_time_series_from_file(filename):
     
    3233max_st = max_sp = max_mom = min_st = min_sp = min_mom = max_t = 0.
    3334count = j = -1
    34 for i in range(47*2):
     35
     36#find maximums for axis for plots
     37for i in range(47*2):  #could alter get_all_swwfiles to get all file extension and get csv files??
    3538    count += 1
    3639    t = str(count/2.).rfind('.')
     
    3841        j += 1
    3942        extra = 'gauges_time_series_%d.csv' %(j)
    40         filename = directory + time_dir1 + sep + extra
     43        filename = directory + time_dir2 + sep + extra
    4144    else:
    4245        extra = 'gauges_time_series_%d.csv' %(j)
    43         filename = directory + time_dir2 + sep + extra
     46#        filename = directory + time_dir2 + sep + extra
     47        filename = boundary_dir + sep + extra
    4448
    4549    t, stage, depth, speed, momentum = get_time_series_from_file(filename)
     
    110114        ii = 0
    111115        extra = 'gauges_time_series_%d.csv' %(j)
    112         filename = directory + time_dir1 + sep + extra
     116        filename = directory + time_dir2 + sep + extra
    113117    else:
    114118        ii = 1
    115119        extra = 'gauges_time_series_%d.csv' %(j)
    116         filename = directory + time_dir2 + sep + extra
     120#        filename = directory + time_dir2 + sep + extra
     121        filename = boundary_dir + sep + extra
    117122       
    118123    t, stage, depth, speed, momentum = get_time_series_from_file(filename)
     
    122127    name = filename[i1:i2-1]
    123128    loc = filename.rfind('1.5') + 4
    124     if filename[loc] is 'o': community = 'onslow'
    125     if filename[loc] is 'e': community = 'exmouth'
     129#    if filename[loc] is 'o': community = 'onslow'
     130#    if filename[loc] is 'e': community = 'exmouth'
    126131
    127132    print 'hello', community, filename
Note: See TracChangeset for help on using the changeset viewer.