phylmd/Interface_surf/calbeta.f

module calbeta_m

  IMPLICIT NONE

contains

  SUBROUTINE calbeta(indice, snow, qsol, vbeta, vcal, vdif)

    ! Author: Z. X. Li (LMD/CNRS)
    ! Date: April 14th, 1994

    ! Calcul de quelques paramètres pour appliquer la couche limite.

    USE indicesol, ONLY: is_lic, is_oce, is_sic, is_ter

    INTEGER, intent(in):: indice
    REAL, intent(in):: snow(:), qsol(:) ! (knon)
    REAL, intent(out):: vbeta(:), vcal(:), vdif(:) ! (knon)

    ! Local:

    REAL, PARAMETER:: tau_gl = 86400. * 5. 
    ! temps de relaxation pour la glace de mer

    REAL, PARAMETER:: max_eau_sol = 150. ! in kg m-2

    ! Pour une épaisseur du sol de 15 cm :
    REAL, PARAMETER:: calsol = 1. / (2.5578E6 * 0.15)
    REAL, PARAMETER:: calsno = 1. / (2.3867E6 * 0.15)
    REAL, PARAMETER:: calice = 1. / (5.1444E6 * 0.15)

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

    select case (indice)
    case(is_oce)
       vbeta = 1.
       vcal = 0.
       vdif = 0.

    case (is_sic)
       vbeta = 1.
       vcal = merge(calsno, calice, snow > 0.)
       vdif = 1. / tau_gl

    case (is_ter)
       vbeta = min(2. * qsol / max_eau_sol, 1.)
       vcal = merge(calsno, calsol, snow > 0.)
       vdif = 0.

    case (is_lic)
       vbeta = 1.
       vcal = merge(calsno, calice, snow > 0.)
       vdif = 0.
    END select

  END SUBROUTINE calbeta

end module calbeta_m
1	guez	101	module calbeta_m
2
3	guez	81	IMPLICIT NONE
4	guez	3
5	guez	101	contains
6	guez	81
7	guez	101	SUBROUTINE calbeta(indice, snow, qsol, vbeta, vcal, vdif)
8	guez	81
9	guez	101	! Author: Z. X. Li (LMD/CNRS)
10			! Date: April 14th, 1994
11	guez	81
12	guez	101	! Calcul de quelques paramètres pour appliquer la couche limite.
13	guez	81
14	guez	101	USE indicesol, ONLY: is_lic, is_oce, is_sic, is_ter
15	guez	81
16	guez	101	INTEGER, intent(in):: indice
17			REAL, intent(in):: snow(:), qsol(:) ! (knon)
18			REAL, intent(out):: vbeta(:), vcal(:), vdif(:) ! (knon)
19	guez	81
20	guez	101	! Local:
21	guez	81
22	guez	101	REAL, PARAMETER:: tau_gl = 86400. * 5.
23			! temps de relaxation pour la glace de mer
24	guez	81
25	guez	101	REAL, PARAMETER:: max_eau_sol = 150. ! in kg m-2
26	guez	81
27	guez	101	! Pour une épaisseur du sol de 15 cm :
28			REAL, PARAMETER:: calsol = 1. / (2.5578E6 * 0.15)
29			REAL, PARAMETER:: calsno = 1. / (2.3867E6 * 0.15)
30			REAL, PARAMETER:: calice = 1. / (5.1444E6 * 0.15)
31	guez	81
32	guez	101	!------------------------------------------------------------
33
34			select case (indice)
35			case(is_oce)
36			vbeta = 1.
37			vcal = 0.
38			vdif = 0.
39
40			case (is_sic)
41			vbeta = 1.
42			vcal = merge(calsno, calice, snow > 0.)
43			vdif = 1. / tau_gl
44
45			case (is_ter)
46			vbeta = min(2. * qsol / max_eau_sol, 1.)
47			vcal = merge(calsno, calsol, snow > 0.)
48			vdif = 0.
49
50			case (is_lic)
51			vbeta = 1.
52			vcal = merge(calsno, calice, snow > 0.)
53			vdif = 0.
54			END select
55
56			END SUBROUTINE calbeta
57
58			end module calbeta_m