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

Last change on this file since 97 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.9 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.
12
13c
14         zsto = dtime
15         zout = dtime * FLOAT(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)
[97]26      CALL histwrite(nid_day,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]27C
28      i = NINT(zout/zsto)
29      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
[97]30      CALL histwrite(nid_day,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]31c
32      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ftsol,zx_tmp_2d)
[97]33      CALL histwrite(nid_day,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]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)
[97]39      CALL histwrite(nid_day,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]40c
41      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
[97]42      CALL histwrite(nid_day,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]43c
44      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
[97]45      CALL histwrite(nid_day,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]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)
[97]55      CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
[3]56     .                                   iim*jjmp1*klev,ndex3d)
57c
58      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
[97]59      CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
[3]60     .                                   iim*jjmp1*klev,ndex3d)
61c
62      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
[97]63      CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
[3]64     .                                   iim*jjmp1*klev,ndex3d)
65c
66      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
[97]67      CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
[3]68     .                                   iim*jjmp1*klev,ndex3d)
69c
70      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
[97]71      CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
[3]72     .                                   iim*jjmp1*klev,ndex3d)
73c
74      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
[97]75      CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
[3]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)
[97]88      CALL histwrite(nid_day,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]89c
90      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
[97]91      CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]92c
93      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
[97]94      CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]95c
96      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
[97]97      CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]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)
[97]103      CALL histwrite(nid_day,"SWnet",itau_w,zx_tmp_3d,
[3]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)
[97]108      CALL histwrite(nid_day,"LWnet",itau_w,zx_tmp_3d,
[3]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)
[97]118      CALL histwrite(nid_day,"dtdyn",itau_w,zx_tmp_3d,
[3]119     .                                   iim*jjmp1*klev,ndex3d)
120c
121c en K/s     
122c      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
[97]123c      CALL histwrite(nid_day,"dtphy",itau_w,zx_tmp_3d,
[3]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)
[97]128      CALL histwrite(nid_day,"dtvdf",itau_w,zx_tmp_3d,
[3]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)
[97]133      CALL histwrite(nid_day,"dtajs",itau_w,zx_tmp_3d,
[3]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)
[97]139      CALL histwrite(nid_day,"dtswr",itau_w,zx_tmp_3d,
[3]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)
[97]145      CALL histwrite(nid_day,"dtlwr",itau_w,zx_tmp_3d,
[3]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)
[97]150c     CALL histwrite(nid_day,"dtec",itau_w,zx_tmp_3d,
[3]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)
[97]155      CALL histwrite(nid_day,"dudyn",itau_w,zx_tmp_3d,
[3]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)
[97]160      CALL histwrite(nid_day,"duvdf",itau_w,zx_tmp_3d,
[3]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)
[97]165      CALL histwrite(nid_day,"duajs",itau_w,zx_tmp_3d,
[3]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)
[97]170      CALL histwrite(nid_day,"dugwo",itau_w,zx_tmp_3d,
[3]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)
[97]175      CALL histwrite(nid_day,"dugwno",itau_w,zx_tmp_3d,
[3]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!!!!!!!!!!!!!!!
[97]181c     zx_tmp_3d=-1.*zx_tmp_3d
182c     CALL histwrite(nid_day,"dvvdf",itau_w,zx_tmp_3d,
[3]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
[97]195      CALL histwrite(nid_day,"trac"//str2,itau_w,zx_tmp_3d,
[3]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,
[97]206c    $      zx_tmp_2d,iim*jjmp1,ndex2d) 
[3]207c     
208c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxv , zx_tmp_2d)
209c      CALL histwrite(nid_day,"tauy_",itau_w,
[97]210c    $      zx_tmp_2d,iim*jjmp1,ndex2d)
[3]211c
212c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
[97]213c     CALL histwrite(nid_day,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]214c
215c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
[97]216c     CALL histwrite(nid_day,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[3]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.