- Timestamp:
- 01/25/11 12:24:31 (13 years ago)
- Location:
- trunk/SRC/Computation
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Computation/e3t_3d.pro
r446 r447 61 61 IF ok[0] NE -1 THEN BEGIN 62 62 ; the bottom of the ocean in 3D index is: 63 bottom = lindgen(nx*ny)+ nx*ny*(temporary(bottom ))[ok ]63 bottom = (lindgen(nx*ny))[ok ] + nx*ny*(temporary(bottom ))[ok ] 64 64 e3t3d[temporary(bottom )] = arr2d[ok ] * (e3t_ps[fstx:lstx, fsty:lsty])[ok ] 65 65 ENDIF 66 66 IF ok2[0] NE -1 THEN BEGIN 67 bottom2 = lindgen(nx*ny)+ nx*ny*(temporary(bottom2))[ok2]67 bottom2 = (lindgen(nx*ny))[ok2] + nx*ny*(temporary(bottom2))[ok2] 68 68 e3t3d[temporary(bottom2)] = arr2d[ok2] * (e3t_ps[fstx:lstx, fsty:lsty])[ok2] 69 69 ENDIF -
trunk/SRC/Computation/e3w_3d.pro
r446 r447 55 55 IF ok[0] NE -1 THEN BEGIN 56 56 ; the bottom of the ocean in 3D index is: 57 bottom = lindgen(nxt*nyt)+ nxt*nyt*(temporary(bottom ))[ok ]57 bottom = (lindgen(nxt*nyt))[ok ] + nxt*nyt*(temporary(bottom ))[ok ] 58 58 e3w_3d[bottom ] = arr2d[ok ] * (e3w_ps[firstxt:lastxt, firstyt:lastyt])[ok ] 59 59 ENDIF 60 60 IF ok2[0] NE -1 THEN BEGIN 61 bottom2 = lindgen(nxt*nyt)+ nxt*nyt*(temporary(bottom2))[ok2]61 bottom2 = (lindgen(nxt*nyt))[ok2] + nxt*nyt*(temporary(bottom2))[ok2] 62 62 e3w_3d[bottom2] = arr2d[ok2] * (e3t_ps[firstxt:lastxt, firstyt:lastyt])[ok2] ; use e3t_ps and not e3w_ps 63 63 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.