surdi.F90

Go to the documentation of this file.

!
! $Id: surdi.F90 2027 2014-04-29 13:38:53Z fairhead $
!
SUBROUTINE surdi

!**** *SURDI*   - INITIALIZE COMMON YOERDI CONTROLLING RADINT

!     PURPOSE.
!     --------
!           INITIALIZE YOERDI, THE COMMON THAT CONTROLS THE
!           RADIATION INTERFACE

!**   INTERFACE.
!     ----------
!        CALL *SURDI* FROM *SURAD*
!              -----        -----

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

!        IMPLICIT ARGUMENTS :
!        --------------------
!        COMMON YOERDI

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

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

!     REFERENCE.
!     ----------
!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS MODEL

!     AUTHOR.
!     -------
!        Original  JEAN-JACQUES MORCRETTE  *ECMWF*
!        Modified   P. Viterbo   99-03-26    Tiling of the land surface
!        Modified   P. Viterbo   24-05-2004  surf library

!     MODIFICATIONS.
!     --------------
!        ORIGINAL : 88-12-15
!        M.Hamrud      01-Oct-2003 CY28 Cleaning
!        JJMorcrette   2004-10-07 Gas concentrations
!     ------------------------------------------------------------------

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

USE yoerdi   , ONLY : rrae     ,&
! & RCARDI   ,RCH4     ,RN2O     ,RO3      ,RCFC11   ,&
 & rcfc12, &   
 & repclc   ,reph2o   ,rsundur  ,&
 & rcco2    ,rcch4    ,rcn2o    ,rccfc11  ,rccfc12


IMPLICIT NONE

!#include "clesphys.h"
REAL(KIND=JPRB) :: ZAIRMWG, ZC11MWG, ZC12MWG, ZCH4MWG, ZCO2MWG, ZN2OMWG, ZO3MWG
REAL(KIND=JPRB) :: ZHOOK_HANDLE

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

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

IF (lhook) CALL dr_hook('SURDI',0,zhook_handle)
rrae = 0.1277e-02_jprb

!* Threshold for computing sunshine duration (W/m2)
rsundur=120._jprb

!*  For sea ice, monthly values are based on Ebert and Curry, 1993, Table 2.
!   We take dry snow albedo as the representative value for non-summer
!   months, and bare sea-ice as the representative value for summer
!   months. The values for Antarctic are shifted six-months.
! All computations brought back to *SUSWN*

!*  Concentration of the various trace gases (IPCC/SACC values for 1990)
!        CO2         CH4        N2O        CFC11       CFC12
!      353ppmv     1.72ppmv   310ppbv     280pptv     484pptv

zairmwg = 28.970_jprb
zco2mwg = 44.011_jprb
zch4mwg = 16.043_jprb
zn2omwg = 44.013_jprb
zo3mwg  = 47.9982_jprb
zc11mwg = 137.3686_jprb
zc12mwg = 120.9140_jprb

!RCARDI  = 353.E-06_JPRB*ZCO2MWG/ZAIRMWG
!RCH4    = 1.72E-06_JPRB*ZCH4MWG/ZAIRMWG
!RN2O    = 310.E-09_JPRB*ZN2OMWG/ZAIRMWG
!RO3     =   1.E-06_JPRB*ZO3MWG /ZAIRMWG
!RCFC11  = 280.E-12_JPRB*ZC11MWG/ZAIRMWG
!RCFC12  = 484.E-12_JPRB*ZC12MWG/ZAIRMWG

!RCARDI  = RCCO2   * ZCO2MWG/ZAIRMWG
!RCH4    = RCCH4   * ZCH4MWG/ZAIRMWG
!RN2O    = RCN2O   * ZN2OMWG/ZAIRMWG
!RO3     = 1.E-06_JPRB*ZO3MWG /ZAIRMWG
!RCFC11  = RCCFC11 * ZC11MWG/ZAIRMWG
!RCFC12  = RCCFC12 * ZC12MWG/ZAIRMWG


repclc=1.e-12_jprb
reph2o=1.e-12_jprb


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

IF (lhook) CALL dr_hook('SURDI',1,zhook_handle)
END SUBROUTINE surdi