Changeset 5501 for LMDZ6/trunk/libf/phylmd

Timestamp:

Jan 21, 2025, 6:46:33 PM (2 days ago)

Author:

fhourdin

Message:

Supression de save

File:

• LMDZ6/trunk/libf/phylmd/lmdz_thermcell_old.f90 (modified) (24 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_thermcell_old.f90

@@ -52 +52 @@
   REAL fraca(ngrid, nlay+1), zw2(ngrid, nlay+1)
 
-  INTEGER, SAVE :: idetr = 3, lev_out = 1
-  !$OMP THREADPRIVATE(idetr,lev_out)
+  INTEGER :: idetr, lev_out
 
   ! local:
   ! ------
 
-  INTEGER, SAVE :: dvdq = 0, flagdq = 0, dqimpl = 1
-  LOGICAL, SAVE :: debut = .TRUE.
-  !$OMP THREADPRIVATE(dvdq,flagdq,debut,dqimpl)
+  INTEGER :: dvdq, flagdq, dqimpl
+  LOGICAL :: debut
+
 
   INTEGER ig, k, l, lmax(klon, klev+1), lmaxa(klon), lmix(klon)
@@ -117 +116 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec, ll
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
+  INTEGER ll
 
 
@@ -126 +122 @@
   ! initialisation:
   ! ---------------
+
+idetr=3
+lev_out=1
 
   sorties = .TRUE.
@@ -174 +173 @@
   ! -----------------------------------------------------------------------
 
-  IF (debut) THEN
-    flagdq = (iflag_thermals-1000)/100
-    dvdq = (iflag_thermals-(1000+flagdq*100))/10
-    IF (flagdq==2) dqimpl = -1
-    IF (flagdq==3) dqimpl = 1
-    debut = .FALSE.
-  END IF
-  PRINT *, 'TH flag th ', iflag_thermals, flagdq, dvdq, dqimpl
+  flagdq = (iflag_thermals-1000)/100
+  dvdq = (iflag_thermals-(1000+flagdq*100))/10
+  IF (flagdq==2) dqimpl = -1
+  IF (flagdq==3) dqimpl = 1
+  !PRINT *, 'TH flag th ', iflag_thermals, flagdq, dvdq, dqimpl
 
   DO l = 2, nlay
@@ -764 +760 @@
 
   INTEGER idetr
-  SAVE idetr
-  DATA idetr/3/
-  !$OMP THREADPRIVATE(idetr)
 
   ! local:
@@ -778 +771 @@
   REAL zmix(klon), fracazmix(klon)
   REAL alpha
-  SAVE alpha
-  DATA alpha/1./
-  !$OMP THREADPRIVATE(alpha)
 
   ! RC
@@ -890 +880 @@
   REAL f_old
   REAL zlevinter(klon)
-  LOGICAL, SAVE :: first = .TRUE.
+  LOGICAL,SAVE :: first = .TRUE.
   !$OMP THREADPRIVATE(first)
   ! data first /.false./
@@ -915 +905 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec, ll
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
-
-
+  INTEGER ll
+
+
+  idetr=3
+  alpha=1.
 
   ! -----------------------------------------------------------------------
@@ -2365 +2354 @@
 
   INTEGER idetr
-  SAVE idetr
-  DATA idetr/3/
-  !$OMP THREADPRIVATE(idetr)
 
   ! local:
@@ -2459 +2445 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec, ll
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
+  INTEGER ll
 
 
@@ -2470 +2453 @@
   ! ---------------
 
+  idetr=3
   sorties = .TRUE.
   IF (ngrid/=klon) THEN
@@ -3294 +3278 @@
 
   INTEGER idetr
-  SAVE idetr
-  DATA idetr/3/
-  !$OMP THREADPRIVATE(idetr)
 
   ! local:
@@ -3375 +3356 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec, ll
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
+  INTEGER ll
 
 
@@ -3385 +3363 @@
   ! ---------------
 
+  idetr=3
   sorties = .TRUE.
   IF (ngrid/=klon) THEN
@@ -4507 +4486 @@
 
   INTEGER idetr
-  SAVE idetr
-  DATA idetr/3/
-  !$OMP THREADPRIVATE(idetr)
 
   ! local:
@@ -4572 +4548 @@
   REAL f(klon), f0(klon)
   REAL zlevinter(klon)
-  LOGICAL first
-  DATA first/.FALSE./
-  SAVE first
-  !$OMP THREADPRIVATE(first)
-  ! RC
 
   CHARACTER *2 str2
@@ -4588 +4559 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec, ll
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
+  INTEGER ll
 
 
@@ -4598 +4566 @@
   ! ---------------
 
+  idetr=3
   sorties = .TRUE.
   IF (ngrid/=klon) THEN
@@ -4612 +4581 @@
   ! print*,'0 OK convect8'
 
+  idetr=3
   DO l = 1, nlay
     DO ig = 1, ngrid
@@ -5337 +5307 @@
 
   INTEGER idetr
-  SAVE idetr
-  DATA idetr/3/
-  !$OMP THREADPRIVATE(idetr)
   ! local:
   ! ------
@@ -5376 +5343 @@
 
   REAL count_time
-  ! integer isplit,nsplit
-  INTEGER isplit, nsplit, ialt
-  PARAMETER (nsplit=10)
-  DATA isplit/0/
-  SAVE isplit
-  !$OMP THREADPRIVATE(isplit)
 
   LOGICAL sorties
@@ -5427 +5388 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
 
 
@@ -5437 +5394 @@
   ! ---------------
 
+  idetr=3
   sorties = .TRUE.
   IF (ngrid/=klon) THEN
@@ -6226 +6184 @@
     ! print*,'15 OK convect8'
 
-    isplit = isplit + 1
 
   END IF