1 |
! 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 |