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

Last change on this file since 105 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.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         if (iflag_trac.eq.1) THEN
107          if (microfi.eq.1) then
108           DO iq=1,nmicro
109         CALL histdef(nid_day, tname(iq), ttext(iq), "n/m2",
110     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
111     .                "ave(X)", zsto,zout)
112           ENDDO
113          endif
114          if (nmicro.lt.nqmax) then
115           DO iq=nmicro+1,nqmax
116         CALL histdef(nid_day, tname(iq), ttext(iq), "ppm",
117     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
118     .                "ave(X)", zsto,zout)
119           ENDDO
120          endif
121         endif
122c
123      ENDIF !lev_histday.GE.2
124c
125c-------------------------------------------------------
126      IF(lev_histday.GE.3) THEN
127c
128cccccccccccccccccc  Radiative transfer
129c
130c 2D
131c
132         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
133     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
134     .                "ave(X)", zsto1,zout)
135c
136         CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
137     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
138     .                "ave(X)", zsto1,zout)
139c
140         CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
141     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
142     .                "ave(X)", zsto1,zout)
143c
144c 3D
145c
146         CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
147     .                iim,jjmp1,nhori, klev,1,klev,nvert,
148     .                32, "ave(X)", zsto1,zout)
149c
150         CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
151     .                iim,jjmp1,nhori, klev,1,klev,nvert,
152     .                32, "ave(X)", zsto1,zout)
153c
154      ENDIF !lev_histday.GE.3
155c
156c-------------------------------------------------------
157      IF(lev_histday.GE.4) THEN
158c
159         CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
160     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
161     .                "ave(X)", zsto,zout)
162c
163         CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
164     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
165     .                "ave(X)", zsto,zout)
166c
167         CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
168     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
169     .                "ave(X)", zsto,zout)
170c
171         CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
172     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
173     .                "ave(X)", zsto,zout)
174c
175         CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
176     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
177     .                "ave(X)", zsto,zout)
178c
179         CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
180     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
181     .                "ave(X)", zsto,zout)
182c
183         CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
184     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
185     .                "ave(X)", zsto,zout)
186c
187      ENDIF !lev_histday.GE.4
188c
189c-------------------------------------------------------
190      IF(lev_histday.GE.5) THEN
191c
192c
193c        call histdef(nid_day, "taux", 
194c    $         "Zonal wind stress", "Pa", 
195c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
196c    $         "ave(X)", zsto,zout)
197c
198c        call histdef(nid_day, "tauy", 
199c    $         "Meridional xind stress", "Pa", 
200c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
201c    $         "ave(X)", zsto,zout)
202c
203c        CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
204c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
205c    .                "ave(X)", zsto,zout)
206c
207c        CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
208c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
209c    .                "ave(X)", zsto,zout)
210c
211      ENDIF !lev_histday.GE.5
212c-------------------------------------------------------
213c
214         CALL histend(nid_day)
215c
216         ndex2d = 0
217         ndex3d = 0
218c
219      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.