1 |
! This is a script for Ferret. |
! This is a script for Ferret. |
2 |
|
|
3 |
! This script compares "coefoz_v2_3.nc" and "coefoz_LMDZ.nc", at a |
! This script compares "coefoz_v2_3.nc" and "coefoz_LMDZ.nc". |
|
! given horizontal position of the LMDZ grid and at a given month. Run |
|
|
! this script from the directory containing "coefoz_LMDZ.nc". |
|
4 |
|
|
5 |
! Reset everything: |
! Reset everything: |
6 |
cancel data_set/all |
cancel data_set/all |
16 |
pplus/reset |
pplus/reset |
17 |
cancel window/all |
cancel window/all |
18 |
|
|
19 |
use/title="48" coefoz_LMDZ |
use coefoz_v2_3.nc |
20 |
use "$HOME/Documents/Utilisation_LMDZ/Input_etat0_lim/Ozone/coefoz_v2_3.nc" |
use coefoz_LMDZ.nc |
21 |
!!use "$HOME/Documents/Utilisation_LMDZ/Results_etat0_lim/45/coefoz_LMDZ.nc" |
|
22 |
|
define axis/z/depth log_plev = log(z[g=p_net[d=COEFOZ_V2_3]]) |
23 |
set region/i=25/j=2/l=9 ! choose horizontal position and month |
define grid/like=p_net_mob/z=log_plev my_grid |
24 |
pplus axlsze .15 |
|
25 |
! First plot coefoz_LMDZ to have longitude in the labels |
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 |
|
|
|
! Mole fraction: |
|
38 |
set window/new |
set window/new |
39 |
plot/symbol/line/vs/d=1/vlog/vlimits=1e4:.001/size=.2/hlog/hlim=1e-8:1e-4/set_up r, pls |
shade/l=76/palette=centered/lev=10c p_net_mob[gz=log_plev@asn] * 1e12 |
|
pplus ylab "pressure (hPa)" |
|
|
go unlabel 6 |
|
|
pplus plot |
|
|
|
|
|
plot/line/symbol/overlay/title="Mobidic regridded in latitude" \ |
|
|
r[d=2, gy=r[d=1]@ave] |
|
|
!!plot/symbol/line/vs/d=3/title="45"/overlay r, pls |
|
40 |
|
|
|
! Production: |
|
41 |
set window/new |
set window/new |
42 |
plot/symbol/line/vs/d=1/vlog/vlimits=1e4:.001/size=.2/set_up p_net, pls |
shade/l=76/lev=1d log(my_thresh[gz=log_plev@asn]) |
43 |
pplus ylab "pressure (hPa)" |
|
44 |
go unlabel 6 |
message |
45 |
pplus plot |
|
46 |
|
set region/z=3 |
47 |
plot/line/symbol/overlay/title="Mobidic regridded in latitude" \ |
palette centered |
48 |
p_net[d=2, gy=p_net[d=1]@ave] |
|
49 |
!!plot/symbol/line/vs/d=3/title="45"/overlay p_net, pls |
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 |