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

Last change on this file since 86 was 3, checked in by slebonnois, 15 years ago

Creation de repertoires:

  • chantiers : pour communiquer sur nos projets de modifs
  • documentation : pour stocker les docs

Ajout de:

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