1 | ; |
---|
2 | ; |
---|
3 | @common |
---|
4 | ;--------------------------------------------------------- |
---|
5 | ; bornes de la grille |
---|
6 | ;---------------------------------------------------------- |
---|
7 | computegrid, 0.5, -89.5, 1, 1, 360, 180 |
---|
8 | ; |
---|
9 | iodir = '/usr1/antiope/smasson/DATA/' |
---|
10 | ; |
---|
11 | ; |
---|
12 | jpk = 33 |
---|
13 | gdept = ncdf_lec('Levitus98_1m_01_12_S.nc', iodir = iodir, var = 'deptht') |
---|
14 | gdepw = .5*(gdept+shift(gdept, 1)) |
---|
15 | gdepw[0] = 0 |
---|
16 | e3t = shift(gdepw, -1)-gdepw |
---|
17 | e3t[jpk-1] = gdept[jpk-1]-gdepw[jpk-1] |
---|
18 | e3w = gdept-shift(gdept, 1) |
---|
19 | e3w[0] = 0 |
---|
20 | ; |
---|
21 | mer = ncdf_lec('Levitus98_1m_01_12_S.nc', iodir = iodir, var = 'vosaline', count = [jpi, jpj, jpk, 1]) |
---|
22 | mer = where(mer GT -50) |
---|
23 | tmask = bytarr(jpi, jpj, jpk) |
---|
24 | tmask[temporary(mer)] = 1b |
---|
25 | ; |
---|
26 | jpt = 1 |
---|
27 | time = 0 |
---|
28 | ;------------------------------------------------------ |
---|
29 | ; lecture de la grille |
---|
30 | ;------------------------------------------------------ |
---|
31 | key_shift = 0 |
---|
32 | key_stride = [1, 1, 1] |
---|
33 | ;------------------------------------------------------------- |
---|
34 | domdef |
---|
35 | ; |
---|
36 | key_periodique = 1 |
---|
37 | triangles = triangule() |
---|
38 | ;---------------------------------------------------------- |
---|
39 | ; bornes pour la lecture des donnees |
---|
40 | ;---------------------------------------------------------- |
---|
41 | jpidta = jpiglo |
---|
42 | jpjdta = jpjglo |
---|
43 | jpkdta = jpkglo |
---|
44 | ixmindta = 0 |
---|
45 | ixmaxdta = jpidta-1 |
---|
46 | iymindta = 0 |
---|
47 | iymaxdta = jpjdta-1 |
---|
48 | izmindta = 0 |
---|
49 | izmaxdta = jpkdta-1 |
---|
50 | |
---|