Changeset 2039 for LMDZ5/trunk/libf/dyn3dmem/guide_loc_mod.F90

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.

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.

File:

• LMDZ5/trunk/libf/dyn3dmem/guide_loc_mod.F90 (modified) (1 diff)

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