source: trunk/SRC/ToBeReviewed/GRILLE/fmask.pro @ 226

Last change on this file since 226 was 226, checked in by pinsard, 17 years ago

corrections of some misspellings in some *.pro
  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 1.5 KB
Line 
1;------------------------------------------------------------
2;------------------------------------------------------------
3;------------------------------------------------------------
4;+
5;
6; @file_comments
7; calculate fmask. Economize memory (we do not need to keep fmask anymore)
8;
9; @categories
10; Utilities
11;
12; @returns
13; A 3d array corresponding to fmask.
14;
15; @examples
16; Used like if fmask was a known array by substituing fmask by (fmask())
17; For example: Instead of write fmask[domainef] we have to write (fmask())[domainef]
18;
19; @history
20; Sebastien Masson (smasson\@lodyc.jussieu.fr)
21;                       11/12/98
22;
23; @version
24; $Id$
25;
26;-
27;------------------------------------------------------------
28;------------------------------------------------------------
29;------------------------------------------------------------
30FUNCTION fmask
31;---------------------------------------------------------
32;
33  compile_opt idl2, strictarrsubs
34;
35@cm_4mesh
36  IF NOT keyword_set(key_forgetold) THEN BEGIN
37@updatenew
38  ENDIF
39;---------------------------------------------------------
40  tempsun = systime(1)          ; For key_performance
41;
42  CASE size(tmask, /n_dimensions) OF
43    2:res = tmask*shift(tmask, -1, 0)*shift(tmask, 0, -1)*shift(tmask, -1, -1)
44    3:res = tmask*shift(tmask, -1, 0, 0)*shift(tmask, 0, -1, 0)*shift(tmask, -1, -1, 0)
45  ENDCASE
46;
47  if NOT keyword_set(key_periodic) then res[jpi-1, *, *] = fmaskredy
48  res[*, jpj-1, *] = fmaskredx
49;
50  if keyword_set(key_performance) THEN print, 'temps fmask', systime(1)-tempsun
51
52  return, res
53end
Note: See TracBrowser for help on using the repository browser.