/[lmdze]/trunk/check_coefoz.jnl
ViewVC logotype

Diff of /trunk/check_coefoz.jnl

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 7 by guez, Mon Mar 31 12:24:17 2008 UTC
# Line 1  Line 1 
1  ! This is a script for Ferret.  ! This is a script for Ferret.
2    
3  ! This script compares "coefoz_v2_3.nc" and "coefoz_LMDZ.nc", at a  ! This script compares "coefoz_v2_3.nc" and "coefoz_LMDZ.nc".
 ! given horizontal position of the LMDZ grid and at a given month. Run  
 ! this script from the directory containing "coefoz_LMDZ.nc".  
4    
5  ! Reset everything:  ! Reset everything:
6  cancel data_set/all  cancel data_set/all
# Line 18  CANCEL VARIABLE/ALL Line 16  CANCEL VARIABLE/ALL
16  pplus/reset  pplus/reset
17  cancel window/all  cancel window/all
18    
19  use/title="48" coefoz_LMDZ  use coefoz_v2_3.nc
20  use "$HOME/Documents/Utilisation_LMDZ/Input_etat0_lim/Ozone/coefoz_v2_3.nc"  use coefoz_LMDZ.nc
21  !!use "$HOME/Documents/Utilisation_LMDZ/Results_etat0_lim/45/coefoz_LMDZ.nc"  
22    define axis/z/depth log_plev = log(z[g=p_net[d=COEFOZ_V2_3]])
23  set region/i=25/j=2/l=9 ! choose horizontal position and month  define grid/like=p_net_mob/z=log_plev my_grid
24  pplus axlsze .15  
25  ! First plot coefoz_LMDZ to have longitude in the labels  let my_var = P_NET_MOB / \
26        P_NET[D=COEFOZ_V2_3, G=P_NET_MOB[D=COEFOZ_LMDZ], GY=@AVE] - 1
27    stat abs(my_var)
28    
29    let MY_THRESH = IF ABS(MY_VAR) GT 0.01 THEN MY_VAR
30    list my_thresh[y=@ngd, z=@ngd, t=@ngd]
31    list sortl(my_thresh[y=@ngd, z=@ngd])
32    ! (gives the value of "l" where the number of points is highest)
33    
34    set window/new
35    shade/l=76/palette=centered/lev=10c \
36       p_net[d=coefoz_v2_3, g=my_grid, gy=@ave, gz=@asn] * 1e12
37    
 ! Mole fraction:  
38  set window/new  set window/new
39  plot/symbol/line/vs/d=1/vlog/vlimits=1e4:.001/size=.2/hlog/hlim=1e-8:1e-4/set_up r, pls  shade/l=76/palette=centered/lev=10c p_net_mob[gz=log_plev@asn] * 1e12
 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  
40    
 ! Production:  
41  set window/new  set window/new
42  plot/symbol/line/vs/d=1/vlog/vlimits=1e4:.001/size=.2/set_up p_net, pls  shade/l=76/lev=1d log(my_thresh[gz=log_plev@asn])
43  pplus ylab "pressure (hPa)"  
44  go unlabel 6  message
45  pplus plot  
46    set region/z=3
47  plot/line/symbol/overlay/title="Mobidic regridded in latitude" \  palette centered
48     p_net[d=2, gy=p_net[d=1]@ave]  
49  !!plot/symbol/line/vs/d=3/title="45"/overlay p_net, pls  set win 1
50    fill/lev=10c P_NET[d=coefoz_v2_3, g=p_net_mob[d=2], gy=@ave] * 1e12
51    
52    set win 2
53    fill/lev=10c P_NET_mob * 1e12
54    
55    set win 3
56    fill/lev=10c P_NET[d=coefoz_v2_3] * 1e12
57    
58    message
59    
60    cancel window 3
61    set axis/modulo time
62    set axis/modulo time1
63    
64    set win 1
65    fill/lev=10c/l=4:17 P_NET[d=coefoz_v2_3] * 1e12
66    
67    set win 2
68    fill/lev=10c/l=120:510 P_NET_mob * 1e12
Legend:
Removed from v.3  
changed lines
  Added in v.7
  ViewVC Help
Powered by ViewVC 1.1.21