/[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 38 by guez, Thu Jan 6 17:52:19 2011 UTC trunk/Sources/phylmd/coefkz2.f revision 233 by guez, Tue Nov 7 10:52:46 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 SUPHEC_M    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 pcfm(klon, klev), pcfh(klon, 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.38  
changed lines
  Added in v.233
  ViewVC Help
Powered by ViewVC 1.1.21