1 |
! This is a script for Ferret. |
2 |
|
3 |
! This script compares column-density from "coefoz_LMDZ.c" and |
4 |
! "histrac.nc", then it compares ozone mole fraction. |
5 |
|
6 |
! Reset everything: |
7 |
cancel data_set/all |
8 |
CANCEL GRID |
9 |
cancel viewport |
10 |
CANCEL EXPRESSION |
11 |
CANCEL LIST/ALL |
12 |
CANCEL MEMORY/all |
13 |
CANCEL SYMBOL/ALL |
14 |
CANCEL REGION/all |
15 |
CANCEL VARIABLE/ALL |
16 |
cancel axis/all |
17 |
pplus/reset |
18 |
cancel window/all |
19 |
set window/new |
20 |
|
21 |
use coefoz_LMDZ.nc |
22 |
use "/usr/local/guez/histrac.nc" |
23 |
|
24 |
contour sigma_mob[d=1, k=45, l=1:90] / 2.69e16 |
25 |
let sigma_o3 = o3 * zmasse |
26 |
contour/over sigma_o3[x=@ave,z=@sum] / 2.14e-5 |
27 |
|
28 |
set window/new |
29 |
define axis/depth/z log_plev= log(z[d=1, g=r_mob]) |
30 |
define axis/depth/z log_PRESNIVS=log(z[g=o3]) - 2. |
31 |
set region/l=2 |
32 |
contour r_mob[d=1, gz=log_plev@asn] * 1e6 |
33 |
contour/over o3[x=@ave, gz=log_presnivs@asn] * 29 / 48 * 1e6 |