Changeset 2418

Timestamp:

Feb 16, 2006, 4:26:43 PM (19 years ago)

Author:

nicholas

Message:

Crude comparison file CCS.py working fine.

File:

• development/momentum_sink/CCS.py (modified) (4 diffs)

development/momentum_sink/CCS.py

@@ -12 +12 @@
     
 
-swwfile_B = project.outputdir + sep  + 'friction_n=10_3135.sww'
-swwfile_F = project.outputdir + sep  + 'friction_n=20_3135.sww'
+swwfile_B = project.outputdir + sep  + 'Buildings_3662.sww'
+
 
 gauge_depth = Numeric.arrayrange(0, 2.5*WidtH, 5)
@@ -27 +27 @@
 quantities = ['stage', 'elevation', 'xmomentum', 'ymomentum']
 
+
+
 f_B = file_function(swwfile_B,
                   quantities = quantities,
@@ -32 +34 @@
                   verbose = True,
                   use_cache = True)
+fric_OK=[]
+Friction = [1,2,3,5,7,10,20,30,50]
+for fric in Friction:
+    swwfile_F = project.outputdir + sep  + 'friction_n=' + repr(fric) + '_3135.sww'
+    f_F = file_function(swwfile_F,
+                      quantities = quantities,
+                      interpolation_points = gauge_locations,
+                      verbose = True,
+                      use_cache = True)
 
-f_F = file_function(swwfile_F,
-                  quantities = quantities,
-                  interpolation_points = gauge_locations,
-                  verbose = True,
-                  use_cache = True)
-
-
-T = Numeric.arrayrange(0,500,20)
-#T = [ 20, 50, 100, 150, 200 ]
-
-for t in T:
-    model_time = []
+    t=500
+    
     stages_B = []
     stages_F = []
@@ -64 +65 @@
         velB = sqrt(uhB*uhB + vhB*vhB) / (wB + 1.e-30) #Absolute velocity
         velF = sqrt(uhF*uhF + vhF*vhF) / (wF + 1.e-30) #Absolute velocity                   
-        #print velB, "buiding velocity"
-        #print velF, "friction velocity"
 
         stages_B.append(wB)
         stages_F.append(wF)
-        #elevations.append(zB)  #Should be constant over time
-        #momenta.append(mB)
-        #velocityB.append(velB)
-        #velocityF.append(velF)
         
-    #print len(stages), "stages"
+        
+    
     diff = abs(Numeric.array(stages_B) - Numeric.array(stages_F))
-    #print diff, "difference"
+
+    if max(diff) < 1.2:
+        print "_-_-_-_-_-_-_-_-_-_-_-_"
+        print fric, " this file is OK"
+        print "-_-_-_-_-_-_-_-_-_-_-_-"
+        fric_OK.append(fric)
+        #print swwfile_F, "this file OK"
+    else:
+        print fric, " this file is definelty NOT OK"
+print fric_OK, "THESE FRICTION VALUES ARE OK"
+
+    
 
 
-    ion()
-    hold(True)
-    subplot(211)
-    
-    plot(gauge_depth, stages_B, '-b',gauge_depth,stages_F,'-r')
-    title('Time_%d_ vs water depth (n=%f, BP=%f)' %(t, fric, porosity))
-    xlabel('gauge length(m)')
-    ylabel('water depths (m)')
-    #savefig('Comp_Time%d_B_vs_F' %t)
-    hold(False)
-#savefig('Final_comp_n=%s' %fric)
-subplot(212)
-subplots_adjust(hspace = 0.4)
-plot(gauge_depth, diff, 'g')
-title('Final Differences in water levels')
-xlabel('gauge length(m)')
-ylabel('water differences (m)')
-savefig('Final_Water_diff_n=%s' %fric)
-print "finished"
-
-