source: LMDZ.3.3/trunk/libf/phylmd/write_histmth.h @ 544

Last change on this file since 544 was 416, checked in by lmdzadmin, 23 years ago

Inclusion initiale
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 18.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, zxqsol,zx_tmp_2d)
43      CALL histwrite(nid_mth,"qsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
44c
45      DO i = 1, klon
46         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
47      ENDDO
48      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
49      CALL histwrite(nid_mth,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
50c
51      DO i = 1, klon
52         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
53      ENDDO
54      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
55      CALL histwrite(nid_mth,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
56c
57      DO i = 1, klon
58         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
59      ENDDO
60      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
61      CALL histwrite(nid_mth,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
62c
63      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
64      CALL histwrite(nid_mth,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
65c
66      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
67      CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
68     .               ndex2d)
69c
70      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
71      CALL histwrite(nid_mth,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
72c
73      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
74      CALL histwrite(nid_mth,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
75c
76      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
77      CALL histwrite(nid_mth,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
78c
79      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
80      CALL histwrite(nid_mth,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
81c
82      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
83      CALL histwrite(nid_mth,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
84c
85      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
86      CALL histwrite(nid_mth,"solldown",itau_w,zx_tmp_2d,iim*jjmp1,
87     .               ndex2d)
88c
89      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
90      CALL histwrite(nid_mth,"tops0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
91c
92      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
93      CALL histwrite(nid_mth,"topl0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
94c
95      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
96      CALL histwrite(nid_mth,"sols0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
97c
98      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
99      CALL histwrite(nid_mth,"soll0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
100c
101      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
102      CALL histwrite(nid_mth,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
103c
104      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
105      CALL histwrite(nid_mth,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
106c
107      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
108      CALL histwrite(nid_mth,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
109c
110c
111c      DO i = 1, klon
112c         zx_tmp_fi2d(i) = fluxu(i,1)
113c      ENDDO
114c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
115c      CALL histwrite(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
116c
117c      DO i = 1, klon
118c         zx_tmp_fi2d(i) = fluxv(i,1)
119c      ENDDO
120c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
121c      CALL histwrite(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
122c
123      DO nsrf = 1, nbsrf
124CYYY
125        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
126        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
127        CALL histwrite(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
128     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
129C
130        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
131        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
132        CALL histwrite(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
133     $      zx_tmp_2d,iim*jjmp1,ndex2d)
134C
135        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
136        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
137        CALL histwrite(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
138     $      zx_tmp_2d,iim*jjmp1,ndex2d)
139C
140        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
141        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
142        CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
143     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
144C
145        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
146        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
147        CALL histwrite(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
148     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
149C     
150        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
151        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
152        CALL histwrite(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
153     $      zx_tmp_2d,iim*jjmp1,ndex2d)
154C
155        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
156        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
157        CALL histwrite(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
158     $      zx_tmp_2d,iim*jjmp1,ndex2d)
159C
160        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
161        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
162        CALL histwrite(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
163     $      zx_tmp_2d,iim*jjmp1,ndex2d)
164c
165      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
166      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
167      CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),itau_w
168     $    ,zx_tmp_2d,iim*jjmp1,ndex2d)
169
170      END DO 
171cXXX      DO i = 1, klon
172cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
173cXXX      ENDDO
174cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
175cXXX      CALL histwrite(nid_mth,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
176c
177      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
178      CALL histwrite(nid_mth,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
179      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
180      CALL histwrite(nid_mth,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
181c
182      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
183      CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
184c
185      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
186      CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
187c
188      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
189      CALL histwrite(nid_mth,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
190c
191      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
192      CALL histwrite(nid_mth,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
193c
194      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
195      CALL histwrite(nid_mth,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
196c
197      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
198      CALL histwrite(nid_mth,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
199c
200      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
201      CALL histwrite(nid_mth,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
202c
203      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
204      CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
205c
206      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
207      CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
208c
209      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
210      CALL histwrite(nid_mth,"uq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
211c
212      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
213      CALL histwrite(nid_mth,"vq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
214cKE43
215      IF (iflag_con .GE. 3) THEN ! sb
216c
217      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
218      CALL histwrite(nid_mth,"cape",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
219c
220      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
221      CALL histwrite(nid_mth,"pbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
222c
223      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
224      CALL histwrite(nid_mth,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
225c
226      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
227      CALL histwrite(nid_mth,"fbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
228c
229c
230      ENDIF
231c34EK
232c
233c Champs 3D:
234C
235      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
236      CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
237     .                                   iim*jjmp1*klev,ndex3d)
238c
239      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
240      CALL histwrite(nid_mth,"ovap",itau_w,zx_tmp_3d,
241     .                                   iim*jjmp1*klev,ndex3d)
242c
243      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
244      CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
245     .                                   iim*jjmp1*klev,ndex3d)
246c
247      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
248      CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
249     .                                   iim*jjmp1*klev,ndex3d)
250c
251      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
252      CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
253     .                                   iim*jjmp1*klev,ndex3d)
254c
255      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
256      CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
257     .                                   iim*jjmp1*klev,ndex3d)
258c
259      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
260      CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
261     .                                   iim*jjmp1*klev,ndex3d)
262c
263      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
264      CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d,
265     .                                   iim*jjmp1*klev,ndex3d)
266c
267      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
268      CALL histwrite(nid_mth,"rhum",itau_w,zx_tmp_3d,
269     .                                   iim*jjmp1*klev,ndex3d)
270c
271      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d)
272      CALL histwrite(nid_mth,"oliq",itau_w,zx_tmp_3d,
273     .                                   iim*jjmp1*klev,ndex3d)
274c
275      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
276      CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d,
277     .                                   iim*jjmp1*klev,ndex3d)
278c
279      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
280      CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d,
281     .                                   iim*jjmp1*klev,ndex3d)
282c
283      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
284      CALL histwrite(nid_mth,"dqdyn",itau_w,zx_tmp_3d,
285     .                                   iim*jjmp1*klev,ndex3d)
286c
287      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_con, zx_tmp_3d)
288      CALL histwrite(nid_mth,"dtcon",itau_w,zx_tmp_3d,
289     .                                   iim*jjmp1*klev,ndex3d)
290c
291      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d)
292      CALL histwrite(nid_mth,"dqcon",itau_w,zx_tmp_3d,
293     .                                   iim*jjmp1*klev,ndex3d)
294c
295      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_lsc, zx_tmp_3d)
296      CALL histwrite(nid_mth,"dtlsc",itau_w,zx_tmp_3d,
297     .                                   iim*jjmp1*klev,ndex3d)
298c
299      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d)
300      CALL histwrite(nid_mth,"dqlsc",itau_w,zx_tmp_3d,
301     .                                   iim*jjmp1*klev,ndex3d)
302c
303      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d)
304      CALL histwrite(nid_mth,"dtvdf",itau_w,zx_tmp_3d,
305     .                                   iim*jjmp1*klev,ndex3d)
306c
307      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
308      CALL histwrite(nid_mth,"dqvdf",itau_w,zx_tmp_3d,
309     .                                   iim*jjmp1*klev,ndex3d)
310c
311      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_eva, zx_tmp_3d)
312      CALL histwrite(nid_mth,"dteva",itau_w,zx_tmp_3d,
313     .                                   iim*jjmp1*klev,ndex3d)
314c
315      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d)
316      CALL histwrite(nid_mth,"dqeva",itau_w,zx_tmp_3d,
317     .                                   iim*jjmp1*klev,ndex3d)
318c
319      zpt_conv = 0.
320      where (ptconv) zpt_conv = 1.
321      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
322      CALL histwrite(nid_mth,"ptconv",itau_w,zx_tmp_3d,
323     .                                   iim*(jjmp1)*klev,ndex3d)
324c
325      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
326      CALL histwrite(nid_mth,"ratqs",itau_w,zx_tmp_3d,
327     .                                   iim*(jjmp1)*klev,ndex3d)
328c
329      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d)
330      CALL histwrite(nid_mth,"dtajs",itau_w,zx_tmp_3d,
331     .                                   iim*jjmp1*klev,ndex3d)
332c
333      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d)
334      CALL histwrite(nid_mth,"dqajs",itau_w,zx_tmp_3d,
335     .                                   iim*jjmp1*klev,ndex3d)
336c
337      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat, zx_tmp_3d)
338      CALL histwrite(nid_mth,"dtswr",itau_w,zx_tmp_3d,
339     .                                   iim*jjmp1*klev,ndex3d)
340c
341      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat0, zx_tmp_3d)
342      CALL histwrite(nid_mth,"dtsw0",itau_w,zx_tmp_3d,
343     .                                   iim*jjmp1*klev,ndex3d)
344c
345      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool, zx_tmp_3d)
346      CALL histwrite(nid_mth,"dtlwr",itau_w,zx_tmp_3d,
347     .                                   iim*jjmp1*klev,ndex3d)
348c
349      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool0, zx_tmp_3d)
350      CALL histwrite(nid_mth,"dtlw0",itau_w,zx_tmp_3d,
351     .                                   iim*jjmp1*klev,ndex3d)
352c
353      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ec, zx_tmp_3d)
354      CALL histwrite(nid_mth,"dtec",itau_w,zx_tmp_3d,
355     .                                   iim*jjmp1*klev,ndex3d)
356c
357      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
358      CALL histwrite(nid_mth,"duvdf",itau_w,zx_tmp_3d,
359     .                                   iim*jjmp1*klev,ndex3d)
360c
361      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
362      CALL histwrite(nid_mth,"dvvdf",itau_w,zx_tmp_3d,
363     .                                   iim*jjmp1*klev,ndex3d)
364c
365      IF (ok_orodr) THEN
366      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d)
367      CALL histwrite(nid_mth,"duoro",itau_w,zx_tmp_3d,
368     .                                   iim*jjmp1*klev,ndex3d)
369c
370      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oro, zx_tmp_3d)
371      CALL histwrite(nid_mth,"dvoro",itau_w,zx_tmp_3d,
372     .                                   iim*jjmp1*klev,ndex3d)
373c
374      ENDIF
375C
376      IF (ok_orolf) THEN
377      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_lif, zx_tmp_3d)
378      CALL histwrite(nid_mth,"dulif",itau_w,zx_tmp_3d,
379     .                                   iim*jjmp1*klev,ndex3d)
380c
381      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_lif, zx_tmp_3d)
382      CALL histwrite(nid_mth,"dvlif",itau_w,zx_tmp_3d,
383     .                                   iim*jjmp1*klev,ndex3d)
384      ENDIF
385C
386      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
387      CALL histwrite(nid_mth,"ozone",itau_w,zx_tmp_3d,
388     .                                   iim*jjmp1*klev,ndex3d)
389c
390      IF (nqmax.GE.3) THEN
391      DO iq=1,nqmax-2
392      IF (iq.LE.99) THEN
393         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d)
394         WRITE(str2,'(i2.2)') iq
395         CALL histwrite(nid_mth,"trac"//str2,itau_w,zx_tmp_3d,
396     .                                   iim*jjmp1*klev,ndex3d)
397      ELSE
398         PRINT*, "Trop de traceurs"
399         CALL abort
400      ENDIF
401      ENDDO
402      ENDIF
403cKE43
404      IF (iflag_con.GE.3) THEN ! (sb)
405c
406      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
407      CALL histwrite(nid_mth,"upwd",itau_w,zx_tmp_3d,
408     .                                   iim*jjmp1*klev,ndex3d)
409c
410      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
411      CALL histwrite(nid_mth,"dnwd",itau_w,zx_tmp_3d,
412     .                                   iim*jjmp1*klev,ndex3d)
413c
414      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
415      CALL histwrite(nid_mth,"dnwd0",itau_w,zx_tmp_3d,
416     .                                   iim*jjmp1*klev,ndex3d)
417c
418      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
419      CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
420     .                                   iim*jjmp1*klev,ndex3d)
421c
422cccIM
423      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
424      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
425      CALL histwrite(nid_mth, "SWupTOA",itau_w,zx_tmp_2d,
426     .                               iim*jjmp1,ndex2d)
427c
428      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
429      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
430      CALL histwrite(nid_mth, "SWupSFC",itau_w,zx_tmp_2d,
431     .                               iim*jjmp1,ndex2d)
432c
433      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
434      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
435      CALL histwrite(nid_mth, "SWdnTOA",itau_w,zx_tmp_2d,
436     .                               iim*jjmp1,ndex2d)
437c
438      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
439      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
440      CALL histwrite(nid_mth, "SWdnSFC",itau_w,zx_tmp_2d,
441     .                               iim*jjmp1,ndex2d)
442c
443     
444      ENDIF
445
446      if (ok_sync) then
447        call histsync(nid_mth)
448      endif
449      ENDIF
450
Note: See TracBrowser for help on using the repository browser.