Changeset 5344

Timestamp:

May 20, 2008, 11:35:54 AM (16 years ago)

Author:

Leharne

Message:

Scaling the validation data to truescale

Location:

anuga_work/development/convergence_okushiri_2008

Files:

• compare_timeseries.py (modified) (5 diffs)
• output_truescale_ch5-7-9.txt (added)

anuga_work/development/convergence_okushiri_2008/compare_timeseries.py

@@ -28 +28 @@
 #-------------------------
 
-finaltime = 22.5
-timestep = 0.05
+finaltime = 450
+timestep = 1
 
-gauge_locations = [[0.000, 1.696]] # Boundary gauge
-gauge_locations += [[4.521, 1.196],  [4.521, 1.696],  [4.521, 2.196]] #Ch 5-7-9
+gauge_locations = [[0.000, 678.4]] # Boundary gauge
+gauge_locations += [[1808.4, 478.4],  [1808.4, 678.4],  [1808.4, 878.4]] #Ch 5-7-9
 gauge_names = ['Boundary', 'ch5', 'ch7', 'ch9']
 
@@ -44 +44 @@
 #-------------------------
 
-print 'Reading', project.boundary_filename
-fid = NetCDFFile(project.boundary_filename, 'r')
+print 'Reading', project_truescale.boundary_filename
+fid = NetCDFFile(project_truescale.boundary_filename, 'r')
 input_time = fid.variables['time'][:]
 validation_data['Boundary'] = fid.variables['stage'][:]
 
 reference_time = []
-fid = open(project.validation_filename)
+fid = open(project_truescale.validation_filename)
 lines = fid.readlines()
 fid.close()
@@ -62 +62 @@
     for j, key in enumerate(gauge_names[1:]):  # Omit boundary gauge
         value = float(fields[1:][j])           # Omit time 
-        validation_data[key].append(value/100) # Convert cm2m
+        validation_data[key].append(value) 
 
 
@@ -82 +82 @@
     sww_filename = sys.argv[1]
 else:    
-    sww_filename = project.output_filename
+    sww_filename = project_truescale.output_filename
     
 f = file_function(sww_filename,
@@ -102 +102 @@
 
 print 'Max runup elevation: ', q
-print 'Max runup elevation (scaled by 400): ', q*400
 print 'Max runup location:  ', loc
-
-#from create_okushiri import gulleys
-#assert is_inside_polygon(loc, gulleys)