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

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

Phasage avec la version de Ionela
IM/LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 25.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
78cIM: 071003
79      zx_tmp_fi2d(1:klon)=evap(1:klon)*86400.
80      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
81      CALL histwrite(nid_mth,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
82c
83      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
84      CALL histwrite(nid_mth,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
85c
86      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
87      CALL histwrite(nid_mth,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
88c
89      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
90      CALL histwrite(nid_mth,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
91c
92      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
93      CALL histwrite(nid_mth,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
94c
95      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
96      CALL histwrite(nid_mth,"LWdnSFC",itau_w,zx_tmp_2d,iim*jjmp1,
97     .               ndex2d)
98cIM: 071003
99      zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
100      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
101      CALL histwrite(nid_mth,"LWupSFC",itau_w,zx_tmp_2d,iim*jjmp1,
102     .               ndex2d)
103cLWupSFC
104      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
105      CALL histwrite(nid_mth,"tops0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
106c
107      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
108      CALL histwrite(nid_mth,"topl0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
109c
110      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
111      CALL histwrite(nid_mth,"sols0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
112c
113      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
114      CALL histwrite(nid_mth,"soll0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
115c
116      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
117      CALL histwrite(nid_mth,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
118c
119      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
120      CALL histwrite(nid_mth,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
121c
122      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
123      CALL histwrite(nid_mth,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
124c
125c
126c      DO i = 1, klon
127c         zx_tmp_fi2d(i) = fluxu(i,1)
128c      ENDDO
129c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
130c      CALL histwrite(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
131c
132c      DO i = 1, klon
133c         zx_tmp_fi2d(i) = fluxv(i,1)
134c      ENDDO
135c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
136c      CALL histwrite(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
137c
138c
139       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
140       CALL histwrite(nid_mth,"ffonte",itau_w,zx_tmp_2d,iim*jjmp1,
141     $                ndex2d)
142c
143       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
144       CALL histwrite(nid_mth,"fqcalving",itau_w,zx_tmp_2d,iim*jjmp1,
145     $                ndex2d)
146c
147
148      DO nsrf = 1, nbsrf
149CYYY
150        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
151        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
152        CALL histwrite(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
153     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
154C
155        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
156        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
157        CALL histwrite(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
158     $      zx_tmp_2d,iim*jjmp1,ndex2d)
159C
160        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
161        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
162        CALL histwrite(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
163     $      zx_tmp_2d,iim*jjmp1,ndex2d)
164C
165        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
166        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
167        CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
168     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
169C
170        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
171        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
172        CALL histwrite(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
173     $      zx_tmp_2d,iim*jjmp1,ndex2d)
174c
175        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
176        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
177        CALL histwrite(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
178     $      zx_tmp_2d,iim*jjmp1,ndex2d)
179C
180        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
181        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
182        CALL histwrite(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
183     $      zx_tmp_2d,iim*jjmp1,ndex2d)
184C
185        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
186        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
187        CALL histwrite(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
188     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
189C     
190        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
191        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
192        CALL histwrite(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
193     $      zx_tmp_2d,iim*jjmp1,ndex2d)
194C
195        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
196        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
197        CALL histwrite(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
198     $      zx_tmp_2d,iim*jjmp1,ndex2d)
199C
200        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
201        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
202        CALL histwrite(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
203     $      zx_tmp_2d,iim*jjmp1,ndex2d)
204c
205      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
206      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
207      CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),itau_w
208     $    ,zx_tmp_2d,iim*jjmp1,ndex2d)
209
210      END DO 
211cXXX      DO i = 1, klon
212cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
213cXXX      ENDDO
214cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
215cXXX      CALL histwrite(nid_mth,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
216c
217      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
218      CALL histwrite(nid_mth,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
219      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
220      CALL histwrite(nid_mth,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
221c
222      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
223      CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
224c
225      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
226      CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
227c
228      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
229      CALL histwrite(nid_mth,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
230c
231      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
232      CALL histwrite(nid_mth,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
233c
234      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
235      CALL histwrite(nid_mth,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
236c
237      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
238      CALL histwrite(nid_mth,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
239c
240      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
241      CALL histwrite(nid_mth,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
242c
243      zx_tmp_fi2d(1:klon) = flwp(1:klon)
244      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
245      CALL histwrite(nid_mth,"lwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
246c
247      zx_tmp_fi2d(1:klon) = fiwp(1:klon)
248      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
249      CALL histwrite(nid_mth,"iwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
250c
251      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
252      CALL histwrite(nid_mth,"lwcon",itau_w,zx_tmp_3d,
253     .               iim*jjmp1*klev,ndex3d)
254c
255      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
256      CALL histwrite(nid_mth,"iwcon",itau_w,zx_tmp_3d,
257     .               iim*jjmp1*klev,ndex3d)
258c
259      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
260      CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
261c
262      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
263      CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
264c
265      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
266      CALL histwrite(nid_mth,"uq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
267c
268      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
269      CALL histwrite(nid_mth,"vq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
270cKE43
271      IF (iflag_con .GE. 3) THEN ! sb
272c
273      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
274      CALL histwrite(nid_mth,"cape",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
275c
276      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
277      CALL histwrite(nid_mth,"pbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
278c
279      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
280      CALL histwrite(nid_mth,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
281c
282      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
283      CALL histwrite(nid_mth,"fbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
284c
285      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
286      CALL histwrite(nid_mth,"upwd",itau_w,zx_tmp_3d,
287     .                                   iim*jjmp1*klev,ndex3d)
288c
289      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
290      CALL histwrite(nid_mth,"dnwd",itau_w,zx_tmp_3d,
291     .                                   iim*jjmp1*klev,ndex3d)
292c
293      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
294      CALL histwrite(nid_mth,"dnwd0",itau_w,zx_tmp_3d,
295     .                                   iim*jjmp1*klev,ndex3d)
296c
297      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
298      CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
299     .                                   iim*jjmp1*klev,ndex3d)
300c
301      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
302      CALL histwrite(nid_mth,"prw",itau_w,zx_tmp_2d,
303     .               iim*jjmp1,ndex2d)
304c
305      ENDIF
306c34EK
307c
308c Champs 3D:
309C
310      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
311      CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
312     .                                   iim*jjmp1*klev,ndex3d)
313c
314      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
315      CALL histwrite(nid_mth,"ovap",itau_w,zx_tmp_3d,
316     .                                   iim*jjmp1*klev,ndex3d)
317cIM: 071003
318      zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
319     .                         (1-qx(1:klon,1:klev,ivap))
320      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
321      CALL histwrite(nid_mth,"wvap",itau_w,zx_tmp_3d,
322     .                                   iim*jjmp1*klev,ndex3d)
323c
324      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
325      CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
326     .                                   iim*jjmp1*klev,ndex3d)
327c
328      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
329      CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
330     .                                   iim*jjmp1*klev,ndex3d)
331c
332      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
333      CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
334     .                                   iim*jjmp1*klev,ndex3d)
335c
336      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
337      CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
338     .                                   iim*jjmp1*klev,ndex3d)
339c
340      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
341      CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
342     .                                   iim*jjmp1*klev,ndex3d)
343c
344      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
345      CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d,
346     .                                   iim*jjmp1*klev,ndex3d)
347c
348      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
349      CALL histwrite(nid_mth,"rnebcon",itau_w,zx_tmp_3d,
350     .                                   iim*jjmp1*klev,ndex3d)
351c
352      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
353      CALL histwrite(nid_mth,"rhum",itau_w,zx_tmp_3d,
354     .                                   iim*jjmp1*klev,ndex3d)
355c
356      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
357      CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d,
358     .                                   iim*jjmp1*klev,ndex3d)
359c
360      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d)
361      CALL histwrite(nid_mth,"oliq",itau_w,zx_tmp_3d,
362     .                                   iim*jjmp1*klev,ndex3d)
363c
364      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
365      CALL histwrite(nid_mth,"dtphy",itau_w,zx_tmp_3d,
366     .                                   iim*jjmp1*klev,ndex3d)
367c
368      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
369      CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d,
370     .                                   iim*jjmp1*klev,ndex3d)
371c
372      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
373      CALL histwrite(nid_mth,"dqdyn",itau_w,zx_tmp_3d,
374     .                                   iim*jjmp1*klev,ndex3d)
375c
376cIM: 101003 : K/30min ==> K/s
377      zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
378      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
379      CALL histwrite(nid_mth,"dtcon",itau_w,zx_tmp_3d,
380     .                                   iim*jjmp1*klev,ndex3d)
381c
382      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d)
383      CALL histwrite(nid_mth,"dqcon",itau_w,zx_tmp_3d,
384     .                                   iim*jjmp1*klev,ndex3d)
385c
386cIM: 101003 : K/30min ==> K/s
387      zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
388      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
389      CALL histwrite(nid_mth,"dtlsc",itau_w,zx_tmp_3d,
390     .                                   iim*jjmp1*klev,ndex3d)
391cIM: 071003
392cIM: 101003 : K/30min ==> K/s
393      zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
394     .                             d_t_eva(1:klon,1:klev))/pdtphys
395      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
396      CALL histwrite(nid_mth,"dtlschr",itau_w,zx_tmp_3d,
397     .                                   iim*jjmp1*klev,ndex3d)
398c
399      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d)
400      CALL histwrite(nid_mth,"dqlsc",itau_w,zx_tmp_3d,
401     .                                   iim*jjmp1*klev,ndex3d)
402c
403cIM: 101003 : K/30min ==> K/s
404      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
405      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
406      CALL histwrite(nid_mth,"dtvdf",itau_w,zx_tmp_3d,
407     .                                   iim*jjmp1*klev,ndex3d)
408c
409      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
410      CALL histwrite(nid_mth,"dqvdf",itau_w,zx_tmp_3d,
411     .                                   iim*jjmp1*klev,ndex3d)
412c
413cIM: 101003 : K/30min ==> K/s
414      zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
415      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
416      CALL histwrite(nid_mth,"dteva",itau_w,zx_tmp_3d,
417     .                                   iim*jjmp1*klev,ndex3d)
418c
419      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d)
420      CALL histwrite(nid_mth,"dqeva",itau_w,zx_tmp_3d,
421     .                                   iim*jjmp1*klev,ndex3d)
422c
423      zpt_conv = 0.
424      where (ptconv) zpt_conv = 1.
425      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
426      CALL histwrite(nid_mth,"ptconv",itau_w,zx_tmp_3d,
427     .                                   iim*(jjmp1)*klev,ndex3d)
428c
429      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
430      CALL histwrite(nid_mth,"ratqs",itau_w,zx_tmp_3d,
431     .                                   iim*(jjmp1)*klev,ndex3d)
432c
433cIM: 101003 : K/30min ==> K/s
434      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
435      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
436      CALL histwrite(nid_mth,"dtajs",itau_w,zx_tmp_3d,
437     .                                   iim*jjmp1*klev,ndex3d)
438c
439      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d)
440      CALL histwrite(nid_mth,"dqajs",itau_w,zx_tmp_3d,
441     .                                   iim*jjmp1*klev,ndex3d)
442c
443cIM: 101003 : K/day ==> K/s
444cIM: LMD_ARMIP3   zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)*pdtphys/RDAY
445      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
446      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
447      CALL histwrite(nid_mth,"dtswr",itau_w,zx_tmp_3d,
448     .                                   iim*jjmp1*klev,ndex3d)
449c
450cIM: 101003 : K/day ==> K/s     
451cIM: LMD_ARMIP3   zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)*pdtphys/RDAY
452      zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
453      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
454      CALL histwrite(nid_mth,"dtsw0",itau_w,zx_tmp_3d,
455     .                                   iim*jjmp1*klev,ndex3d)
456c
457cIM: 101003 : K/day ==> K/s     
458cIM: LMD_ARMIP3     zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)*pdtphys/RDAY
459      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
460      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
461      CALL histwrite(nid_mth,"dtlwr",itau_w,zx_tmp_3d,
462     .                                   iim*jjmp1*klev,ndex3d)
463c
464cIM: 101003 : K/day ==> K/s     
465cIM: LMD_ARMIP3     zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)*pdtphys/RDAY
466      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
467      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
468      CALL histwrite(nid_mth,"dtlw0",itau_w,zx_tmp_3d,
469     .                                   iim*jjmp1*klev,ndex3d)
470c
471cIM: 101003 : deja en K/s     
472      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
473      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
474      CALL histwrite(nid_mth,"dtec",itau_w,zx_tmp_3d,
475     .                                   iim*jjmp1*klev,ndex3d)
476c
477      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
478      CALL histwrite(nid_mth,"duvdf",itau_w,zx_tmp_3d,
479     .                                   iim*jjmp1*klev,ndex3d)
480c
481      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
482      CALL histwrite(nid_mth,"dvvdf",itau_w,zx_tmp_3d,
483     .                                   iim*jjmp1*klev,ndex3d)
484c
485      IF (ok_orodr) THEN
486      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d)
487      CALL histwrite(nid_mth,"duoro",itau_w,zx_tmp_3d,
488     .                                   iim*jjmp1*klev,ndex3d)
489c
490      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oro, zx_tmp_3d)
491      CALL histwrite(nid_mth,"dvoro",itau_w,zx_tmp_3d,
492     .                                   iim*jjmp1*klev,ndex3d)
493c
494      ENDIF
495C
496      IF (ok_orolf) THEN
497      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_lif, zx_tmp_3d)
498      CALL histwrite(nid_mth,"dulif",itau_w,zx_tmp_3d,
499     .                                   iim*jjmp1*klev,ndex3d)
500c
501      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_lif, zx_tmp_3d)
502      CALL histwrite(nid_mth,"dvlif",itau_w,zx_tmp_3d,
503     .                                   iim*jjmp1*klev,ndex3d)
504      ENDIF
505C
506      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
507      CALL histwrite(nid_mth,"ozone",itau_w,zx_tmp_3d,
508     .                                   iim*jjmp1*klev,ndex3d)
509c
510      IF (nqmax.GE.3) THEN
511      DO iq=1,nqmax-2
512      IF (iq.LE.99) THEN
513         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d)
514         WRITE(str2,'(i2.2)') iq
515         CALL histwrite(nid_mth,"trac"//str2,itau_w,zx_tmp_3d,
516     .                                   iim*jjmp1*klev,ndex3d)
517      ELSE
518         PRINT*, "Trop de traceurs"
519         CALL abort
520      ENDIF
521      ENDDO
522      ENDIF
523c
524c     zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
525      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 2 )
526      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
527      CALL histwrite(nid_mth, "SWupTOA",itau_w,zx_tmp_2d,
528     .                               iim*jjmp1,ndex2d)
529c
530c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
531      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
532      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
533      CALL histwrite(nid_mth, "SWupSFC",itau_w,zx_tmp_2d,
534     .                               iim*jjmp1,ndex2d)
535c
536c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
537      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 2 )
538      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
539      CALL histwrite(nid_mth, "SWdnTOA",itau_w,zx_tmp_2d,
540     .                               iim*jjmp1,ndex2d)
541c
542c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
543      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
544      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
545      CALL histwrite(nid_mth, "SWdnSFC",itau_w,zx_tmp_2d,
546     .                               iim*jjmp1,ndex2d)
547c
548c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, klevp1)
549      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 2 )
550      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
551      CALL histwrite(nid_mth, "SWupTOAclr",itau_w,zx_tmp_2d,
552     .                               iim*jjmp1,ndex2d)
553c
554c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
555      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
556      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
557      CALL histwrite(nid_mth, "SWupSFCclr",itau_w,zx_tmp_2d,
558     .                               iim*jjmp1,ndex2d)
559c
560c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, klevp1)
561      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 2 )
562      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
563      CALL histwrite(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_2d,
564     .                               iim*jjmp1,ndex2d)
565c
566c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
567      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
568      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
569      CALL histwrite(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_2d,
570     .                               iim*jjmp1,ndex2d)
571c
572c  Champs interpolles sur des niveaux de pression
573
574      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
575      CALL histwrite(nid_mth,"u850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
576
577      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
578      CALL histwrite(nid_mth,"v850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
579
580      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
581      CALL histwrite(nid_mth,"u500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
582
583      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
584      CALL histwrite(nid_mth,"v500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
585
586      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
587      CALL histwrite(nid_mth,"u200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
588
589      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
590      CALL histwrite(nid_mth,"v200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
591
592      CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
593      CALL histwrite(nid_mth,"phi500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
594
595      CALL gr_fi_ecrit(1, klon,iim,jjmp1, w500,zx_tmp_2d)
596      CALL histwrite(nid_mth,"w500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
597c
598      if (ok_sync) then
599        call histsync(nid_mth)
600      endif
601      ENDIF
Note: See TracBrowser for help on using the repository browser.