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

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

Parent Directory | Revision Log | View Patch Patch

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

 Legend:



Removed from v.233
 


changed lines


 
Added in v.251
 Legend:



Removed from v.233
 


changed lines


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

	ViewVC Help
Powered by ViewVC 1.1.21