source: LMDZ4/trunk/libf/phylmd/ini_histmthNMC.h @ 1352

Last change on this file since 1352 was 1352, checked in by musat, 14 years ago

Add 3 output files for standard pressure levels AR5 exercice and flags
to manage their computation and output frequencies
histhfNMC.nc with 3 standard pressure levels
histdayNMC.nc with 8 (or may have 17) standard pressure levels
histmthNMC.nc with 17 standard pressure levels
Add 3 flags in the physiq.def file: freq_calNMC(3), freq_outNMC(3) and lev_histdayNMC
freq_calNMC(3) : computation frequency of variables on standard pressure levels

and by default has fallowing values (in fact physics' time step dtime)

freq_calNMC(1)=900.
freq_calNMC(2)=900.
freq_calNMC(3)=900.
freq_outNMC(3) : output frequency of variables on standard pressure levels

with following default values

freq_out(1) = 2592000. (30 days)
freq_out(2) = 86400. (1 day)
freq_out(3) = 21600. (6 hours)
lev_histdayNMC is 8 by default but may be switched to 17 (if we need more levels for a particular run)
IM
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.8 KB
Line 
1!
2! $Header$
3!
4c$OMP MASTER
5       zout = freq_outNMC(1)
6c
7       idayref = day_ref
8       CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
9c
10cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
11cym         DO i = 1, iim
12cym            zx_lon(i,1) = rlon(i+1)
13cym            zx_lon(i,jjmp1) = rlon(i+1)
14cym         ENDDO
15         DO ll=1,klev
16            znivsig(ll)=float(ll)
17         ENDDO
18cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
19cym         CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
20cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
21cym     .                 nhori, nid_mthnmc)
22
23         CALL histbeg_phy("histmthNMC",itau_phy, zjulian, dtime, 
24     .                 nhori, nid_mthnmc)
25c
26         CALL histvert(nid_mthnmc, "plev", "pressure", "Pa",
27     .                 nlevSTD, rlevSTD, nvert)
28ccc
29ccc Champs 3D interpolles sur des niveaux de pression du NMC
30ccc
31c
32c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh
33c             sur les niveaux du NMC et on somme & moyenne
34c             toutes les freq_moyNMC secondes par des routines undefSTD et
35c             moy_undefSTD pour eliminer les valeurs "undef"
36c             de la moyenne mensuelle
37c ======> le "inst(X)" ci-dessous est par consequence factice !
38c
39c
40          CALL histdef(nid_mthnmc, "tnondef",
41     .                 "Valeurs non-definies","-",
42     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
43     .                "inst(X)", zout,zout)
44c
45          CALL histdef(nid_mthnmc, "ta",
46     .                 "Air temperature","K",
47     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
48     .                "inst(X)", zout,zout)
49c
50         CALL histdef(nid_mthnmc, "zg",
51     .                "Geopotential height", "m",
52     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
53     .                "inst(X)", zout,zout)
54c
55          CALL histdef(nid_mthnmc, "hus",
56     .                 "Specific humidity","1",
57     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
58     .                "inst(X)", zout,zout)
59c
60         CALL histdef(nid_mthnmc, "hur",
61     .                 "Relative humidity", "%",
62     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
63     .                "inst(X)", zout,zout)
64c
65          CALL histdef(nid_mthnmc, "ua",
66     .                 "Eastward wind","m s-1",
67     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
68     .                "inst(X)", zout,zout)
69c
70          CALL histdef(nid_mthnmc, "va",
71     .                 "Northward wind","m s-1",
72     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
73     .                "inst(X)", zout,zout)
74c
75          CALL histdef(nid_mthnmc, "wap",
76     .                 "Lagrangian tendency of air pressure","Pa s-1",
77     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
78     .                "inst(X)", zout,zout)
79c
80          CALL histdef(nid_mthnmc, "psbg",
81     .         "Pressure sfce below ground","%",
82     .         iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
83     .         "inst(X)", zout,zout)
84c
85          CALL histdef(nid_mthnmc, "uv",
86     .         "uv ",
87     .         "m2/s2",iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
88     .         "inst(X)", zout,zout)
89c
90          CALL histdef(nid_mthnmc, "vq",
91     .         "vq ",
92     .         "m/s * (kg/kg)",iim,jj_nb,nhori, 
93     .          nlevSTD,1,nlevSTD, nvert, 32,
94     .         "inst(X)", zout,zout)
95c
96          CALL histdef(nid_mthnmc, "vT",
97     .         "vT ", 
98     .         "mK/s",iim,jj_nb,nhori, 
99     .          nlevSTD,1,nlevSTD, nvert, 32,
100     .         "inst(X)", zout,zout)
101c
102          CALL histdef(nid_mthnmc, "wq",
103     .         "wq ", 
104     .         "(Pa/s)*(kg/kg)",iim,jj_nb,nhori,
105     .          nlevSTD,1,nlevSTD, nvert, 32,
106     .         "inst(X)", zout,zout)
107c
108          CALL histdef(nid_mthnmc, "vphi",
109     .         "vphi ", 
110     .         "m2/s",iim,jj_nb,nhori, 
111     .          nlevSTD,1,nlevSTD, nvert, 32,
112     .         "inst(X)", zout,zout)
113c
114          CALL histdef(nid_mthnmc, "wT",
115     .         "wT ", 
116     .         "K*Pa/s",iim,jj_nb,nhori,
117     .          nlevSTD,1,nlevSTD, nvert, 32,
118     .         "inst(X)", zout,zout)
119c
120          CALL histdef(nid_mthnmc, "uxu",
121     .         "u2 ", 
122     .         "m2/s2",iim,jj_nb,nhori,
123     .          nlevSTD,1,nlevSTD, nvert, 32,
124     .         "inst(X)", zout,zout)
125c
126          CALL histdef(nid_mthnmc, "vxv",
127     .         "v2 ", 
128     .         "m2/s2",iim,jj_nb,nhori,
129     .          nlevSTD,1,nlevSTD, nvert, 32,
130     .         "inst(X)", zout,zout)
131c
132          CALL histdef(nid_mthnmc, "TxT",
133     .         "T2 ", 
134     .         "K2",iim,jj_nb,nhori,
135     .          nlevSTD,1,nlevSTD, nvert, 32,
136     .         "inst(X)", zout,zout)
137c
138
139         CALL histend(nid_mthnmc)
140c
141c$OMP END MASTER
Note: See TracBrowser for help on using the repository browser.