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

Last change on this file since 436 was 433, checked in by lmdzadmin, 22 years ago

Convergence avec la version de Ionela dec 2002

YOMCST.? : suppression RI0 (IM)
albedo.F : facteur 1.2 sur le nouveau calcul (IM)
clesphys.h : rajout de différentes ctes (concentration des gaz) (IM)
clmain.F : separation des flux LW, SW (JLD)

remplace qsurf par yqsol (IM)

conf_phys.F90 : rajout de différentes ctes (gaz + orbite) (IM)
convect3.F : DPINV+SIGD*0.5*(EVAP(1)+EVAP(2)) (SBL)
cv3_routines.F:
cvparam3.h : compatibilite avec conema3 TEMPORAIRE (FH)
phyetat0.F : lecture de co2_ppm et solaire pour tests de coherence
phyredem.F : co2_ppm et solaire passé en common
physiq.F : separation flux LW, SW

rajout diagnostiques (slp, w500)
suppression iflag_con = 4
clwcon0=qcondc (FH)
position dU "ENDIF ! ok_cvl"

radlwsw.F : passage des concentrations gaz dans un common (IM)

PEMIS(i) = 1.0 (JLD pour cohérence ORCHIDEE)

stdlevvar.F90 :
suphec.F : suppression init. des ctes orbitales (IM)

