libf/phylmd/calbeta.f

      SUBROUTINE calbeta(dtime,indice,knon,snow,qsol,
     .                    vbeta,vcal,vdif)
      use dimens_m
      use indicesol
      use dimphy
      use iniprint
      use YOMCST
      IMPLICIT none
c======================================================================
c Auteur(s): Z.X. Li (LMD/CNRS) (adaptation du GCM du LMD)
c date: 19940414
c======================================================================
c
c Calculer quelques parametres pour appliquer la couche limite
c ------------------------------------------------------------
      REAL tau_gl ! temps de relaxation pour la glace de mer
ccc      PARAMETER (tau_gl=86400.0*30.0)
      PARAMETER (tau_gl=86400.0*5.0)
      REAL mx_eau_sol
      PARAMETER (mx_eau_sol=150.0)
c
      REAL calsol, calsno, calice ! epaisseur du sol: 0.15 m
      PARAMETER (calsol=1.0/(2.5578E+06*0.15))
      PARAMETER (calsno=1.0/(2.3867E+06*0.15))
      PARAMETER (calice=1.0/(5.1444E+06*0.15))
C
      INTEGER i
c
      REAL dtime
      REAL snow(klon), qsol(klon)
      INTEGER indice, knon
C
      REAL vbeta(klon)
      REAL vcal(klon)
      REAL vdif(klon)
C

      IF (indice.EQ.is_oce) THEN
      DO i = 1, knon
          vcal(i)   = 0.0
          vbeta(i)  = 1.0
          vdif(i) = 0.0
      ENDDO
      ENDIF
c
      IF (indice.EQ.is_sic) THEN
      DO i = 1, knon
          vcal(i) = calice
          IF (snow(i) .GT. 0.0) vcal(i) = calsno
          vbeta(i)  = 1.0
          vdif(i) = 1.0/tau_gl
ccc          vdif(i) = calice/tau_gl ! c'etait une erreur
      ENDDO
      ENDIF
c
      IF (indice.EQ.is_ter) THEN
      DO i = 1, knon
          vcal(i) = calsol
          IF (snow(i) .GT. 0.0) vcal(i) = calsno
          vbeta(i)  = MIN(2.0*qsol(i)/mx_eau_sol, 1.0)
          vdif(i) = 0.0
      ENDDO
      ENDIF
c
      IF (indice.EQ.is_lic) THEN
      DO i = 1, knon
          vcal(i) = calice
          IF (snow(i) .GT. 0.0) vcal(i) = calsno
          vbeta(i)  = 1.0
          vdif(i) = 0.0
      ENDDO
      ENDIF
c
      RETURN
      END
1	guez	3	SUBROUTINE calbeta(dtime,indice,knon,snow,qsol,
2			. vbeta,vcal,vdif)
3			use dimens_m
4			use indicesol
5			use dimphy
6			use iniprint
7			use YOMCST
8			IMPLICIT none
9			c======================================================================
10			c Auteur(s): Z.X. Li (LMD/CNRS) (adaptation du GCM du LMD)
11			c date: 19940414
12			c======================================================================
13			c
14			c Calculer quelques parametres pour appliquer la couche limite
15			c ------------------------------------------------------------
16			REAL tau_gl ! temps de relaxation pour la glace de mer
17			ccc PARAMETER (tau_gl=86400.0*30.0)
18			PARAMETER (tau_gl=86400.0*5.0)
19			REAL mx_eau_sol
20			PARAMETER (mx_eau_sol=150.0)
21			c
22			REAL calsol, calsno, calice ! epaisseur du sol: 0.15 m
23			PARAMETER (calsol=1.0/(2.5578E+06*0.15))
24			PARAMETER (calsno=1.0/(2.3867E+06*0.15))
25			PARAMETER (calice=1.0/(5.1444E+06*0.15))
26			C
27			INTEGER i
28			c
29			REAL dtime
30			REAL snow(klon), qsol(klon)
31			INTEGER indice, knon
32			C
33			REAL vbeta(klon)
34			REAL vcal(klon)
35			REAL vdif(klon)
36			C
37
38			IF (indice.EQ.is_oce) THEN
39			DO i = 1, knon
40			vcal(i) = 0.0
41			vbeta(i) = 1.0
42			vdif(i) = 0.0
43			ENDDO
44			ENDIF
45			c
46			IF (indice.EQ.is_sic) THEN
47			DO i = 1, knon
48			vcal(i) = calice
49			IF (snow(i) .GT. 0.0) vcal(i) = calsno
50			vbeta(i) = 1.0
51			vdif(i) = 1.0/tau_gl
52			ccc vdif(i) = calice/tau_gl ! c'etait une erreur
53			ENDDO
54			ENDIF
55			c
56			IF (indice.EQ.is_ter) THEN
57			DO i = 1, knon
58			vcal(i) = calsol
59			IF (snow(i) .GT. 0.0) vcal(i) = calsno
60			vbeta(i) = MIN(2.0*qsol(i)/mx_eau_sol, 1.0)
61			vdif(i) = 0.0
62			ENDDO
63			ENDIF
64			c
65			IF (indice.EQ.is_lic) THEN
66			DO i = 1, knon
67			vcal(i) = calice
68			IF (snow(i) .GT. 0.0) vcal(i) = calsno
69			vbeta(i) = 1.0
70			vdif(i) = 0.0
71			ENDDO
72			ENDIF
73			c
74			RETURN
75			END