Changeset 9078

Timestamp:

Apr 10, 2014, 3:01:05 PM (11 years ago)

Author:

davies

Message:

Forcing 'height' to be non-negative for plot_utils.Make_Geotif
(Tsunami algorithm would not necessarily do this automatically)

File:

• trunk/anuga_core/source/anuga/utilities/plot_utils.py (modified) (2 diffs)

trunk/anuga_core/source/anuga/utilities/plot_utils.py

@@ -738 +738 @@
                 if(output_quantity=='depth'):
                     gridq=p2.height[myTS,:][gridqInd]
+                    gridq=gridq*(gridq>=0.) # Force positive depth (tsunami alg)
                 if(output_quantity=='velocity'):
                     gridq=p2.vel[myTS,:][gridqInd]
@@ -751 +752 @@
                 if(output_quantity=='depth'):
                     gridq=p2.height.max(axis=0)[gridqInd]
+                    gridq=gridq*(gridq>=0.) # Force positive depth (tsunami alg)
                 if(output_quantity=='velocity'):
                     gridq=p2.vel.max(axis=0)[gridqInd]