--- trunk/Scripts/check_coefoz.jnl	2008/02/27 13:16:39	3
+++ trunk/Scripts/check_coefoz.jnl	2008/03/31 12:24:17	7
@@ -1,8 +1,6 @@
 ! 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".
+! This script compares "coefoz_v2_3.nc" and "coefoz_LMDZ.nc".
 
 ! Reset everything:
 cancel data_set/all
@@ -18,32 +16,53 @@
 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
+use coefoz_v2_3.nc
+use coefoz_LMDZ.nc
+
+define axis/z/depth log_plev = log(z[g=p_net[d=COEFOZ_V2_3]])
+define grid/like=p_net_mob/z=log_plev my_grid
+
+let my_var = P_NET_MOB / \
+    P_NET[D=COEFOZ_V2_3, G=P_NET_MOB[D=COEFOZ_LMDZ], GY=@AVE] - 1
+stat abs(my_var)
+
+let MY_THRESH = IF ABS(MY_VAR) GT 0.01 THEN MY_VAR
+list my_thresh[y=@ngd, z=@ngd, t=@ngd]
+list sortl(my_thresh[y=@ngd, z=@ngd])
+! (gives the value of "l" where the number of points is highest)
+
+set window/new
+shade/l=76/palette=centered/lev=10c \
+   p_net[d=coefoz_v2_3, g=my_grid, gy=@ave, gz=@asn] * 1e12
 
-! 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
+shade/l=76/palette=centered/lev=10c p_net_mob[gz=log_plev@asn] * 1e12
 
-! 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
+shade/l=76/lev=1d log(my_thresh[gz=log_plev@asn])
+
+message
+
+set region/z=3
+palette centered
+
+set win 1
+fill/lev=10c P_NET[d=coefoz_v2_3, g=p_net_mob[d=2], gy=@ave] * 1e12
+
+set win 2
+fill/lev=10c P_NET_mob * 1e12
+
+set win 3
+fill/lev=10c P_NET[d=coefoz_v2_3] * 1e12
+
+message
+
+cancel window 3
+set axis/modulo time
+set axis/modulo time1
+
+set win 1
+fill/lev=10c/l=4:17 P_NET[d=coefoz_v2_3] * 1e12
+
+set win 2
+fill/lev=10c/l=120:510 P_NET_mob * 1e12