source: trunk/libf/phyvenus/ini_histmth.h @ 98

Last change on this file since 98 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: 9.7 KB
Line 
1!
2! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ini_histmth.h,v 1.4 2005/03/09 12:30:16 fairhead Exp $
3!
4      IF (ok_mensuel) THEN
5c
6         zsto = dtime
7         zout = dtime * ecrit_mth
8c
9c zsto1: pour des valeurs "instantanees" mensuelles
10         zsto1 = dtime * ecrit_mth
11c zsto2: pour des flux radiatifs calcules avec un autre pas de temps
12         zsto2 = dtime * radpas
13         PRINT*,' zsto,zsto1,zsto2,zout=',zsto, zsto1, zsto2,zout
14c
15         idayref = day_ref
16         CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
17c
18         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
19         DO i = 1, iim
20            zx_lon(i,1) = rlond(i+1)
21            zx_lon(i,jjmp1) = rlond(i+1)
22         ENDDO
23         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)
24         
25c VENUS: regardee a l'envers!!!!!!!!!!!!!!!
26c    d'ou le -1.* devant zx_lon et zx_lat       
27         CALL histbeg("histmth.nc", iim,-1.*zx_lon(:,1), jjmp1,
28     .                 -1.*zx_lat(1,:),
29     .                 1,iim,1,jjmp1, itau_phy, zjulian, zsto, 
30     .                 nhori, nid_mth)
31         write(*,*)'Mensuel ', itau_phy, zjulian
32         CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
33     .                 klev, presnivs, nvert)
34
35c-------------------------------------------------------
36      IF(lev_histmth.GE.1) THEN
37c
38ccccccccccccc 2D fields, basics
39c
40         CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
41     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
42     .                "once",  zsto,zout)
43c
44         CALL histdef(nid_mth, "aire", "Grid area", "-",
45     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
46     .                "once",  zsto,zout)
47c
48         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
49     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
50     .                "ave(X)", zsto,zout)
51c
52         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
53     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
54     .                "ave(X)", zsto,zout)
55c
56c        CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
57c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
58c    .                "ave(X)", zsto,zout)
59c
60c        CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
61c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
62c    .                "ave(X)", zsto,zout)
63c
64         CALL histdef(nid_mth, "cdragh", "Drag coef on T", "-",
65     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
66     .                "ave(X)", zsto,zout)
67c
68         CALL histdef(nid_mth, "cdragm", "Drag coef on U", "-",
69     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
70     .                "ave(X)", zsto,zout)
71c
72      ENDIF !lev_histmth.GE.1
73c
74c-------------------------------------------------------
75      IF(lev_histmth.GE.2) THEN
76c
77ccccccccccccc 3D fields, basics
78c
79         CALL histdef(nid_mth, "temp", "Air temperature", "K",
80     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
81     .                "ave(X)", zsto,zout)
82c
83         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
84     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
85     .                "ave(X)", zsto,zout)
86c
87         CALL histdef(nid_mth, "geop", "Geopotential height", "m",
88     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
89     .                "ave(X)", zsto,zout)
90c
91         CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
92     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
93     .                "ave(X)", zsto,zout)
94c
95         CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
96     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
97     .                "ave(X)", zsto,zout)
98c
99         CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
100     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
101     .                "ave(X)", zsto,zout)
102c
103         CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2",
104     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
105     .                "ave(X)", zsto,zout)
106c
107         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
108     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
109     .                "ave(X)", zsto,zout)
110c
111c        CALL histdef(nid_mth, "mang", "Angular momentum", "kg m2/s",
112c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
113c    .                "ave(X)", zsto,zout)
114c
115         CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s",
116     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
117     .                "ave(X)", zsto,zout)
118c
119c plusieurs traceurs
120          if (iflag_trac.eq.1) THEN
121            DO iq=1,nqmax
122             IF (iq.LE.99) THEN
123          WRITE(str2,'(i2.2)') iq
124          CALL histdef(nid_mth, tnom(iq), ttext(iq), "ppm",
125     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
126     .                "ave(X)", zsto,zout)
127             ELSE
128          PRINT*, "Trop de traceurs"
129          CALL abort
130             ENDIF
131            ENDDO
132          endif
133c
134         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
135     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
136     .                "ave(X)", zsto2,zout)
137c
138      ENDIF !lev_histmth.GE.2
139c
140c-------------------------------------------------------
141      IF(lev_histmth.GE.3) THEN
142c
143cccccccccccccccccc  Radiative transfer
144c
145c 2D
146c
147         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
148     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
149     .                "ave(X)", zsto2,zout)
150c
151         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
152     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
153     .                "ave(X)", zsto2,zout)
154c
155         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
156     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
157     .                "ave(X)", zsto2,zout)
158c
159c 3D
160c
161         CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2",
162     .                iim,jjmp1,nhori, klev,1,klev,nvert,
163     .                32, "ave(X)", zsto2,zout)
164c
165         CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2",
166     .                iim,jjmp1,nhori, klev,1,klev,nvert,
167     .                32, "ave(X)", zsto2,zout)
168c
169         CALL histdef(nid_mth, "fluxvdf", "PBL net flux","W/m2",
170     .                iim,jjmp1,nhori, klev,1,klev,nvert,
171     .                32, "ave(X)", zsto,zout)
172c
173         CALL histdef(nid_mth, "fluxdyn", "Dyn. net flux","W/m2",
174     .                iim,jjmp1,nhori, klev,1,klev,nvert,
175     .                32, "ave(X)", zsto,zout)
176c
177         CALL histdef(nid_mth, "fluxajs", "Dry adj. net flux","W/m2",
178     .                iim,jjmp1,nhori, klev,1,klev,nvert,
179     .                32, "ave(X)", zsto,zout)
180c
181c        CALL histdef(nid_mth, "fluxec", "Cin. net flux","W/m2",
182c    .                iim,jjmp1,nhori, klev,1,klev,nvert,
183c    .                32, "ave(X)", zsto,zout)
184c
185      ENDIF !lev_histmth.GE.3
186c
187c-------------------------------------------------------
188      IF(lev_histmth.GE.4) THEN
189c
190         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
191     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
192     .                "ave(X)", zsto,zout)
193c
194c        CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
195c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
196c    .                "ave(X)", zsto,zout)
197c
198         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
199     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
200     .                "ave(X)", zsto,zout)
201c
202         CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
203     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
204     .                "ave(X)", zsto,zout)
205c
206         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
207     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
208     .                "ave(X)", zsto,zout)
209c
210         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
211     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
212     .                "ave(X)", zsto,zout)
213c
214c        CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
215c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
216c    .                "ave(X)", zsto,zout)
217c
218         CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2",
219     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
220     .                "ave(X)", zsto,zout)
221c
222         CALL histdef(nid_mth, "dugwo", "GW oro dU", "m/s2",
223     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
224     .                "ave(X)", zsto,zout)
225c
226         CALL histdef(nid_mth, "dugwno", "GW non-oro dU", "m/s2",
227     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
228     .                "ave(X)", zsto,zout)
229c
230c        CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
231c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
232c    .                "ave(X)", zsto,zout)
233c
234      ENDIF !lev_histmth.GE.4
235c
236c-------------------------------------------------------
237      IF(lev_histmth.GE.5) THEN
238c
239c        call histdef(nid_mth, "taux", 
240c    $         "Zonal wind stress", "Pa", 
241c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
242c    $         "ave(X)", zsto,zout)
243c
244c        call histdef(nid_mth, "tauy", 
245c    $         "Meridional xind stress", "Pa", 
246c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
247c    $         "ave(X)", zsto,zout)
248c
249c        CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
250c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
251c    .                "ave(X)", zsto,zout)
252c
253c        CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
254c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
255c    .                "ave(X)", zsto,zout)
256c
257      ENDIF !lev_histmth.GE.5
258c-------------------------------------------------------
259c
260         CALL histend(nid_mth)
261c
262         ndex2d = 0
263         ndex3d = 0
264c
265      ENDIF ! fin de test sur ok_mensuel
Note: See TracBrowser for help on using the repository browser.