[131] | 1 | ;+ |
---|
[232] | 2 | ; |
---|
[58] | 3 | ; definition of the domain on ORCA2 but with i/j index instead of |
---|
| 4 | ; longitude/latitude |
---|
| 5 | ; |
---|
[131] | 6 | ; different possibilities are proposed |
---|
[237] | 7 | ; |
---|
[131] | 8 | ; @history |
---|
| 9 | ; Sebastien Masson |
---|
[58] | 10 | ; |
---|
[237] | 11 | ; @version |
---|
| 12 | ; $Id$ |
---|
| 13 | ; |
---|
[325] | 14 | ; @todo |
---|
| 15 | ; find the right way to avoid |
---|
| 16 | ; IDLDOC: unknown tag "file_comments" in file tst_initorca05_index.pro |
---|
| 17 | ; |
---|
[131] | 18 | ;- |
---|
| 19 | ; |
---|
[58] | 20 | ; 1) definition with a file |
---|
| 21 | ; |
---|
| 22 | file = 'Levitus98_1m_01_12_Temperature_Pot_ORCA05.nc' |
---|
| 23 | initncdf, file, xminmesh = 1, xmaxmesh = -1 $ |
---|
| 24 | , yminmesh = 0, ymaxmesh = -1, /xyindex $ |
---|
[77] | 25 | , /lookalldir |
---|
[58] | 26 | ; |
---|
| 27 | ; 2) definition with the meshmask |
---|
| 28 | ; |
---|
| 29 | ; file = 'micromeshmaskORCA05.nc' |
---|
| 30 | ; initncdf, file, xaxisname = 'glamt', yaxisname = 'gphit' $ |
---|
| 31 | ; , /xyindex, zaxisname = 'gdept' $ |
---|
| 32 | ; , yminmesh = 0, ymaxmesh = -1, xminmesh = 1, xmaxmesh = -1 $ |
---|
[77] | 33 | ; , /lookalldir |
---|
[58] | 34 | ; |
---|
| 35 | ; 3) by calling directly computegrid |
---|
| 36 | ; |
---|
| 37 | ; if we don't know jpiglo and jpjglo, we call ncdf_meshread with /getdimensions |
---|
| 38 | ; ncdf_meshread, 'micromeshmaskORCA05.nc', /getdimensions $ |
---|
| 39 | ; , iodirectory = [iodir, homedir, !path], /recursive |
---|
| 40 | ;--------------------------------------------------------- |
---|
| 41 | ; define the grid parameters |
---|
| 42 | ;--------------------------------------------------------- |
---|
| 43 | ; computegrid, 0, 0, 1, 1, jpiglo, jpjglo, xminmesh = 1, xmaxmesh = -1 $ |
---|
| 44 | ; , yminmesh = 0, ymaxmesh = -1, onearth = 0 $ |
---|
[77] | 45 | ; , /lookalldir |
---|