sucld.F90

Go to the documentation of this file.

SUBROUTINE sucld ( KLEV , PETA )

!**** *SUCLD*   - INITIALIZE COMMON YOECLD CONTROLLING *CLOUD*

!     PURPOSE.
!     --------
!           INITIALIZE YOECLD

!**   INTERFACE.
!     ----------
!        CALL *SUCLD* FROM *SUPHEC*
!              -----        ------

!        EXPLICIT ARGUMENTS :
!        --------------------
!        NONE

!        IMPLICIT ARGUMENTS :
!        --------------------
!        COMMON YOECLD

!     METHOD.
!     -------
!        SEE DOCUMENTATION

!     EXTERNALS.
!     ----------
!        NONE

!     REFERENCE.
!     ----------
!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE
!     "INTEGRATED FORECASTING SYSTEM"

!     AUTHOR.
!     -------
!        JEAN-JACQUES MORCRETTE  *ECMWF*

!     MODIFICATIONS.
!     --------------
!        ORIGINAL : 89-12-15
!        M.Hamrud      01-Oct-2003 CY28 Cleaning

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

USE parkind1  ,ONLY : jpim     ,jprb
USE yomhook   ,ONLY : lhook,   dr_hook

USE yoecld   , ONLY : ceta     ,ranva    ,ranvb    ,ranvh    ,&
 & rcca     ,rccb     ,rccc     ,rcfct    ,rclwmr   ,&
 & rcscal   ,retahb   ,retamb   ,rloia    ,rloib    ,&
 & rloic    ,rloid    ,rlonia   ,rlonib   ,rrhh     ,&
 & rrhl     ,rrhm     ,rgammas  ,repscr   ,repsec   ,&
 & lomega  

IMPLICIT NONE

INTEGER(KIND=JPIM),INTENT(IN)    :: KLEV 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PETA(klev) 
!      ----------------------------------------------------------------

INTEGER(KIND=JPIM) :: JK
REAL(KIND=JPRB) :: ZHOOK_HANDLE

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

!*       1.    SET VALUES
!              ----------

IF (lhook) CALL dr_hook('SUCLD',0,zhook_handle)
ranva  = 2._jprb
ranvb  = 0.3_jprb
ranvh  = 0.4_jprb
rcca   = 0.125_jprb
rccb   = 1.5_jprb
rccc   = 0.8_jprb
rcfct  = 0.400_jprb
rcscal = 1.0e+11_jprb

retahb = 0.45_jprb
retamb = 0.80_jprb

rloia  = 1.0e+02_jprb
rloib  =-10.00_jprb
rloic  =-0.9_jprb
rloid  = 5.0_jprb

rlonia = -0.1_jprb
rlonib = -10.0_jprb

rrhh   = 0.9_jprb
rrhm   = 0.8_jprb
rrhl   = 0.70_jprb

rgammas= 0.05_jprb
rclwmr = 1.e-04_jprb
lomega =.true.

repsec = 1.0e-12_jprb
repscr = 1.0e-12_jprb

DO jk=1,klev
  ceta(jk)=peta(jk)
ENDDO

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

IF (lhook) CALL dr_hook('SUCLD',1,zhook_handle)
END SUBROUTINE sucld