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