Sources/phylmd/coefkzmin.f

module coefkzmin_m

  IMPLICIT NONE

contains

  SUBROUTINE coefkzmin(ngrid, ypaprs, ypplay, yu, yv, yt, yq, ycoefm, km, kn)

    ! From LMDZ4/libf/phylmd/coefkzmin.F, version 1.1.1.1 2004/05/19 12:53:08

    ! Entrées modifiées en attendant une version où les zlev et zlay
    ! soient disponibles.

    USE dimphy, ONLY: klev, klon
    USE suphec_m, ONLY: rd, rg, rkappa

    integer, intent(in):: ngrid
    REAL, intent(in):: ypaprs(klon, klev+1), ypplay(klon, klev)
    REAL, intent(in):: yu(klon, klev), yv(klon, klev) ! wind, in m s-1
    REAL, intent(in):: yt(klon, klev) ! temperature, in K
    REAL, intent(in):: yq(klon, klev)
    REAL, intent(in):: ycoefm(:) ! (ngrid) drag coefficient

    REAL, intent(inout):: km(klon, klev)
    ! coefficient de diffusion turbulente de quantité de mouvement (au
    ! bas de chaque couche) (en sortie : la valeur à la fin du pas de
    ! temps), m2 s-1

    REAL, intent(inout):: kn(klon, klev)
    ! coefficient de diffusion turbulente des scalaires (au bas de
    ! chaque couche) (en sortie : la valeur à la fin du pas de temps), m2 s-1

    ! Local:

    real ustar(ngrid) ! u*
    real zlay(ngrid, klev) ! in m
    integer i, k
    real pblhmin(ngrid)
    real, parameter:: coriol = 1e-4

    REAL zlev(ngrid, 2: klev)
    ! altitude at level (interface between layer with same index), in m

    REAL teta(ngrid, klev)
    ! température potentielle au centre de chaque couche (la valeur au
    ! debut du pas de temps)

    real, PARAMETER:: kap = 0.4

    !---------------------------------------------------------------------

    ! Debut de la partie qui doit etre incluse a terme dans clmain.

    do i = 1, ngrid
       zlay(i, 1) = RD * yt(i, 1) * 2 / (ypaprs(i, 1) + ypplay(i, 1)) &
            * (ypaprs(i, 1) - ypplay(i, 1)) / RG
    enddo

    do k = 2, klev
       do i = 1, ngrid
          zlay(i, k) = zlay(i, k-1) + RD * 0.5 * (yt(i, k - 1) + yt(i, k)) &
               / ypaprs(i, k) * (ypplay(i, k - 1) - ypplay(i, k)) / RG
       enddo
    enddo

    do k=1, klev
       do i = 1, ngrid
          ! Attention : on passe la temperature potentielle virtuelle
          ! pour le calcul de K.
          teta(i, k) = yt(i, k) * (ypaprs(i, 1) / ypplay(i, k))**rkappa &
               * (1. + 0.61 * yq(i, k))
       enddo
    enddo

    forall (k = 2: klev) zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k-1))
    ustar = SQRT(ycoefm * (yu(:ngrid, 1)**2 + yv(:ngrid, 1)**2))

    ! Fin de la partie qui doit être incluse à terme dans clmain

    ! Cette routine est ecrite pour avoir en entree ustar, teta et zlev
    ! Ici, on a inclus le calcul de ces trois variables dans la routine
    ! coefkzmin en attendant une nouvelle version de la couche limite
    ! ou ces variables seront disponibles.

    ! Debut de la routine coefkzmin proprement dite

    pblhmin = 0.07 * ustar / coriol

    do k = 2, klev
       do i = 1, ngrid
          if (teta(i, 2) > teta(i, 1)) then
             kn(i, k) = kap * zlev(i, k) * ustar(i) &
                  * (max(1. - zlev(i, k) / pblhmin(i), 0.))**2
          else
             kn(i, k) = 0. ! min n'est utilisé que pour les SL stables
          endif
          km(i, k) = kn(i, k)
       enddo
    enddo

  end SUBROUTINE coefkzmin

