suclop15.F90

Go to the documentation of this file.

!OPTIONS XOPT(NOEVAL)
SUBROUTINE suclop15

!**** *SUCLOP15*  - INITIALIZE COMMON YOMCLOP15
!****               FROZEN VERSION (CYCLE 15) OF SUCLOP

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

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

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

!        IMPLICIT ARGUMENTS :
!        --------------------
!        COMMON YOMCLOP15

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

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

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

!     AUTHOR.
!     -------
!        96-11: Ph. Dandin. Meteo-France
!        ORIGINAL : J.-J. MORCRETTE         *ECMWF*

!     MODIFICATIONS.
!     --------------
!        M.Hamrud      01-Oct-2003 CY28 Cleaning
!        P. Marquet    14-Feb-2006 REFFWIA15 + NAMCLOP15 introduced

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

USE parkind1  ,ONLY : jpim     ,jprb
USE yomhook   ,ONLY : lhook,   dr_hook
 
! Ce qui concerne NULNAM commente par MPL le 15.04.09
!USE YOMLUN   , ONLY : NULNAM, NULOUT
USE yomlun   , ONLY : nulout
USE yomclop15, ONLY : ryfwca15 ,ryfwcb15 ,ryfwcc15 ,ryfwcd15 ,&
 & ryfwce15 ,ryfwcf15 ,rebcua15 ,rebcub15 ,rebcuc15 ,&
 & rebcud15 ,rebcue15 ,rebcuf15 ,rebcug15 ,rebcuh15 ,&
 & reffia15 ,reffib15 ,rtiw15   ,rriw15   ,reffwia15 

IMPLICIT NONE

#include "posnam.intfb.h"
#include "namclop15.h"


!*       1.    Set default values.
!              -------------------

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

! SW : 2 spectral intervals

REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook('SUCLOP15',0,zhook_handle)
rebcua15(1)= 3.448e-03_jprb
rebcua15(2)= 3.448e-03_jprb
rebcub15(1)= 2.431_jprb
rebcub15(2)= 2.431_jprb
rebcuc15(1)= 0.99999_jprb
rebcuc15(2)= 0.975634_jprb
rebcud15(1)= 0._jprb
rebcud15(2)= 2.487e-04_jprb
rebcue15(1)= 0.7661_jprb
rebcue15(2)= 0.7866_jprb
rebcuf15(1)= 5.851e-04_jprb
rebcuf15(2)= 5.937e-04_jprb

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

rebcug15= 1.07677_jprb
rebcuh15= 0.00267_jprb

! Ice particle Effective Radius as a function of LWC

reffia15= 40._jprb
reffib15=  0._jprb

! Water and Ice particle Effective Radius in the RADLSW15 formulae

reffwia15= 10._jprb

!* Water cloud properties - from Fouquart (1987)

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

ryfwca15(1)= 0._jprb
ryfwca15(2)= 0._jprb
ryfwcb15(1)= 1.5_jprb
ryfwcb15(2)= 1.5_jprb
ryfwcc15(1)= 0.9999_jprb
ryfwcc15(2)= 0.9988_jprb
ryfwcd15(1)= 5.000e-04_jprb
ryfwcd15(2)= 2.500e-03_jprb
ryfwce15(1)= 0.5_jprb
ryfwce15(2)= 0.05_jprb
ryfwcf15(1)= 0.865_jprb
ryfwcf15(2)= 0.910_jprb

!* Liquid/Solid water transition

rtiw15= 263._jprb
rriw15= 20._jprb


!*       2.    Modify default values.
!              ----------------------

! Ce qui concerne NAMCLOP15 commente par MPL le 15.04.09
!CALL POSNAM(NULNAM,'NAMCLOP15')
!READ       (NULNAM, NAMCLOP15)


!*       3.    Print final values.
!              -------------------

WRITE(unit=nulout,fmt='('' - - - - - - - - -'')')
WRITE(unit=nulout,fmt='('' COMMON YOMCLOP15 '')')
WRITE(unit=nulout,fmt='('' - - - - - - - - -'')')
WRITE(unit=nulout,fmt='( '' REFFWIA15 = '',E11.4 )') &
                          & reffwia15

IF (lhook) CALL dr_hook('SUCLOP15',1,zhook_handle)

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

END SUBROUTINE suclop15