source: trunk/preproc_mask.pro @ 2

Last change on this file since 2 was 2, checked in by pinsard, 18 years ago

initial import from /usr/work/fvi/OPA/geomag/
File size: 1.4 KB
Line 
1;------------------------------------------------------------
2;------------------------------------------------------------
3;------------------------------------------------------------
4;+
5; NAME: preproc_mask.pro
6;
7; PURPOSE: Extends land-sea mask by one point over ocean
8;
9; CATEGORY: Function
10;
11; CALLING SEQUENCE: preproc_mask, mask
12;
13; INPUTS:
14;          mask : mask field composed of 0 and 1 values
15;                 with 0 over land points       
16;
17; KEYWORD PARAMETERS: None
18;
19; OUTPUTS:
20;          mask : mask field extended
21;                 
22;
23; COMMON BLOCKS:
24;          common_interp
25;
26; SIDE EFFECTS:
27;
28; RESTRICTIONS:
29;
30; EXAMPLE:
31;
32; MODIFICATION HISTORY: 11/99 A. Jouzeau
33;
34;-
35;------------------------------------------------------------
36;------------------------------------------------------------
37;------------------------------------------------------------
38FUNCTION preproc_mask, mask
39@common_interp
40;
41printf, 40, ''
42printf, 40, 'Entering preproc_mask'
43zero = fltarr(jpiatm+2, jpjatm+2)
44zero(*, *) = 1
45zero[1:jpiatm, 1:jpjatm] = mask
46mask = zero
47;
48mask = mask*shift(mask,  1, 0)*shift(mask, 0,  1)*shift(mask,  1,  1) $
49           *shift(mask, -1, 0)*shift(mask, 0, -1)*shift(mask, -1, -1) $
50           *shift(mask, -1, 1)*shift(mask, 1, -1)
51z = fltarr(jpiatm, jpjatm)
52z = mask[1:jpiatm, 1:jpjatm]
53printf, 40, ''
54printf, 40, 'Extension done'
55;
56return, z
57END
Note: See TracBrowser for help on using the repository browser.