source: trunk/libf/phyvenus/ini_histday.h @ 6

Last change on this file since 6 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: 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         z4julian = zjulian
12         CALL ymds2ju(annee_ref, 1, idayref, zero, z4julian)
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         z4_lon = zx_lon
21         z4_lat = zx_lat
22         zpresnivs=presnivs
23         CALL histbeg("histday", iim,z4_lon(:,1), jjmp1,z4_lat(1,:),
24     .                 1,iim,1,jjmp1, itau_phy, z4julian, zsto, 
25     .                 nhori, nid_day)
26         write(*,*)'Journee ', itau_phy, zjulian
27         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
28     .                 klev, zpresnivs, 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.