source: trunk/libf/phytitan/ini_histday.h @ 107

Last change on this file since 107 was 106, checked in by slebonnois, 14 years ago

SL : mise a jour de phytitan pour etre conforme aux sources actuelles
utilisees sur gnome.
File size: 9.7 KB
Line 
1!
2! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ini_histday.h,v 1.2 2005/01/27 10:06:12 fairhead Exp $
3!
4      IF (ok_journe) THEN
5c
6         zsto = dtime
7         zout = dtime * FLOAT(ecrit_day)
8c zsto1: pour des flux radiatifs calcules tous les radpas appels physiq
9         zsto1= dtime * FLOAT(radpas)
10c
11         idayref = day_ref
12         CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
13c
14         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
15         DO i = 1, iim
16            zx_lon(i,1) = rlond(i+1)
17            zx_lon(i,jjmp1) = rlond(i+1)
18         ENDDO
19         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)
20         CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
21     .                 1,iim,1,jjmp1, itau_phy, zjulian, zsto, 
22     .                 nhori, nid_day)
23         write(*,*)'Journee ', itau_phy, zjulian
24         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
25     .                 klev, presnivs, nvert)
26
27c-------------------------------------------------------
28      IF(lev_histday.GE.1) THEN
29c
30ccccccccccccc 2D fields, invariables
31c
32         CALL histdef(nid_day, "phis", "Surface geop. height", "-",
33     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
34     .                "once",  zsto,zout)
35c
36         CALL histdef(nid_day, "aire", "Grid area", "-",
37     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
38     .                "once",  zsto,zout)
39c
40ccccccc axe Ls
41         CALL histdef(nid_day, "ls", "Solar longitude", "degrees",
42     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
43     .                "ave(X)", zsto,zout)
44c
45ccccccccccccc 2D fields, variables
46c
47         CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
48     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
49     .                "ave(X)", zsto,zout)
50c
51         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
52     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
53     .                "ave(X)", zsto,zout)
54c
55c        CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
56c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
57c    .                "ave(X)", zsto,zout)
58c
59c        CALL histdef(nid_day, "ve", "Merid energy transport", "-",
60c     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
61c     .                "ave(X)", zsto,zout)
62c
63      ENDIF !lev_histday.GE.1
64c
65c-------------------------------------------------------
66      IF(lev_histday.GE.2) THEN
67c
68ccccccccccccc 3D fields, basics
69c
70         CALL histdef(nid_day, "temp", "Air temperature", "K",
71     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
72     .                "ave(X)", zsto,zout)
73c
74         CALL histdef(nid_day, "pres", "Air pressure", "Pa",
75     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
76     .                "ave(X)", zsto,zout)
77c
78         CALL histdef(nid_day, "geop", "Geopotential height", "m",
79     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
80     .                "ave(X)", zsto,zout)
81c
82         CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
83     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
84     .                "ave(X)", zsto,zout)
85c
86         CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
87     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
88     .                "ave(X)", zsto,zout)
89c
90         CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
91     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
92     .                "ave(X)", zsto,zout)
93c
94         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
95     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
96     .                "ave(X)", zsto1,zout)
97c
98         CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2",
99     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
100     .                "ave(X)", zsto,zout)
101c
102         CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2",
103     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
104     .                "ave(X)", zsto,zout)
105c
106      ENDIF !lev_histday.GE.2
107c
108c-------------------------------------------------------
109      IF(lev_histday.GE.3) THEN
110c
111cccccccccccccccccc  Tracers
112c
113         if (iflag_trac.eq.1) THEN
114          if (microfi.eq.1) then
115           DO iq=1,nmicro
116         CALL histdef(nid_day, tname(iq), ttext(iq), "n/m2",
117     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
118     .                "ave(X)", zsto,zout)
119           ENDDO
120          endif
121          if (nmicro.lt.nqmax) then
122           DO iq=nmicro+1,nqmax
123         CALL histdef(nid_day, tname(iq), ttext(iq), "ppm",
124     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
125     .                "ave(X)", zsto,zout)
126           ENDDO
127          endif
128         endif
129c
130cccccccccccccccccc  Radiative transfer
131c
132c 2D
133c
134         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
135     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
136     .                "ave(X)", zsto1,zout)
137c
138         CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
139     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
140     .                "ave(X)", zsto1,zout)
141c
142         CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
143     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
144     .                "ave(X)", zsto1,zout)
145c
146c 3D
147c
148         CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
149     .                iim,jjmp1,nhori, klev,1,klev,nvert,
150     .                32, "ave(X)", zsto1,zout)
151c
152         CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
153     .                iim,jjmp1,nhori, klev,1,klev,nvert,
154     .                32, "ave(X)", zsto1,zout)
155c
156         DO k=7,NSPECV,10
157           write(str1,'(i2.2)') k
158         CALL histdef(nid_day,"thv"//str1,"Haze Opa Vis",
159     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
160     .                "ave(X)",zsto1,zout)
161         ENDDO
162c
163         DO k=7,NSPECV,10
164           write(str1,'(i2.2)') k
165         CALL histdef(nid_day,"khv"//str1,"Haze ext Vis ",
166     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
167     .                "ave(X)",zsto1,zout)
168         ENDDO
169c
170         DO k=7,NSPECV,10
171           write(str1,'(i2.2)') k
172         CALL histdef(nid_day,"tgv"//str1,"Haze Opa Vis",
173     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
174     .                "ave(X)",zsto1,zout)
175         ENDDO
176c
177         DO k=7,NSPECV,10
178           write(str1,'(i2.2)') k
179         CALL histdef(nid_day,"kgv"//str1,"Haze ext Vis ",
180     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
181     .                "ave(X)",zsto1,zout)
182         ENDDO
183c
184         DO k=8,NSPECI,10
185           write(str1,'(i2.2)') k
186         CALL histdef(nid_day,"thi"//str1,"Haze Opa IR",
187     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
188     .                "ave(X)",zsto1,zout)
189         ENDDO
190c
191         DO k=8,NSPECI,10
192           write(str1,'(i2.2)') k
193         CALL histdef(nid_day,"khi"//str1,"Haze ext IR ",
194     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
195     .                "ave(X)",zsto1,zout)
196         ENDDO
197c
198         DO k=8,NSPECI,10
199           write(str1,'(i2.2)') k
200         CALL histdef(nid_day,"tgi"//str1,"Haze Opa IR",
201     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
202     .                "ave(X)",zsto1,zout)
203         ENDDO
204c
205         DO k=8,NSPECI,10
206           write(str1,'(i2.2)') k
207         CALL histdef(nid_day,"kgi"//str1,"Haze ext IR ",
208     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
209     .                "ave(X)",zsto1,zout)
210         ENDDO
211c
212      ENDIF !lev_histday.GE.3
213c
214c-------------------------------------------------------
215      IF(lev_histday.GE.4) THEN
216c
217         CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
218     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
219     .                "ave(X)", zsto,zout)
220c
221         CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
222     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
223     .                "ave(X)", zsto,zout)
224c
225         CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
226     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
227     .                "ave(X)", zsto,zout)
228c
229         CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
230     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
231     .                "ave(X)", zsto,zout)
232c
233         CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
234     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
235     .                "ave(X)", zsto,zout)
236c
237         CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
238     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
239     .                "ave(X)", zsto,zout)
240c
241c        CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
242c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
243c    .                "ave(X)", zsto,zout)
244c
245      ENDIF !lev_histday.GE.4
246c
247c-------------------------------------------------------
248      IF(lev_histday.GE.5) THEN
249c
250c
251c        call histdef(nid_day, "taux", 
252c    $         "Zonal wind stress", "Pa", 
253c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
254c    $         "ave(X)", zsto,zout)
255c
256c        call histdef(nid_day, "tauy", 
257c    $         "Meridional xind stress", "Pa", 
258c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
259c    $         "ave(X)", zsto,zout)
260c
261c        CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
262c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
263c    .                "ave(X)", zsto,zout)
264c
265c        CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
266c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
267c    .                "ave(X)", zsto,zout)
268c
269      ENDIF !lev_histday.GE.5
270c-------------------------------------------------------
271c
272         CALL histend(nid_day)
273c
274         ndex2d = 0
275         ndex3d = 0
276c
277      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.