Ignore:
Timestamp:
May 8, 2014, 2:11:19 AM (11 years ago)
Author:
fhourdin
Message:
  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.

  1. Correction du guidage dans dyn3dmem.
  1. 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.

  1. Bug fixing for nudging in dyn3dmem.
  1. 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.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/dyn3d/leapfrog.F

    r2038 r2039  
    384384         endif
    385385
     386! Appel a geopot ajoute le 2014/05/08 pour garantir la convergence numerique
     387! avec dyn3dmem
     388         CALL geopot  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )
     389
    386390!           rdaym_ini  = itau * dtvr / daysec
    387391!           rdayvrai   = rdaym_ini  + day_ini
Note: See TracChangeset for help on using the changeset viewer.