Changeset 3948

Timestamp:

Nov 8, 2006, 6:14:08 PM (18 years ago)

Author:

ole

Message:

Work in progress

File:

• anuga_core/source/anuga/shallow_water/data_manager.py (modified) (3 diffs)

anuga_core/source/anuga/shallow_water/data_manager.py

@@ -380 +380 @@
 
         # FIXME (HACK)
-        truncation = self.domain.number_of_full_nodes        
-        #print len(z), len(Z), truncation
+        if len(z) <> len(Z):
+            truncation = self.domain.number_of_full_nodes        
+            Z = Z[:truncation]
+            print len(z), len(Z), truncation, len(self.domain.get_nodes())
         
-        z[:] = Z[:truncation].astype(self.precision)
+        z[:] = Z.astype(self.precision)
 
         #FIXME: Backwards compatibility
         z = fid.variables['z']
-        z[:] = Z[:truncation].astype(self.precision)
+        z[:] = Z.astype(self.precision)
         ################################
 
@@ -501 +503 @@
                 # HACK
                 truncation = self.domain.number_of_full_nodes
-                A = A[:truncation]
+                
 
                 #FIXME: Make this general (see below)
@@ -508 +510 @@
                     #print z[:]
                     #print A-z[:]
+
+                    # HACK
+                    if len(z) <> len(A):                
+                        A = A[:truncation]
+
+                    
                     A = choose(A-z[:] >= self.minimum_storable_height,
                                (z[:], A))