Changeset 234 for inundation/ga/storm_surge/pyvolution/quantity.py

Timestamp:

Aug 27, 2004, 4:33:24 PM (21 years ago)

Author:

ole

Message:

Played with protection against mnimal heights and
balanced limiting. Looks good.
The demo show_balanced_limiters now runs!

File:

• inundation/ga/storm_surge/pyvolution/quantity.py (modified) (3 diffs)

inundation/ga/storm_surge/pyvolution/quantity.py

@@ -183 +183 @@
                 x2, y2 = self.domain.centroids[k2] #V2 centroid
 
-                #if k < 4:
-                #    print 'k, k0, k1, k2 = %d, %d, %d, %d' %(k, k0, k1, k2)
-                #    #print 'k = %d: (%f, %f), (%f, %f), (%f, %f)'\
-                #    #      %(k, x0, y0, x1, y1, x2, y2)
-                #    print 'k = %d: (%f, %f, %f)' %(k, q0, q1, q2)   
-                
                 #Gradient
                 a[k], b[k] = gradient(x0, y0, x1, y1, x2, y2, q0, q1, q2)
@@ -220 +214 @@
             dq[:,i] = qv[:,i] - qc
         
-        #Find min and max at this and neighbour's centroids
+        #Find min and max of this and neighbour's centroid values
         qmax = zeros(qc.shape, Float)
         qmin = zeros(qc.shape, Float)        
@@ -275 +269 @@
         
         a, b = self.compute_gradients()
-        ##for k in range(4):
-        ##    print 'Gradients %d: %16.12f, %16.12f' %(k, a[k], b[k])
 
         V = self.domain.get_vertex_coordinates()