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

Diff of /trunk/phylmd/Interface_surf/coefkz2.f

Parent Directory | Revision Log | View Patch Patch

-trunk/libf/phylmd/coefkz2.f90
revision 57 by guez,
Mon Jan 30 12:54:02 2012 UTC
+trunk/Sources/phylmd/coefkz2.f
revision 248 by guez,
Fri Jan  5 16:40:13 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, knon, 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
-   ! t--------input-R- temperature (K)
+     ! ou une forte inversion est presente
+     ! On peut dire qu'il represente la convection peu profonde
-   ! pcfm-----output-R- coefficients a calculer (vitesse)
-   ! pcfh-----output-R- coefficients a calculer (chaleur et humidite)
+     use dimens_m
+     use indicesol
-   ! Arguments:
+     use dimphy
+     use conf_gcm_m
-   INTEGER knon, nsrf
+     use SUPHEC_M
-   REAL paprs(klon, klev+1), pplay(klon, klev)
-   REAL t(klon, klev)
+     ! Arguments:
+     ! nsrf-----input-I- indicateur de la nature du sol
-   REAL pcfm(klon, klev), pcfh(klon, klev)
+     ! knon-----input-I- nombre de points a traiter
+     ! paprs----input-R- pression a chaque intercouche (en Pa)
-   ! Quelques constantes et options:
+     ! pplay----input-R- pression au milieu de chaque couche (en Pa)
+     ! t--------input-R- temperature (K)
-   REAL prandtl
-   PARAMETER (prandtl=0.4)
+     ! pcfm-----output-R- coefficients a calculer (vitesse)
-   REAL kstable
+     ! pcfh-----output-R- coefficients a calculer (chaleur et humidite)
-   PARAMETER (kstable=0.002)
-   REAL mixlen ! constante controlant longueur de melange
+     ! Arguments:
-   PARAMETER (mixlen=35.0)
-   REAL seuil ! au-dela l'inversion est consideree trop faible
+     INTEGER knon
-   PARAMETER (seuil=-0.02)
+     integer, intent(in):: nsrf
+     REAL, intent(in):: paprs(klon, klev+1), pplay(klon, klev)
-   ! Variables locales:
+     REAL, intent(in):: t(klon, klev)
-   INTEGER i, k, invb(knon)
+     REAL, intent(out):: pcfm(:, 2:), pcfh(:, 2:) ! (knon, 2:klev)
-   REAL zl2(knon)
-   REAL zdthmin(knon), zdthdp
+     ! Quelques constantes et options:
-   !----------------------------------------------------------
+     REAL prandtl
+     PARAMETER (prandtl=0.4)
-   ! Initialiser les sorties
+     REAL kstable
-   DO k = 1, klev
+     PARAMETER (kstable=0.002)
-      DO i = 1, knon
+     REAL mixlen ! constante controlant longueur de melange
-         pcfm(i, k) = 0.0
+     PARAMETER (mixlen=35.0)
-         pcfh(i, k) = 0.0
+     REAL seuil ! au-dela l'inversion est consideree trop faible
-      ENDDO
+     PARAMETER (seuil=-0.02)
-   ENDDO
+     ! Variables locales:
-   ! Chercher la zone d'inversion forte
+     INTEGER i, k, invb(knon)
-   DO i = 1, knon
+     REAL zl2(knon)
-      invb(i) = klev
+     REAL zdthmin(knon), zdthdp
-      zdthmin(i)=0.0
-   ENDDO
+     !----------------------------------------------------------
-   DO k = 2, klev/ 2 - 1
-      DO i = 1, knon
+     ! Initialiser les sorties
-         zdthdp = (t(i, k) - t(i, k + 1)) / (pplay(i, k) - pplay(i, k + 1)) &
+     DO k = 2, klev
-              - RD * 0.5 * (t(i, k) + t(i, k + 1)) / RCPD / paprs(i, k + 1)
+        DO i = 1, knon
-         zdthdp = zdthdp * 100.
+           pcfm(i, k) = 0.0
-         IF (pplay(i, k) > 0.8 * paprs(i, 1) .AND. zdthdp < zdthmin(i)) THEN
+           pcfh(i, k) = 0.0
-            zdthmin(i) = zdthdp
+        ENDDO
-            invb(i) = k
+     ENDDO
-         ENDIF
-      ENDDO
+     ! Chercher la zone d'inversion forte
-   ENDDO
+     DO i = 1, knon
-   ! Introduire une diffusion:
+        invb(i) = klev
-   DO k = 2, klev
+        zdthmin(i)=0.0
-      DO i = 1, knon
+     ENDDO
-         ! si on est sur ocean et s'il n'y a pas d'inversion ou si
+     DO k = 2, klev/ 2 - 1
-         ! l'inversion est trop faible:
+        DO i = 1, knon
-         IF ((nsrf.EQ.is_oce) .AND.  &
+           zdthdp = (t(i, k) - t(i, k + 1)) / (pplay(i, k) - pplay(i, k + 1)) &
-              ((invb(i).EQ.klev) .OR. (zdthmin(i) > seuil))) THEN
+                - RD * 0.5 * (t(i, k) + t(i, k + 1)) / RCPD / paprs(i, k + 1)
-            zl2(i)=(mixlen*MAX(0.0, (paprs(i, k)-paprs(i, klev+1)) &
+           zdthdp = zdthdp * 100.
-                 /(paprs(i, 2)-paprs(i, klev+1))))**2
+           IF (pplay(i, k) > 0.8 * paprs(i, 1) .AND. zdthdp < zdthmin(i)) THEN
-            pcfm(i, k)= zl2(i)* kstable
+              zdthmin(i) = zdthdp
-            pcfh(i, k) = pcfm(i, k) /prandtl ! h et m different
+              invb(i) = k
-         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.57
 


changed lines


 
Added in v.248
 Legend:



Removed from v.57
 


changed lines


 
Added in v.248
-Removed from v.57
+Added in v.248

	ViewVC Help
Powered by ViewVC 1.1.21