source: LMDZ4/trunk/libf/phylmd/ini_bilKP_ave.h @ 5442

Last change on this file since 5442 was 1403, checked in by Laurent Fairhead, 14 years ago

Merged LMDZ4V5.0-dev branch changes r1292:r1399 to trunk.

Validation:
Validation consisted in compiling the HEAD revision of the trunk,
LMDZ4V5.0-dev branch and the merged sources and running different
configurations on local and SX8 machines comparing results.

Local machine: bench configuration, 32x24x11, gfortran

  • IPSLCM5A configuration (comparison between trunk and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent
  • MH07 configuration, new physics package (comparison between LMDZ4V5.0-dev branch and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent

SX8 machine (brodie), 96x95x39 on 4 processors:

  • IPSLCM5A configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent
  • MH07 configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent

Changes to the makegcm and create_make_gcm scripts to take into account
main programs in F90 files

Fusion de la branche LMDZ4V5.0-dev (r1292:r1399) au tronc principal

Validation:
La validation a consisté à compiler la HEAD de le trunk et de la banche
LMDZ4V5.0-dev et les sources fusionnées et de faire tourner le modéle selon
différentes configurations en local et sur SX8 et de comparer les résultats

En local: 32x24x11, config bench/gfortran

  • pour une config IPSLCM5A (comparaison tronc/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux (à part sur RN et Pb)
    • fichiers histoire égaux
  • pour une config nlle physique (MH07) (comparaison LMDZ4v5.0-dev/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux
    • fichiers histoire équivalents

Sur brodie, 96x95x39 sur 4 proc:

  • pour une config IPSLCM5A:
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc
  • pour une config MH07
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc

Changement sur makegcm et create_make-gcm pour pouvoir prendre en compte des
programmes principaux en *F90
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 8.8 KB
Line 
1c
2c $Id: ini_bilKP_ave.h 1403 2010-07-01 09:02:53Z fhourdin $
3c
4      IF (ok_journe) THEN
5c
6         zsto = dtime
7         zout = ecrit_day
8         typeval=tave
9c
10         idayref = day_ref
11         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
12c
13cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
14cym         DO i = 1, iim
15cym            zx_lon(i,1) = rlon(i+1)
16cym            zx_lon(i,jjmp1) = rlon(i+1)
17cym         ENDDO
18         DO ll=1,klev
19            znivsig(ll)=REAL(ll)
20         ENDDO
21cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
22cym         write(*,*)'zx_lon = ',zx_lon(:,1)
23cym         write(*,*)'zx_lat = ',zx_lat(1,:)
24cym         CALL histbeg("histbilKP_ave", iim,zx_lon(:,1), jjmp1,
25cym     .                zx_lat(1,:),
26cym     .                1,iim,1,jjmp1, itau_phy, zjulian, dtime,
27cym     .                nhori, nid_bilKPave)
28         CALL histbeg_phy("histbilKP_ave", itau_phy, zjulian, dtime,
29     .                nhori, nid_bilKPave)
30
31         write(*,*)'Journee ', itau_phy, zjulian
32         CALL histvert(nid_bilKPave, "presnivs",
33     .                "Vertical levels","mb",
34     .                 klev, presnivs/100., nvert)
35c
36c
37c Champs 3D:
38c
39         CALL histdef(nid_bilKPave,"ue",
40     .   "Zonal energy transport","-",
41     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 
42     .                typeval, zsto,zout)
43c
44         CALL histdef(nid_bilKPave,"ve",
45     .   "Merid energy transport","-",
46     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 
47     .                typeval, zsto,zout)
48c
49         CALL histdef(nid_bilKPave,"uq",
50     .   "Zonal humidity transport","-",
51     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 
52     .                typeval, zsto,zout)
53c
54         CALL histdef(nid_bilKPave,"vq",
55     .   "Merid humidity transport","-",
56     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 
57     .                typeval, zsto,zout)
58c
59c Champs 3D:
60c
61         CALL histdef(nid_bilKPave,"temp",
62     .   "Air temperature","K",
63     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
64     .                typeval, zsto,zout)
65c
66         CALL histdef(nid_bilKPave,"ovap",
67     .   "Specific humidity","Kg/Kg",
68     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
69     .                typeval, zsto,zout)
70c
71         CALL histdef(nid_bilKPave,"geop",
72     .   "Geopotential height","m",
73     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
74     .                typeval, zsto,zout)
75c
76         CALL histdef(nid_bilKPave,"vitu",
77     .   "Zonal wind","m/s",
78     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
79     .                typeval, zsto,zout)
80c
81         CALL histdef(nid_bilKPave,"vitv",
82     .   "Meridional wind","m/s",
83     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
84     .                typeval, zsto,zout)
85c
86         CALL histdef(nid_bilKPave, "vitw", 
87     .   "Vertical wind", "m/s",
88     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
89     .                typeval, zsto,zout)
90c
91         CALL histdef(nid_bilKPave, "pres", 
92     .   "Inter-Layer Air pressure",
93     .                "Pa",
94     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
95     .                typeval, zsto,zout)
96c
97         CALL histdef(nid_bilKPave, "play", 
98     .   "Mean-Layer Air pressure",
99     .                "Pa",
100     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
101     .                typeval, zsto,zout)
102c
103         CALL histdef(nid_bilKPave, "oliq", 
104     .   "Liquid water content", 
105     .                "kg/kg",
106     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
107     .                typeval, zsto,zout)
108c
109         CALL histdef(nid_bilKPave, "dtdyn", 
110     .   "Dynamics dT", "K/s",
111     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
112     .                typeval, zsto,zout)
113c
114         CALL histdef(nid_bilKPave, "dqdyn", 
115     .   "Dynamics dQ", "Kg/Kg/s",
116     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
117     .                typeval, zsto,zout)
118c
119         CALL histdef(nid_bilKPave, "dtcon", 
120     .   "Convection dT", "K/s",
121     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
122     .                typeval, zsto,zout)
123c
124         CALL histdef(nid_bilKPave, "ducon", 
125     .   "Convection du", "m/s2",
126     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
127     .                typeval, zsto,zout)
128c
129         CALL histdef(nid_bilKPave, "dvcon", 
130     .   "Convection dv", "m/s2",
131     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
132     .                typeval, zsto,zout)
133c
134         CALL histdef(nid_bilKPave,"dqcon",
135     .   "Convection dQ","Kg/Kg/s",
136     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
137     .                typeval, zsto,zout)
138c
139         CALL histdef(nid_bilKPave, "dtlsc", 
140     .   "Condensation dT", "K/s",
141     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
142     .                typeval, zsto,zout)
143c
144         CALL histdef(nid_bilKPave,"dqlsc",
145     .   "Condensation dQ","Kg/Kg/s",
146     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
147     .                typeval, zsto,zout)
148c
149         CALL histdef(nid_bilKPave,"dtvdf",
150     .   "Boundary-layer dT","K/s",
151     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
152     .                typeval, zsto,zout)
153c
154         CALL histdef(nid_bilKPave, "dqvdf", 
155     .   "Boundary-layer dQ", 
156     .               "Kg/Kg/s",
157     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
158     .                typeval, zsto,zout)
159c
160         CALL histdef(nid_bilKPave,"dtajs",
161     .   "Ajustement sec dT","K/s",
162     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
163     .                typeval, zsto,zout)
164c
165         CALL histdef(nid_bilKPave, "dqajs",
166     .   "Ajustement sec dQ", 
167     .               "Kg/Kg/s",
168     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
169     .                typeval, zsto,zout)
170c
171         CALL histdef(nid_bilKPave, "dteva", 
172     .   "Reevaporation dT", "K/s",
173     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
174     .                typeval, zsto,zout)
175c
176         CALL histdef(nid_bilKPave,"dqeva",
177     .   "Reevaporation dQ",
178     .                "Kg/Kg/s",
179     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
180     .                typeval, zsto,zout)
181
182c
183         CALL histdef(nid_bilKPave, "dtswr", 
184     .   "SW radiation dT", "K/s",
185     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
186     .                typeval, zsto,zout)
187c
188         CALL histdef(nid_bilKPave, "dtsw0", 
189     .   "SW radiation dT", "K/s",
190     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
191     .                typeval, zsto,zout)
192c
193         CALL histdef(nid_bilKPave, "dtlwr", 
194     .   "LW radiation dT", "K/s",
195     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
196     .                typeval, zsto,zout)
197c
198         CALL histdef(nid_bilKPave, "dtlw0", 
199     .   "LW radiation dT", "K/s",
200     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
201     .                typeval, zsto,zout)
202c
203         CALL histdef(nid_bilKPave,"duvdf",
204     .   "Boundary-layer dU","m/s2",
205     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
206     .                typeval, zsto,zout)
207c
208         CALL histdef(nid_bilKPave,"dvvdf",
209     .   "Boundary-layer dV","m/s2",
210     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
211     .                typeval, zsto,zout)
212c
213         IF (ok_orodr) THEN
214         IF (ok_orolf) THEN
215         CALL histdef(nid_bilKPave, "duoli",
216     .   "Orography dU", "m/s2",
217     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
218     .                typeval, zsto,zout)
219c
220         CALL histdef(nid_bilKPave, "dvoli", 
221     .   "Orography dV", "m/s2",
222     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
223     .                typeval, zsto,zout)
224c
225         ENDIF
226         ENDIF
227C
228         CALL histdef(nid_bilKPave, "duphy",
229     .   "Physiq dU","-",
230     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
231     .                typeval, zsto,zout)
232C
233         CALL histdef(nid_bilKPave, "dvphy",
234     .   "Physiq dV","-",
235     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
236     .                typeval, zsto,zout)
237C
238         CALL histdef(nid_bilKPave, "dtphy",
239     .   "Physiq dT","-",
240     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
241     .                typeval, zsto,zout)
242C
243         CALL histdef(nid_bilKPave, "dqphy",
244     .   "Physiq dQ","-",
245     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
246     .                typeval, zsto,zout)
247C
248         CALL histdef(nid_bilKPave, "dqlphy",
249     .   "Physiq dQl","-",
250     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
251     .                typeval, zsto,zout)
252C
253C
254         CALL histend(nid_bilKPave)
255c
256         ndex2d = 0
257         ndex3d = 0
258c
259      ENDIF ! fin de test sur ok_journe
Note: See TracBrowser for help on using the repository browser.