suclop.F90

Go to the documentation of this file.

SUBROUTINE suclop

!**** *SUCLOP*  - INITIALIZE COMMON YOECLOP

!     PURPOSE.
!     --------
!           INITIALIZE YOMCLOP, WITH CLOUD OPTICAL PARAMETERS

!**   INTERFACE.
!     ----------
!        *CALL*  SUCLOP
!     FROM *SUECRAD*

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

!        IMPLICIT ARGUMENTS :
!        --------------------
!        COMMON YOECLOP

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

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

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

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

!     MODIFICATIONS.
!     --------------
!        ORIGINAL : 92-02-29
!        M.Hamrud      01-Oct-2003 CY28 Cleaning
!     ------------------------------------------------------------------

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

USE yoeclop  , ONLY : ryfwca   ,ryfwcb   ,ryfwcc   ,ryfwcd   ,&
 & ryfwce   ,ryfwcf   ,rebcua   ,rebcub   ,rebcuc   ,&
 & rebcud   ,rebcue   ,rebcuf   ,rebcug   ,rebcuh   ,&
 & reffia   ,reffib   ,rtiw     ,rriw  

IMPLICIT NONE

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

!* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992

! SW : 2 spectral intervals

REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook('SUCLOP',0,zhook_handle)
rebcua(1)= 3.448e-03_jprb
rebcua(2)= 3.448e-03_jprb
rebcub(1)= 2.431_jprb
rebcub(2)= 2.431_jprb
rebcuc(1)= 0.99999_jprb
rebcuc(2)= 0.975634_jprb
rebcud(1)= 0._jprb
rebcud(2)= 2.487e-04_jprb
rebcue(1)= 0.7661_jprb
rebcue(2)= 0.7866_jprb
rebcuf(1)= 5.851e-04_jprb
rebcuf(2)= 5.937e-04_jprb

! LW : spectrally averaged with reference Planck function at 257 K

rebcug= 1.07677_jprb
rebcuh= 0.00267_jprb

! Ice particle Effective Radius as a function of LWC

reffia= 40._jprb
reffib=  0._jprb

!* Water cloud properties - from Fouquart (1987)

! SW : 2 spectral intervals: parameters as a function of Reff

ryfwca(1)= 0._jprb
ryfwca(2)= 0._jprb
ryfwcb(1)= 1.5_jprb
ryfwcb(2)= 1.5_jprb
ryfwcc(1)= 0.9999_jprb
ryfwcc(2)= 0.9988_jprb
ryfwcd(1)= 5.000e-04_jprb
ryfwcd(2)= 2.500e-03_jprb
ryfwce(1)= 0.5_jprb
ryfwce(2)= 0.05_jprb
ryfwcf(1)= 0.865_jprb
ryfwcf(2)= 0.910_jprb

!* Liquid/Solid water transition

rtiw= 263._jprb
rriw= 20._jprb
IF (lhook) CALL dr_hook('SUCLOP',1,zhook_handle)

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

END SUBROUTINE suclop