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	guez	11	! 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]