Changeset 2039

Timestamp:

May 8, 2014, 2:11:19 AM (10 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.

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.

Location:

LMDZ5/trunk/libf

Files:

• dyn3d/leapfrog.F (modified) (1 diff)
• dyn3dmem/guide_loc_mod.F90 (modified) (1 diff)
• dyn3dpar/leapfrog_p.F (modified) (1 diff)
• phylmd/cv3_routines.F90 (modified) (1 diff)

LMDZ5/trunk/libf/dyn3d/leapfrog.F

@@ -384 +384 @@
          endif
 
+! Appel a geopot ajoute le 2014/05/08 pour garantir la convergence numerique
+! avec dyn3dmem
+         CALL geopot  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )
+
 !           rdaym_ini  = itau * dtvr / daysec
 !           rdayvrai   = rdaym_ini  + day_ini

LMDZ5/trunk/libf/dyn3dmem/guide_loc_mod.F90

@@ -2287 +2287 @@
     CASE("u","ua")
 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-        DO l=1,llm ; print*,'l field-glo',vsize,llm,l ; field_glo(:,2:jjm,l)=field_glo(:,2:jjm,l)/cu(:,2:jjm) ; ENDDO
-        field_glo(:,1,:)=0. ; field_glo(:,jjp1,:)=0.
+        DO l=1,llm 
+            field_glo(:,2:jjm,l)=field_glo(:,2:jjm,l)/cu(:,2:jjm)
+            field_glo(:,1,l)=0. ; field_glo(:,jjp1,l)=0.
+        ENDDO
     CASE("v","va")
 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-        DO l=1,llm ; field_glo(:,:,l)=field_glo(:,:,l)/cv(:,:) ; ENDDO
+        DO l=1,llm
+           field_glo(:,:,l)=field_glo(:,:,l)/cv(:,:)
+        ENDDO
     END SELECT
 

LMDZ5/trunk/libf/dyn3dpar/leapfrog_p.F

@@ -717 +717 @@
            CALL exner_milieu_p( ip1jmp1, ps, p, pks, pk, pkf )
          endif
+! Appel a geopot ajoute le 2014/05/08 pour garantir la convergence numerique
+! avec dyn3dmem
+      CALL geopot  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )
 c$OMP BARRIER
            jD_cur = jD_ref + day_ini - day_ref

LMDZ5/trunk/libf/phylmd/cv3_routines.F90

@@ -129 +129 @@
   END IF ! (first)
 
-  print*,'tau=',tau
+! print*,'tau=',tau
   beta = 1.0 - delt/tau
   alpha1 = 1.5E-3