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

Last change on this file since 102 was 97, checked in by slebonnois, 14 years ago

Serie de modifs SL pour homogeneisation des phytitan et phyvenus
Ca touche aussi aux liens phy/dyn (surtout a propos de clesphy0),
a verifier avec les autres, donc...
File size: 7.4 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.
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)
24      CALL histwrite(nid_day,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
25C
26      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
27      CALL histwrite(nid_day,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
28c
29ccccccc axe Ls
30      do j=1,jjmp1
31       do i=1,iim
32        zx_tmp_2d(i,j)=zls*180./RPI      ! zls est en radians !!
33       enddo
34      enddo
35      CALL histwrite(nid_day,"ls",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
36c
37ccccccccccccc 2D fields, variables
38c
39      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ftsol,zx_tmp_2d)
40      CALL histwrite(nid_day,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
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)
46      CALL histwrite(nid_day,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
47c
48c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
49c     CALL histwrite(nid_day,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
50c
51c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
52c     CALL histwrite(nid_day,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
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)
62      CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
63     .                                   iim*jjmp1*klev,ndex3d)
64c
65      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
66      CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
67     .                                   iim*jjmp1*klev,ndex3d)
68c
69      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
70      CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
71     .                                   iim*jjmp1*klev,ndex3d)
72c
73      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
74      CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
75     .                                   iim*jjmp1*klev,ndex3d)
76c
77      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
78      CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
79     .                                   iim*jjmp1*klev,ndex3d)
80c
81      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
82      CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
83     .                                   iim*jjmp1*klev,ndex3d)
84c
85      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
86      CALL histwrite(nid_day,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
87c
88      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_dyn, zx_tmp_3d)
89      CALL histwrite(nid_day,"dudyn",itau_w,zx_tmp_3d,
90     .                                   iim*jjmp1*klev,ndex3d)
91c
92      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
93      CALL histwrite(nid_day,"duvdf",itau_w,zx_tmp_3d,
94     .                                   iim*jjmp1*klev,ndex3d)
95c
96         if (iflag_trac.eq.1) THEN
97          if (microfi.eq.1) then
98           DO iq=1,nmicro
99       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qaer(1,1,iq), zx_tmp_3d)
100       CALL histwrite(nid_day,tname(iq),itau_w,zx_tmp_3d,
101     .                                   iim*jjmp1*klev,ndex3d)
102           ENDDO
103          endif
104          if (nmicro.lt.nqmax) then
105           DO iq=nmicro+1,nqmax
106       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,tr_seri(1,1,iq),zx_tmp_3d)
107       CALL histwrite(nid_day,tname(iq),itau_w,zx_tmp_3d,
108     .                                   iim*jjmp1*klev,ndex3d)
109           ENDDO
110          endif
111         endif
112c
113      ENDIF !lev_histday.GE.2
114c
115c-------------------------------------------------------
116      IF(lev_histday.GE.3) THEN
117c
118cccccccccccccccccc  Radiative transfer
119c
120c 2D
121c
122      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
123      CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
124c
125      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
126      CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
127c
128      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
129      CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
130c
131c 3D
132c
133      zx_tmp_fi3d(1:klon,1:klev)=swnet(1:klon,1:klev)
134      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
135      CALL histwrite(nid_day,"SWnet",itau_w,zx_tmp_3d,
136     .                                   iim*jjmp1*klev,ndex3d)
137c
138      zx_tmp_fi3d(1:klon,1:klev)=lwnet(1:klon,1:klev)
139      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
140      CALL histwrite(nid_day,"LWnet",itau_w,zx_tmp_3d,
141     .                                   iim*jjmp1*klev,ndex3d)
142c
143      ENDIF !lev_histday.GE.3
144c
145c-------------------------------------------------------
146      IF(lev_histday.GE.4) THEN
147c
148      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
149      CALL histwrite(nid_day,"dtdyn",itau_w,zx_tmp_3d,
150     .                                   iim*jjmp1*klev,ndex3d)
151c
152      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
153      CALL histwrite(nid_day,"dtphy",itau_w,zx_tmp_3d,
154     .                                   iim*jjmp1*klev,ndex3d)
155c K/s
156      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)
157      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
158      CALL histwrite(nid_day,"dtvdf",itau_w,zx_tmp_3d,
159     .                                   iim*jjmp1*klev,ndex3d)
160c
161c K/s
162      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)
163      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
164      CALL histwrite(nid_day,"dtajs",itau_w,zx_tmp_3d,
165     .                                   iim*jjmp1*klev,ndex3d)
166c
167c K/s
168      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)
169      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
170      CALL histwrite(nid_day,"dtswr",itau_w,zx_tmp_3d,
171     .                                   iim*jjmp1*klev,ndex3d)
172c
173c K/s     
174      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)
175      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
176      CALL histwrite(nid_day,"dtlwr",itau_w,zx_tmp_3d,
177     .                                   iim*jjmp1*klev,ndex3d)
178c K/s     
179      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
180      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
181      CALL histwrite(nid_day,"dtec",itau_w,zx_tmp_3d,
182     .                                   iim*jjmp1*klev,ndex3d)
183c
184      ENDIF !lev_histday.GE.4
185c
186c-------------------------------------------------------
187      IF(lev_histday.GE.5) THEN
188c
189c
190c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxu , zx_tmp_2d)
191c      CALL histwrite(nid_day,"taux_",itau_w,
192c    $      zx_tmp_2d,iim*jjmp1,ndex2d) 
193c     
194c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxv , zx_tmp_2d)
195c      CALL histwrite(nid_day,"tauy_",itau_w,
196c    $      zx_tmp_2d,iim*jjmp1,ndex2d)
197c
198c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
199c     CALL histwrite(nid_day,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
200c
201c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
202c     CALL histwrite(nid_day,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
203c
204      ENDIF !lev_histday.GE.5
205c-------------------------------------------------------
206c
207      if (ok_sync) then
208        call histsync(nid_day)
209      endif
210
211      ENDIF
Note: See TracBrowser for help on using the repository browser.