Changeset 217

Timestamp:

02/21/07 15:59:11 (17 years ago)

Author:

smasson

Message:

minor bugfix in computegrid

Location:

trunk/SRC

Files:

• Documentation/xmldoc/savesaxo.sh (modified) (1 diff)
• Grid/computegrid.pro (modified) (1 diff)

trunk/SRC/Documentation/xmldoc/savesaxo.sh

@@ -96 +96 @@
 dooldtest=0
 doidlwavedoc=0
-doidldoc=0svn 
+doidldoc=0
 doTARsrc=0
 doTARdata=0

trunk/SRC/Grid/computegrid.pro

@@ -1014 +1014 @@
 ; default mask eq 1
   if NOT keyword_set(mask) then tmask = -1 ELSE tmask = mask
-  IF keyword_set(romsh) AND tmask[0] NE -1 THEN tmask = tmask[0:jpiglo-1, 0:jpjglo-1]
-  IF n_elements(mask) EQ jpiglo*jpjglo AND jpkglo GT 1 THEN BEGIN
-    tmask = tmask[*]#replicate(1, jpkglo)
-    tmask = reform(tmask, jpiglo, jpjglo, jpkglo, /overwrite)
-  ENDIF  
-  IF jpiglo EQ 1 OR jpjglo EQ 1 THEN tmask = reform(tmask, jpiglo, jpjglo, jpkglo, /overwrite)
 ;
   if tmask[0] NE -1 then BEGIN
+    IF keyword_set(romsh) THEN tmask = tmask[0:jpiglo-1, 0:jpjglo-1]
+    IF n_elements(mask) EQ jpiglo*jpjglo AND jpkglo GT 1 THEN BEGIN
+      tmask = tmask[*]#replicate(1, jpkglo)
+      tmask = reform(tmask, jpiglo, jpjglo, jpkglo, /overwrite)
+    ENDIF  
+    IF jpiglo EQ 1 OR jpjglo EQ 1 THEN tmask = reform(tmask, jpiglo, jpjglo, jpkglo, /overwrite)
     tmask = byte(tmask[ixminmesh:ixmaxmesh, iyminmesh:iymaxmesh, izminmesh:izmaxmesh])
     tmask = reform(tmask, jpi, jpj, jpk, /over)