surdi15.F90

Go to the documentation of this file.

!OPTIONS XOPT(NOEVAL)
SUBROUTINE surdi15

!**** *SURDI15*   - INITIALIZE COMMON YOMRDI15 CONTROLLING RADINT
!****               FROZEN VERSION (CYCLE 15) OF SURDI

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

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

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

!        IMPLICIT ARGUMENTS :
!        --------------------
!        COMMON YOMRDI15

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

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

!     REFERENCE.
!     ----------
!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE
!     "IN CORE MODEL"

!     AUTHOR.
!     -------
!        96-11: Ph. Dandin. Meteo-France
!        ORIGINAL : 88-12-15 BY JEAN-JACQUES MORCRETTE  *ECMWF*

!     MODIFICATIONS.
!     --------------
!        M.Hamrud      01-Oct-2003 CY28 Cleaning
!        A.Alias       05-12-2005 greenhouse gases variables (M.Deque)
!        A.Alias       13-06-2006 RI0 value can be changed via namscen.h

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

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
USE yomlun_ifsaux   , ONLY : nulout

USE yomcst          , ONLY : ri0
USE yomrdi15        , ONLY : rsdtsn15 ,rrae15   ,rmu0015  ,ralbice15,&
 & ralbsea15,ralbsnm15,ralbsno15,rcardi15 ,remiss15 ,&
 & rsnowal15,rvlbdc15 ,rch415   ,rn2o15   ,rcfc1115 ,&
 & rcfc1215 ,repalb15 ,repclc15 ,reph2o15  
USE yomrdu15        , ONLY : repsec15

IMPLICIT NONE

REAL(KIND=JPRB) :: XCARDI, XCFC11, XCFC12, XCH4, XN2O
REAL(KIND=JPRB) :: ZAIRMWG, ZC11MWG, ZC12MWG, ZCH4MWG, ZCO2MWG, ZN2OMWG, ZSUPSAT
REAL(KIND=JPRB) :: ZHOOK_HANDLE

#include "namscen.h"
#include "posnam.intfb.h"
!      ----------------------------------------------------------------

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

IF (lhook) CALL dr_hook('SURDI15',0,zhook_handle)
rsdtsn15  = 5.0_jprb

rrae15 = 0.1277e-02_jprb
rmu0015 = rrae15/sqrt(rrae15*(rrae15+2.0_jprb))

ralbice15 = 0.55_jprb
ralbsea15 = 0.07_jprb
ralbsno15 = 0.80_jprb
ralbsnm15 = 0.40_jprb
rsnowal15 = 0.01_jprb
!*  Concentration of the various trace gases (IPCC/SACC values for 1990)
!        CO2         CH4        N2O        CFC11       CFC12
!      353ppmv     1.72ppmv   310ppbv     280pptv     484pptv

xcardi  = 353.e-06_jprb
xch4    = 1.72e-06_jprb
xn2o    = 310.e-09_jprb
xcfc11  = 280.e-12_jprb
xcfc12  = 484.e-12_jprb

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

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

WRITE(nulout,'( &
 & '' CO2   = '',E14.7,'' CH4   = '',E14.7 &
 & ,'' N2O   = '',E14.7,'' CFC11 = '',E14.7 &
 & ,'' CFC12 = '',E14.7,'' RI0 = '',E14.7 &
 & )') xcardi,xch4,xn2o,xcfc11,xcfc12,ri0  

rcardi15  = xcardi*zco2mwg/zairmwg
rch415    = xch4*zch4mwg/zairmwg
rn2o15    = xn2o*zn2omwg/zairmwg
rcfc1115  = xcfc11*zc11mwg/zairmwg
rcfc1215  = xcfc12*zc12mwg/zairmwg
remiss15  = 0.996_jprb
!ZSUPSAT = 0.01_JPRB
rvlbdc15  = 0.5_jprb

repsec15=1.e-12_jprb
repclc15=1.e-12_jprb
reph2o15=1.e-12_jprb
repalb15=1.e-12_jprb

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

IF (lhook) CALL dr_hook('SURDI15',1,zhook_handle)
END SUBROUTINE surdi15