source: trunk/LMDZ.VENUS/libf/phyvenus/ini_histday.h @ 777

Last change on this file since 777 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: 8.1 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)
8         zsto1= dtime * FLOAT(ecrit_day)
9c
10         idayref = day_ref
11         CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
12c
13         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
14         DO i = 1, iim
15            zx_lon(i,1) = rlond(i+1)
16            zx_lon(i,jjmp1) = rlond(i+1)
17         ENDDO
18         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)
19
20c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
21c    d'ou le -1.* devant zx_lon et zx_lat       
22         CALL histbeg("histday.nc", iim,-1.*zx_lon(:,1), jjmp1,
23     .                 -1.*zx_lat(1,:),
24     .                 1,iim,1,jjmp1, itau_phy, zjulian, zsto, 
25     .                 nhori, nid_day)
26         write(*,*)'Journee ', itau_phy, zjulian
27         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
28     .                 klev, presnivs, nvert)
29
30c-------------------------------------------------------
31      IF(lev_histday.GE.1) THEN
32c
33ccccccccccccc 2D fields, basics
34c
35         CALL histdef(nid_day, "phis", "Surface geop. height", "-",
36     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
37     .                "once",  zsto,zout)
38c
39         CALL histdef(nid_day, "aire", "Grid area", "-",
40     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
41     .                "once",  zsto,zout)
42c
43         CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
44     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
45     .                "ave(X)", zsto,zout)
46c
47         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
48     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
49     .                "ave(X)", zsto,zout)
50c
51         CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
52     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
53     .                "ave(X)", zsto,zout)
54c
55         CALL histdef(nid_day, "ve", "Merid energy transport", "-",
56     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
57     .                "ave(X)", zsto,zout)
58c
59      ENDIF !lev_histday.GE.1
60c
61c-------------------------------------------------------
62      IF(lev_histday.GE.2) THEN
63c
64ccccccccccccc 3D fields, basics
65c
66         CALL histdef(nid_day, "temp", "Air temperature", "K",
67     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
68     .                "ave(X)", zsto,zout)
69c
70         CALL histdef(nid_day, "pres", "Air pressure", "Pa",
71     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
72     .                "ave(X)", zsto,zout)
73c
74         CALL histdef(nid_day, "geop", "Geopotential height", "m",
75     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
76     .                "ave(X)", zsto,zout)
77c
78         CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
79     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
80     .                "ave(X)", zsto,zout)
81c
82         CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
83     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
84     .                "ave(X)", zsto,zout)
85c
86         CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
87     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
88     .                "ave(X)", zsto,zout)
89c
90      ENDIF !lev_histday.GE.2
91c
92c-------------------------------------------------------
93      IF(lev_histday.GE.3) THEN
94c
95cccccccccccccccccc  Radiative transfer
96c
97c 2D
98c
99         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
100     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
101     .                "ave(X)", zsto1,zout)
102c
103         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
104     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
105     .                "ave(X)", zsto1,zout)
106c
107         CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
108     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
109     .                "ave(X)", zsto1,zout)
110c
111         CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
112     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
113     .                "ave(X)", zsto1,zout)
114c
115c 3D
116c
117         CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
118     .                iim,jjmp1,nhori, klev,1,klev,nvert,
119     .                32, "ave(X)", zsto1,zout)
120c
121         CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
122     .                iim,jjmp1,nhori, klev,1,klev,nvert,
123     .                32, "ave(X)", zsto1,zout)
124c
125      ENDIF !lev_histday.GE.3
126c
127c-------------------------------------------------------
128      IF(lev_histday.GE.4) THEN
129c
130         CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
131     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
132     .                "ave(X)", zsto,zout)
133c
134c         CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
135c     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
136c     .                "ave(X)", zsto,zout)
137c
138         CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
139     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
140     .                "ave(X)", zsto,zout)
141c
142         CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
143     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
144     .                "ave(X)", zsto,zout)
145c
146         CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
147     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
148     .                "ave(X)", zsto,zout)
149c
150         CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
151     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
152     .                "ave(X)", zsto,zout)
153c
154c         CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
155c     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
156c     .                "ave(X)", zsto,zout)
157c
158         CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2",
159     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
160     .                "ave(X)", zsto,zout)
161c
162         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
163     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
164     .                "ave(X)", zsto,zout)
165c
166         CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2",
167     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
168     .                "ave(X)", zsto,zout)
169c
170         CALL histdef(nid_mth, "dugwo", "GW oro dU", "m/s2",
171     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
172     .                "ave(X)", zsto,zout)
173c
174         CALL histdef(nid_mth, "dugwno", "GW non-oro dU", "m/s2",
175     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
176     .                "ave(X)", zsto,zout)
177c
178c        CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
179c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
180c    .                "ave(X)", zsto,zout)
181c
182      ENDIF !lev_histday.GE.4
183c
184c-------------------------------------------------------
185      IF(lev_histday.GE.5) THEN
186c
187         if (nqmax.GE.3) THEN
188           DO iq=1,nqmax-2
189            IF (iq.LE.99) THEN
190         WRITE(str2,'(i2.2)') iq
191         CALL histdef(nid_day, "trac"//str2, "Tracer No."//str2, "-",
192     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
193     .                "ave(X)", zsto,zout)
194            ELSE
195         PRINT*, "Trop de traceurs"
196         CALL abort
197            ENDIF
198           ENDDO
199         endif
200c
201c        call histdef(nid_day, "taux", 
202c    $         "Zonal wind stress", "Pa", 
203c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
204c    $         "ave(X)", zsto,zout)
205c
206c        call histdef(nid_day, "tauy", 
207c    $         "Meridional xind stress", "Pa", 
208c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
209c    $         "ave(X)", zsto,zout)
210c
211c        CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
212c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
213c    .                "ave(X)", zsto,zout)
214c
215c        CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
216c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
217c    .                "ave(X)", zsto,zout)
218c
219      ENDIF !lev_histday.GE.5
220c-------------------------------------------------------
221c
222         CALL histend(nid_day)
223c
224         ndex2d = 0
225         ndex3d = 0
226c
227      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.