Changeset 5144 for LMDZ6/branches/Amaury_dev/libf/phylmd/stratocu_if.F90

Timestamp:

Jul 29, 2024, 11:01:04 PM (3 months ago)

Author:

abarral

Message:

Put YOMCST.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/stratocu_if.F90 (modified) (1 diff)

LMDZ6/branches/Amaury_dev/libf/phylmd/stratocu_if.F90

@@ -1 @@
-  SUBROUTINE stratocu_if(klon,klev,pctsrf,paprs, pplay,t &
-,seuil_inversion,weak_inversion,dthmin)
+SUBROUTINE stratocu_if(klon, klev, pctsrf, paprs, pplay, t &
+        , seuil_inversion, weak_inversion, dthmin)
 
   USE indice_sol_mod
+  USE lmdz_yomcst
 
-IMPLICIT NONE
+  IMPLICIT NONE
 
-!======================================================================
-! 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
+  !======================================================================
+  ! 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
 
-! Arguments:
-! klon-----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)
+  ! Arguments:
+  ! klon-----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)
 
-! weak_inversion-----logical
-!======================================================================
+  ! weak_inversion-----logical
+  !======================================================================
 
-! Arguments:
+  ! Arguments:
 
-    INTEGER, INTENT(IN)                       :: klon,klev
-    REAL, DIMENSION(klon, klev+1), INTENT(IN) ::  paprs
-    REAL, DIMENSION(klon, klev), INTENT(IN)   ::  pplay
-    REAL, DIMENSION(klon, 4), INTENT(IN)   ::  pctsrf
-    REAL, DIMENSION(klon, klev), INTENT(IN)   :: t
-    
-    REAL, DIMENSION(klon), INTENT(OUT)  :: weak_inversion
+  INTEGER, INTENT(IN) :: klon, klev
+  REAL, DIMENSION(klon, klev + 1), INTENT(IN) :: paprs
+  REAL, DIMENSION(klon, klev), INTENT(IN) :: pplay
+  REAL, DIMENSION(klon, 4), INTENT(IN) :: pctsrf
+  REAL, DIMENSION(klon, klev), INTENT(IN) :: t
 
-! Quelques constantes et options:
+  REAL, DIMENSION(klon), INTENT(OUT) :: weak_inversion
 
-    REAL seuil_inversion ! au-dela l'inversion est consideree trop faible
-!    PARAMETER (seuil=-0.1)
+  ! Quelques constantes et options:
 
-! Variables locales:
+  REAL seuil_inversion ! au-dela l'inversion est consideree trop faible
+  !    PARAMETER (seuil=-0.1)
 
-    INTEGER i, k, invb(klon)
-    REAL zl2(klon)
-    REAL dthmin(klon), zdthdp
+  ! Variables locales:
 
-    INCLUDE "YOMCST.h"
+  INTEGER i, k, invb(klon)
+  REAL zl2(klon)
+  REAL dthmin(klon), zdthdp
 
-! Chercher la zone d'inversion forte
+  ! Chercher la zone d'inversion forte
 
+  DO i = 1, klon
+    invb(i) = klev
+    dthmin(i) = 0.0
+  ENDDO
+  DO k = 2, klev / 2 - 1
     DO i = 1, klon
-       invb(i) = klev
-       dthmin(i)=0.0
+      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.0
+      IF (pplay(i, k)>0.8 * paprs(i, 1) .AND. &
+              zdthdp<dthmin(i)) THEN
+        dthmin(i) = zdthdp
+        invb(i) = k
+      ENDIF
     ENDDO
-    DO k = 2, klev/2-1
-       DO i = 1, klon
-          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.0
-          IF (pplay(i,k)>0.8*paprs(i,1) .AND. &
-               zdthdp<dthmin(i) ) THEN
-             dthmin(i) = zdthdp
-             invb(i) = k
-          ENDIF
-       ENDDO
-    ENDDO
+  ENDDO
 
-! Introduire une diffusion:
+  ! Introduire une diffusion:
 
-    DO i = 1, klon
-       IF ( (pctsrf(i,is_oce) < 0.5) .OR. &
-          (invb(i) == klev) .OR. (dthmin(i) > seuil_inversion) ) THEN 
-          weak_inversion(i)=1.
-       ELSE
-          weak_inversion(i)=0.
-       ENDIF
-    ENDDO
+  DO i = 1, klon
+    IF ((pctsrf(i, is_oce) < 0.5) .OR. &
+            (invb(i) == klev) .OR. (dthmin(i) > seuil_inversion)) THEN
+      weak_inversion(i) = 1.
+    ELSE
+      weak_inversion(i) = 0.
+    ENDIF
+  ENDDO
 
-  END SUBROUTINE stratocu_if
+END SUBROUTINE stratocu_if