Changeset 1696

Timestamp:

Aug 8, 2005, 5:59:09 PM (19 years ago)

Author:

ole

Message:

Played with lwru2 validation (momentum of BC)

File:

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

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

@@ -59 +59 @@
 
     fid.createVariable('xmomentum', Float, ('number_of_timesteps',))
-    fid.variables['xmomentum'][:] = 0.0 #0.04
+    #fid.variables['xmomentum'][:] = 0.005
+    #fid.variables['xmomentum'][:-1] = Q[1:]-Q[:-1]/0.05 #Too high
+    #fid.variables['xmomentum'][:-1] = Q[1:]-Q[:-1]/0.5   #Too low
+    #fid.variables['xmomentum'][:-1] = Q[1:]-Q[:-1]/0.1   #Too high
+    #fid.variables['xmomentum'][:-1] = Q[1:]-Q[:-1]/0.25   #Too high
+    #fid.variables['xmomentum'][:-1] = Q[1:]-Q[:-1]/0.35 #Too low
+    #fid.variables['xmomentum'][:-1] = Q[1:]-Q[:-1]/0.3  #Too low
+    #fid.variables['xmomentum'][-1] = 0.0
+    
+
+    fid.variables['xmomentum'][:] = 0.003
+    
     fid.createVariable('ymomentum', Float, ('number_of_timesteps',))
     fid.variables['ymomentum'][:] = 0.0
@@ -83 +94 @@
     
 try:
-    raise Exception
+    #raise Exception
     print 'Read pickled domain'
     fid = open(picklefile)
@@ -157 +168 @@
 import time
 t0 = time.time()
-domain.visualise = True
+domain.visualise = False
 
 for t in domain.evolve(yieldstep = 0.05, finaltime = 22.5):