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	guez	12	! 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