--- trunk/Sources/phylmd/coefkz2.f	2017/11/07 10:52:46	233
+++ trunk/Sources/phylmd/coefkz2.f	2017/11/07 12:20:42	234
@@ -1,99 +1,99 @@
 module coefkz2_m
-  
+
   IMPLICIT none
 
 contains
 
   SUBROUTINE coefkz2(nsrf, knon, paprs, pplay, t, pcfm, pcfh)
 
-  ! J'introduit un peu de diffusion sauf dans les endroits
-  ! ou une forte inversion est presente
-  ! On peut dire qu'il represente la convection peu profonde
-
-  use dimens_m
-  use indicesol
-  use dimphy
-  use conf_gcm_m
-  use SUPHEC_M
-
-  ! Arguments:
-  ! nsrf-----input-I- indicateur de la nature du sol
-  ! knon-----input-I- nombre de points a traiter
-  ! paprs----input-R- pression a chaque intercouche (en Pa)
-  ! pplay----input-R- pression au milieu de chaque couche (en Pa)
-  ! t--------input-R- temperature (K)
-
-  ! pcfm-----output-R- coefficients a calculer (vitesse)
-  ! pcfh-----output-R- coefficients a calculer (chaleur et humidite)
-
-  ! Arguments:
-
-  INTEGER knon, nsrf
-  REAL paprs(klon, klev+1), pplay(klon, klev)
-  REAL t(klon, klev)
-
-  REAL pcfm(klon, klev), pcfh(klon, klev)
-
-  ! Quelques constantes et options:
-
-  REAL prandtl
-  PARAMETER (prandtl=0.4)
-  REAL kstable
-  PARAMETER (kstable=0.002)
-  REAL mixlen ! constante controlant longueur de melange
-  PARAMETER (mixlen=35.0)
-  REAL seuil ! au-dela l'inversion est consideree trop faible
-  PARAMETER (seuil=-0.02)
-
-  ! Variables locales:
-
-  INTEGER i, k, invb(knon)
-  REAL zl2(knon)
-  REAL zdthmin(knon), zdthdp
-
-  !----------------------------------------------------------
-
-  ! Initialiser les sorties
-  DO k = 1, klev
-     DO i = 1, knon
-        pcfm(i, k) = 0.0
-        pcfh(i, k) = 0.0
-     ENDDO
-  ENDDO
-
-  ! Chercher la zone d'inversion forte
-
-  DO i = 1, knon
-     invb(i) = klev
-     zdthmin(i)=0.0
-  ENDDO
-  DO k = 2, klev/ 2 - 1
-     DO i = 1, knon
-        zdthdp = (t(i, k) - t(i, k + 1)) / (pplay(i, k) - pplay(i, k + 1)) &
-             - RD * 0.5 * (t(i, k) + t(i, k + 1)) / RCPD / paprs(i, k + 1)
-        zdthdp = zdthdp * 100.
-        IF (pplay(i, k) > 0.8 * paprs(i, 1) .AND. zdthdp < zdthmin(i)) THEN
-           zdthmin(i) = zdthdp
-           invb(i) = k
-        ENDIF
-     ENDDO
-  ENDDO
-
-  ! Introduire une diffusion:
-  DO k = 2, klev
-     DO i = 1, knon
-        ! si on est sur ocean et s'il n'y a pas d'inversion ou si
-        ! l'inversion est trop faible:
-        IF ((nsrf.EQ.is_oce) .AND.  &  
-             ((invb(i).EQ.klev) .OR. (zdthmin(i) > seuil))) THEN
-           zl2(i)=(mixlen*MAX(0.0, (paprs(i, k)-paprs(i, klev+1)) &
-                /(paprs(i, 2)-paprs(i, klev+1))))**2
-           pcfm(i, k)= zl2(i)* kstable
-           pcfh(i, k) = pcfm(i, k) /prandtl ! h et m different
-        ENDIF
-     ENDDO
-  ENDDO
+    ! J'introduit un peu de diffusion sauf dans les endroits
+    ! ou une forte inversion est presente
+    ! On peut dire qu'il represente la convection peu profonde
+
+    use dimens_m
+    use indicesol
+    use dimphy
+    use conf_gcm_m
+    use SUPHEC_M
+
+    ! Arguments:
+    ! nsrf-----input-I- indicateur de la nature du sol
+    ! knon-----input-I- nombre de points a traiter
+    ! paprs----input-R- pression a chaque intercouche (en Pa)
+    ! pplay----input-R- pression au milieu de chaque couche (en Pa)
+    ! t--------input-R- temperature (K)
+
+    ! pcfm-----output-R- coefficients a calculer (vitesse)
+    ! pcfh-----output-R- coefficients a calculer (chaleur et humidite)
+
+    ! Arguments:
+
+    INTEGER knon, nsrf
+    REAL paprs(klon, klev+1), pplay(klon, klev)
+    REAL t(klon, klev)
+
+    REAL, intent(out):: pcfm(:, :), pcfh(:, :) ! (knon, klev)
+
+    ! Quelques constantes et options:
+
+    REAL prandtl
+    PARAMETER (prandtl=0.4)
+    REAL kstable
+    PARAMETER (kstable=0.002)
+    REAL mixlen ! constante controlant longueur de melange
+    PARAMETER (mixlen=35.0)
+    REAL seuil ! au-dela l'inversion est consideree trop faible
+    PARAMETER (seuil=-0.02)
+
+    ! Variables locales:
+
+    INTEGER i, k, invb(knon)
+    REAL zl2(knon)
+    REAL zdthmin(knon), zdthdp
+
+    !----------------------------------------------------------
+
+    ! Initialiser les sorties
+    DO k = 1, klev
+       DO i = 1, knon
+          pcfm(i, k) = 0.0
+          pcfh(i, k) = 0.0
+       ENDDO
+    ENDDO
+
+    ! Chercher la zone d'inversion forte
+
+    DO i = 1, knon
+       invb(i) = klev
+       zdthmin(i)=0.0
+    ENDDO
+    DO k = 2, klev/ 2 - 1
+       DO i = 1, knon
+          zdthdp = (t(i, k) - t(i, k + 1)) / (pplay(i, k) - pplay(i, k + 1)) &
+               - RD * 0.5 * (t(i, k) + t(i, k + 1)) / RCPD / paprs(i, k + 1)
+          zdthdp = zdthdp * 100.
+          IF (pplay(i, k) > 0.8 * paprs(i, 1) .AND. zdthdp < zdthmin(i)) THEN
+             zdthmin(i) = zdthdp
+             invb(i) = k
+          ENDIF
+       ENDDO
+    ENDDO
+
+    ! Introduire une diffusion:
+    DO k = 2, klev
+       DO i = 1, knon
+          ! si on est sur ocean et s'il n'y a pas d'inversion ou si
+          ! l'inversion est trop faible:
+          IF ((nsrf.EQ.is_oce) .AND.  &  
+               ((invb(i).EQ.klev) .OR. (zdthmin(i) > seuil))) THEN
+             zl2(i)=(mixlen*MAX(0.0, (paprs(i, k)-paprs(i, klev+1)) &
+                  /(paprs(i, 2)-paprs(i, klev+1))))**2
+             pcfm(i, k)= zl2(i)* kstable
+             pcfh(i, k) = pcfm(i, k) /prandtl ! h et m different
+          ENDIF
+       ENDDO
+    ENDDO
 
-END SUBROUTINE coefkz2
+  END SUBROUTINE coefkz2
 
 end module coefkz2_m