trunk/phylmd/calbeta.f

SUBROUTINE calbeta(dtime, indice, knon, snow, qsol, vbeta, vcal, vdif)
  USE dimens_m
  USE indicesol
  USE dimphy
  USE conf_gcm_m
  USE suphec_m
  IMPLICIT NONE
  ! ======================================================================
  ! Auteur(s): Z.X. Li (LMD/CNRS) (adaptation du GCM du LMD)
  ! date: 19940414
  ! ======================================================================

  ! Calculer quelques parametres pour appliquer la couche limite
  ! ------------------------------------------------------------
  REAL tau_gl ! temps de relaxation pour la glace de mer
  ! cc      PARAMETER (tau_gl=86400.0*30.0)
  PARAMETER (tau_gl=86400.0*5.0)
  REAL mx_eau_sol
  PARAMETER (mx_eau_sol=150.0)

  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))

  INTEGER i

  REAL dtime
  REAL snow(klon), qsol(klon)
  INTEGER indice, knon

  REAL vbeta(klon)
  REAL vcal(klon)
  REAL vdif(klon)


  IF (indice==is_oce) THEN
    DO i = 1, knon
      vcal(i) = 0.0
      vbeta(i) = 1.0
      vdif(i) = 0.0
    END DO
  END IF

  IF (indice==is_sic) THEN
    DO i = 1, knon
      vcal(i) = calice
      IF (snow(i)>0.0) vcal(i) = calsno
      vbeta(i) = 1.0
      vdif(i) = 1.0/tau_gl
      ! cc          vdif(i) = calice/tau_gl ! c'etait une erreur
    END DO
  END IF

  IF (indice==is_ter) THEN
    DO i = 1, knon
      vcal(i) = calsol
      IF (snow(i)>0.0) vcal(i) = calsno
      vbeta(i) = min(2.0*qsol(i)/mx_eau_sol, 1.0)
      vdif(i) = 0.0
    END DO
  END IF

  IF (indice==is_lic) THEN
    DO i = 1, knon
      vcal(i) = calice
      IF (snow(i)>0.0) vcal(i) = calsno
      vbeta(i) = 1.0
      vdif(i) = 0.0
    END DO
  END IF

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