source: trunk/LMDZ.TITAN/libf/phytitan/write_histmth.h @ 201

Last change on this file since 201 was 175, checked in by slebonnois, 14 years ago

S.LEBONNOIS:

  • Revision majeure de la physique Titan => ajout des nuages version 10 bins (Jeremie Burgalat) Cette version reste a tester mais avec clouds=0, on reste sur l'ancienne.
  • Quelques ajouts dans la doc.
File size: 18.2 KB
Line 
1      IF (ok_mensuel) THEN
2c
3      ndex2d = 0
4      ndex3d = 0
5      zx_tmp_2d = 0.
6      zx_tmp_3d = 0.
7      zx_tmp_fi2d=0.
8      zx_tmp_fi3d=0.
9c
10         zsto = dtime
11         zout = dtime * FLOAT(ecrit_mth)
12         itau_w = itau_phy + itap
13c
14c-------------------------------------------------------
15      IF(lev_histmth.GE.1) THEN
16c
17ccccccccccccc 2D fields, invariables
18c
19      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
20      CALL histwrite(nid_mth,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
21C
22      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
23      CALL histwrite(nid_mth,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
24c
25ccccccc axe Ls
26      do j=1,jjmp1
27       do i=1,iim
28        zx_tmp_2d(i,j)=zls*180./RPI      ! zls est en radians !!
29       enddo
30      enddo
31      CALL histwrite(nid_mth,"ls",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
32c
33ccccccccccccc 2D fields, variables
34c
35      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ftsol,zx_tmp_2d)
36      CALL histwrite(nid_mth,"tsol",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
44c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
45c     CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
46c
47c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
48c     CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
49c
50c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
51c     CALL histwrite(nid_mth,"cdragh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
52c
53c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
54c     CALL histwrite(nid_mth,"cdragm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
55c
56      ENDIF !lev_histmth.GE.1
57c
58c-------------------------------------------------------
59      IF(lev_histmth.GE.2) THEN
60c
61ccccccccccccc 3D fields, basics
62c
63      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
64      CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
65     .                                   iim*jjmp1*klev,ndex3d)
66c
67      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
68      CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
69     .                                   iim*jjmp1*klev,ndex3d)
70c
71      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
72      CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
73     .                                   iim*jjmp1*klev,ndex3d)
74c
75      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
76      CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
77     .                                   iim*jjmp1*klev,ndex3d)
78c
79      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
80      CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
81     .                                   iim*jjmp1*klev,ndex3d)
82c
83      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
84      CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
85     .                                   iim*jjmp1*klev,ndex3d)
86c
87c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ycoefh, zx_tmp_3d)
88c     CALL histwrite(nid_mth,"Kz",itau_w,zx_tmp_3d,
89c    .                                   iim*jjmp1*klev,ndex3d)
90c
91      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
92      CALL histwrite(nid_mth,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
93c
94cccccccccccccccccc  Tracers
95c
96         if (iflag_trac.eq.1) THEN
97          if (microfi.ge.1) then
98c          DO iq=1,nmicro
99c      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qaer(1,1,iq), zx_tmp_3d)
100c      CALL histwrite(nid_mth,tname(iq),itau_w,zx_tmp_3d,
101c    .                                   iim*jjmp1*klev,ndex3d)
102c          ENDDO
103c    -------   NB AER TOT
104               do i=1,klon
105                 do j=1,klev
106                   zx_tmp_fi3d(i,j)= SUM(qaer(i,j,1:nrad))
107                 enddo
108               enddo
109       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
110       CALL histwrite(nid_mth,"qaer",itau_w,zx_tmp_3d,
111     .                                   iim*jjmp1*klev,ndex3d)
112c    -------   NB NOY TOT
113               do i=1,klon
114                 do j=1,klev
115                   zx_tmp_fi3d(i,j)= SUM(qaer(i,j,nrad+1:2*nrad))
116                 enddo
117               enddo
118       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
119       CALL histwrite(nid_mth,"qnoy",itau_w,zx_tmp_3d,
120     .                                   iim*jjmp1*klev,ndex3d)
121c    -------   V GLA1 TOT
122               do i=1,klon
123                 do j=1,klev
124                   zx_tmp_fi3d(i,j)= SUM(qaer(i,j,2*nrad+1:3*nrad))
125                 enddo
126               enddo
127       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
128       CALL histwrite(nid_mth,"qgl1",itau_w,zx_tmp_3d,
129     .                                   iim*jjmp1*klev,ndex3d)
130c    -------   V GLA2 TOT
131               do i=1,klon
132                 do j=1,klev
133                   zx_tmp_fi3d(i,j)= SUM(qaer(i,j,3*nrad+1:4*nrad))
134                 enddo
135               enddo
136       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
137       CALL histwrite(nid_mth,"qgl2",itau_w,zx_tmp_3d,
138     .                                   iim*jjmp1*klev,ndex3d)
139c    -------   V GLA3 TOT
140               do i=1,klon
141                 do j=1,klev
142                   zx_tmp_fi3d(i,j)= SUM(qaer(i,j,4*nrad+1:5*nrad))
143                 enddo
144               enddo
145       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
146       CALL histwrite(nid_mth,"qgl3",itau_w,zx_tmp_3d,
147     .                                   iim*jjmp1*klev,ndex3d)
148c --------------
149c ----- SATURATION ESP NUAGES
150             if (clouds.eq.1) then
151
152       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, satch4,zx_tmp_3d)
153       CALL histwrite(nid_mth,"ch4sat", itau_w, zx_tmp_3d,
154     .                                   iim*jjmp1*klev,ndex3d)
155
156       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, satc2h6,zx_tmp_3d)
157       CALL histwrite(nid_mth,"c2h6sat", itau_w, zx_tmp_3d,
158     .                                   iim*jjmp1*klev,ndex3d)
159
160       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, satc2h2,zx_tmp_3d)
161       CALL histwrite(nid_mth,"c2h2sat", itau_w, zx_tmp_3d,
162     .                                   iim*jjmp1*klev,ndex3d)
163c --------------
164c ----- RESERVOIR DE SURFACE
165       CALL gr_fi_ecrit(1, klon,iim,jjmp1,reservoir,zx_tmp_2d)
166       CALL histwrite(nid_mth,"reserv",itau_w,zx_tmp_2d,
167     .                        iim*jjmp1,ndex2d)
168c --------------
169c ----- PRECIPITATIONS
170c       -----  CH4
171       CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,1),zx_tmp_2d)
172       CALL histwrite(nid_mth,"prech4",itau_w,zx_tmp_2d,
173     .                        iim*jjmp1,ndex2d)
174c       -----  C2H6
175       CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,2),zx_tmp_2d)
176       CALL histwrite(nid_mth,"prec2h6",itau_w,zx_tmp_2d,
177     .                        iim*jjmp1,ndex2d)
178c       -----  C2H2
179       CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,3),zx_tmp_2d)
180       CALL histwrite(nid_mth,"prec2h2",itau_w,zx_tmp_2d,
181     .                        iim*jjmp1,ndex2d)
182c
183c --------------
184c ----- FLUX GLACE
185       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,1),zx_tmp_3d)
186       CALL histwrite(nid_mth,"flxgl1", itau_w, zx_tmp_3d,
187     .                                   iim*jjmp1*klev,ndex3d)
188       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,2),zx_tmp_3d)
189       CALL histwrite(nid_mth,"flxgl2", itau_w, zx_tmp_3d,
190     .                                   iim*jjmp1*klev,ndex3d)
191       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,3),zx_tmp_3d)
192       CALL histwrite(nid_mth,"flxgl3", itau_w, zx_tmp_3d,
193     .                                   iim*jjmp1*klev,ndex3d)
194c
195c --------------
196c ----- RAYON MOYEN GOUTTE
197       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rmcloud,zx_tmp_3d)
198       CALL histwrite(nid_mth,"rcldbar", itau_w, zx_tmp_3d,
199     .                                   iim*jjmp1*klev,ndex3d)
200c
201             endif
202          endif
203c
204c --------------
205c ----- TRACEURS CHIMIQUES
206          if (nmicro.lt.nqmax) then
207           DO iq=nmicro+1,nqmax
208       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,tr_seri(1,1,iq),zx_tmp_3d)
209       CALL histwrite(nid_mth,tname(iq),itau_w,zx_tmp_3d,
210     .                                   iim*jjmp1*klev,ndex3d)
211           ENDDO
212c Condensation:
213c          DO iq=nmicro+1,nqmax
214c      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,d_tr_mph(1,1,iq),zx_tmp_3d)
215c      CALL histwrite(nid_mth,"c_"//tname(iq),itau_w,zx_tmp_3d,
216c    .                                   iim*jjmp1*klev,ndex3d)
217c          ENDDO
218          endif
219         endif
220c
221      ENDIF !lev_histmth.GE.2
222c
223c-------------------------------------------------------
224      IF(lev_histmth.GE.3) THEN
225c
226cccccccccccccccccc  Radiative transfer
227c
228c 2D
229c
230      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
231      CALL histwrite(nid_mth,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
232c
233      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
234      CALL histwrite(nid_mth,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
235c
236      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
237      CALL histwrite(nid_mth,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
238c
239c 3D
240c
241      zx_tmp_fi3d(1:klon,1:klev)=swnet(1:klon,1:klev)
242      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
243      CALL histwrite(nid_mth,"SWnet",itau_w,zx_tmp_3d,
244     .                                   iim*jjmp1*klev,ndex3d)
245c
246      zx_tmp_fi3d(1:klon,1:klev)=lwnet(1:klon,1:klev)
247      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
248      CALL histwrite(nid_mth,"LWnet",itau_w,zx_tmp_3d,
249     .                                   iim*jjmp1*klev,ndex3d)
250c
251c --------------
252c ----- OPACITE BRUME
253       do k=7,NSPECV,10
254         do i=1,klon
255         do l=1,klev
256           t_tauhvd(i,l)=TAUHVD(i,klev-l+1,k)
257         enddo
258         enddo
259         write(str1,'(i2.2)') k
260      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
261      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
262      CALL histwrite(nid_mth,"thv"//str1,itau_w,zx_tmp_3d,
263     .                                   iim*jjmp1*klev,ndex3d)
264       enddo      ! fin boucle NSPECV
265
266       do k=8,NSPECI,10
267         do i=1,klon
268         do l=1,klev
269           t_tauhvd(i,l)=TAUHID(i,klev-l+1,k)
270         enddo
271         enddo
272         write(str1,'(i2.2)') k
273      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
274      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
275      CALL histwrite(nid_mth,"thi"//str1,itau_w,zx_tmp_3d,
276     .                                   iim*jjmp1*klev,ndex3d)
277       enddo      ! fin boucle NSPECI
278c
279c --------------
280c ----- EXTINCTION BRUME
281       do k=7,NSPECV,10
282         do i=1,klon
283         do l=1,klev
284          if(l.ne.klev)
285     s     t_khvd(i,l)=TAUHVD(i,klev-l+1,k)
286     s                -TAUHVD(i,klev-l+1-1,k)
287          if(l.eq.klev)
288     s     t_khvd(i,l)=TAUHVD(i,klev-l+1,k)
289
290         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
291         enddo
292         enddo
293         write(str1,'(i2.2)') k
294      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
295      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
296      CALL histwrite(nid_mth,"khv"//str1,itau_w,zx_tmp_3d,
297     .                                   iim*jjmp1*klev,ndex3d)
298       enddo      ! fin boucle NSPECV
299
300       do k=8,NSPECI,10
301         do i=1,klon
302         do l=1,klev
303          if(l.ne.klev)
304     s     t_khvd(i,l)=TAUHID(i,klev-l+1,k)
305     s                -TAUHID(i,klev-l+1-1,k)
306          if(l.eq.klev)
307     s     t_khvd(i,l)=TAUHID(i,klev-l+1,k)
308
309         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
310         enddo
311         enddo
312         write(str1,'(i2.2)') k
313      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
314      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
315      CALL histwrite(nid_mth,"khi"//str1,itau_w,zx_tmp_3d,
316     .                                   iim*jjmp1*klev,ndex3d)
317       enddo      ! fin boucle NSPECI
318c
319c --------------
320c ----- OPACITE GAZ
321       do k=7,NSPECV,10
322         do i=1,klon
323         do l=1,klev
324           t_tauhvd(i,l)=TAUGVD(i,klev-l+1,k)
325         enddo
326         enddo
327         write(str1,'(i2.2)') k
328      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
329      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
330      CALL histwrite(nid_mth,"tgv"//str1,itau_w,zx_tmp_3d,
331     .                                   iim*jjmp1*klev,ndex3d)
332       enddo      ! fin boucle NSPECV
333
334       do k=8,NSPECI,10
335         do i=1,klon
336         do l=1,klev
337           t_tauhvd(i,l)=TAUGID(i,klev-l+1,k)
338         enddo
339         enddo
340         write(str1,'(i2.2)') k
341      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
342      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
343      CALL histwrite(nid_mth,"tgi"//str1,itau_w,zx_tmp_3d,
344     .                                   iim*jjmp1*klev,ndex3d)
345       enddo      ! fin boucle NSPECI
346c
347c --------------
348c ----- EXTINCTION GAZ
349       do k=7,NSPECV,10
350         do i=1,klon
351         do l=1,klev
352          if(l.ne.klev)
353     s     t_khvd(i,l)=TAUGVD(i,klev-l+1,k)
354     s                -TAUGVD(i,klev-l+1-1,k)
355          if(l.eq.klev)
356     s     t_khvd(i,l)=TAUGVD(i,klev-l+1,k)
357
358         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
359         enddo
360         enddo
361         write(str1,'(i2.2)') k
362      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
363      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
364      CALL histwrite(nid_mth,"kgv"//str1,itau_w,zx_tmp_3d,
365     .                                   iim*jjmp1*klev,ndex3d)
366       enddo      ! fin boucle NSPECV
367
368       do k=8,NSPECI,10
369         do i=1,klon
370         do l=1,klev
371          if(l.ne.klev)
372     s     t_khvd(i,l)=TAUGID(i,klev-l+1,k)
373     s                -TAUGID(i,klev-l+1-1,k)
374
375          if(l.eq.klev)
376     s     t_khvd(i,l)=TAUGID(i,klev-l+1,k)
377
378         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
379         enddo
380         enddo
381         write(str1,'(i2.2)') k
382      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
383      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
384      CALL histwrite(nid_mth,"kgi"//str1,itau_w,zx_tmp_3d,
385     .                                   iim*jjmp1*klev,ndex3d)
386       enddo      ! fin boucle NSPECI
387
388c --------------
389c ----- OPACITE NUAGES (ATTENTION PROXY)
390         if (clouds.eq.1) then
391           zx_tmp_fi3d(1:klon,1:klev)=occcld(1:klon,1:klev)
392           CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
393           CALL histwrite(nid_mth,"tcld",itau_w,zx_tmp_3d,
394     .                                   iim*jjmp1*klev,ndex3d)
395c --------------
396c ----- EXTINCTION NUAGES (ATTENTION PROXY)
397           do i=1,klon
398             t_kcld(i,klev)=occcld(i,klev)
399     .       /(zzlev(i,klev+1)-zzlev(i,klev))
400             do j=klev-1,1,-1
401               t_kcld(i,j)=(occcld(i,j)-occcld(i,j+1))
402     .         /(zzlev(i,j+1)-zzlev(i,j))
403             enddo
404           enddo
405           zx_tmp_fi3d(1:klon,1:klev)=t_kcld(1:klon,1:klev)
406           CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
407           CALL histwrite(nid_mth,"kcld",itau_w,zx_tmp_3d,
408     .                                   iim*jjmp1*klev,ndex3d)
409c
410c --------------
411c ----- OCCURENCE NUAGES
412           do k=1,12
413             write(str1,'(i2.2)') k
414             zx_tmp_fi3d(1:klon,1:klev)=occcld_m(1:klon,1:klev,k)
415             CALL histwrite(nid_mth,"occcld"//str1,itau_w,zx_tmp_3d,
416     .                                   iim*jjmp1*klev,ndex3d)
417           enddo
418c
419        endif 
420c
421      ENDIF !lev_histmth.GE.3
422c
423c-------------------------------------------------------
424      IF(lev_histmth.GE.4) THEN
425c
426      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
427      CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d,
428     .                                   iim*jjmp1*klev,ndex3d)
429c
430      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
431      CALL histwrite(nid_mth,"dtphy",itau_w,zx_tmp_3d,
432     .                                   iim*jjmp1*klev,ndex3d)
433c K/s
434      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)
435      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
436      CALL histwrite(nid_mth,"dtvdf",itau_w,zx_tmp_3d,
437     .                                   iim*jjmp1*klev,ndex3d)
438c
439c K/s
440      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)
441      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
442      CALL histwrite(nid_mth,"dtajs",itau_w,zx_tmp_3d,
443     .                                   iim*jjmp1*klev,ndex3d)
444c
445c K/s
446      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)
447      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
448      CALL histwrite(nid_mth,"dtswr",itau_w,zx_tmp_3d,
449     .                                   iim*jjmp1*klev,ndex3d)
450c
451c K/s     
452      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)
453      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
454      CALL histwrite(nid_mth,"dtlwr",itau_w,zx_tmp_3d,
455     .                                   iim*jjmp1*klev,ndex3d)
456c K/s     
457c     zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
458c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
459c     CALL histwrite(nid_mth,"dtec",itau_w,zx_tmp_3d,
460c    .                                   iim*jjmp1*klev,ndex3d)
461c
462      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
463      CALL histwrite(nid_mth,"duvdf",itau_w,zx_tmp_3d,
464     .                                   iim*jjmp1*klev,ndex3d)
465c
466      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_dyn, zx_tmp_3d)
467      CALL histwrite(nid_mth,"dudyn",itau_w,zx_tmp_3d,
468     .                                   iim*jjmp1*klev,ndex3d)
469c
470c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
471c     CALL histwrite(nid_mth,"dvvdf",itau_w,zx_tmp_3d,
472c    .                                   iim*jjmp1*klev,ndex3d)
473c
474      ENDIF !lev_histmth.GE.4
475c
476c-------------------------------------------------------
477      IF(lev_histmth.GE.5) THEN
478c
479c
480c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxu , zx_tmp_2d)
481c      CALL histwrite(nid_mth,"taux_",itau_w,
482c    $      zx_tmp_2d,iim*jjmp1,ndex2d) 
483c     
484c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxv , zx_tmp_2d)
485c      CALL histwrite(nid_mth,"tauy_",itau_w,
486c    $      zx_tmp_2d,iim*jjmp1,ndex2d)
487c
488c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
489c     CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
490c
491c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
492c     CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
493c
494      ENDIF !lev_histmth.GE.5
495c-------------------------------------------------------
496c
497      if (ok_sync) then
498        call histsync(nid_mth)
499      endif
500
501      ENDIF
Note: See TracBrowser for help on using the repository browser.