surdi.F90

Go to the documentation of this file.

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