Changeset 1674

Timestamp:

Aug 2, 2005, 6:04:27 PM (19 years ago)

Author:

ole

Message:

Introduced file boundary to lwru2.py

File:

• inundation/ga/storm_surge/validation/LWRU2/lwru2.py (modified) (3 diffs)

inundation/ga/storm_surge/validation/LWRU2/lwru2.py

@@ -8 +8 @@
 
 """
+
+def prepare_timeboundary(filename):
+    """Converting benchmark 2 time series to NetCDF sww file.
+    """
+
+
+    print 'Preparing time boundary from %s' %filename
+    from Numeric import array
+
+    fid = open(filename)
+
+    #Skip first line
+    line = fid.readline()
+
+    lines = fid.readlines()
+    fid.close()
+
+    N = len(lines)
+
+    T = zeros(N, Float)       #Time
+    Q = zeros(N, Float)  #Values
+
+    for i, line in enumerate(lines):
+        fields = line.split()
+        
+        T[i] = float(fields[0])
+        Q[i] = float(fields[1])
+
+    #Create sww file
+    from Scientific.IO.NetCDF import NetCDFFile
+
+    fid = NetCDFFile(filename[-4:] + '.sww', 'w')
+
+    fid.institution = 'Geoscience Australia'
+    fid.description = 'Input wave for Benchmark 2'
+    fid.starttime = 0.0
+    fid.createDimension('number_of_timesteps', len(T))
+    fid.createVariable('time', Float, ('number_of_timesteps',))
+    fid.variables['time'][:] = T
+            
+    fid.createVariable('stage', Float, ('number_of_timesteps',))
+    fid.variables['stage'][:] = Q[:]
+
+    fid.createVariable('xmomentum', Float, ('number_of_timesteps',))
+    fid.variables['xmomentum'][:] = 0.0
+
+    fid.createVariable('ymomentum', Float, ('number_of_timesteps',))
+    fid.variables['ymomentum'][:] = 0.0                
+
+
+
+    fid.close()
+
 
 
@@ -24 +77 @@
 except:
 
+    #Preparing time boundary
+    prepare_timeboundary('Benchmark_2_input.txt')
+    
     #Preparing points
     from pyvolution.data_manager import xya2pts
-    xya2pts('bathymetry.txt', verbose = True)
+    xya2pts('Benchmark_2_Bathymetry.txt', verbose = True)
 
 
@@ -75 +131 @@
 print 'Boundaries'
 Br = Reflective_boundary(domain)
-Bf = File_boundary('input_wave.txt', domain)
+Bf = File_boundary('Benchmark_2_input.sww', domain, verbose=True)
 
 #Set boundary conditions