Changeset 217 for trunk/SRC/Grid/computegrid.pro
- Timestamp:
- 02/21/07 15:59:11 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Grid/computegrid.pro
r213 r217 1014 1014 ; default mask eq 1 1015 1015 if NOT keyword_set(mask) then tmask = -1 ELSE tmask = mask 1016 IF keyword_set(romsh) AND tmask[0] NE -1 THEN tmask = tmask[0:jpiglo-1, 0:jpjglo-1]1017 IF n_elements(mask) EQ jpiglo*jpjglo AND jpkglo GT 1 THEN BEGIN1018 tmask = tmask[*]#replicate(1, jpkglo)1019 tmask = reform(tmask, jpiglo, jpjglo, jpkglo, /overwrite)1020 ENDIF1021 IF jpiglo EQ 1 OR jpjglo EQ 1 THEN tmask = reform(tmask, jpiglo, jpjglo, jpkglo, /overwrite)1022 1016 ; 1023 1017 if tmask[0] NE -1 then BEGIN 1018 IF keyword_set(romsh) THEN tmask = tmask[0:jpiglo-1, 0:jpjglo-1] 1019 IF n_elements(mask) EQ jpiglo*jpjglo AND jpkglo GT 1 THEN BEGIN 1020 tmask = tmask[*]#replicate(1, jpkglo) 1021 tmask = reform(tmask, jpiglo, jpjglo, jpkglo, /overwrite) 1022 ENDIF 1023 IF jpiglo EQ 1 OR jpjglo EQ 1 THEN tmask = reform(tmask, jpiglo, jpjglo, jpkglo, /overwrite) 1024 1024 tmask = byte(tmask[ixminmesh:ixmaxmesh, iyminmesh:iymaxmesh, izminmesh:izmaxmesh]) 1025 1025 tmask = reform(tmask, jpi, jpj, jpk, /over)
Note: See TracChangeset
for help on using the changeset viewer.