source: LMDZ6/branches/Amaury_dev/libf/phylmd/calbeta.F90 @ 5136

Last change on this file since 5136 was 5134, checked in by abarral, 4 months ago

Replace academic.h, alpale.h, comdissip.h, comdissipn.h, comdissnew.h by modules
Remove unused clesph0.h
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.6 KB
RevLine 
[5099]1
[782]2! $Header$
[5099]3
[782]4SUBROUTINE calbeta(dtime,indice,knon,snow,qsol, &
5     vbeta,vcal,vdif)
6
7  USE dimphy
[1785]8  USE indice_sol_mod
[3780]9
[5134]10  IMPLICIT NONE
[3780]11
[4593]12  INCLUDE "flux_arp.h"
[3780]13
[782]14!======================================================================
15! Auteur(s): Z.X. Li (LMD/CNRS) (adaptation du GCM au LMD)
16! date: 19940414
17!======================================================================
[5099]18
[782]19! Calculer quelques parametres pour appliquer la couche limite
[1785]20! ------------------------------------------------------------
[782]21! Variables d'entrees
22!****************************************************************************************
23  REAL, INTENT(IN)                   :: dtime
24  INTEGER, INTENT(IN)                :: indice
25  INTEGER, INTENT(IN)                :: knon
26  REAL, DIMENSION(klon), INTENT(IN)  :: snow
27  REAL, DIMENSION(klon), INTENT(IN)  :: qsol
28
29 
30! Variables de sorties
31!****************************************************************************************
32  REAL, DIMENSION(klon), INTENT(OUT) :: vbeta
33  REAL, DIMENSION(klon), INTENT(OUT) :: vcal
34  REAL, DIMENSION(klon), INTENT(OUT) :: vdif
35
36! Variables locales
37!****************************************************************************************
38  REAL, PARAMETER :: tau_gl=86400.0*5.0 ! temps de relaxation pour la glace de mer
39!cc      PARAMETER (tau_gl=86400.0*30.0)
40  REAL, PARAMETER :: mx_eau_sol=150.0
41  REAL, PARAMETER :: calsol=1.0/(2.5578E+06*0.15)
42  REAL, PARAMETER :: calsno=1.0/(2.3867E+06*0.15)
43  REAL, PARAMETER :: calice=1.0/(5.1444E+06*0.15)
44 
45  INTEGER         :: i
46
47!**************************************************************************************** 
48   
49  vbeta(:) = 0.0
50  vcal(:) = 0.0
51  vdif(:) = 0.0
52 
[5082]53  IF (indice==is_oce) THEN
[782]54     DO i = 1, knon
55        vcal(i)   = 0.0
56        vbeta(i)  = 1.0
57        vdif(i) = 0.0
58     ENDDO
59  ENDIF
60 
[5082]61  IF (indice==is_sic) THEN
[782]62     DO i = 1, knon
63        vcal(i) = calice
[5082]64        IF (snow(i) > 0.0) vcal(i) = calsno
[782]65        vbeta(i)  = 1.0
66        vdif(i) = 1.0/tau_gl
67!          vdif(i) = calice/tau_gl ! c'etait une erreur
68     ENDDO
69  ENDIF
70 
[5082]71  IF (indice==is_ter) THEN
[782]72     DO i = 1, knon
73        vcal(i) = calsol
[5082]74        IF (snow(i) > 0.0) vcal(i) = calsno
[782]75        vbeta(i)  = MIN(2.0*qsol(i)/mx_eau_sol, 1.0)
76        vdif(i) = 0.0
77     ENDDO
78  ENDIF
79 
[5082]80  IF (indice==is_lic) THEN
[782]81     DO i = 1, knon
82        vcal(i) = calice
[5082]83        IF (snow(i) > 0.0) vcal(i) = calsno
[782]84        vbeta(i)  = 1.0
85        vdif(i) = 0.0
86     ENDDO
87  ENDIF
[3780]88
89  ! EV: when beta is prescribed for 1D cases:
[5082]90  IF (knon==1 .AND. ok_prescr_beta) THEN
[3780]91     DO i = 1, knon
92          vbeta(i)=betaevap
93      ENDDO
94  ENDIF
[782]95 
96END SUBROUTINE calbeta
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
Note: See TracBrowser for help on using the repository browser.