source: trunk/libf/phyvenus/write_histday.h @ 86

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