lmdze/trunk/press_field.jnl

! This is a script for Ferret.

! This script computes the pressure field.

! 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 start
let eta = z[gz=1:10:1] ! 1: llm + 1: 1
let/title="pressure"/units="hPa" \
   pressure = (reshape(ap, eta) + reshape(bp, eta) * ps) / 100.

! Now we can see or plot what we want. For example:
list pressure[x=0, y=0]
stat pressure[k=7]
1	! This is a script for Ferret.
2
3	! This script computes the pressure field.
4
5	! Reset everything:
6	cancel data_set/all
7	CANCEL GRID
8	cancel viewport
9	CANCEL EXPRESSION
10	CANCEL LIST/ALL
11	CANCEL MEMORY/all
12	CANCEL SYMBOL/ALL
13	CANCEL REGION/all
14	CANCEL VARIABLE/ALL
15	cancel axis/all
16	pplus/reset
17	cancel window/all
18	set window/new
19
20	use start
21	let eta = z[gz=1:10:1] ! 1: llm + 1: 1
22	let/title="pressure"/units="hPa" \
23	pressure = (reshape(ap, eta) + reshape(bp, eta) * ps) / 100.
24
25	! Now we can see or plot what we want. For example:
26	list pressure[x=0, y=0]
27	stat pressure[k=7]