lmdze/trunk/comp_o3.jnl

! This is a script for Ferret.

! This script compares column-density from "coefoz_LMDZ.c" and
! "histrac.nc", then it compares ozone mole fraction.

! Reset everything:
cancel data_set/all
CANCEL GRID
cancel viewport
CANCEL EXPRESSION
CANCEL LIST/ALL
CANCEL MEMORY/all
CANCEL SYMBOL/ALL
CANCEL REGION/all
CANCEL VARIABLE/ALL
cancel axis/all
pplus/reset
cancel window/all
set window/new

use coefoz_LMDZ.nc
use "/usr/local/guez/histrac.nc"

contour sigma_mob[d=1, k=45, l=1:90] / 2.69e16
let sigma_o3 = o3 * zmasse
contour/over sigma_o3[x=@ave,z=@sum] / 2.14e-5

set window/new
define axis/depth/z log_plev= log(z[d=1, g=r_mob])
define axis/depth/z log_PRESNIVS=log(z[g=o3]) - 2.
set region/l=2
contour r_mob[d=1, gz=log_plev@asn] * 1e6
contour/over o3[x=@ave, gz=log_presnivs@asn] * 29 / 48 * 1e6
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