[2] | 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 | |
---|