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

Last change on this file since 902 was 808, checked in by slebonnois, 13 years ago

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