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

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

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