source: trunk/LMDZ.TITAN/libf/phytitan/write_histday.h @ 324

Last change on this file since 324 was 306, checked in by slebonnois, 14 years ago

SLebonnois: correction de bugs dans la physique Titan:

  • effg.F : Z doit etre en km, donc conversion
  • optc*_1pt_2.F : On utilise cfffv11 et plus optfrac Du coup, les fichiers input testag* ne sont plus necessaires.
  • phytrac.F : passage de la tendance aerosols en intensif dans tous les cas


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)
[306]114
115             if (clouds.eq.1) then
116
[175]117c    -------   NB NOY TOT
118               do i=1,klon
119                 do j=1,klev
120                   zx_tmp_fi3d(i,j)= SUM(qaer(i,j,nrad+1:2*nrad))
121                 enddo
122               enddo
123       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
124       CALL histwrite(nid_day,"qnoy",itau_w,zx_tmp_3d,
125     .                                   iim*jjmp1*klev,ndex3d)
126c    -------   V GLA1 TOT
127               do i=1,klon
128                 do j=1,klev
129                   zx_tmp_fi3d(i,j)= SUM(qaer(i,j,2*nrad+1:3*nrad))
130                 enddo
131               enddo
132       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
133       CALL histwrite(nid_day,"qgl1",itau_w,zx_tmp_3d,
134     .                                   iim*jjmp1*klev,ndex3d)
135c    -------   V GLA2 TOT
136               do i=1,klon
137                 do j=1,klev
138                   zx_tmp_fi3d(i,j)= SUM(qaer(i,j,3*nrad+1:4*nrad))
139                 enddo
140               enddo
141       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
142       CALL histwrite(nid_day,"qgl2",itau_w,zx_tmp_3d,
143     .                                   iim*jjmp1*klev,ndex3d)
144c    -------   V GLA3 TOT
145               do i=1,klon
146                 do j=1,klev
147                   zx_tmp_fi3d(i,j)= SUM(qaer(i,j,4*nrad+1:5*nrad))
148                 enddo
149               enddo
150       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
151       CALL histwrite(nid_day,"qgl3",itau_w,zx_tmp_3d,
152     .                                   iim*jjmp1*klev,ndex3d)
153c --------------
154c ----- SATURATION ESP NUAGES
155       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, satch4,zx_tmp_3d)
156       CALL histwrite(nid_day,"ch4sat", itau_w, zx_tmp_3d,
157     .                                   iim*jjmp1*klev,ndex3d)
158
159       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, satc2h6,zx_tmp_3d)
160       CALL histwrite(nid_day,"c2h6sat", itau_w, zx_tmp_3d,
161     .                                   iim*jjmp1*klev,ndex3d)
162
163       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, satc2h2,zx_tmp_3d)
164       CALL histwrite(nid_day,"c2h2sat", itau_w, zx_tmp_3d,
165     .                                   iim*jjmp1*klev,ndex3d)
166c --------------
167c ----- RESERVOIR DE SURFACE
168       CALL gr_fi_ecrit(1, klon,iim,jjmp1,reservoir,zx_tmp_2d)
169       CALL histwrite(nid_day,"reserv",itau_w,zx_tmp_2d,
170     .                        iim*jjmp1,ndex2d)
171c --------------
172c ----- PRECIPITATIONS
173c       -----  CH4
174       CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,1),zx_tmp_2d)
175       CALL histwrite(nid_day,"prech4",itau_w,zx_tmp_2d,
176     .                        iim*jjmp1,ndex2d)
177c       -----  C2H6
178       CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,2),zx_tmp_2d)
179       CALL histwrite(nid_day,"prec2h6",itau_w,zx_tmp_2d,
180     .                        iim*jjmp1,ndex2d)
181c       -----  C2H2
182       CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,3),zx_tmp_2d)
183       CALL histwrite(nid_day,"prec2h2",itau_w,zx_tmp_2d,
184     .                        iim*jjmp1,ndex2d)
185c
186c --------------
187c ----- FLUX GLACE
188       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,1),zx_tmp_3d)
189       CALL histwrite(nid_day,"flxgl1", itau_w, zx_tmp_3d,
190     .                                   iim*jjmp1*klev,ndex3d)
191       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,2),zx_tmp_3d)
192       CALL histwrite(nid_day,"flxgl2", itau_w, zx_tmp_3d,
193     .                                   iim*jjmp1*klev,ndex3d)
194       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,3),zx_tmp_3d)
195       CALL histwrite(nid_day,"flxgl3", itau_w, zx_tmp_3d,
196     .                                   iim*jjmp1*klev,ndex3d)
197c
198c --------------
199c ----- RAYON MOYEN GOUTTE
200       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rmcloud,zx_tmp_3d)
201       CALL histwrite(nid_day,"rcldbar", itau_w, zx_tmp_3d,
202     .                                   iim*jjmp1*klev,ndex3d)
203c
204             endif
[3]205          endif
[175]206c
207c --------------
208c ----- TRACEURS CHIMIQUES
[3]209          if (nmicro.lt.nqmax) then
210           DO iq=nmicro+1,nqmax
211       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,tr_seri(1,1,iq),zx_tmp_3d)
[110]212       CALL histwrite(nid_day,tname(iq),itau_w,zx_tmp_3d,
[3]213     .                                   iim*jjmp1*klev,ndex3d)
214           ENDDO
215          endif
216         endif
217c
[175]218      ENDIF !lev_histday.GE.2
219c
220c-------------------------------------------------------
221      IF(lev_histday.GE.3) THEN
222c
[3]223cccccccccccccccccc  Radiative transfer
224c
225c 2D
226c
227      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
[97]228      CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]229c
230      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
[97]231      CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]232c
233      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
[97]234      CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]235c
236c 3D
237c
238      zx_tmp_fi3d(1:klon,1:klev)=swnet(1:klon,1:klev)
239      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[97]240      CALL histwrite(nid_day,"SWnet",itau_w,zx_tmp_3d,
[3]241     .                                   iim*jjmp1*klev,ndex3d)
242c
243      zx_tmp_fi3d(1:klon,1:klev)=lwnet(1:klon,1:klev)
244      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[97]245      CALL histwrite(nid_day,"LWnet",itau_w,zx_tmp_3d,
[3]246     .                                   iim*jjmp1*klev,ndex3d)
247c
[175]248c --------------
249c ----- OPACITE BRUME
[106]250       do k=7,NSPECV,10
251         do i=1,klon
252         do l=1,klev
253           t_tauhvd(i,l)=TAUHVD(i,klev-l+1,k)
254         enddo
255         enddo
256         write(str1,'(i2.2)') k
257      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
258      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[110]259      CALL histwrite(nid_day,"thv"//str1,itau_w,zx_tmp_3d,
[106]260     .                                   iim*jjmp1*klev,ndex3d)
261       enddo      ! fin boucle NSPECV
262
[175]263       do k=8,NSPECI,10
264         do i=1,klon
265         do l=1,klev
266           t_tauhvd(i,l)=TAUHID(i,klev-l+1,k)
267         enddo
268         enddo
269         write(str1,'(i2.2)') k
270      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
271      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
272      CALL histwrite(nid_day,"thi"//str1,itau_w,zx_tmp_3d,
273     .                                   iim*jjmp1*klev,ndex3d)
274       enddo      ! fin boucle NSPECI
275c
276c --------------
277c ----- EXTINCTION BRUME
[106]278       do k=7,NSPECV,10
279         do i=1,klon
280         do l=1,klev
[175]281          if(l.ne.klev)
282     s     t_khvd(i,l)=TAUHVD(i,klev-l+1,k)
283     s                -TAUHVD(i,klev-l+1-1,k)
[106]284          if(l.eq.klev)
[175]285     s     t_khvd(i,l)=TAUHVD(i,klev-l+1,k)
[106]286
287         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
288         enddo
289         enddo
290         write(str1,'(i2.2)') k
291      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
292      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[110]293      CALL histwrite(nid_day,"khv"//str1,itau_w,zx_tmp_3d,
[106]294     .                                   iim*jjmp1*klev,ndex3d)
295       enddo      ! fin boucle NSPECV
296
[175]297       do k=8,NSPECI,10
[106]298         do i=1,klon
299         do l=1,klev
[175]300          if(l.ne.klev)
301     s     t_khvd(i,l)=TAUHID(i,klev-l+1,k)
302     s                -TAUHID(i,klev-l+1-1,k)
303          if(l.eq.klev)
304     s     t_khvd(i,l)=TAUHID(i,klev-l+1,k)
305
306         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
[106]307         enddo
308         enddo
309         write(str1,'(i2.2)') k
[175]310      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
[106]311      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[175]312      CALL histwrite(nid_day,"khi"//str1,itau_w,zx_tmp_3d,
[106]313     .                                   iim*jjmp1*klev,ndex3d)
[175]314       enddo      ! fin boucle NSPECI
315c
316c --------------
317c ----- OPACITE GAZ
[106]318       do k=7,NSPECV,10
319         do i=1,klon
320         do l=1,klev
[175]321           t_tauhvd(i,l)=TAUGVD(i,klev-l+1,k)
[106]322         enddo
323         enddo
324         write(str1,'(i2.2)') k
[175]325      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
[106]326      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[175]327      CALL histwrite(nid_day,"tgv"//str1,itau_w,zx_tmp_3d,
[106]328     .                                   iim*jjmp1*klev,ndex3d)
329       enddo      ! fin boucle NSPECV
330
331       do k=8,NSPECI,10
332         do i=1,klon
333         do l=1,klev
[175]334           t_tauhvd(i,l)=TAUGID(i,klev-l+1,k)
[106]335         enddo
336         enddo
337         write(str1,'(i2.2)') k
338      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
339      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[175]340      CALL histwrite(nid_day,"tgi"//str1,itau_w,zx_tmp_3d,
[106]341     .                                   iim*jjmp1*klev,ndex3d)
342       enddo      ! fin boucle NSPECI
[175]343c
344c --------------
345c ----- EXTINCTION GAZ
346       do k=7,NSPECV,10
[106]347         do i=1,klon
348         do l=1,klev
[175]349          if(l.ne.klev)
350     s     t_khvd(i,l)=TAUGVD(i,klev-l+1,k)
351     s                -TAUGVD(i,klev-l+1-1,k)
[106]352          if(l.eq.klev)
[175]353     s     t_khvd(i,l)=TAUGVD(i,klev-l+1,k)
[106]354
355         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
356         enddo
357         enddo
358         write(str1,'(i2.2)') k
359      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
360      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[175]361      CALL histwrite(nid_day,"kgv"//str1,itau_w,zx_tmp_3d,
[106]362     .                                   iim*jjmp1*klev,ndex3d)
[175]363       enddo      ! fin boucle NSPECV
[106]364
365       do k=8,NSPECI,10
366         do i=1,klon
367         do l=1,klev
[175]368          if(l.ne.klev)
369     s     t_khvd(i,l)=TAUGID(i,klev-l+1,k)
370     s                -TAUGID(i,klev-l+1-1,k)
[106]371
372          if(l.eq.klev)
[175]373     s     t_khvd(i,l)=TAUGID(i,klev-l+1,k)
[106]374
375         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
376         enddo
377         enddo
378         write(str1,'(i2.2)') k
379      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
380      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[110]381      CALL histwrite(nid_day,"kgi"//str1,itau_w,zx_tmp_3d,
[106]382     .                                   iim*jjmp1*klev,ndex3d)
383       enddo      ! fin boucle NSPECI
384
[175]385c --------------
386c ----- OPACITE NUAGES (ATTENTION PROXY)
387         if (clouds.eq.1) then
388           zx_tmp_fi3d(1:klon,1:klev)=occcld(1:klon,1:klev)
389           CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
390           CALL histwrite(nid_day,"tcld",itau_w,zx_tmp_3d,
391     .                                   iim*jjmp1*klev,ndex3d)
392c --------------
393c ----- EXTINCTION NUAGES (ATTENTION PROXY)
394           do i=1,klon
395             t_kcld(i,klev)=occcld(i,klev)
396     .       /(zzlev(i,klev+1)-zzlev(i,klev))
397             do j=klev-1,1,-1
398               t_kcld(i,j)=(occcld(i,j)-occcld(i,j+1))
399     .         /(zzlev(i,j+1)-zzlev(i,j))
400             enddo
401           enddo
402           zx_tmp_fi3d(1:klon,1:klev)=t_kcld(1:klon,1:klev)
403           CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
404           CALL histwrite(nid_day,"kcld",itau_w,zx_tmp_3d,
405     .                                   iim*jjmp1*klev,ndex3d)
406        endif 
407c
[3]408      ENDIF !lev_histday.GE.3
409c
410c-------------------------------------------------------
411      IF(lev_histday.GE.4) THEN
412c
413      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
[97]414      CALL histwrite(nid_day,"dtdyn",itau_w,zx_tmp_3d,
[3]415     .                                   iim*jjmp1*klev,ndex3d)
416c
417      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
[97]418      CALL histwrite(nid_day,"dtphy",itau_w,zx_tmp_3d,
[3]419     .                                   iim*jjmp1*klev,ndex3d)
420c K/s
421      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)
422      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[97]423      CALL histwrite(nid_day,"dtvdf",itau_w,zx_tmp_3d,
[3]424     .                                   iim*jjmp1*klev,ndex3d)
425c
426c K/s
427      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)
428      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[97]429      CALL histwrite(nid_day,"dtajs",itau_w,zx_tmp_3d,
[3]430     .                                   iim*jjmp1*klev,ndex3d)
431c
432c K/s
433      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)
434      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[97]435      CALL histwrite(nid_day,"dtswr",itau_w,zx_tmp_3d,
[3]436     .                                   iim*jjmp1*klev,ndex3d)
437c
438c K/s     
439      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)
440      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
[97]441      CALL histwrite(nid_day,"dtlwr",itau_w,zx_tmp_3d,
[3]442     .                                   iim*jjmp1*klev,ndex3d)
443c K/s     
[106]444c     zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
445c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
446c     CALL histwrite(nid_day,"dtec",itau_w,zx_tmp_3d,
447c    .                                   iim*jjmp1*klev,ndex3d)
[3]448c
449      ENDIF !lev_histday.GE.4
450c
451c-------------------------------------------------------
452      IF(lev_histday.GE.5) THEN
453c
454c
455c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxu , zx_tmp_2d)
456c      CALL histwrite(nid_day,"taux_",itau_w,
[97]457c    $      zx_tmp_2d,iim*jjmp1,ndex2d) 
[3]458c     
459c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxv , zx_tmp_2d)
460c      CALL histwrite(nid_day,"tauy_",itau_w,
[97]461c    $      zx_tmp_2d,iim*jjmp1,ndex2d)
[3]462c
463c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
[97]464c     CALL histwrite(nid_day,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]465c
466c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
[97]467c     CALL histwrite(nid_day,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]468c
469      ENDIF !lev_histday.GE.5
470c-------------------------------------------------------
471c
472      if (ok_sync) then
473        call histsync(nid_day)
474      endif
475
476      ENDIF
Note: See TracBrowser for help on using the repository browser.