nouvelles E/S (ini_hist..., write_hist...)

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 20.8 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) = fsollw( 1 : klon, nsrf)
146        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
147        CALL histwrite(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
148     $      zx_tmp_2d,iim*jjmp1,ndex2d)
149cIM cf. JLD
150        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
151        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
152        CALL histwrite(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
153     $      zx_tmp_2d,iim*jjmp1,ndex2d)
154C
155        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
156        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
157        CALL histwrite(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
158     $      zx_tmp_2d,iim*jjmp1,ndex2d)
159C
160        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
161        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
162        CALL histwrite(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
163     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
164C     
165        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
166        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
167        CALL histwrite(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
168     $      zx_tmp_2d,iim*jjmp1,ndex2d)
169C
170        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
171        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
172        CALL histwrite(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
173     $      zx_tmp_2d,iim*jjmp1,ndex2d)
174C
175        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
176        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
177        CALL histwrite(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
178     $      zx_tmp_2d,iim*jjmp1,ndex2d)
179c
180      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
181      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
182      CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),itau_w
183     $    ,zx_tmp_2d,iim*jjmp1,ndex2d)
184
185      END DO 
186cXXX      DO i = 1, klon
187cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
188cXXX      ENDDO
189cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
190cXXX      CALL histwrite(nid_mth,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
191c
192      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
193      CALL histwrite(nid_mth,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
194      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
195      CALL histwrite(nid_mth,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
196c
197      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
198      CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
199c
200      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
201      CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
202c
203      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
204      CALL histwrite(nid_mth,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
205c
206      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
207      CALL histwrite(nid_mth,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
208c
209      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
210      CALL histwrite(nid_mth,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
211c
212      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
213      CALL histwrite(nid_mth,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
214c
215      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
216      CALL histwrite(nid_mth,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
217c
218      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
219      CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
220c
221      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
222      CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
223c
224      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
225      CALL histwrite(nid_mth,"uq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
226c
227      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
228      CALL histwrite(nid_mth,"vq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
229cKE43
230      IF (iflag_con .GE. 3) THEN ! sb
231c
232      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
233      CALL histwrite(nid_mth,"cape",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
234c
235      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
236      CALL histwrite(nid_mth,"pbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
237c
238      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
239      CALL histwrite(nid_mth,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
240c
241      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
242      CALL histwrite(nid_mth,"fbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
243c
244c
245      ENDIF
246c34EK
247c
248c Champs 3D:
249C
250      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
251      CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
252     .                                   iim*jjmp1*klev,ndex3d)
253c
254      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
255      CALL histwrite(nid_mth,"ovap",itau_w,zx_tmp_3d,
256     .                                   iim*jjmp1*klev,ndex3d)
257c
258      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
259      CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
260     .                                   iim*jjmp1*klev,ndex3d)
261c
262      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
263      CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
264     .                                   iim*jjmp1*klev,ndex3d)
265c
266      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
267      CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
268     .                                   iim*jjmp1*klev,ndex3d)
269c
270      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
271      CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
272     .                                   iim*jjmp1*klev,ndex3d)
273c
274      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
275      CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
276     .                                   iim*jjmp1*klev,ndex3d)
277c
278      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
279      CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d,
280     .                                   iim*jjmp1*klev,ndex3d)
281c
282      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
283      CALL histwrite(nid_mth,"rnebcon",itau_w,zx_tmp_3d,
284     .                                   iim*jjmp1*klev,ndex3d)
285c
286      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
287      CALL histwrite(nid_mth,"rhum",itau_w,zx_tmp_3d,
288     .                                   iim*jjmp1*klev,ndex3d)
289c
290      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d)
291      CALL histwrite(nid_mth,"oliq",itau_w,zx_tmp_3d,
292     .                                   iim*jjmp1*klev,ndex3d)
293c
294      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
295      CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d,
296     .                                   iim*jjmp1*klev,ndex3d)
297c
298      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
299      CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d,
300     .                                   iim*jjmp1*klev,ndex3d)
301c
302      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
303      CALL histwrite(nid_mth,"dqdyn",itau_w,zx_tmp_3d,
304     .                                   iim*jjmp1*klev,ndex3d)
305c
306      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_con, zx_tmp_3d)
307      CALL histwrite(nid_mth,"dtcon",itau_w,zx_tmp_3d,
308     .                                   iim*jjmp1*klev,ndex3d)
309c
310      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d)
311      CALL histwrite(nid_mth,"dqcon",itau_w,zx_tmp_3d,
312     .                                   iim*jjmp1*klev,ndex3d)
313c
314      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_lsc, zx_tmp_3d)
315      CALL histwrite(nid_mth,"dtlsc",itau_w,zx_tmp_3d,
316     .                                   iim*jjmp1*klev,ndex3d)
317c
318      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d)
319      CALL histwrite(nid_mth,"dqlsc",itau_w,zx_tmp_3d,
320     .                                   iim*jjmp1*klev,ndex3d)
321c
322      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d)
323      CALL histwrite(nid_mth,"dtvdf",itau_w,zx_tmp_3d,
324     .                                   iim*jjmp1*klev,ndex3d)
325c
326      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
327      CALL histwrite(nid_mth,"dqvdf",itau_w,zx_tmp_3d,
328     .                                   iim*jjmp1*klev,ndex3d)
329c
330      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_eva, zx_tmp_3d)
331      CALL histwrite(nid_mth,"dteva",itau_w,zx_tmp_3d,
332     .                                   iim*jjmp1*klev,ndex3d)
333c
334      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d)
335      CALL histwrite(nid_mth,"dqeva",itau_w,zx_tmp_3d,
336     .                                   iim*jjmp1*klev,ndex3d)
337c
338      zpt_conv = 0.
339      where (ptconv) zpt_conv = 1.
340      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
341      CALL histwrite(nid_mth,"ptconv",itau_w,zx_tmp_3d,
342     .                                   iim*(jjmp1)*klev,ndex3d)
343c
344      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
345      CALL histwrite(nid_mth,"ratqs",itau_w,zx_tmp_3d,
346     .                                   iim*(jjmp1)*klev,ndex3d)
347c
348      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d)
349      CALL histwrite(nid_mth,"dtajs",itau_w,zx_tmp_3d,
350     .                                   iim*jjmp1*klev,ndex3d)
351c
352      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d)
353      CALL histwrite(nid_mth,"dqajs",itau_w,zx_tmp_3d,
354     .                                   iim*jjmp1*klev,ndex3d)
355c
356      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat, zx_tmp_3d)
357      CALL histwrite(nid_mth,"dtswr",itau_w,zx_tmp_3d,
358     .                                   iim*jjmp1*klev,ndex3d)
359c
360      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat0, zx_tmp_3d)
361      CALL histwrite(nid_mth,"dtsw0",itau_w,zx_tmp_3d,
362     .                                   iim*jjmp1*klev,ndex3d)
363c
364      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool, zx_tmp_3d)
365      CALL histwrite(nid_mth,"dtlwr",itau_w,zx_tmp_3d,
366     .                                   iim*jjmp1*klev,ndex3d)
367c
368      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool0, zx_tmp_3d)
369      CALL histwrite(nid_mth,"dtlw0",itau_w,zx_tmp_3d,
370     .                                   iim*jjmp1*klev,ndex3d)
371c
372      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ec, zx_tmp_3d)
373      CALL histwrite(nid_mth,"dtec",itau_w,zx_tmp_3d,
374     .                                   iim*jjmp1*klev,ndex3d)
375c
376      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
377      CALL histwrite(nid_mth,"duvdf",itau_w,zx_tmp_3d,
378     .                                   iim*jjmp1*klev,ndex3d)
379c
380      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
381      CALL histwrite(nid_mth,"dvvdf",itau_w,zx_tmp_3d,
382     .                                   iim*jjmp1*klev,ndex3d)
383c
384      IF (ok_orodr) THEN
385      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d)
386      CALL histwrite(nid_mth,"duoro",itau_w,zx_tmp_3d,
387     .                                   iim*jjmp1*klev,ndex3d)
388c
389      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oro, zx_tmp_3d)
390      CALL histwrite(nid_mth,"dvoro",itau_w,zx_tmp_3d,
391     .                                   iim*jjmp1*klev,ndex3d)
392c
393      ENDIF
394C
395      IF (ok_orolf) THEN
396      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_lif, zx_tmp_3d)
397      CALL histwrite(nid_mth,"dulif",itau_w,zx_tmp_3d,
398     .                                   iim*jjmp1*klev,ndex3d)
399c
400      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_lif, zx_tmp_3d)
401      CALL histwrite(nid_mth,"dvlif",itau_w,zx_tmp_3d,
402     .                                   iim*jjmp1*klev,ndex3d)
403      ENDIF
404C
405      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
406      CALL histwrite(nid_mth,"ozone",itau_w,zx_tmp_3d,
407     .                                   iim*jjmp1*klev,ndex3d)
408c
409      IF (nqmax.GE.3) THEN
410      DO iq=1,nqmax-2
411      IF (iq.LE.99) THEN
412         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d)
413         WRITE(str2,'(i2.2)') iq
414         CALL histwrite(nid_mth,"trac"//str2,itau_w,zx_tmp_3d,
415     .                                   iim*jjmp1*klev,ndex3d)
416      ELSE
417         PRINT*, "Trop de traceurs"
418         CALL abort
419      ENDIF
420      ENDDO
421      ENDIF
422cKE43
423      IF (iflag_con.GE.3) THEN ! (sb)
424c
425      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
426      CALL histwrite(nid_mth,"upwd",itau_w,zx_tmp_3d,
427     .                                   iim*jjmp1*klev,ndex3d)
428c
429      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
430      CALL histwrite(nid_mth,"dnwd",itau_w,zx_tmp_3d,
431     .                                   iim*jjmp1*klev,ndex3d)
432c
433      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
434      CALL histwrite(nid_mth,"dnwd0",itau_w,zx_tmp_3d,
435     .                                   iim*jjmp1*klev,ndex3d)
436c
437      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
438      CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
439     .                                   iim*jjmp1*klev,ndex3d)
440c
441cccIM
442      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
443      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
444      CALL histwrite(nid_mth, "SWupTOA",itau_w,zx_tmp_2d,
445     .                               iim*jjmp1,ndex2d)
446c
447      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
448      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
449      CALL histwrite(nid_mth, "SWupSFC",itau_w,zx_tmp_2d,
450     .                               iim*jjmp1,ndex2d)
451c
452      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
453      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
454      CALL histwrite(nid_mth, "SWdnTOA",itau_w,zx_tmp_2d,
455     .                               iim*jjmp1,ndex2d)
456c
457      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
458      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
459      CALL histwrite(nid_mth, "SWdnSFC",itau_w,zx_tmp_2d,
460     .                               iim*jjmp1,ndex2d)
461c
462cccIM clear sky
463      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, klevp1)
464      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
465      CALL histwrite(nid_mth, "SWupTOAclr",itau_w,zx_tmp_2d,
466     .                               iim*jjmp1,ndex2d)
467c
468      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
469      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
470      CALL histwrite(nid_mth, "SWupSFCclr",itau_w,zx_tmp_2d,
471     .                               iim*jjmp1,ndex2d)
472c
473      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, klevp1)
474      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
475      CALL histwrite(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_2d,
476     .                               iim*jjmp1,ndex2d)
477c
478      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
479      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
480      CALL histwrite(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_2d,
481     .                               iim*jjmp1,ndex2d)
482c
483c  Champs interpolles sur des niveaux de pression
484
485      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
486      CALL histwrite(nid_mth,"u850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
487
488      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
489      CALL histwrite(nid_mth,"v850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
490
491      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
492      CALL histwrite(nid_mth,"u500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
493
494      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
495      CALL histwrite(nid_mth,"v500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
496
497      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
498      CALL histwrite(nid_mth,"u200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
499
500      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
501      CALL histwrite(nid_mth,"v200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
502
503      CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
504      CALL histwrite(nid_mth,"phi500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
505
506      CALL gr_fi_ecrit(1, klon,iim,jjmp1, w500,zx_tmp_2d)
507      CALL histwrite(nid_mth,"w500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
508
509
510     
511       ENDIF
512
513      if (ok_sync) then
514        call histsync(nid_mth)
515      endif
516      ENDIF
517
Note: See TracBrowser for help on using the repository browser.