Changeset 8251 for trunk/anuga_work/development/2010-projects/anuga_1d/base/generic_domain.py

Timestamp:

Nov 18, 2011, 7:55:23 PM (12 years ago)

Author:

steve

Message:

Tracking down NaN in 1d anuga code

File:

• trunk/anuga_work/development/2010-projects/anuga_1d/base/generic_domain.py (modified) (3 diffs)

trunk/anuga_work/development/2010-projects/anuga_1d/base/generic_domain.py

@@ -1059 +1059 @@
         """
 
+        import numpy as np
+        Stage      = self.quantities['stage']
+        print 'w_ex_update', np.any(np.isnan(Stage.explicit_update))
 
         # Compute fluxes across each element edge
         self.compute_fluxes()
 
+        print 'w_ex_update', np.any(np.isnan(Stage.explicit_update))
+
         # Update timestep to fit yieldstep and finaltime
-        self.update_timestep(yieldstep, finaltime)
+        self.update_timestep(yieldstep, finaltime)     
 
         # Update conserved quantities
@@ -1390 +1395 @@
 
 
+        Stage      = self.quantities['stage']
+        Xmom       = self.quantities['xmomentum']
+        Bed        = self.quantities['elevation']
+        Height     = self.quantities['height']
+        Velocity   = self.quantities['velocity']
+
+        #Arrays
+        w_C   = Stage.centroid_values
+        uh_C  = Xmom.centroid_values
+        z_C   = Bed.centroid_values
+        h_C   = Height.centroid_values
+        u_C   = Velocity.centroid_values
+
+        import numpy as np
+        print '== before forcing ======='
+        print 'w_C', np.any(np.isnan(w_C))
+        print 'uh_C', np.any(np.isnan(uh_C))
+        print 'z_C', np.any(np.isnan(z_C))
+        print 'h_C', np.any(np.isnan(h_C))
+        print 'u_C', np.any(np.isnan(u_C))
+        print 'w_ex_update', np.any(np.isnan(Stage.explicit_update))
+
         timestep = self.timestep
+
+
 
         #Compute forcing terms
         self.compute_forcing_terms()
+
+        print '==after forcing ======='
+        print 'w_C', np.any(np.isnan(w_C))
+        print 'uh_C', np.any(np.isnan(uh_C))
+        print 'z_C', np.any(np.isnan(z_C))
+        print 'h_C', np.any(np.isnan(h_C))
+        print 'u_C', np.any(np.isnan(u_C))
+        print 'w_ex_update', np.any(np.isnan(Stage.explicit_update))
 
         #Update conserved_quantities
@@ -1400 +1437 @@
             Q.update(timestep)
 
-
+        print '==after quantity update  ======='
+        print 'w_C', np.any(np.isnan(w_C))
+        print 'uh_C', np.any(np.isnan(uh_C))
+        print 'z_C', np.any(np.isnan(z_C))
+        print 'h_C', np.any(np.isnan(h_C))
+        print 'u_C', np.any(np.isnan(u_C))
+        print 'w_ex_update', np.any(np.isnan(Stage.explicit_update))
 
 if __name__ == "__main__":