1 | ; |
---|
2 | ; definition of the domain on ORCA2 but with i/j index instead of |
---|
3 | ; longitude/latitude |
---|
4 | ; |
---|
5 | ; deferent possibilities are proposed |
---|
6 | ; |
---|
7 | ; 1) definition with a file |
---|
8 | ; |
---|
9 | file = 'Levitus98_1m_01_12_Temperature_Pot_ORCA05.nc' |
---|
10 | initncdf, file, xminmesh = 1, xmaxmesh = -1 $ |
---|
11 | , yminmesh = 0, ymaxmesh = -1, /xyindex, stride = [2, 1, 1] $ |
---|
12 | , /lookalldir |
---|
13 | ; |
---|
14 | ; 2) definition with the meshmask |
---|
15 | ; |
---|
16 | ; file = 'micromeshmaskORCA05.nc' |
---|
17 | ; initncdf, file, xaxisname = 'glamt', yaxisname = 'gphit' $ |
---|
18 | ; , /xyindex, zaxisname = 'gdept', stride = [2, 1, 1] $ |
---|
19 | ; , yminmesh = 0, ymaxmesh = -1, xminmesh = 1, xmaxmesh = -1 $ |
---|
20 | ; , /lookalldir |
---|
21 | ; |
---|
22 | ; 3) by calling directly computegrid |
---|
23 | ; |
---|
24 | ; if we don't know jpiglo and jpjglo, we call ncdf_meshread with /getdimensions |
---|
25 | ; ncdf_meshread, 'micromeshmaskORCA05.nc', /getdimensions $ |
---|
26 | ; , iodirectory = [iodir, homedir, !path], /recursive |
---|
27 | ;--------------------------------------------------------- |
---|
28 | ; define the grid parameters |
---|
29 | ;--------------------------------------------------------- |
---|
30 | ; computegrid, 0, 0, 1, 1, jpiglo, jpjglo, xminmesh = 1, xmaxmesh = -1 $ |
---|
31 | ; , yminmesh = 0, ymaxmesh = -1, onearth = 0, stride = [2, 1, 1] $ |
---|
32 | ; , /lookalldir |
---|