end module coefkzmin_m
1	guez	47	module coefkzmin_m
2	guez	3
3	guez	47	IMPLICIT NONE
4	guez	3
5	guez	47	contains
6	guez	3
7	guez	47	SUBROUTINE coefkzmin(ngrid, ypaprs, ypplay, yu, yv, yt, yq, ycoefm, km, kn)
8	guez	3
9	guez	47	! From LMDZ4/libf/phylmd/coefkzmin.F, version 1.1.1.1 2004/05/19 12:53:08
10	guez	3
11	guez	47	! Entrées modifiées en attendant une version où les zlev et zlay
12			! soient disponibles.
13	guez	3
14	guez	47	USE dimphy, ONLY: klev, klon
15			USE suphec_m, ONLY: rd, rg, rkappa
16	guez	3
17	guez	47	integer, intent(in):: ngrid
18			REAL, intent(in):: ypaprs(klon, klev+1), ypplay(klon, klev)
19			REAL, intent(in):: yu(klon, klev), yv(klon, klev) ! wind, in m s-1
20			REAL, intent(in):: yt(klon, klev) ! temperature, in K
21			REAL, intent(in):: yq(klon, klev)
22	guez	70	REAL, intent(in):: ycoefm(:) ! (ngrid) drag coefficient
23	guez	3
24	guez	47	REAL, intent(inout):: km(klon, klev)
25			! coefficient de diffusion turbulente de quantité de mouvement (au
26			! bas de chaque couche) (en sortie : la valeur à la fin du pas de
27			! temps), m2 s-1
28	guez	3
29	guez	47	REAL, intent(inout):: kn(klon, klev)
30			! coefficient de diffusion turbulente des scalaires (au bas de
31			! chaque couche) (en sortie : la valeur à la fin du pas de temps), m2 s-1
32	guez	3
33	guez	47	! Local:
34	guez	3
35	guez	47	real ustar(ngrid) ! u*
36			real zlay(ngrid, klev) ! in m
37			integer i, k
38			real pblhmin(ngrid)
39			real, parameter:: coriol = 1e-4
40	guez	3
41	guez	47	REAL zlev(ngrid, 2: klev)
42			! altitude at level (interface between layer with same index), in m
43	guez	3
44	guez	47	REAL teta(ngrid, klev)
45			! température potentielle au centre de chaque couche (la valeur au
46			! debut du pas de temps)
47	guez	3
48	guez	47	real, PARAMETER:: kap = 0.4
49	guez	3
50	guez	47	!---------------------------------------------------------------------
51	guez	3
52	guez	47	! Debut de la partie qui doit etre incluse a terme dans clmain.
53	guez	3
54	guez	47	do i = 1, ngrid
55			zlay(i, 1) = RD * yt(i, 1) * 2 / (ypaprs(i, 1) + ypplay(i, 1)) &
56			* (ypaprs(i, 1) - ypplay(i, 1)) / RG
57			enddo
58	guez	3
59	guez	47	do k = 2, klev
60			do i = 1, ngrid
61			zlay(i, k) = zlay(i, k-1) + RD * 0.5 * (yt(i, k - 1) + yt(i, k)) &
62			/ ypaprs(i, k) * (ypplay(i, k - 1) - ypplay(i, k)) / RG
63			enddo
64			enddo
65	guez	3
66	guez	47	do k=1, klev
67			do i = 1, ngrid
68			! Attention : on passe la temperature potentielle virtuelle
69			! pour le calcul de K.
70			teta(i, k) = yt(i, k) * (ypaprs(i, 1) / ypplay(i, k))**rkappa &
71			* (1. + 0.61 * yq(i, k))
72			enddo
73			enddo
74	guez	3
75	guez	47	forall (k = 2: klev) zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k-1))
76	guez	70	ustar = SQRT(ycoefm * (yu(:ngrid, 1)2 + yv(:ngrid, 1)2))
77	guez	3
78	guez	47	! Fin de la partie qui doit être incluse à terme dans clmain
79	guez	3
80	guez	47	! Cette routine est ecrite pour avoir en entree ustar, teta et zlev
81			! Ici, on a inclus le calcul de ces trois variables dans la routine
82			! coefkzmin en attendant une nouvelle version de la couche limite
83			! ou ces variables seront disponibles.
84	guez	3
85	guez	47	! Debut de la routine coefkzmin proprement dite
86
87			pblhmin = 0.07 * ustar / coriol
88
89			do k = 2, klev
90			do i = 1, ngrid
91			if (teta(i, 2) > teta(i, 1)) then
92			kn(i, k) = kap * zlev(i, k) * ustar(i) &
93			* (max(1. - zlev(i, k) / pblhmin(i), 0.))**2
94			else
95			kn(i, k) = 0. ! min n'est utilisé que pour les SL stables
96			endif
97			km(i, k) = kn(i, k)
98			enddo
99			enddo
100
101			end SUBROUTINE coefkzmin
102
103			end module coefkzmin_m