Changeset 454 for trunk/SRC/Computation/e3t_3d.pro

Timestamp:

04/07/11 10:03:36 (13 years ago)

Author:

smasson

Message:

bugfix to keep double precision of vertical scale factors in all cases

File:

• trunk/SRC/Computation/e3t_3d.pro (modified) (1 diff)

trunk/SRC/Computation/e3t_3d.pro

@@ -44 +44 @@
     keyword_set(e1) AND keyword_set(e2):arr2d = e1t[fstx:lstx, fsty:lsty] * e2t[fstx:lstx, fsty:lsty]
     keyword_set(e1)                    :arr2d = e1t[fstx:lstx, fsty:lsty]
-    keyword_set(e2):arr2d =                             e2t[fstx:lstx, fsty:lsty]
-    ELSE:arr2d = replicate(1., nx*ny)
+                        keyword_set(e2):arr2d =                             e2t[fstx:lstx, fsty:lsty]
+    ELSE:arr2d = replicate(1.d, nx*ny)
   ENDCASE
   e3t3d = arr2d[*] # e3t[firstzt:lastzt]