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

Last change on this file since 5305 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
RevLine 
[644]1c
[1403]2c $Id: ini_bilKP_ave.h 1403 2010-07-01 09:02:53Z abarral $
[644]3c
4      IF (ok_journe) THEN
5c
6         zsto = dtime
[684]7         zout = ecrit_day
[644]8         typeval=tave
9c
10         idayref = day_ref
11         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
12c
[766]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
[644]18         DO ll=1,klev
[1403]19            znivsig(ll)=REAL(ll)
[644]20         ENDDO
[766]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,
[644]29     .                nhori, nid_bilKPave)
[766]30
[644]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","-",
[766]41     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 
[644]42     .                typeval, zsto,zout)
43c
44         CALL histdef(nid_bilKPave,"ve",
45     .   "Merid energy transport","-",
[766]46     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 
[644]47     .                typeval, zsto,zout)
48c
49         CALL histdef(nid_bilKPave,"uq",
50     .   "Zonal humidity transport","-",
[766]51     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 
[644]52     .                typeval, zsto,zout)
53c
54         CALL histdef(nid_bilKPave,"vq",
55     .   "Merid humidity transport","-",
[766]56     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 
[644]57     .                typeval, zsto,zout)
58c
59c Champs 3D:
60c
61         CALL histdef(nid_bilKPave,"temp",
62     .   "Air temperature","K",
[766]63     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]64     .                typeval, zsto,zout)
65c
66         CALL histdef(nid_bilKPave,"ovap",
67     .   "Specific humidity","Kg/Kg",
[766]68     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]69     .                typeval, zsto,zout)
70c
71         CALL histdef(nid_bilKPave,"geop",
72     .   "Geopotential height","m",
[766]73     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]74     .                typeval, zsto,zout)
75c
76         CALL histdef(nid_bilKPave,"vitu",
77     .   "Zonal wind","m/s",
[766]78     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]79     .                typeval, zsto,zout)
80c
81         CALL histdef(nid_bilKPave,"vitv",
82     .   "Meridional wind","m/s",
[766]83     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]84     .                typeval, zsto,zout)
85c
86         CALL histdef(nid_bilKPave, "vitw", 
87     .   "Vertical wind", "m/s",
[766]88     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]89     .                typeval, zsto,zout)
90c
91         CALL histdef(nid_bilKPave, "pres", 
92     .   "Inter-Layer Air pressure",
93     .                "Pa",
[766]94     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]95     .                typeval, zsto,zout)
96c
97         CALL histdef(nid_bilKPave, "play", 
98     .   "Mean-Layer Air pressure",
99     .                "Pa",
[766]100     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]101     .                typeval, zsto,zout)
102c
103         CALL histdef(nid_bilKPave, "oliq", 
104     .   "Liquid water content", 
105     .                "kg/kg",
[766]106     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]107     .                typeval, zsto,zout)
108c
109         CALL histdef(nid_bilKPave, "dtdyn", 
110     .   "Dynamics dT", "K/s",
[766]111     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]112     .                typeval, zsto,zout)
113c
114         CALL histdef(nid_bilKPave, "dqdyn", 
115     .   "Dynamics dQ", "Kg/Kg/s",
[766]116     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]117     .                typeval, zsto,zout)
118c
119         CALL histdef(nid_bilKPave, "dtcon", 
120     .   "Convection dT", "K/s",
[766]121     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]122     .                typeval, zsto,zout)
123c
124         CALL histdef(nid_bilKPave, "ducon", 
125     .   "Convection du", "m/s2",
[766]126     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]127     .                typeval, zsto,zout)
128c
129         CALL histdef(nid_bilKPave, "dvcon", 
130     .   "Convection dv", "m/s2",
[766]131     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]132     .                typeval, zsto,zout)
133c
134         CALL histdef(nid_bilKPave,"dqcon",
135     .   "Convection dQ","Kg/Kg/s",
[766]136     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]137     .                typeval, zsto,zout)
138c
139         CALL histdef(nid_bilKPave, "dtlsc", 
140     .   "Condensation dT", "K/s",
[766]141     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]142     .                typeval, zsto,zout)
143c
144         CALL histdef(nid_bilKPave,"dqlsc",
145     .   "Condensation dQ","Kg/Kg/s",
[766]146     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]147     .                typeval, zsto,zout)
148c
149         CALL histdef(nid_bilKPave,"dtvdf",
150     .   "Boundary-layer dT","K/s",
[766]151     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]152     .                typeval, zsto,zout)
153c
154         CALL histdef(nid_bilKPave, "dqvdf", 
155     .   "Boundary-layer dQ", 
156     .               "Kg/Kg/s",
[766]157     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]158     .                typeval, zsto,zout)
159c
160         CALL histdef(nid_bilKPave,"dtajs",
161     .   "Ajustement sec dT","K/s",
[766]162     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]163     .                typeval, zsto,zout)
164c
165         CALL histdef(nid_bilKPave, "dqajs",
166     .   "Ajustement sec dQ", 
167     .               "Kg/Kg/s",
[766]168     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]169     .                typeval, zsto,zout)
170c
171         CALL histdef(nid_bilKPave, "dteva", 
172     .   "Reevaporation dT", "K/s",
[766]173     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]174     .                typeval, zsto,zout)
175c
176         CALL histdef(nid_bilKPave,"dqeva",
177     .   "Reevaporation dQ",
178     .                "Kg/Kg/s",
[766]179     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]180     .                typeval, zsto,zout)
181
182c
183         CALL histdef(nid_bilKPave, "dtswr", 
184     .   "SW radiation dT", "K/s",
[766]185     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]186     .                typeval, zsto,zout)
187c
188         CALL histdef(nid_bilKPave, "dtsw0", 
189     .   "SW radiation dT", "K/s",
[766]190     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]191     .                typeval, zsto,zout)
192c
193         CALL histdef(nid_bilKPave, "dtlwr", 
194     .   "LW radiation dT", "K/s",
[766]195     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]196     .                typeval, zsto,zout)
197c
198         CALL histdef(nid_bilKPave, "dtlw0", 
199     .   "LW radiation dT", "K/s",
[766]200     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]201     .                typeval, zsto,zout)
202c
203         CALL histdef(nid_bilKPave,"duvdf",
204     .   "Boundary-layer dU","m/s2",
[766]205     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]206     .                typeval, zsto,zout)
207c
208         CALL histdef(nid_bilKPave,"dvvdf",
209     .   "Boundary-layer dV","m/s2",
[766]210     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]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",
[766]217     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]218     .                typeval, zsto,zout)
219c
220         CALL histdef(nid_bilKPave, "dvoli", 
221     .   "Orography dV", "m/s2",
[766]222     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]223     .                typeval, zsto,zout)
224c
225         ENDIF
226         ENDIF
227C
228         CALL histdef(nid_bilKPave, "duphy",
229     .   "Physiq dU","-",
[766]230     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]231     .                typeval, zsto,zout)
232C
233         CALL histdef(nid_bilKPave, "dvphy",
234     .   "Physiq dV","-",
[766]235     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]236     .                typeval, zsto,zout)
237C
238         CALL histdef(nid_bilKPave, "dtphy",
239     .   "Physiq dT","-",
[766]240     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]241     .                typeval, zsto,zout)
242C
243         CALL histdef(nid_bilKPave, "dqphy",
244     .   "Physiq dQ","-",
[766]245     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]246     .                typeval, zsto,zout)
247C
248         CALL histdef(nid_bilKPave, "dqlphy",
249     .   "Physiq dQl","-",
[766]250     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
[644]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.