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

Timestamp:

01/25/11 12:24:31 (13 years ago)

Author:

smasson

Message:

bugfix related to changeset:445

File:

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

trunk/SRC/Computation/e3w_3d.pro

@@ -55 +55 @@
     IF ok[0]  NE -1 THEN BEGIN
 ; the bottom of the ocean in 3D index is:
-      bottom  = lindgen(nxt*nyt) + nxt*nyt*(temporary(bottom ))[ok ]
+      bottom  = (lindgen(nxt*nyt))[ok ] + nxt*nyt*(temporary(bottom ))[ok ]
       e3w_3d[bottom ] = arr2d[ok ] * (e3w_ps[firstxt:lastxt, firstyt:lastyt])[ok ]
     ENDIF
     IF ok2[0] NE -1 THEN BEGIN
-      bottom2 = lindgen(nxt*nyt) + nxt*nyt*(temporary(bottom2))[ok2]
+      bottom2 = (lindgen(nxt*nyt))[ok2] + nxt*nyt*(temporary(bottom2))[ok2]
       e3w_3d[bottom2] = arr2d[ok2] * (e3t_ps[firstxt:lastxt, firstyt:lastyt])[ok2] ; use e3t_ps and not e3w_ps
     ENDIF