/[lmdze]/trunk/phylmd/coefkz2.f

Diff of /trunk/phylmd/coefkz2.f

Parent Directory | Revision Log | View Patch Patch

-trunk/phylmd/coefkz2.f90
revision 76 by guez,
Fri Nov 15 18:45:49 2013 UTC
+trunk/Sources/phylmd/coefkz2.f
revision 251 by guez,
Mon Jan  8 14:12:02 2018 UTC
 Line 1
- SUBROUTINE coefkz2(nsrf, knon, paprs, pplay, t, pcfm, pcfh)
+ module coefkz2_m
-   ! 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
    IMPLICIT none
-   ! Arguments:
+ contains
-   ! nsrf-----input-I- indicateur de la nature du sol
-   ! knon-----input-I- nombre de points a traiter
+   SUBROUTINE coefkz2(nsrf, paprs, pplay, t, pcfm, pcfh)
-   ! paprs----input-R- pression a chaque intercouche (en Pa)
-   ! pplay----input-R- pression au milieu de chaque couche (en Pa)
+     ! J'introduit un peu de diffusion sauf dans les endroits o\`u une
-   ! t--------input-R- temperature (K)
+     ! forte inversion est pr\'esente. On peut dire que la diffusion
+     ! repr\'esente la convection peu profonde.
-   ! pcfm-----output-R- coefficients a calculer (vitesse)
-   ! pcfh-----output-R- coefficients a calculer (chaleur et humidite)
+     use indicesol, only: is_oce
+     use dimphy, only: klev
-   ! Arguments:
+     use SUPHEC_M, only: RCPD, rd
-   INTEGER knon, nsrf
+     integer, intent(in):: nsrf ! indicateur de la nature du sol
-   REAL paprs(klon, klev+1), pplay(klon, klev)
-   REAL t(klon, klev)
+     REAL, intent(in):: paprs(:, :) ! (knon, klev+1)
+     ! pression a chaque intercouche (en Pa)
-   REAL pcfm(klon, klev), pcfh(klon, klev)
+     REAL, intent(in):: pplay(:, :) ! (knon, klev)
-   ! Quelques constantes et options:
+     ! pression au milieu de chaque couche (en Pa)
-   REAL prandtl
+     REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K)
-   PARAMETER (prandtl=0.4)
-   REAL kstable
+     REAL, intent(out):: pcfm(:, 2:) ! (knon, 2:klev) coefficient vitesse
-   PARAMETER (kstable=0.002)
-   REAL mixlen ! constante controlant longueur de melange
+     REAL, intent(out):: pcfh(:, 2:) ! (knon, 2:klev)
-   PARAMETER (mixlen=35.0)
+     ! coefficient chaleur et humidite)
-   REAL seuil ! au-dela l'inversion est consideree trop faible
-   PARAMETER (seuil=-0.02)
+     ! Local:
-   ! Variables locales:
+     ! Quelques constantes et options:
-   INTEGER i, k, invb(knon)
+     REAL prandtl
-   REAL zl2(knon)
+     PARAMETER (prandtl=0.4)
-   REAL zdthmin(knon), zdthdp
+     REAL kstable
+     PARAMETER (kstable=0.002)
-   !----------------------------------------------------------
+     REAL mixlen ! constante controlant longueur de melange
+     PARAMETER (mixlen=35.0)
-   ! Initialiser les sorties
+     REAL seuil ! au-dela l'inversion est consideree trop faible
-   DO k = 1, klev
+     PARAMETER (seuil=-0.02)
-      DO i = 1, knon
-         pcfm(i, k) = 0.0
+     INTEGER knon ! nombre de points a traiter
-         pcfh(i, k) = 0.0
+     INTEGER i, k
-      ENDDO
+     INTEGER invb(size(paprs, 1)) ! (knon)
-   ENDDO
+     REAL zl2(size(paprs, 1)) ! (knon)
+     REAL zdthmin(size(paprs, 1)) ! (knon)
-   ! Chercher la zone d'inversion forte
+     real zdthdp
-   DO i = 1, knon
+     !----------------------------------------------------------
-      invb(i) = klev
-      zdthmin(i)=0.0
+     knon = size(paprs, 1)
-   ENDDO
-   DO k = 2, klev/ 2 - 1
+     ! Initialiser les sorties
-      DO i = 1, knon
+     DO k = 2, klev
-         zdthdp = (t(i, k) - t(i, k + 1)) / (pplay(i, k) - pplay(i, k + 1)) &
+        DO i = 1, knon
-              - RD * 0.5 * (t(i, k) + t(i, k + 1)) / RCPD / paprs(i, k + 1)
+           pcfm(i, k) = 0.0
-         zdthdp = zdthdp * 100.
+           pcfh(i, k) = 0.0
-         IF (pplay(i, k) > 0.8 * paprs(i, 1) .AND. zdthdp < zdthmin(i)) THEN
+        ENDDO
-            zdthmin(i) = zdthdp
+     ENDDO
-            invb(i) = k
-         ENDIF
+     ! Chercher la zone d'inversion forte
-      ENDDO
-   ENDDO
+     DO i = 1, knon
+        invb(i) = klev
-   ! Introduire une diffusion:
+        zdthmin(i)=0.0
-   DO k = 2, klev
+     ENDDO
-      DO i = 1, knon
+     DO k = 2, klev/ 2 - 1
-         ! si on est sur ocean et s'il n'y a pas d'inversion ou si
+        DO i = 1, knon
-         ! l'inversion est trop faible:
+           zdthdp = (t(i, k) - t(i, k + 1)) / (pplay(i, k) - pplay(i, k + 1)) &
-         IF ((nsrf.EQ.is_oce) .AND.  &
+                - RD * 0.5 * (t(i, k) + t(i, k + 1)) / RCPD / paprs(i, k + 1)
-              ((invb(i).EQ.klev) .OR. (zdthmin(i) > seuil))) THEN
+           zdthdp = zdthdp * 100.
-            zl2(i)=(mixlen*MAX(0.0, (paprs(i, k)-paprs(i, klev+1)) &
+           IF (pplay(i, k) > 0.8 * paprs(i, 1) .AND. zdthdp < zdthmin(i)) THEN
-                 /(paprs(i, 2)-paprs(i, klev+1))))**2
+              zdthmin(i) = zdthdp
-            pcfm(i, k)= zl2(i)* kstable
+              invb(i) = k
-            pcfh(i, k) = pcfm(i, k) /prandtl ! h et m different
+           ENDIF
-         ENDIF
+        ENDDO
-      ENDDO
+     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

 Legend:



Removed from v.76
 


changed lines


 
Added in v.251
 Legend:



Removed from v.76
 


changed lines


 
Added in v.251
-Removed from v.76
+Added in v.251

	ViewVC Help
Powered by ViewVC 1.1.21