/[lmdze]/trunk/Sources/phylmd/coefkz2.f
ViewVC logotype

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

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/libf/phylmd/coefkz2.f revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC trunk/Sources/phylmd/coefkz2.f revision 234 by guez, Tue Nov 7 12:20:42 2017 UTC
# Line 1  Line 1 
1        SUBROUTINE coefkz2(nsrf, knon, paprs, pplay,t,  module coefkz2_m
2       .                  pcfm, pcfh)  
3        use dimens_m    IMPLICIT none
4        use indicesol  
5        use dimphy  contains
6        use iniprint  
7        use YOMCST    SUBROUTINE coefkz2(nsrf, knon, paprs, pplay, t, pcfm, pcfh)
8        IMPLICIT none  
9  c======================================================================      ! J'introduit un peu de diffusion sauf dans les endroits
10  c J'introduit un peu de diffusion sauf dans les endroits      ! ou une forte inversion est presente
11  c ou une forte inversion est presente      ! On peut dire qu'il represente la convection peu profonde
12  c On peut dire qu'il represente la convection peu profonde  
13  c      use dimens_m
14  c Arguments:      use indicesol
15  c nsrf-----input-I- indicateur de la nature du sol      use dimphy
16  c knon-----input-I- nombre de points a traiter      use conf_gcm_m
17  c paprs----input-R- pression a chaque intercouche (en Pa)      use SUPHEC_M
18  c pplay----input-R- pression au milieu de chaque couche (en Pa)  
19  c t--------input-R- temperature (K)      ! Arguments:
20  c      ! nsrf-----input-I- indicateur de la nature du sol
21  c pcfm-----output-R- coefficients a calculer (vitesse)      ! knon-----input-I- nombre de points a traiter
22  c pcfh-----output-R- coefficients a calculer (chaleur et humidite)      ! paprs----input-R- pression a chaque intercouche (en Pa)
23  c======================================================================      ! pplay----input-R- pression au milieu de chaque couche (en Pa)
24  c      ! t--------input-R- temperature (K)
25  c Arguments:  
26  c      ! pcfm-----output-R- coefficients a calculer (vitesse)
27        INTEGER knon, nsrf      ! pcfh-----output-R- coefficients a calculer (chaleur et humidite)
28        REAL paprs(klon,klev+1), pplay(klon,klev)  
29        REAL t(klon,klev)      ! Arguments:
30  c  
31        REAL pcfm(klon,klev), pcfh(klon,klev)      INTEGER knon, nsrf
32  c      REAL paprs(klon, klev+1), pplay(klon, klev)
33  c Quelques constantes et options:      REAL t(klon, klev)
34  c  
35        REAL prandtl      REAL, intent(out):: pcfm(:, :), pcfh(:, :) ! (knon, klev)
36        PARAMETER (prandtl=0.4)  
37        REAL kstable      ! Quelques constantes et options:
38        PARAMETER (kstable=0.002)  
39  ccc      PARAMETER (kstable=0.001)      REAL prandtl
40        REAL mixlen ! constante controlant longueur de melange      PARAMETER (prandtl=0.4)
41        PARAMETER (mixlen=35.0)      REAL kstable
42        REAL seuil ! au-dela l'inversion est consideree trop faible      PARAMETER (kstable=0.002)
43        PARAMETER (seuil=-0.02)      REAL mixlen ! constante controlant longueur de melange
44  ccc      PARAMETER (seuil=-0.04)      PARAMETER (mixlen=35.0)
45  ccc      PARAMETER (seuil=-0.06)      REAL seuil ! au-dela l'inversion est consideree trop faible
46  ccc      PARAMETER (seuil=-0.09)      PARAMETER (seuil=-0.02)
47  c  
48  c Variables locales:      ! Variables locales:
49  c  
50        INTEGER i, k, invb(knon)      INTEGER i, k, invb(knon)
51        REAL zl2(knon)      REAL zl2(knon)
52        REAL zdthmin(knon), zdthdp      REAL zdthmin(knon), zdthdp
53  c  
54  c Initialiser les sorties      !----------------------------------------------------------
55  c  
56        DO k = 1, klev      ! Initialiser les sorties
57        DO i = 1, knon      DO k = 1, klev
58           pcfm(i,k) = 0.0         DO i = 1, knon
59           pcfh(i,k) = 0.0            pcfm(i, k) = 0.0
60        ENDDO            pcfh(i, k) = 0.0
61        ENDDO         ENDDO
62  c      ENDDO
63  c Chercher la zone d'inversion forte  
64  c      ! Chercher la zone d'inversion forte
65        DO i = 1, knon  
66           invb(i) = klev      DO i = 1, knon
67           zdthmin(i)=0.0         invb(i) = klev
68        ENDDO         zdthmin(i)=0.0
69        DO k = 2, klev/2-1      ENDDO
70        DO i = 1, knon      DO k = 2, klev/ 2 - 1
71           zdthdp = (t(i,k)-t(i,k+1))/(pplay(i,k)-pplay(i,k+1))         DO i = 1, knon
72       .          - RD * 0.5*(t(i,k)+t(i,k+1))/RCPD/paprs(i,k+1)            zdthdp = (t(i, k) - t(i, k + 1)) / (pplay(i, k) - pplay(i, k + 1)) &
73           zdthdp = zdthdp * 100.0                 - RD * 0.5 * (t(i, k) + t(i, k + 1)) / RCPD / paprs(i, k + 1)
74           IF (pplay(i,k).GT.0.8*paprs(i,1) .AND.            zdthdp = zdthdp * 100.
75       .       zdthdp.LT.zdthmin(i) ) THEN            IF (pplay(i, k) > 0.8 * paprs(i, 1) .AND. zdthdp < zdthmin(i)) THEN
76              zdthmin(i) = zdthdp               zdthmin(i) = zdthdp
77              invb(i) = k               invb(i) = k
78           ENDIF            ENDIF
79        ENDDO         ENDDO
80        ENDDO      ENDDO
81  c  
82  c Introduire une diffusion:      ! Introduire une diffusion:
83  c      DO k = 2, klev
84        DO k = 2, klev         DO i = 1, knon
85        DO i = 1, knon            ! si on est sur ocean et s'il n'y a pas d'inversion ou si
86  cIM cf FH/GK   IF ( (nsrf.NE.is_oce) .OR.  ! si ce n'est pas sur l'ocean            ! l'inversion est trop faible:
87  cIM cf FH/GK  .     (invb(i).EQ.klev) .OR. ! s'il n'y a pas d'inversion            IF ((nsrf.EQ.is_oce) .AND.  &  
88        !IM cf JLD/ GKtest TERkz2                 ((invb(i).EQ.klev) .OR. (zdthmin(i) > seuil))) THEN
89        ! IF ( (nsrf.EQ.is_ter) .OR.  ! si on est sur la terre               zl2(i)=(mixlen*MAX(0.0, (paprs(i, k)-paprs(i, klev+1)) &
90        ! fin GKtest                    /(paprs(i, 2)-paprs(i, klev+1))))**2
91        IF ( (nsrf.EQ.is_oce) .AND.  ! si on est sur ocean et si               pcfm(i, k)= zl2(i)* kstable
92       .     ( (invb(i).EQ.klev) .OR.      ! s'il n'y a pas d'inversion               pcfh(i, k) = pcfm(i, k) /prandtl ! h et m different
93       .     (zdthmin(i).GT.seuil) ) )THEN ! si l'inversion est trop faible            ENDIF
94           zl2(i)=(mixlen*MAX(0.0,(paprs(i,k)-paprs(i,klev+1))         ENDDO
95       .                       /(paprs(i,2)-paprs(i,klev+1)) ))**2      ENDDO
96           pcfm(i,k)= zl2(i)* kstable  
97           pcfh(i,k) = pcfm(i,k) /prandtl ! h et m different    END SUBROUTINE coefkz2
98        ENDIF  
99        ENDDO  end module coefkz2_m
       ENDDO  
 c  
       RETURN  
       END  
Legend:
Removed from v.3  
changed lines
  Added in v.234
  ViewVC Help
Powered by ViewVC 1.1.21