;------------------------------------------------------------
;------------------------------------------------------------
;------------------------------------------------------------
;+
; NAME:vmask
;
; PURPOSE:calcule vmask
;
; CATEGORY:fonction economisatrice de memoire (plus besion de garder
; vmask)  
;
; CALLING SEQUENCE:res=vmask()
;
; OUTPUTS:un tableau 3d correspondant a vmask
;
; EXAMPLE:s''utilise comme si vmask etait un tableau connu en
; replacant ds la syntaxe vmask par vmask(). par ex:
; au lieu de taper vmask[domainef] il faut taper vmask()[domainef]
;
; MODIFICATION HISTORY: Sebastien Masson (smasson@lodyc.jussieu.fr)
;                       11/12/98
;                       20/9/99 cas jpk=1 (merci jpblod@ipsl.jussieu.fr)
;-
;------------------------------------------------------------
;------------------------------------------------------------
;------------------------------------------------------------
FUNCTION vmask
;
  compile_opt idl2, strictarrsubs
;
@common
  tempsun = systime(1)          ; pour key_performance
;
  CASE size(tmask, /n_dimensions) OF
    2:res = tmask*shift(tmask, 0, -1)
    3:res = tmask*shift(tmask, 0, -1, 0)
  ENDCASE
;
  res[*, jpj-1, *] = vmaskred
  if keyword_set(key_performance) THEN print, 'temps vmask', systime(1)-tempsun 
;
  return, res
end