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

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

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