1 |
! This is a script for Ferret. |
2 |
|
3 |
! This script compares "coefoz_v2_3.nc" and "coefoz_LMDZ.nc". |
4 |
|
5 |
! Reset everything: |
6 |
cancel data_set/all |
7 |
cancel axis/all |
8 |
cancel viewport |
9 |
CANCEL EXPRESSION |
10 |
CANCEL GRID |
11 |
CANCEL LIST/ALL |
12 |
CANCEL MEMORY/all |
13 |
CANCEL SYMBOL/ALL |
14 |
CANCEL REGION/all |
15 |
CANCEL VARIABLE/ALL |
16 |
pplus/reset |
17 |
cancel window/all |
18 |
|
19 |
use coefoz_v2_3.nc |
20 |
use coefoz_LMDZ.nc |
21 |
|
22 |
define axis/z/depth log_plev = log(z[g=p_net[d=COEFOZ_V2_3]]) |
23 |
define grid/like=p_net_mob/z=log_plev my_grid |
24 |
|
25 |
let my_var = P_NET_MOB / \ |
26 |
P_NET[D=COEFOZ_V2_3, G=P_NET_MOB[D=COEFOZ_LMDZ], GY=@AVE] - 1 |
27 |
stat abs(my_var) |
28 |
|
29 |
let MY_THRESH = IF ABS(MY_VAR) GT 0.01 THEN MY_VAR |
30 |
list my_thresh[y=@ngd, z=@ngd, t=@ngd] |
31 |
list sortl(my_thresh[y=@ngd, z=@ngd]) |
32 |
! (gives the value of "l" where the number of points is highest) |
33 |
|
34 |
set window/new |
35 |
shade/l=76/palette=centered/lev=10c \ |
36 |
p_net[d=coefoz_v2_3, g=my_grid, gy=@ave, gz=@asn] * 1e12 |
37 |
|
38 |
set window/new |
39 |
shade/l=76/palette=centered/lev=10c p_net_mob[gz=log_plev@asn] * 1e12 |
40 |
|
41 |
set window/new |
42 |
shade/l=76/lev=1d log(my_thresh[gz=log_plev@asn]) |
43 |
|
44 |
message |
45 |
|
46 |
set region/z=3 |
47 |
palette centered |
48 |
|
49 |
set win 1 |
50 |
fill/lev=10c P_NET[d=coefoz_v2_3, g=p_net_mob[d=2], gy=@ave] * 1e12 |
51 |
|
52 |
set win 2 |
53 |
fill/lev=10c P_NET_mob * 1e12 |
54 |
|
55 |
set win 3 |
56 |
fill/lev=10c P_NET[d=coefoz_v2_3] * 1e12 |
57 |
|
58 |
message |
59 |
|
60 |
cancel window 3 |
61 |
set axis/modulo time |
62 |
set axis/modulo time1 |
63 |
|
64 |
set win 1 |
65 |
fill/lev=10c/l=4:17 P_NET[d=coefoz_v2_3] * 1e12 |
66 |
|
67 |
set win 2 |
68 |
fill/lev=10c/l=120:510 P_NET_mob * 1e12 |