Index: /LMDZ6/trunk/libf/phylmd/lmdz_thermcell_old.f90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/lmdz_thermcell_old.f90	(revision 5500)
+++ /LMDZ6/trunk/libf/phylmd/lmdz_thermcell_old.f90	(revision 5501)
@@ -52,13 +52,12 @@
   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,8 +116,5 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec, ll
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
+  INTEGER ll
 
 
@@ -126,4 +122,7 @@
   ! initialisation:
   ! ---------------
+
+idetr=3
+lev_out=1
 
   sorties = .TRUE.
@@ -174,12 +173,9 @@
   ! -----------------------------------------------------------------------
 
-  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,7 +760,4 @@
 
   INTEGER idetr
-  SAVE idetr
-  DATA idetr/3/
-  !$OMP THREADPRIVATE(idetr)
 
   ! local:
@@ -778,7 +771,4 @@
   REAL zmix(klon), fracazmix(klon)
   REAL alpha
-  SAVE alpha
-  DATA alpha/1./
-  !$OMP THREADPRIVATE(alpha)
 
   ! RC
@@ -890,5 +880,5 @@
   REAL f_old
   REAL zlevinter(klon)
-  LOGICAL, SAVE :: first = .TRUE.
+  LOGICAL,SAVE :: first = .TRUE.
   !$OMP THREADPRIVATE(first)
   ! data first /.false./
@@ -915,10 +905,9 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec, ll
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
-
-
+  INTEGER ll
+
+
+  idetr=3
+  alpha=1.
 
   ! -----------------------------------------------------------------------
@@ -2365,7 +2354,4 @@
 
   INTEGER idetr
-  SAVE idetr
-  DATA idetr/3/
-  !$OMP THREADPRIVATE(idetr)
 
   ! local:
@@ -2459,8 +2445,5 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec, ll
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
+  INTEGER ll
 
 
@@ -2470,4 +2453,5 @@
   ! ---------------
 
+  idetr=3
   sorties = .TRUE.
   IF (ngrid/=klon) THEN
@@ -3294,7 +3278,4 @@
 
   INTEGER idetr
-  SAVE idetr
-  DATA idetr/3/
-  !$OMP THREADPRIVATE(idetr)
 
   ! local:
@@ -3375,8 +3356,5 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec, ll
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
+  INTEGER ll
 
 
@@ -3385,4 +3363,5 @@
   ! ---------------
 
+  idetr=3
   sorties = .TRUE.
   IF (ngrid/=klon) THEN
@@ -4507,7 +4486,4 @@
 
   INTEGER idetr
-  SAVE idetr
-  DATA idetr/3/
-  !$OMP THREADPRIVATE(idetr)
 
   ! local:
@@ -4572,9 +4548,4 @@
   REAL f(klon), f0(klon)
   REAL zlevinter(klon)
-  LOGICAL first
-  DATA first/.FALSE./
-  SAVE first
-  !$OMP THREADPRIVATE(first)
-  ! RC
 
   CHARACTER *2 str2
@@ -4588,8 +4559,5 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec, ll
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
+  INTEGER ll
 
 
@@ -4598,4 +4566,5 @@
   ! ---------------
 
+  idetr=3
   sorties = .TRUE.
   IF (ngrid/=klon) THEN
@@ -4612,4 +4581,5 @@
   ! print*,'0 OK convect8'
 
+  idetr=3
   DO l = 1, nlay
     DO ig = 1, ngrid
@@ -5337,7 +5307,4 @@
 
   INTEGER idetr
-  SAVE idetr
-  DATA idetr/3/
-  !$OMP THREADPRIVATE(idetr)
   ! local:
   ! ------
@@ -5376,10 +5343,4 @@
 
   REAL count_time
-  ! integer isplit,nsplit
-  INTEGER isplit, nsplit, ialt
-  PARAMETER (nsplit=10)
-  DATA isplit/0/
-  SAVE isplit
-  !$OMP THREADPRIVATE(isplit)
 
   LOGICAL sorties
@@ -5427,8 +5388,4 @@
   EXTERNAL scopy
 
-  INTEGER ncorrec
-  SAVE ncorrec
-  DATA ncorrec/0/
-  !$OMP THREADPRIVATE(ncorrec)
 
 
@@ -5437,4 +5394,5 @@
   ! ---------------
 
+  idetr=3
   sorties = .TRUE.
   IF (ngrid/=klon) THEN
@@ -6226,5 +6184,4 @@
     ! print*,'15 OK convect8'
 
-    isplit = isplit + 1
 
   END IF