source: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_bilKP_ave.h @ 715

Last change on this file since 715 was 644, checked in by Laurent Fairhead, 19 years ago

Synchronisation avec tous les diagnostiques de Ionela IM
Inclusion du slab ocean IM
LF

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 8.6 KB
Line 
1c
2c $Header$
3c
4      IF (ok_journe) THEN
5c
6         zsto = dtime
7         zout = dtime * FLOAT(ecrit_day)
8         typeval=tave
9c
10         idayref = day_ref
11         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
12c
13         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
14         DO i = 1, iim
15            zx_lon(i,1) = rlon(i+1)
16            zx_lon(i,jjmp1) = rlon(i+1)
17         ENDDO
18         DO ll=1,klev
19            znivsig(ll)=float(ll)
20         ENDDO
21         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
22         write(*,*)'zx_lon = ',zx_lon(:,1)
23         write(*,*)'zx_lat = ',zx_lat(1,:)
24         CALL histbeg("histbilKP_ave", iim,zx_lon(:,1), jjmp1,
25     .                zx_lat(1,:),
26     .                1,iim,1,jjmp1, itau_phy, zjulian, dtime,
27     .                nhori, nid_bilKPave)
28         write(*,*)'Journee ', itau_phy, zjulian
29         CALL histvert(nid_bilKPave, "presnivs",
30     .                "Vertical levels","mb",
31     .                 klev, presnivs/100., nvert)
32c
33c
34c Champs 3D:
35c
36         CALL histdef(nid_bilKPave,"ue",
37     .   "Zonal energy transport","-",
38     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
39     .                typeval, zsto,zout)
40c
41         CALL histdef(nid_bilKPave,"ve",
42     .   "Merid energy transport","-",
43     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
44     .                typeval, zsto,zout)
45c
46         CALL histdef(nid_bilKPave,"uq",
47     .   "Zonal humidity transport","-",
48     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
49     .                typeval, zsto,zout)
50c
51         CALL histdef(nid_bilKPave,"vq",
52     .   "Merid humidity transport","-",
53     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32, 
54     .                typeval, zsto,zout)
55c
56c Champs 3D:
57c
58         CALL histdef(nid_bilKPave,"temp",
59     .   "Air temperature","K",
60     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
61     .                typeval, zsto,zout)
62c
63         CALL histdef(nid_bilKPave,"ovap",
64     .   "Specific humidity","Kg/Kg",
65     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
66     .                typeval, zsto,zout)
67c
68         CALL histdef(nid_bilKPave,"geop",
69     .   "Geopotential height","m",
70     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
71     .                typeval, zsto,zout)
72c
73         CALL histdef(nid_bilKPave,"vitu",
74     .   "Zonal wind","m/s",
75     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
76     .                typeval, zsto,zout)
77c
78         CALL histdef(nid_bilKPave,"vitv",
79     .   "Meridional wind","m/s",
80     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
81     .                typeval, zsto,zout)
82c
83         CALL histdef(nid_bilKPave, "vitw", 
84     .   "Vertical wind", "m/s",
85     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
86     .                typeval, zsto,zout)
87c
88         CALL histdef(nid_bilKPave, "pres", 
89     .   "Inter-Layer Air pressure",
90     .                "Pa",
91     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
92     .                typeval, zsto,zout)
93c
94         CALL histdef(nid_bilKPave, "play", 
95     .   "Mean-Layer Air pressure",
96     .                "Pa",
97     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
98     .                typeval, zsto,zout)
99c
100         CALL histdef(nid_bilKPave, "oliq", 
101     .   "Liquid water content", 
102     .                "kg/kg",
103     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
104     .                typeval, zsto,zout)
105c
106         CALL histdef(nid_bilKPave, "dtdyn", 
107     .   "Dynamics dT", "K/s",
108     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
109     .                typeval, zsto,zout)
110c
111         CALL histdef(nid_bilKPave, "dqdyn", 
112     .   "Dynamics dQ", "Kg/Kg/s",
113     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
114     .                typeval, zsto,zout)
115c
116         CALL histdef(nid_bilKPave, "dtcon", 
117     .   "Convection dT", "K/s",
118     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
119     .                typeval, zsto,zout)
120c
121         CALL histdef(nid_bilKPave, "ducon", 
122     .   "Convection du", "m/s2",
123     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
124     .                typeval, zsto,zout)
125c
126         CALL histdef(nid_bilKPave, "dvcon", 
127     .   "Convection dv", "m/s2",
128     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
129     .                typeval, zsto,zout)
130c
131         CALL histdef(nid_bilKPave,"dqcon",
132     .   "Convection dQ","Kg/Kg/s",
133     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
134     .                typeval, zsto,zout)
135c
136         CALL histdef(nid_bilKPave, "dtlsc", 
137     .   "Condensation dT", "K/s",
138     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
139     .                typeval, zsto,zout)
140c
141         CALL histdef(nid_bilKPave,"dqlsc",
142     .   "Condensation dQ","Kg/Kg/s",
143     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
144     .                typeval, zsto,zout)
145c
146         CALL histdef(nid_bilKPave,"dtvdf",
147     .   "Boundary-layer dT","K/s",
148     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
149     .                typeval, zsto,zout)
150c
151         CALL histdef(nid_bilKPave, "dqvdf", 
152     .   "Boundary-layer dQ", 
153     .               "Kg/Kg/s",
154     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
155     .                typeval, zsto,zout)
156c
157         CALL histdef(nid_bilKPave,"dtajs",
158     .   "Ajustement sec dT","K/s",
159     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
160     .                typeval, zsto,zout)
161c
162         CALL histdef(nid_bilKPave, "dqajs",
163     .   "Ajustement sec dQ", 
164     .               "Kg/Kg/s",
165     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
166     .                typeval, zsto,zout)
167c
168         CALL histdef(nid_bilKPave, "dteva", 
169     .   "Reevaporation dT", "K/s",
170     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
171     .                typeval, zsto,zout)
172c
173         CALL histdef(nid_bilKPave,"dqeva",
174     .   "Reevaporation dQ",
175     .                "Kg/Kg/s",
176     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
177     .                typeval, zsto,zout)
178
179c
180         CALL histdef(nid_bilKPave, "dtswr", 
181     .   "SW radiation dT", "K/s",
182     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
183     .                typeval, zsto,zout)
184c
185         CALL histdef(nid_bilKPave, "dtsw0", 
186     .   "SW radiation dT", "K/s",
187     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
188     .                typeval, zsto,zout)
189c
190         CALL histdef(nid_bilKPave, "dtlwr", 
191     .   "LW radiation dT", "K/s",
192     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
193     .                typeval, zsto,zout)
194c
195         CALL histdef(nid_bilKPave, "dtlw0", 
196     .   "LW radiation dT", "K/s",
197     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
198     .                typeval, zsto,zout)
199c
200         CALL histdef(nid_bilKPave,"duvdf",
201     .   "Boundary-layer dU","m/s2",
202     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
203     .                typeval, zsto,zout)
204c
205         CALL histdef(nid_bilKPave,"dvvdf",
206     .   "Boundary-layer dV","m/s2",
207     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
208     .                typeval, zsto,zout)
209c
210         IF (ok_orodr) THEN
211         IF (ok_orolf) THEN
212         CALL histdef(nid_bilKPave, "duoli",
213     .   "Orography dU", "m/s2",
214     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
215     .                typeval, zsto,zout)
216c
217         CALL histdef(nid_bilKPave, "dvoli", 
218     .   "Orography dV", "m/s2",
219     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
220     .                typeval, zsto,zout)
221c
222         ENDIF
223         ENDIF
224C
225         CALL histdef(nid_bilKPave, "duphy",
226     .   "Physiq dU","-",
227     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
228     .                typeval, zsto,zout)
229C
230         CALL histdef(nid_bilKPave, "dvphy",
231     .   "Physiq dV","-",
232     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
233     .                typeval, zsto,zout)
234C
235         CALL histdef(nid_bilKPave, "dtphy",
236     .   "Physiq dT","-",
237     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
238     .                typeval, zsto,zout)
239C
240         CALL histdef(nid_bilKPave, "dqphy",
241     .   "Physiq dQ","-",
242     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
243     .                typeval, zsto,zout)
244C
245         CALL histdef(nid_bilKPave, "dqlphy",
246     .   "Physiq dQl","-",
247     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
248     .                typeval, zsto,zout)
249C
250C
251         CALL histend(nid_bilKPave)
252c
253         ndex2d = 0
254         ndex3d = 0
255c
256      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.