1 | ;------------------------------------------------------------ |
---|
2 | ; MODIFICATION HISTORY: Sebastien Masson (smasson@lodyc.jussieu.fr) |
---|
3 | ; 26/5/98 |
---|
4 | ;------------------------------------------------------------ |
---|
5 | ;---------------------------------------------------------- |
---|
6 | ; bornes de la grille |
---|
7 | ;---------------------------------------------------------- |
---|
8 | ixminmesh = 1 |
---|
9 | ixmaxmesh = 770 |
---|
10 | iyminmesh = 2 |
---|
11 | iymaxmesh = 1294 |
---|
12 | jpk = 1 |
---|
13 | ; |
---|
14 | jpt = 1 |
---|
15 | time = 0 |
---|
16 | ;-------------------------------------------------------- |
---|
17 | ; definition par defaut du fichier que l''on traite |
---|
18 | ;---------------------------------------------------------- |
---|
19 | ;------------------------------------------------------ |
---|
20 | ; |
---|
21 | iodir= '/usr1/antiope/smasson/DATA/' |
---|
22 | prefix='Clipper' |
---|
23 | ; |
---|
24 | ;------------------------------------------------------- |
---|
25 | ; lecture de la grille |
---|
26 | ;------------------------------------------------------ |
---|
27 | .r /usr1/com/smasson/IDL_NICO/read_coordinates |
---|
28 | read_coordinates, iodir,'coordinates.new' |
---|
29 | tmask = replicate(1b, jpi, jpj, jpk) |
---|
30 | gdept = 5 |
---|
31 | gdepw = 0 |
---|
32 | e3t = 10 |
---|
33 | e3w = 10 |
---|
34 | ;------------------------------------------------------------- |
---|
35 | ;irregularmask |
---|
36 | ;------------------------------------------------------------- |
---|
37 | domdef, grille='T' |
---|
38 | ;---------------------------------------------------------- |
---|
39 | ;---------------------------------------------------------- |
---|
40 | .r /usr1/com/smasson/IDL_NICO/read_dimg_body |
---|
41 | .r /usr1/com/smasson/IDL_NICO/read_dimg_head |
---|
42 | read_dimg_body, iodir, 'ATL6-V5_S_y1958m02d14.dimg', sss |
---|
43 | sss = sss[ixminmesh:ixmaxmesh, iyminmesh:iymaxmesh] |
---|
44 | tmask[where(sss EQ 0)] = 0 |
---|
45 | key_periodique = 0 |
---|
46 | triangles = triangule() |
---|
47 | ; |
---|
48 | vargrid = 'T' |
---|
49 | varname = 'sss' |
---|
50 | vardate= '0' |
---|
51 | varexp = 'Clipper' |
---|
52 | varunit = 'psu' |
---|
53 | ; |
---|
54 | ; |
---|
55 | premierxu = premierxt & premierxv = premierxt & premierxf = premierxt |
---|
56 | premieryu = premieryt & premieryv = premieryt & premieryf = premieryt |
---|
57 | dernierxu = dernierxt & dernierxv = dernierxt & dernierxf = dernierxt |
---|
58 | dernieryu = dernieryt & dernieryv = dernieryt & dernieryf = dernieryt |
---|