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

Last change on this file since 524 was 474, checked in by slebonnois, 13 years ago

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