source: LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/ini_histmthNMC.h @ 1844

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