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