Changeset 13 for trunk/ToBeReviewed/GRILLE/fmask.pro

Timestamp:

04/28/06 14:18:03 (18 years ago)

Author:

pinsard

Message:

upgrade of GRILLE/Utilities according to cerbere.lodyc.jussieu.fr: /usr/home/smasson/SAXO_RD/ : files

File:

• trunk/ToBeReviewed/GRILLE/fmask.pro (copied) (copied from trunk/GRILLE/fmask.pro) (1 diff)

trunk/ToBeReviewed/GRILLE/fmask.pro

@@ -25 +25 @@
 ;------------------------------------------------------------
 FUNCTION fmask
-@common
-   tempsun = systime(1)         ; pour key_performance
-   if jpk EQ 1 then begin
-      res = tmask*shift(tmask, -1, 0)*shift(tmask, 0, -1)*shift(tmask, -1, -1)
-      if NOT keyword_set(key_periodique) then res[jpi-1, *] = fmaskredy
-      res[*, jpj-1] = fmaskredx
-   ENDIF ELSE BEGIN
-      res = tmask*shift(tmask, -1, 0, 0)*shift(tmask, 0, -1, 0)*shift(tmask, -1, -1, 0)
-      if NOT keyword_set(key_periodique) then res[jpi-1, *, *] = fmaskredy
-      res[*, jpj-1, *] = fmaskredx
-   ENDELSE
-   if keyword_set(key_performance) THEN print, 'temps fmask', systime(1)-tempsun 
-
-   return, res
+;---------------------------------------------------------
+@cm_4mesh
+  IF NOT keyword_set(key_forgetold) THEN BEGIN
+@updatenew
+  ENDIF
+;---------------------------------------------------------
+  tempsun = systime(1)          ; pour key_performance
+; 
+  CASE size(tmask, /n_dimensions) OF
+    2:res = tmask*shift(tmask, -1, 0)*shift(tmask, 0, -1)*shift(tmask, -1, -1)
+    3:res = tmask*shift(tmask, -1, 0, 0)*shift(tmask, 0, -1, 0)*shift(tmask, -1, -1, 0)
+  ENDCASE
+;
+  if NOT keyword_set(key_periodic) then res[jpi-1, *, *] = fmaskredy
+  res[*, jpj-1, *] = fmaskredx
+;
+  if keyword_set(key_performance) THEN print, 'temps fmask', systime(1)-tempsun 
+  
+  return, res
 end