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

Last change on this file since 456 was 455, checked in by lmdzadmin, 22 years ago

Ajout slp, ffonte, fqcalving IM/FH/JLD
  • 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",itap,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",itap,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",itap,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",itap,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",itap,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",itap,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",itap,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",itap,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.