source: trunk/LMDZ.VENUS/libf/phyvenus/ini_histmth.h @ 799

Last change on this file since 799 was 101, checked in by slebonnois, 14 years ago

SL: modifications pour arriver a compiler le gcm VENUS !
Ca marche !
A noter: modifs de makelmdz
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, tname(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.