Changeset 1959

Timestamp:

Jul 3, 2018, 2:43:42 PM (6 years ago)

Author:

emillour

Message:

Common (parallel) dynamics:
Add some initializations of global arrays. This is not a bug fix, strictly
speaking, as results are unchanged. But in some cases computations extend
into halos and these cases are identified as errors with recent versions
of ifort (2018) with the -init=arrays,snan debugging option.
EM

Location:

trunk/LMDZ.COMMON/libf/dyn3dpar

Files:

• advtrac_p.F90 (modified) (1 diff)
• integrd_p.F (modified) (1 diff)
• leapfrog_p.F (modified) (1 diff)

trunk/LMDZ.COMMON/libf/dyn3dpar/advtrac_p.F90

@@ -47 +47 @@
   !-------------------------------------------------------------
 
-  REAL pbaruc(ip1jmp1,llm),pbarvc(ip1jm,llm)
-  REAL massem(ip1jmp1,llm),zdp(ip1jmp1)
+  REAL,SAVE :: pbaruc(ip1jmp1,llm)=0
+  REAL,SAVE :: pbarvc(ip1jm,llm)=0
+  REAL,SAVE :: massem(ip1jmp1,llm)
+  REAL zdp(ip1jmp1)
   REAL,SAVE::pbarug(ip1jmp1,llm),pbarvg(ip1jm,llm),wg(ip1jmp1,llm) 
   REAL (kind=kind(1.d0)) :: t_initial, t_final, tps_cpu
-  INTEGER iadvtr
   INTEGER ij,l,iq,iiq
   REAL zdpmin, zdpmax
-  SAVE iadvtr, massem, pbaruc, pbarvc
-  DATA iadvtr/0/
+  INTEGER,SAVE :: iadvtr=0
   !$OMP THREADPRIVATE(iadvtr)
   !----------------------------------------------------------

trunk/LMDZ.COMMON/libf/dyn3dpar/integrd_p.F

@@ -76 +76 @@
       EXTERNAL SSUM
       INTEGER ijb,ije,jjb,jje
-      REAL,SAVE :: ps(ip1jmp1)
+      REAL,SAVE :: ps(ip1jmp1)=0
       LOGICAL :: checksum
       INTEGER :: stop_it

trunk/LMDZ.COMMON/libf/dyn3dpar/leapfrog_p.F

@@ -292 +292 @@
         dpfi(:)     =0.
         dq(:,:,:)   =0.
+        dp(:)=0
+        pbaru(:,:)=0
+        pbarv(:,:)=0
 
       CALL pression ( ip1jmp1, ap, bp, ps, p       )