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