Changeset 5397

Timestamp:

Jun 10, 2008, 4:08:35 PM (15 years ago)

Author:

Leharne

Message:

Output timeseries at gauge locations to a new text file

File:

• anuga_work/development/convergence_okushiri_2008/compare_timeseries.py (modified) (2 diffs)

anuga_work/development/convergence_okushiri_2008/compare_timeseries.py

@@ -98 +98 @@
 from anuga.utilities.polygon import is_inside_polygon
 
-q = get_maximum_inundation_elevation(sww_filename)
-loc = get_maximum_inundation_location(sww_filename)
+q = get_maximum_inundation_elevation(sww_filename, tolerance=.4)
+loc = get_maximum_inundation_location(sww_filename, tolerance=.4)
 
 print 'Max runup elevation: ', q
@@ -152 +152 @@
 
 
+    # Write timeseries at gauge locations to file
+    filename=name+'_truescale_model.txt'
+    fid=open(filename,'w')
+    fid.write('Time (s), Stage (m)\n')
+    for i, t in enumerate(reference_time):
+        v=model[i]
+        fid.write('%.2f, %.2f\n' %(t, v))
 
+    fid.close()        
+    
 
     if plotting is True: