Calcul de bils_ech seulement si iflag_ener_conserv >=1
IM
- This line, and those below, will be ignored--
M ener_conserv.F90
M phys_output_write_mod.F90
12 lines of code changed in 2 files:
Just some printing improvement.
5 lines of code changed in 1 file:
Added a check. nstra can be 0 for small values of llm. If nstra
computed by sugwd_strato is 0 then there is an array index out of
bounds later in gwprofil_strato. The check does not cost much because
sugwd_strato is called only once per run.
1 lines of code changed in 1 file:
Forgot this in revision 2021.
1 lines of code changed in 1 file:
Nouvelles options des thermiques en route vers les stratocumulus
+ sorties de la discretisation verticale dans le vieux disvert_old
New options in thermals for stratocumulus
244 lines of code changed in 2 files:
Modification de la nouvelle discretisation verticale pour l'option
vert_sampling=strato_custom
Inclusion d'une version strato_correct qui correige une erreur dans strato
Modified verticale discretisation vertical for stratospheric versions.
The new version can be controled by parameter in .def files :
! Al the parameter are given in km assuming a given scalehigh
vert_scale=7. ! scale hight
vert_dzmin=0.02 ! width of first layer
vert_dzlow=1. ! dz in the low atmosphere
vert_z0low=8. ! height at which resolution recches dzlow
vert_dzmid=3. ! dz in the mid atmsophere
vert_z0mid=70. ! height at which resolution recches dzmid
vert_h_mid=20. ! width of the transition
vert_dzhig=11. ! dz in the high atmsophere
vert_z0hig=80. ! height at which resolution recches dz
vert_h_hig=20. ! width of the transition
133 lines of code changed in 1 file:
Initialisation Rayleigh scattering
S. Baek, MPL
1 lines of code changed in 1 file:
Pour compilation avec ifort
--------------
For ifort compilation
S. Baek
21 lines of code changed in 5 files:
Ajout d'un diagnostique pour la conservation de l'energie.
Added diagnostic for energie conservation.
31 lines of code changed in 4 files:
Corrections associees a la conservation de l'eau par la convection
(Jean-Yves Grandpeix)
Les difficultés sont liées aux niveaux de condensation très bas.
Finalement, j'ai opté pour une valeur minimale du delta_P menant à la
condensation (j'ai pris 1 hPa, soit 10m de soulèvement). Ça évite les
divisions par zero. Et comme les calculs sont faux, de toute façon, si le
niveau de condensation est dans la premiere couche, ce n'est pas trop
grave.
Avec ce cv3_vertmix et dans le cas twpice, j'ai à nouveau conservation
(si ok_conserv_q=y) avec la version svn 2029 (sauf une petite violation par
les thermiques, due vraisemblablement à la valeur de dqimpl (= -1 par
defaut pour les flags distribués avec les cas 1D)).
Corrections associated with numerical conservation of water in the convection.
71 lines of code changed in 1 file:
Added a "strato_custom" mode in disvert (for development purposes).
This is triggered by setting "vert_sampling=strato_custom" in a def file.
More work is needed to clean disvert (and merge it with disvert_noterre).
DC+EM
150 lines of code changed in 3 files:
1. Inclusion d'un appel supplementaire a geopot dans leapfrog dans
les versions dyn3d et dyn3dpar pour garantir la convergence
des 3 versions dyn3d/dyn3dpar/dyn3dmem.
La convergence fonctionne avec un compilateur gfortran 4.6.3 / openmpi
distribué sous ubuntu avec la nouvelle physique (NPv3.2) et le guidage.
options : gfortran -fdefault-real-8 -DNC_DOUBLE -O3
Une seule exception : si on guide avec dyn3dpar et openMP.
2. Correction du guidage dans dyn3dmem.
3. un print supprime dans cv3_toutines.F90
Rm : On ne devrait perdre la convergence numérique avec les précédentes
versions que pour la nouvelle physique puisque le géopotentiel
n'est utilisé dans la physique que par les theermiques et sisvat.
=====================================================================
1. Added call to geopot in leapfrog in dyn3d and dyn3dpar in order to
insure numerical convergence with dyn3dmem.
The convergence dyn3d/dyn3dpar/dyn3dmem is satistied with
gfortran 4.6.3 / openmpi, new physics (NPv3.2) and nudging.
options : gfortran -fdefault-real-8 -DNC_DOUBLE -O3
Only one exception : with nudging & dynd3dpar & openMP.
2. Bug fixing for nudging in dyn3dmem.
3. print supressed in cv3_routines.F90
Rm : Numerical convergence with previous releases should be lost for new
physics only since the geopotential is used only by thermals and
sisvat in the physics.
15 lines of code changed in 4 files:
Possibilité d'imposer le nombre de pas de temps de la simulation avec nday<0
If nday<0, nday is the total number of time step of the simulation
34 lines of code changed in 6 files:
PVteta computed by PVtheta was not used. Also there were a couple of
problems with PVtheta:
-- PVtheta calls tetalevel in phylmd, and interpolates at
Earth-specific values of potential temperature.
-- PVtheta calls tourabs, and the computation of rot in tourabs should
be modified (it is not correct when there is a zoom).
-- Even when there is no zoom, the computation of rot in tourabs
should probably be modified: directly combine vcov and ucov, as in
tourpot, instead of dividing by cv and cu.
9 lines of code changed in 11 files:
Correction de la version openMP du guidage
Bug fixing for openMP
83 lines of code changed in 1 file:
Correction pour compilation mixte MPI/OpenMP
Bug fixing for mixt MPI/OpenMP compilation.
3 lines of code changed in 1 file:
Nettoyage. Passage a gfortran par défaut pour ciclad.
0 lines of code changed in 1 file: