1 |
guez |
4 |
! This is a script for Ferret. |
2 |
|
|
|
3 |
|
|
! This script visualizes data from the file "limit.nc" or |
4 |
|
|
! "startphy.nc". |
5 |
|
|
! It reshapes NetCDF variables from the "physics" grid to the "dynamics" grid. |
6 |
|
|
|
7 |
|
|
! Reset everything: |
8 |
|
|
cancel data_set/all |
9 |
|
|
CANCEL GRID |
10 |
|
|
cancel viewport |
11 |
|
|
CANCEL EXPRESSION |
12 |
|
|
CANCEL LIST/ALL |
13 |
|
|
CANCEL MEMORY/all |
14 |
|
|
CANCEL SYMBOL/ALL |
15 |
|
|
CANCEL REGION/all |
16 |
|
|
CANCEL VARIABLE/ALL |
17 |
|
|
cancel axis/all |
18 |
|
|
pplus/reset |
19 |
|
|
cancel window/all |
20 |
|
|
set window/new |
21 |
|
|
|
22 |
|
|
use limit |
23 |
|
|
let iim = 32 |
24 |
|
|
let klon = `FOCE, RETURN=IEND` |
25 |
|
|
let jjm = (klon - 2) / iim + 1 |
26 |
|
|
define axis/x=1:`iim`:1 long_index |
27 |
|
|
define axis/y=1:`jjm-1`:1 lat_index |
28 |
|
|
let dyn_grid= x[gx=long_index] + y[gy=lat_index] |
29 |
|
|
use startphy |
30 |
|
|
shade reshape(sst[i=2:`klon-1`, d=1], dyn_grid), \ |
31 |
|
|
reshape(longitude[i=2:`klon-1`], dyn_grid), \ |
32 |
|
|
reshape(latitude[i=2:`klon-1`], dyn_grid) |
33 |
|
|
go land |