source: trunk/libf/phytitan/write_histday.h @ 6

Last change on this file since 6 was 3, checked in by slebonnois, 15 years ago

Creation de repertoires:

  • chantiers : pour communiquer sur nos projets de modifs
  • documentation : pour stocker les docs

Ajout de:

  • libf/phytitan : physique de Titan
  • libf/chimtitan: chimie de Titan
  • libf/phyvenus : physique de Venus
File size: 8.3 KB
Line 
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      z4_tmp_2d = 0.
10      zx_tmp_3d = 0.
11      z4_tmp_3d = 0.
12      zx_tmp_fi2d=0.
13      zx_tmp_fi3d=0.
14c
15         zsto = dtime
16         zout = dtime * FLOAT(ecrit_day)
17         itau_w = itau_phy + itap
18
19c
20c-------------------------------------------------------
21      IF(lev_histday.GE.1) THEN
22c
23ccccccccccccc 2D fields, invariables
24c
25      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
26      z4_tmp_2d=zx_tmp_2d
27      CALL histwrite(nid_day,"phis",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
28C
29      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
30      z4_tmp_2d=zx_tmp_2d
31      CALL histwrite(nid_day,"aire",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
32c
33ccccccc axe Ls
34      do j=1,jjmp1
35       do i=1,iim
36        z4_tmp_2d(i,j)=zls*180./RPI      ! zls est en radians !!
37       enddo
38      enddo
39      CALL histwrite(nid_day,"ls",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
40c
41ccccccccccccc 2D fields, variables
42c
43      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ftsol,zx_tmp_2d)
44      z4_tmp_2d=zx_tmp_2d
45      CALL histwrite(nid_day,"tsol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
46c
47      DO i = 1, klon
48         zx_tmp_fi2d(i) = paprs(i,1)
49      ENDDO
50      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
51      z4_tmp_2d=zx_tmp_2d
52      CALL histwrite(nid_day,"psol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
53c
54c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
55c     z4_tmp_2d=zx_tmp_2d
56c     CALL histwrite(nid_day,"ue",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
57c
58c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
59c     z4_tmp_2d=zx_tmp_2d
60c     CALL histwrite(nid_day,"ve",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
61c
62      ENDIF !lev_histday.GE.1
63c
64c-------------------------------------------------------
65      IF(lev_histday.GE.2) THEN
66c
67ccccccccccccc 3D fields, basics
68c
69      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
70      z4_tmp_3d=zx_tmp_3d
71      CALL histwrite(nid_day,"temp",itau_w,z4_tmp_3d,
72     .                                   iim*jjmp1*klev,ndex3d)
73c
74      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
75      z4_tmp_3d=zx_tmp_3d
76      CALL histwrite(nid_day,"pres",itau_w,z4_tmp_3d,
77     .                                   iim*jjmp1*klev,ndex3d)
78c
79      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
80      z4_tmp_3d=zx_tmp_3d
81      CALL histwrite(nid_day,"geop",itau_w,z4_tmp_3d,
82     .                                   iim*jjmp1*klev,ndex3d)
83c
84      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
85      z4_tmp_3d=zx_tmp_3d
86      CALL histwrite(nid_day,"vitu",itau_w,z4_tmp_3d,
87     .                                   iim*jjmp1*klev,ndex3d)
88c
89      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
90      z4_tmp_3d=zx_tmp_3d
91      CALL histwrite(nid_day,"vitv",itau_w,z4_tmp_3d,
92     .                                   iim*jjmp1*klev,ndex3d)
93c
94      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
95      z4_tmp_3d=zx_tmp_3d
96      CALL histwrite(nid_day,"vitw",itau_w,z4_tmp_3d,
97     .                                   iim*jjmp1*klev,ndex3d)
98c
99      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
100      z4_tmp_2d=zx_tmp_2d
101      CALL histwrite(nid_day,"tops",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
102c
103      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_dyn, zx_tmp_3d)
104      z4_tmp_3d=zx_tmp_3d
105      CALL histwrite(nid_day,"dudyn",itau_w,z4_tmp_3d,
106     .                                   iim*jjmp1*klev,ndex3d)
107c
108      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
109      z4_tmp_3d=zx_tmp_3d
110      CALL histwrite(nid_day,"duvdf",itau_w,z4_tmp_3d,
111     .                                   iim*jjmp1*klev,ndex3d)
112c
113         if (iflag_trac.eq.1) THEN
114          if (microfi.eq.1) then
115           DO iq=1,nmicro
116       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qaer(1,1,iq), zx_tmp_3d)
117       z4_tmp_3d=zx_tmp_3d
118       CALL histwrite(nid_day,tname(iq),itau_w,z4_tmp_3d,
119     .                                   iim*jjmp1*klev,ndex3d)
120           ENDDO
121          endif
122          if (nmicro.lt.nqmax) then
123           DO iq=nmicro+1,nqmax
124       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,tr_seri(1,1,iq),zx_tmp_3d)
125       z4_tmp_3d=zx_tmp_3d
126       CALL histwrite(nid_day,tname(iq),itau_w,z4_tmp_3d,
127     .                                   iim*jjmp1*klev,ndex3d)
128           ENDDO
129          endif
130         endif
131c
132      ENDIF !lev_histday.GE.2
133c
134c-------------------------------------------------------
135      IF(lev_histday.GE.3) THEN
136c
137cccccccccccccccccc  Radiative transfer
138c
139c 2D
140c
141      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
142      z4_tmp_2d=zx_tmp_2d
143      CALL histwrite(nid_day,"topl",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
144c
145      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
146      z4_tmp_2d=zx_tmp_2d
147      CALL histwrite(nid_day,"sols",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
148c
149      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
150      z4_tmp_2d=zx_tmp_2d
151      CALL histwrite(nid_day,"soll",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
152c
153c 3D
154c
155      zx_tmp_fi3d(1:klon,1:klev)=swnet(1:klon,1:klev)
156      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
157      z4_tmp_3d=zx_tmp_3d
158      CALL histwrite(nid_day,"SWnet",itau_w,z4_tmp_3d,
159     .                                   iim*jjmp1*klev,ndex3d)
160c
161      zx_tmp_fi3d(1:klon,1:klev)=lwnet(1:klon,1:klev)
162      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
163      z4_tmp_3d=zx_tmp_3d
164      CALL histwrite(nid_day,"LWnet",itau_w,z4_tmp_3d,
165     .                                   iim*jjmp1*klev,ndex3d)
166c
167      ENDIF !lev_histday.GE.3
168c
169c-------------------------------------------------------
170      IF(lev_histday.GE.4) THEN
171c
172      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
173      z4_tmp_3d=zx_tmp_3d
174      CALL histwrite(nid_day,"dtdyn",itau_w,z4_tmp_3d,
175     .                                   iim*jjmp1*klev,ndex3d)
176c
177      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
178      z4_tmp_3d=zx_tmp_3d
179      CALL histwrite(nid_day,"dtphy",itau_w,z4_tmp_3d,
180     .                                   iim*jjmp1*klev,ndex3d)
181c K/s
182      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)
183      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
184      z4_tmp_3d=zx_tmp_3d
185      CALL histwrite(nid_day,"dtvdf",itau_w,z4_tmp_3d,
186     .                                   iim*jjmp1*klev,ndex3d)
187c
188c K/s
189      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)
190      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
191      z4_tmp_3d=zx_tmp_3d
192      CALL histwrite(nid_day,"dtajs",itau_w,z4_tmp_3d,
193     .                                   iim*jjmp1*klev,ndex3d)
194c
195c K/s
196      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)
197      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
198      z4_tmp_3d=zx_tmp_3d
199      CALL histwrite(nid_day,"dtswr",itau_w,z4_tmp_3d,
200     .                                   iim*jjmp1*klev,ndex3d)
201c
202c K/s     
203      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)
204      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
205      z4_tmp_3d=zx_tmp_3d
206      CALL histwrite(nid_day,"dtlwr",itau_w,z4_tmp_3d,
207     .                                   iim*jjmp1*klev,ndex3d)
208c K/s     
209      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
210      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
211      z4_tmp_3d=zx_tmp_3d
212      CALL histwrite(nid_day,"dtec",itau_w,z4_tmp_3d,
213     .                                   iim*jjmp1*klev,ndex3d)
214c
215      ENDIF !lev_histday.GE.4
216c
217c-------------------------------------------------------
218      IF(lev_histday.GE.5) THEN
219c
220c
221c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxu , zx_tmp_2d)
222c      z4_tmp_2d=zx_tmp_2d
223c      CALL histwrite(nid_day,"taux_",itau_w,
224c    $      z4_tmp_2d,iim*jjmp1,ndex2d) 
225c     
226c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxv , zx_tmp_2d)
227c      z4_tmp_2d=zx_tmp_2d
228c      CALL histwrite(nid_day,"tauy_",itau_w,
229c    $      z4_tmp_2d,iim*jjmp1,ndex2d)
230c
231c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
232c      z4_tmp_2d=zx_tmp_2d
233c     CALL histwrite(nid_day,"cdrm",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
234c
235c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
236c      z4_tmp_2d=zx_tmp_2d
237c     CALL histwrite(nid_day,"cdrh",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
238c
239      ENDIF !lev_histday.GE.5
240c-------------------------------------------------------
241c
242      if (ok_sync) then
243        call histsync(nid_day)
244      endif
245
246      ENDIF
Note: See TracBrowser for help on using the repository browser.