Changeset 450 for trunk/SRC/Computation/e3w_3d.pro

Timestamp:

01/31/11 16:41:25 (13 years ago)

Author:

smasson

Message:

improve the way to take into account continents in e3t/w

File:

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

trunk/SRC/Computation/e3w_3d.pro

@@ -48 +48 @@
 ; level of the bottom of the ocean
     bottom = total(tmask[firstxt:lastxt, firstyt:lastyt, *], 3)
+    sea = where(bottom NE 0)
     bottom2 = long(temporary(bottom)) - firstzw
     bottom  = bottom2 - 1L
-    ok  = where(bottom  GE 0 AND bottom  LT nzw)
-    ok2 = where(bottom2 GE 0 AND bottom2 LT nzw)
+    ok  = inter(sea, where(bottom  GE 0 AND bottom  LT nzw))
+    ok2 = inter(sea, where(bottom2 GE 0 AND bottom2 LT nzw))
 ; apply e3w_ps to e3w_3D at the bottom of the ocean
     IF ok[0]  NE -1 THEN BEGIN