Parent Directory | Revision Log
This revision is not in working order. Pending some moving of files. Important changes. In the program "etat0_lim": ozone coefficients from Mobidic are regridded in time instead of pressure ; consequences in "etat0". In the program "gcm", ozone coefficients from Mobidic are read once per day only for the current day and regridded in pressure ; consequences in "o3_chem_m", "regr_pr_coefoz", "phytrac" and "regr_pr_comb_coefoz_m". NetCDF95 is a library and does not export NetCDF. New variables "nag_gl_options", "nag_fcalls_options" and "nag_cross_options" in "nag_tools.mk". "check_coefoz.jnl" rewritten entirely for new version of "coefoz_LMDZ.nc". Target "obj_etat0_lim" moved from "GNUmakefile" to "nag_rules.mk". Added some "intent" attributes in "calfis", "clmain", "clqh", "cltrac", "cltracrn", "cvltr", "ini_undefSTD", "moy_undefSTD", "nflxtr", "phystokenc", "phytrac", "readsulfate", "readsulfate_preind" and "undefSTD". In "dynetat0", "dynredem0" and "gcm", "phis" has rank 2 instead of 1. "phis" has assumed shape in "dynredem0". Added module containing "dynredem0". Changed some calls with NetCDF Fortran 77 interface to calls with NetCDF95 interface. Replaced calls to "ssum" by calls to "sum" in "inigeom". In "make.sh", new option "-c" to change compiler. In "aaam_bud", argument "rjour" deleted. In "physiq": renamed some variables; deleted variable "xjour". In "phytrac": renamed some variables; new argument "lmt_pas".
1 | ! This is a script for Ferret. |
2 | |
3 | ! This script compares "coefoz_v2_3.nc" and "coefoz_LMDZ.nc". |
4 | |
5 | ! Reset everything: |
6 | cancel data_set/all |
7 | cancel axis/all |
8 | cancel viewport |
9 | CANCEL EXPRESSION |
10 | CANCEL GRID |
11 | CANCEL LIST/ALL |
12 | CANCEL MEMORY/all |
13 | CANCEL SYMBOL/ALL |
14 | CANCEL REGION/all |
15 | CANCEL VARIABLE/ALL |
16 | pplus/reset |
17 | cancel window/all |
18 | |
19 | use coefoz_v2_3.nc |
20 | use coefoz_LMDZ.nc |
21 | |
22 | define axis/z/depth log_plev = log(z[g=p_net[d=COEFOZ_V2_3]]) |
23 | define grid/like=p_net_mob/z=log_plev my_grid |
24 | |
25 | 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 | |
38 | set window/new |
39 | shade/l=76/palette=centered/lev=10c p_net_mob[gz=log_plev@asn] * 1e12 |
40 | |
41 | set window/new |
42 | shade/l=76/lev=1d log(my_thresh[gz=log_plev@asn]) |
43 | |
44 | message |
45 | |
46 | set region/z=3 |
47 | palette centered |
48 | |
49 | 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 |
ViewVC Help | |
Powered by ViewVC 1.1.21 |