Changeset 5397


Ignore:
Timestamp:
Jun 10, 2008, 4:08:35 PM (14 years ago)
Author:
Leharne
Message:

Output timeseries at gauge locations to a new text file

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/convergence_okushiri_2008/compare_timeseries.py

    r5345 r5397  
    9898from anuga.utilities.polygon import is_inside_polygon
    9999
    100 q = get_maximum_inundation_elevation(sww_filename)
    101 loc = get_maximum_inundation_location(sww_filename)
     100q = get_maximum_inundation_elevation(sww_filename, tolerance=.4)
     101loc = get_maximum_inundation_location(sww_filename, tolerance=.4)
    102102
    103103print 'Max runup elevation: ', q
     
    152152
    153153
     154    # Write timeseries at gauge locations to file
     155    filename=name+'_truescale_model.txt'
     156    fid=open(filename,'w')
     157    fid.write('Time (s), Stage (m)\n')
     158    for i, t in enumerate(reference_time):
     159        v=model[i]
     160        fid.write('%.2f, %.2f\n' %(t, v))
    154161
     162    fid.close()       
     163   
    155164
    156165    if plotting is True:
Note: See TracChangeset for help on using the changeset viewer.