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

Last change on this file since 1661 was 1530, checked in by emillour, 9 years ago

Venus and Titan GCMs:
Updates in the physics to keep up with updates in LMDZ5 (up to
LMDZ5 trunk, rev 2350) concerning dynamics/physics separation:

  • Adapted makelmdz and makelmdz_fcm script to stop if trying to compile 1d model or newstart or start2archive in parallel.
  • got rid of references to "dimensions.h" in physics. Within physics packages, use nbp_lon (=iim), nbp_lat (=jjmp1) and nbp_lev (=llm) from module mod_grid_phy_lmdz (in phy_common) instead. Only partially done for Titan, because of many hard-coded commons; a necessary first step will be to clean these up (using modules).

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