source: trunk/LMDZ.VENUS/libf/phyvenus/write_histday.h @ 808

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