| 1 | ; + |
|---|
| 2 | ; @file_comments |
|---|
| 3 | ; initialisation on ORCA2 grid |
|---|
| 4 | ; |
|---|
| 5 | ; @history |
|---|
| 6 | ; F. Vivier 2003 |
|---|
| 7 | ; F. Pinsard adaptation to new SAXO |
|---|
| 8 | ; according to SAXO ./SRC/Tests/tst_initorca2.pro |
|---|
| 9 | ; |
|---|
| 10 | ; @todo |
|---|
| 11 | ; parametrisation for ORCA025 |
|---|
| 12 | ; |
|---|
| 13 | ; zminmesh and izmaxmesh may be set to 0,0 because we only work with |
|---|
| 14 | ; surface problem (no z) whatever the resolution of ORCA. this will safe lot |
|---|
| 15 | ; of memory |
|---|
| 16 | ; |
|---|
| 17 | ; use ncdf_getaxis to recover meshmask dimension |
|---|
| 18 | ; on x, retrieve -1 and -1 (bande de recouvrement haut et bas) |
|---|
| 19 | ; on y, for ORCA2, retrieve -2 (bande de recouvrement ++ et ++) |
|---|
| 20 | ; |
|---|
| 21 | ; reminder : ORCA2 = 182 x 149 and ORCA025 1442 x 1021 |
|---|
| 22 | ; |
|---|
| 23 | ; for ORCA025, user /keepcont in triangule call. |
|---|
| 24 | ; |
|---|
| 25 | ; @version |
|---|
| 26 | ; $Id$ |
|---|
| 27 | ;- |
|---|
| 28 | ; |
|---|
| 29 | @cm_4mesh |
|---|
| 30 | @cm_4cal |
|---|
| 31 | @cm_4data |
|---|
| 32 | ; |
|---|
| 33 | ;--------------------------------------------------------- |
|---|
| 34 | ; boundaries of the read grid regarding to the original grid |
|---|
| 35 | ;---------------------------------------------------------- |
|---|
| 36 | ; |
|---|
| 37 | ixminmesh =1 |
|---|
| 38 | ixmaxmesh =180 |
|---|
| 39 | ; |
|---|
| 40 | iyminmesh =0 |
|---|
| 41 | iymaxmesh =147 |
|---|
| 42 | ; |
|---|
| 43 | izminmesh =-1 |
|---|
| 44 | izmaxmesh =-1 |
|---|
| 45 | ; |
|---|
| 46 | jpt = 1 |
|---|
| 47 | time = 0 |
|---|
| 48 | ;------------------------------------------------------ |
|---|
| 49 | ; read the grid |
|---|
| 50 | ;------------------------------------------------------ |
|---|
| 51 | iodir = '/usr/work/sur/fvi/OPA/ORCA2/' |
|---|
| 52 | |
|---|
| 53 | key_stride = [1, 1, 1] |
|---|
| 54 | ncdf_meshlec, 'meshmask_bab.nc', glamboundary = [20, 380] |
|---|
| 55 | ;------------------------------------------------------------- |
|---|
| 56 | domdef |
|---|
| 57 | ; |
|---|
| 58 | triangles_list = triangule() |
|---|
| 59 | ;---------------------------------------------------------- |
|---|
| 60 | ; boundaries of the data regarding to the original grid |
|---|
| 61 | ;---------------------------------------------------------- |
|---|
| 62 | jpidta = jpiglo |
|---|
| 63 | jpjdta = jpjglo |
|---|
| 64 | jpkdta = jpkglo |
|---|
| 65 | ixmindta = 0 |
|---|
| 66 | ixmaxdta = jpidta-1 |
|---|
| 67 | iymindta = 0 |
|---|
| 68 | iymaxdta = jpjdta-1 |
|---|
| 69 | izmindta = 0 |
|---|
| 70 | izmaxdta = jpkdta-1 |
|---|
| 71 | ;---------------------------------------------------------- |
|---|
| 72 | @updateold |
|---|
