lmdze/trunk/check_coefoz.jnl

! This is a script for Ferret.

! This script compares "coefoz_v2_3.nc" and "coefoz_LMDZ.nc", at a
! given horizontal position of the LMDZ grid and at a given month. Run
! this script from the directory containing "coefoz_LMDZ.nc".

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

use/title="48" coefoz_LMDZ
use "$HOME/Documents/Utilisation_LMDZ/Input_etat0_lim/Ozone/coefoz_v2_3.nc"
!!use "$HOME/Documents/Utilisation_LMDZ/Results_etat0_lim/45/coefoz_LMDZ.nc"

set region/i=25/j=2/l=9 ! choose horizontal position and month
pplus axlsze .15
! First plot coefoz_LMDZ to have longitude in the labels

! Mole fraction:
set window/new
plot/symbol/line/vs/d=1/vlog/vlimits=1e4:.001/size=.2/hlog/hlim=1e-8:1e-4/set_up r, pls
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

! Production:
set window/new
plot/symbol/line/vs/d=1/vlog/vlimits=1e4:.001/size=.2/set_up p_net, pls
pplus ylab "pressure (hPa)"
go unlabel 6
pplus plot

plot/line/symbol/overlay/title="Mobidic regridded in latitude" \
   p_net[d=2, gy=p_net[d=1]@ave]
!!plot/symbol/line/vs/d=3/title="45"/overlay p_net, pls
1	guez	3	! This is a script for Ferret.
2
3			! This script compares "coefoz_v2_3.nc" and "coefoz_LMDZ.nc", at a
4			! given horizontal position of the LMDZ grid and at a given month. Run
5			! this script from the directory containing "coefoz_LMDZ.nc".
6
7			! Reset everything:
8			cancel data_set/all
9			cancel axis/all
10			cancel viewport
11			CANCEL EXPRESSION
12			CANCEL GRID
13			CANCEL LIST/ALL
14			CANCEL MEMORY/all
15			CANCEL SYMBOL/ALL
16			CANCEL REGION/all
17			CANCEL VARIABLE/ALL
18			pplus/reset
19			cancel window/all
20
21			use/title="48" coefoz_LMDZ
22			use "$HOME/Documents/Utilisation_LMDZ/Input_etat0_lim/Ozone/coefoz_v2_3.nc"
23			!!use "$HOME/Documents/Utilisation_LMDZ/Results_etat0_lim/45/coefoz_LMDZ.nc"
24
25			set region/i=25/j=2/l=9 ! choose horizontal position and month
26			pplus axlsze .15
27			! First plot coefoz_LMDZ to have longitude in the labels
28
29			! Mole fraction:
30			set window/new
31			plot/symbol/line/vs/d=1/vlog/vlimits=1e4:.001/size=.2/hlog/hlim=1e-8:1e-4/set_up r, pls
32			pplus ylab "pressure (hPa)"
33			go unlabel 6
34			pplus plot
35
36			plot/line/symbol/overlay/title="Mobidic regridded in latitude" \
37			r[d=2, gy=r[d=1]@ave]
38			!!plot/symbol/line/vs/d=3/title="45"/overlay r, pls
39
40			! Production:
41			set window/new
42			plot/symbol/line/vs/d=1/vlog/vlimits=1e4:.001/size=.2/set_up p_net, pls
43			pplus ylab "pressure (hPa)"
44			go unlabel 6
45			pplus plot
46
47			plot/line/symbol/overlay/title="Mobidic regridded in latitude" \
48			p_net[d=2, gy=p_net[d=1]@ave]
49			!!plot/symbol/line/vs/d=3/title="45"/overlay p_net, pls