Changeset 450 for trunk/SRC/Computation/e3t_3d.pro
- Timestamp:
- 01/31/11 16:41:25 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Computation/e3t_3d.pro
r447 r450 44 44 keyword_set(e1) AND keyword_set(e2):arr2d = e1t[fstx:lstx, fsty:lsty] * e2t[fstx:lstx, fsty:lsty] 45 45 keyword_set(e1) :arr2d = e1t[fstx:lstx, fsty:lsty] 46 46 keyword_set(e2):arr2d = e2t[fstx:lstx, fsty:lsty] 47 47 ELSE:arr2d = replicate(1., nx*ny) 48 48 ENDCASE … … 54 54 ; level of the bottom of the ocean 55 55 bottom = total(tmask[fstx:lstx, fsty:lsty, *], 3) 56 sea = where(bottom NE 0) 56 57 bottom2 = long(temporary(bottom)) - firstzt 57 58 bottom = bottom2 - 1L 58 ok = where(bottom GE 0 AND bottom LT nzt)59 ok2 = where(bottom2 GE 0 AND bottom2 LT nzt)59 ok = inter(sea, where(bottom GE 0 AND bottom LT nzt)) 60 ok2 = inter(sea, where(bottom2 GE 0 AND bottom2 LT nzt)) 60 61 ; apply e3t_ps to e3t_3D at the bottom of the ocean 61 62 IF ok[0] NE -1 THEN BEGIN
Note: See TracChangeset
for help on using the changeset viewer.