source: LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histmth.h @ 450

Last change on this file since 450 was 450, checked in by lmdzadmin, 21 years ago

Un ENDIF de trop IM
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 21.0 KB
Line 
1      IF (ok_mensuel) THEN
2c
3      ndex2d = 0
4      ndex3d = 0
5c
6c Champs 2D:
7c
8         zsto = dtime
9         zout = dtime * ecrit_mth
10         itau_w = itau_phy + itap
11
12      i = NINT(zout/zsto)
13      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
14      CALL histwrite(nid_mth,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
15C
16      i = NINT(zout/zsto)
17      CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
18      CALL histwrite(nid_mth,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
19
20      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
21      CALL histwrite(nid_mth,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
22c
23cccIM
24      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
25      CALL histwrite(nid_mth,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
26c
27      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
28      CALL histwrite(nid_mth,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
29c
30      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
31      CALL histwrite(nid_mth,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
32C
33      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
34      CALL histwrite(nid_mth,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
35C
36      DO i = 1, klon
37         zx_tmp_fi2d(i) = paprs(i,1)
38      ENDDO
39      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
40      CALL histwrite(nid_mth,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
41c
42      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
43      CALL histwrite(nid_mth,"qsurf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
44c
45      if (.not. ok_veget) then
46          CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
47          CALL histwrite(nid_mth,"qsol",itau_w,zx_tmp_2d,iim*jjmp1
48     $        ,ndex2d)
49      endif
50c
51      DO i = 1, klon
52         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
53      ENDDO
54      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
55      CALL histwrite(nid_mth,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
56c
57      DO i = 1, klon
58         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
59      ENDDO
60      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
61      CALL histwrite(nid_mth,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
62c
63      DO i = 1, klon
64         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
65      ENDDO
66      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
67      CALL histwrite(nid_mth,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
68c
69      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
70      CALL histwrite(nid_mth,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
71c
72      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
73      CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
74     .               ndex2d)
75c
76      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
77      CALL histwrite(nid_mth,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
78c
79      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
80      CALL histwrite(nid_mth,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
81c
82      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
83      CALL histwrite(nid_mth,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
84c
85      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
86      CALL histwrite(nid_mth,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
87c
88      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
89      CALL histwrite(nid_mth,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
90c
91      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
92      CALL histwrite(nid_mth,"solldown",itau_w,zx_tmp_2d,iim*jjmp1,
93     .               ndex2d)
94c
95      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
96      CALL histwrite(nid_mth,"tops0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
97c
98      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
99      CALL histwrite(nid_mth,"topl0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
100c
101      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
102      CALL histwrite(nid_mth,"sols0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
103c
104      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
105      CALL histwrite(nid_mth,"soll0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
106c
107      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
108      CALL histwrite(nid_mth,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
109c
110      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
111      CALL histwrite(nid_mth,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
112c
113      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
114      CALL histwrite(nid_mth,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
115c
116c
117c      DO i = 1, klon
118c         zx_tmp_fi2d(i) = fluxu(i,1)
119c      ENDDO
120c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
121c      CALL histwrite(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
122c
123c      DO i = 1, klon
124c         zx_tmp_fi2d(i) = fluxv(i,1)
125c      ENDDO
126c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
127c      CALL histwrite(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
128c
129      DO nsrf = 1, nbsrf
130CYYY
131        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
132        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
133        CALL histwrite(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
134     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
135C
136        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
137        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
138        CALL histwrite(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
139     $      zx_tmp_2d,iim*jjmp1,ndex2d)
140C
141        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
142        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
143        CALL histwrite(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
144     $      zx_tmp_2d,iim*jjmp1,ndex2d)
145C
146        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
147        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
148        CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
149     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
150C
151        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
152        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
153        CALL histwrite(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
154     $      zx_tmp_2d,iim*jjmp1,ndex2d)
155cIM cf. JLD
156        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
157        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
158        CALL histwrite(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
159     $      zx_tmp_2d,iim*jjmp1,ndex2d)
160C
161        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
162        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
163        CALL histwrite(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
164     $      zx_tmp_2d,iim*jjmp1,ndex2d)
165C
166        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
167        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
168        CALL histwrite(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
169     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
170C     
171        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
172        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
173        CALL histwrite(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
174     $      zx_tmp_2d,iim*jjmp1,ndex2d)
175C
176        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
177        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
178        CALL histwrite(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
179     $      zx_tmp_2d,iim*jjmp1,ndex2d)
180C
181        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
182        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
183        CALL histwrite(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
184     $      zx_tmp_2d,iim*jjmp1,ndex2d)
185c
186      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
187      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
188      CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),itau_w
189     $    ,zx_tmp_2d,iim*jjmp1,ndex2d)
190
191      END DO 
192cXXX      DO i = 1, klon
193cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
194cXXX      ENDDO
195cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
196cXXX      CALL histwrite(nid_mth,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
197c
198      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
199      CALL histwrite(nid_mth,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
200      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
201      CALL histwrite(nid_mth,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
202c
203      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
204      CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
205c
206      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
207      CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
208c
209      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
210      CALL histwrite(nid_mth,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
211c
212      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
213      CALL histwrite(nid_mth,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
214c
215      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
216      CALL histwrite(nid_mth,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
217c
218      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
219      CALL histwrite(nid_mth,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
220c
221      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
222      CALL histwrite(nid_mth,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
223c
224      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
225      CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
226c
227      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
228      CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
229c
230      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
231      CALL histwrite(nid_mth,"uq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
232c
233      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
234      CALL histwrite(nid_mth,"vq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
235cKE43
236      IF (iflag_con .GE. 3) THEN ! sb
237c
238      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
239      CALL histwrite(nid_mth,"cape",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
240c
241      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
242      CALL histwrite(nid_mth,"pbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
243c
244      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
245      CALL histwrite(nid_mth,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
246c
247      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
248      CALL histwrite(nid_mth,"fbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
249c
250c
251      ENDIF
252c34EK
253c
254c Champs 3D:
255C
256      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
257      CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
258     .                                   iim*jjmp1*klev,ndex3d)
259c
260      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
261      CALL histwrite(nid_mth,"ovap",itau_w,zx_tmp_3d,
262     .                                   iim*jjmp1*klev,ndex3d)
263c
264      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
265      CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
266     .                                   iim*jjmp1*klev,ndex3d)
267c
268      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
269      CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
270     .                                   iim*jjmp1*klev,ndex3d)
271c
272      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
273      CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
274     .                                   iim*jjmp1*klev,ndex3d)
275c
276      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
277      CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
278     .                                   iim*jjmp1*klev,ndex3d)
279c
280      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
281      CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
282     .                                   iim*jjmp1*klev,ndex3d)
283c
284      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
285      CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d,
286     .                                   iim*jjmp1*klev,ndex3d)
287c
288      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
289      CALL histwrite(nid_mth,"rnebcon",itau_w,zx_tmp_3d,
290     .                                   iim*jjmp1*klev,ndex3d)
291c
292      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
293      CALL histwrite(nid_mth,"rhum",itau_w,zx_tmp_3d,
294     .                                   iim*jjmp1*klev,ndex3d)
295c
296      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d)
297      CALL histwrite(nid_mth,"oliq",itau_w,zx_tmp_3d,
298     .                                   iim*jjmp1*klev,ndex3d)
299c
300      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
301      CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d,
302     .                                   iim*jjmp1*klev,ndex3d)
303c
304      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
305      CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d,
306     .                                   iim*jjmp1*klev,ndex3d)
307c
308      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
309      CALL histwrite(nid_mth,"dqdyn",itau_w,zx_tmp_3d,
310     .                                   iim*jjmp1*klev,ndex3d)
311c
312      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_con, zx_tmp_3d)
313      CALL histwrite(nid_mth,"dtcon",itau_w,zx_tmp_3d,
314     .                                   iim*jjmp1*klev,ndex3d)
315c
316      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d)
317      CALL histwrite(nid_mth,"dqcon",itau_w,zx_tmp_3d,
318     .                                   iim*jjmp1*klev,ndex3d)
319c
320      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_lsc, zx_tmp_3d)
321      CALL histwrite(nid_mth,"dtlsc",itau_w,zx_tmp_3d,
322     .                                   iim*jjmp1*klev,ndex3d)
323c
324      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d)
325      CALL histwrite(nid_mth,"dqlsc",itau_w,zx_tmp_3d,
326     .                                   iim*jjmp1*klev,ndex3d)
327c
328      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d)
329      CALL histwrite(nid_mth,"dtvdf",itau_w,zx_tmp_3d,
330     .                                   iim*jjmp1*klev,ndex3d)
331c
332      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
333      CALL histwrite(nid_mth,"dqvdf",itau_w,zx_tmp_3d,
334     .                                   iim*jjmp1*klev,ndex3d)
335c
336      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_eva, zx_tmp_3d)
337      CALL histwrite(nid_mth,"dteva",itau_w,zx_tmp_3d,
338     .                                   iim*jjmp1*klev,ndex3d)
339c
340      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d)
341      CALL histwrite(nid_mth,"dqeva",itau_w,zx_tmp_3d,
342     .                                   iim*jjmp1*klev,ndex3d)
343c
344      zpt_conv = 0.
345      where (ptconv) zpt_conv = 1.
346      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
347      CALL histwrite(nid_mth,"ptconv",itau_w,zx_tmp_3d,
348     .                                   iim*(jjmp1)*klev,ndex3d)
349c
350      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
351      CALL histwrite(nid_mth,"ratqs",itau_w,zx_tmp_3d,
352     .                                   iim*(jjmp1)*klev,ndex3d)
353c
354      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d)
355      CALL histwrite(nid_mth,"dtajs",itau_w,zx_tmp_3d,
356     .                                   iim*jjmp1*klev,ndex3d)
357c
358      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d)
359      CALL histwrite(nid_mth,"dqajs",itau_w,zx_tmp_3d,
360     .                                   iim*jjmp1*klev,ndex3d)
361c
362      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat, zx_tmp_3d)
363      CALL histwrite(nid_mth,"dtswr",itau_w,zx_tmp_3d,
364     .                                   iim*jjmp1*klev,ndex3d)
365c
366      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat0, zx_tmp_3d)
367      CALL histwrite(nid_mth,"dtsw0",itau_w,zx_tmp_3d,
368     .                                   iim*jjmp1*klev,ndex3d)
369c
370      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool, zx_tmp_3d)
371      CALL histwrite(nid_mth,"dtlwr",itau_w,zx_tmp_3d,
372     .                                   iim*jjmp1*klev,ndex3d)
373c
374      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool0, zx_tmp_3d)
375      CALL histwrite(nid_mth,"dtlw0",itau_w,zx_tmp_3d,
376     .                                   iim*jjmp1*klev,ndex3d)
377c
378      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ec, zx_tmp_3d)
379      CALL histwrite(nid_mth,"dtec",itau_w,zx_tmp_3d,
380     .                                   iim*jjmp1*klev,ndex3d)
381c
382      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
383      CALL histwrite(nid_mth,"duvdf",itau_w,zx_tmp_3d,
384     .                                   iim*jjmp1*klev,ndex3d)
385c
386      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
387      CALL histwrite(nid_mth,"dvvdf",itau_w,zx_tmp_3d,
388     .                                   iim*jjmp1*klev,ndex3d)
389c
390      IF (ok_orodr) THEN
391      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d)
392      CALL histwrite(nid_mth,"duoro",itau_w,zx_tmp_3d,
393     .                                   iim*jjmp1*klev,ndex3d)
394c
395      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oro, zx_tmp_3d)
396      CALL histwrite(nid_mth,"dvoro",itau_w,zx_tmp_3d,
397     .                                   iim*jjmp1*klev,ndex3d)
398c
399      ENDIF
400C
401      IF (ok_orolf) THEN
402      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_lif, zx_tmp_3d)
403      CALL histwrite(nid_mth,"dulif",itau_w,zx_tmp_3d,
404     .                                   iim*jjmp1*klev,ndex3d)
405c
406      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_lif, zx_tmp_3d)
407      CALL histwrite(nid_mth,"dvlif",itau_w,zx_tmp_3d,
408     .                                   iim*jjmp1*klev,ndex3d)
409      ENDIF
410C
411      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
412      CALL histwrite(nid_mth,"ozone",itau_w,zx_tmp_3d,
413     .                                   iim*jjmp1*klev,ndex3d)
414c
415      IF (nqmax.GE.3) THEN
416      DO iq=1,nqmax-2
417      IF (iq.LE.99) THEN
418         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d)
419         WRITE(str2,'(i2.2)') iq
420         CALL histwrite(nid_mth,"trac"//str2,itau_w,zx_tmp_3d,
421     .                                   iim*jjmp1*klev,ndex3d)
422      ELSE
423         PRINT*, "Trop de traceurs"
424         CALL abort
425      ENDIF
426      ENDDO
427      ENDIF
428cKE43
429      IF (iflag_con.GE.3) THEN ! (sb)
430c
431      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
432      CALL histwrite(nid_mth,"upwd",itau_w,zx_tmp_3d,
433     .                                   iim*jjmp1*klev,ndex3d)
434c
435      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
436      CALL histwrite(nid_mth,"dnwd",itau_w,zx_tmp_3d,
437     .                                   iim*jjmp1*klev,ndex3d)
438c
439      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
440      CALL histwrite(nid_mth,"dnwd0",itau_w,zx_tmp_3d,
441     .                                   iim*jjmp1*klev,ndex3d)
442c
443      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
444      CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
445     .                                   iim*jjmp1*klev,ndex3d)
446c
447cccIM
448      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
449      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
450      CALL histwrite(nid_mth, "SWupTOA",itau_w,zx_tmp_2d,
451     .                               iim*jjmp1,ndex2d)
452c
453      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
454      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
455      CALL histwrite(nid_mth, "SWupSFC",itau_w,zx_tmp_2d,
456     .                               iim*jjmp1,ndex2d)
457c
458      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
459      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
460      CALL histwrite(nid_mth, "SWdnTOA",itau_w,zx_tmp_2d,
461     .                               iim*jjmp1,ndex2d)
462c
463      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
464      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
465      CALL histwrite(nid_mth, "SWdnSFC",itau_w,zx_tmp_2d,
466     .                               iim*jjmp1,ndex2d)
467c
468cccIM clear sky
469      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, klevp1)
470      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
471      CALL histwrite(nid_mth, "SWupTOAclr",itau_w,zx_tmp_2d,
472     .                               iim*jjmp1,ndex2d)
473c
474      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
475      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
476      CALL histwrite(nid_mth, "SWupSFCclr",itau_w,zx_tmp_2d,
477     .                               iim*jjmp1,ndex2d)
478c
479      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, klevp1)
480      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
481      CALL histwrite(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_2d,
482     .                               iim*jjmp1,ndex2d)
483c
484      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
485      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
486      CALL histwrite(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_2d,
487     .                               iim*jjmp1,ndex2d)
488c
489      ENDIF
490c
491c  Champs interpolles sur des niveaux de pression
492
493      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
494      CALL histwrite(nid_mth,"u850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
495
496      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
497      CALL histwrite(nid_mth,"v850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
498
499      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
500      CALL histwrite(nid_mth,"u500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
501
502      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
503      CALL histwrite(nid_mth,"v500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
504
505      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
506      CALL histwrite(nid_mth,"u200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
507
508      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
509      CALL histwrite(nid_mth,"v200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
510
511      CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
512      CALL histwrite(nid_mth,"phi500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
513
514      CALL gr_fi_ecrit(1, klon,iim,jjmp1, w500,zx_tmp_2d)
515      CALL histwrite(nid_mth,"w500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
516
517
518     
519
520      if (ok_sync) then
521        call histsync(nid_mth)
522      endif
523       ENDIF
524
Note: See TracBrowser for help on using the repository browser.