source: LMDZ5/branches/IPSLCM6.0.11/libf/phymar/surdi.F90

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

!     MODIFICATIONS.
!     --------------
!        ORIGINAL : 88-12-15
!     ------------------------------------------------------------------

#include "tsmbkind.h"

USE YOERDI   , ONLY : RRAE     ,RALBSEAD ,&
            &RALBICEVS_AR,RALBICENI_AR,RALBICEVS_AN,RALBICENI_AN,&
            &RALBSFO  ,REMISD   ,REMISL   ,REMISN   ,REMISS   ,&
            &RCARDI   ,RCH4     ,RN2O     ,RO3      ,RCFC11   ,&
            &RCFC12   ,REPALB   ,REPCLC   ,REPH2O   ,RSUNDUR

IMPLICIT NONE


!     LOCAL REAL SCALARS
REAL_B :: ZAIRMWG, ZC11MWG, ZC12MWG, ZCH4MWG, ZCO2MWG, ZN2OMWG, ZO3MWG
INTEGER_M :: JM, IM

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

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

RRAE = 0.1277E-02_JPRB

!* Threshold for computing sunshine duration (W/m2)
RSUNDUR=120._JPRB

!*  Ocean surface albedo for diffuse radiation (Taylor et al., 1997)      
RALBSEAD = 0.06_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.
!*  Sea ice surf. albedo for visible rad. (snow covered; Ebert and Curry, 1993)
RALBICEVS_AR(1:12) = (/0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB,&
                      &0.975_JPRB,0.876_JPRB,0.778_JPRB,0.778_JPRB,&
                      &0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB/)
!*  Sea ice surf. albedo for near-IR rad. (snow covered; Ebert and Curry, 1993)
RALBICENI_AR(1:12) = (/0.664_JPRB,0.664_JPRB,0.664_JPRB,0.664_JPRB,&
                      &0.664_JPRB,0.476_JPRB,0.288_JPRB,0.288_JPRB,&
                      &0.664_JPRB,0.664_JPRB,0.664_JPRB,0.664_JPRB/)
DO JM=1,12
  IM=MOD(JM+5,12)+1
  RALBICEVS_AN(JM)=RALBICEVS_AR(IM)
  RALBICENI_AN(JM)=RALBICENI_AR(IM)
ENDDO
!*  Snow albedo in the presence of high vegetation
RALBSFO=0.15_JPRB

!- sea surface emissivity and other surfaces outside the window region     
REMISS  = 0.99_JPRB
!- snow window emissivity      
REMISN  = 0.98_JPRB
!- land window emissivity      
REMISL  = 0.96_JPRB
!- desert window emissivity (lower bound when dry)      
REMISD  = 0.93_JPRB

!*  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

REPCLC=1.E-12_JPRB
REPH2O=1.E-12_JPRB
REPALB=1.E-12_JPRB

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

RETURN
END SUBROUTINE SURDI