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