| 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 | |
|---|
