trunk/Scripts/phys_grid.jnl

! This is a script for Ferret.

! This script visualizes data from the file "limit.nc" or
! "startphy.nc".
! It reshapes NetCDF variables from the "physics" grid to the "dynamics" grid.

! Reset everything:
cancel data_set/all
CANCEL GRID
cancel viewport
CANCEL EXPRESSION
CANCEL LIST/ALL
CANCEL MEMORY/all
CANCEL SYMBOL/ALL
CANCEL REGION/all
CANCEL VARIABLE/ALL
cancel axis/all
pplus/reset
cancel window/all
set window/new

use limit
let iim = 32
let klon = `FOCE, RETURN=IEND`
let jjm = (klon - 2) / iim + 1
define axis/x=1:`iim`:1 long_index
define axis/y=1:`jjm-1`:1 lat_index
let dyn_grid= x[gx=long_index] + y[gy=lat_index]
use startphy
shade reshape(sst[i=2:`klon-1`, d=1], dyn_grid), \
   reshape(longitude[i=2:`klon-1`], dyn_grid), \
   reshape(latitude[i=2:`klon-1`], dyn_grid)
go land
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