yoesw.F90

Go to the documentation of this file.

MODULE yoesw

USE parkind1  ,ONLY : jpim     ,jprb

IMPLICIT NONE

SAVE

!     ------------------------------------------------------------------
!*    ** *YOESW* - COEFFICIENTS FOR SHORTWAVE RADIATION TRANSFER
!     ------------------------------------------------------------------

REAL(KIND=JPRB) :: apad(6,3,7)
REAL(KIND=JPRB) :: bpad(6,3,7)
REAL(KIND=JPRB) :: rray(6,6)
REAL(KIND=JPRB), ALLOCATABLE :: rsun(:)
!$OMP THREADPRIVATE(RSUN)
REAL(KIND=JPRB) :: rpdh1
REAL(KIND=JPRB) :: rpdu1
REAL(KIND=JPRB) :: rpnh
REAL(KIND=JPRB) :: rpnu
REAL(KIND=JPRB) :: rswce(6)
REAL(KIND=JPRB) :: rswcp(6)
REAL(KIND=JPRB) :: rtdh2o
REAL(KIND=JPRB) :: rtdumg
REAL(KIND=JPRB) :: rth2o
REAL(KIND=JPRB) :: rtumg
REAL(KIND=JPRB) :: d(6,3)
REAL(KIND=JPRB) :: rexpo3(6,2,7)
INTEGER(KIND=JPIM) :: nexpo3(6)

REAL(KIND=JPRB) :: ryfwca(6)
REAL(KIND=JPRB) :: ryfwcb(6)
REAL(KIND=JPRB) :: ryfwcc(6)
REAL(KIND=JPRB) :: ryfwcd(6)
REAL(KIND=JPRB) :: ryfwce(6)
REAL(KIND=JPRB) :: ryfwcf(6)

REAL(KIND=JPRB) :: rebcua(6)
REAL(KIND=JPRB) :: rebcub(6)
REAL(KIND=JPRB) :: rebcuc(6)
REAL(KIND=JPRB) :: rebcud(6)
REAL(KIND=JPRB) :: rebcue(6)
REAL(KIND=JPRB) :: rebcuf(6)
REAL(KIND=JPRB) :: rebcug(16)
REAL(KIND=JPRB) :: rebcuh(16)
REAL(KIND=JPRB) :: rebcui(6)
REAL(KIND=JPRB) :: rebcuj(6)

REAL(KIND=JPRB) :: raswca(6)
REAL(KIND=JPRB) :: raswcb(6)
REAL(KIND=JPRB) :: raswcc(6)
REAL(KIND=JPRB) :: raswcd(6)
REAL(KIND=JPRB) :: raswce(6)
REAL(KIND=JPRB) :: raswcf(6)

REAL(KIND=JPRB) :: rfueta(16,3),rfuetb(16,4), rfuetc(16,4)
REAL(KIND=JPRB) :: rfulio(16,3)
REAL(KIND=JPRB) :: rhsavi(16,3)
REAL(KIND=JPRB) :: rlilia(16,5),rlilib(16,4)

REAL(KIND=JPRB) :: rflaa0(6)
REAL(KIND=JPRB) :: rflaa1(6)
REAL(KIND=JPRB) :: rflbb0(6)
REAL(KIND=JPRB) :: rflbb1(6)
REAL(KIND=JPRB) :: rflbb2(6)
REAL(KIND=JPRB) :: rflbb3(6)
REAL(KIND=JPRB) :: rflcc0(6)
REAL(KIND=JPRB) :: rflcc1(6)
REAL(KIND=JPRB) :: rflcc2(6)
REAL(KIND=JPRB) :: rflcc3(6)

REAL(KIND=JPRB) :: rfuaa0(6)
REAL(KIND=JPRB) :: rfuaa1(6)
REAL(KIND=JPRB) :: rfubb0(6)
REAL(KIND=JPRB) :: rfubb1(6)
REAL(KIND=JPRB) :: rfubb2(6)
REAL(KIND=JPRB) :: rfubb3(6)
REAL(KIND=JPRB) :: rfucc0(6)
REAL(KIND=JPRB) :: rfucc1(6)
REAL(KIND=JPRB) :: rfucc2(6)
REAL(KIND=JPRB) :: rfucc3(6)
REAL(KIND=JPRB) :: rfldd0(6)
REAL(KIND=JPRB) :: rfldd1(6)
REAL(KIND=JPRB) :: rfldd2(6)
REAL(KIND=JPRB) :: rfldd3(6)

REAL(KIND=JPRB) :: rsushe(6)
REAL(KIND=JPRB) :: rsushf(6)
REAL(KIND=JPRB) :: rsushh(6)
REAL(KIND=JPRB) :: rsushk(6)
REAL(KIND=JPRB) :: rsusha(6)
REAL(KIND=JPRB) :: rsushg(6)
REAL(KIND=JPRB) :: rsushfa(4)
REAL(KIND=JPRB) :: rsushc
REAL(KIND=JPRB) :: rsushd

REAL(KIND=JPRB) :: reffia
REAL(KIND=JPRB) :: reffib
REAL(KIND=JPRB) :: rtiw
REAL(KIND=JPRB) :: rriw
REAL(KIND=JPRB) :: rroma(6)
REAL(KIND=JPRB) :: rromb(6)
REAL(KIND=JPRB) :: rrasy(6)

REAL(KIND=JPRB) :: rhsra(6)
REAL(KIND=JPRB) :: rhsrb(6)
REAL(KIND=JPRB) :: rhsrc(6)
REAL(KIND=JPRB) :: rhsrd(6)
REAL(KIND=JPRB) :: rhsre(6)
REAL(KIND=JPRB) :: rhsrf(6)
REAL(KIND=JPRB) :: rhsrta
REAL(KIND=JPRB) :: rhsrtb

REAL(KIND=JPRB) :: rtaua(6,6)
REAL(KIND=JPRB) :: rpiza(6,6)
REAL(KIND=JPRB) :: rcga(6,6)
REAL(KIND=JPRB) :: raer(6,6)

REAL(KIND=JPRB) :: rtweight(8)
REAL(KIND=JPRB) :: rweigs(6)
REAL(KIND=JPRB) :: rweigv(6)

REAL(KIND=JPRB) :: rweight(6,8)
INTEGER(KIND=JPIM) :: nmpsrtm(14), ntyps

REAL(KIND=JPRB) :: radjust

!        * E.C.M.W.F. PHYSICS PACKAGE *

!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14

!  NAME     TYPE     PURPOSE
!  ----  :  ----   : ---------------------------------------------------
!  APAD  :  REAL     PADE APPROXIMANTS NUMERATOR
!  BPAD  :  REAL     PADE APPROXIMANTS DENOMINATOR
!  D     :  REAL     TRANSMISSION LIMIT FOR INFINITE ABSORBER AMOUNT
!  RRAY  :  REAL     RAYLEIGH SCATTERING COEFFICIENTS
!  RSUN  :  REAL     SOLAR FRACTION IN SPECTRAL INTERVALS
!  RPDH1 :  1 + EXPONENT PRESSURE DEPENDENCE H2O
!  RPDU1 :  1 + EXPONENT PRESSURE DEPENDENCE UNIFORMLY MIXED GASES
!  RPNH  :  REFERENCE PRESSURE FACTOR FOR H2O
!  RPNU  :  REFERENCE PRESSURE FACTOR FOR UNIFORMLY MIXED GASES
!  RSWCE :  E-TYPE, H2O CONTINUUM ABSORPTION COEFFICIENT 
!  RSWCP :  P-TYPE, H2O CONTINUUM ABSORPTION COEFFICIENT 
!  RTDH2O:  EXPONENT TEMPERATURE DEPENDENCE H2O
!  RTDUMG:  EXPONENT TEMPERATURE DEPENDENCE UNIFORMLY MIXED GASES
!  RTH2O :  REFERENCE TEMPERATURE H2O
!  RTUMG :  REFERENCE TEMPERATURE UNIFORMLY MIXED GASES
!     -----------------------------------------------------------------

!        * E.C.M.W.F. PHYSICS PACKAGE *

!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14

!  NAME     TYPE     PURPOSE
!  ----  :  ----   : ---------------------------------------------------
!*    FOUQUART (1987) WATER CLOUD OPTICAL PROPERTIES

! RYFWCA :  REAL   : C1 IN OPTICAL THICKNESS FORMULA
! RYFWCB :  REAL   : C2 IN OPTICAL THICKNESS FORMULA
! RYFWCC :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
! RYFWCD :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
! RYFWCE :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
! RYFWCF :  REAL   : ASSYMETRY FACTOR

!*    SLINGO (1989) WATER CLOUD OPTICAL PROPERTIES

! RASWCA :  REAL   : C1 IN OPTICAL THICKNESS FORMULA
! RASWCB :  REAL   : C2 IN OPTICAL THICKNESS FORMULA
! RASWCC :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
! RASWCD :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
! RASWCE :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
! RASWCF :  REAL   : ASSYMETRY FACTOR

!*   SAVIJARVI (1998) WATER CLOUD OPTICAL PROPERTIES (RRTM)

! RHSAVI : REAL    : MASS ABSORPTION COEFFICIENTS (POLYNOMIAL DEVELOPM)

!*   LINDNER,LI (2000) WATER CLOUD OPTICAL PROPERTIES (RRTM)

! RLILIA : REAL    : MASS ABSORPTION COEFFICIENTS (POLYNOMIAL DEVELOPM)
! RLILIB : REAL    : 1-SSA COEFFICIENTS  (POLYNOMIAL DEVELOPM)

!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM EBERT-CURRY (1992)

! REBCUA :  REAL   : C1 IN OPTICAL THICKNESS FORMULA
! REBCUB :  REAL   : C2 IN OPTICAL THICKNESS FORMULA
! REBCUC :  REAL   : 1-C3  IN SINGLE SCATTERING ALBEDO FORMULA
! REBCUD :  REAL   : C4 IN SINGLE SCATTERING ALBEDO FORMULA
! REBCUE :  REAL   : C5 IN ASSYMETRY FACTOR FORMULA
! REBCUF :  REAL   : C6 IN ASSYMETRY FACTOR FORMULA
! REBCUG :  REAL   : C7 IN MASS ABSORPTION COEFFICIENT FORMULA
! REBCUH :  REAL   : C8 IN MASS ABSORPTION COEFFICIENT FORMULA
! REBCUI :  REAL   : C7 IN MASS ABSORPTION COEFFICIENT SPECTRAL FORMULA
! REBCUJ :  REAL   : C8 IN MASS ABSORPTION COEFFICIENT SPECTRAL FORMULA

!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM SUN-SHINE (1995)

! RSHSUE :  REAL   : E IN SINGLE SCATTERING ALBEDO FORMULA
! RSHSUF :  REAL   : F IN SINGLE SCATTERING ALBEDO FORMULA
! RSHSUH :  REAL   : H IN ASSYMETRY FACTOR FORMULA
! RSHSUK :  REAL   : K IN ASSYMETRY FACTOR FORMULA
! RSHSUA :  REAL   : ALPHA IN SSA CORRECTION FACTOR FORMULA
! RSHSUG :  REAL   : GAMMA IN ASSYMETRY CORRECTION FACTOR FORMULA
! RSHSUFA:  REAL   : COEFFICIENTS IN TEMPERATURE CORRECTION FACTOR

! REFFIA :  REAL   : C9  IN EFFECTIVE RADIUS FORMULA
! REFFIB :  REAL   : C10 IN EFFECTIVE RADIUS FORMULA

!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM FU-LIOU (1993)

! RFULIO :  REAL   : COEFFICIENTS IN EXPRESSION FOR LW EXTINCTION COEFF.
! RFLAA  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW EXTINCTION COEFF.
! RFLBB  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW SINGLE SCATT.ALB.
! RFLCC  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR
! RFLDD  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR

!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM FU (1996) & FU ET AL. (1998)

! RFUETA :  REAL   : COEFFICIENTS IN EXPRESSION FOR LW EXTINCTION COEFF.
! RFUAA  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW EXTINCTION COEFF.
! RFUBB  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW SINGLE SCATT.ALB.
! RFUCC  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR

!*    TRANSITION BETWEEN LIQUID AND SOLID WATER

! RTIW   :  REAL   : TEMPERATURE THRESHOLD
! RRIW   :  REAL   : TRANSITION RANGE

!*    RAIN OPTICAL PROPERTIES FROM SAVIJARVI (1996)

! RROMA  :  REAL   : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO
! RROMB  :  REAL   : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO
! RRASY  :  REAL   : COEFFICIENTS FOR ASSYMETRY FACTOR
! RHSRA  :  REAL   : COEFFICIENTS FOR OPTICAL THICKNESS
! RHSRB  :  REAL   : COEFFICIENTS FOR OPTICAL THICKNESS
! RHSRC  :  REAL   : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO
! RHSRD  :  REAL   : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO
! RHSRE  :  REAL   : COEFFICIENTS FOR ASSYMETRY FACTOR 
! RHSRF  :  REAL   : COEFFICIENTS FOR ASSYMETRY FACTOR
! RHSRTA :  REAL   : COEFFICIENTS FOR OPTICAL THICKNESS
! RHSRTB :  REAL   : COEFFICIENTS FOR OPTICAL THICKNESS
!     -----------------------------------------------------------------

!        * E.C.M.W.F. PHYSICS PACKAGE *

!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14

!  NAME     TYPE     PURPOSE
!  ----  :  ----   : -------
!  RTAUA :  REAL     S.W. NORMALIZED OPTICAL THICKNESS AT 0.55 MICRON
!  RPIZA :  REAL     S.W. SINGLE SCATTERING ALBEDO
!  RCGA  :  REAL     S.W. ASSYMETRY FACTOR
!  RAER  :  REAL     L.W. ABSORPTION COEFFICIENTS
!     -----------------------------------------------------------------

!        * E.C.M.W.F. PHYSICS PACKAGE *

!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14

!  NAME     TYPE     PURPOSE
!  ----  :  ----   : -------
!RTWEIGHT:  REAL     S.W. INTEGRATED WEIGHT 
! NMPSRTM: INTEGER  : Indices for mapping SW[1:6] albedo into SRTM[1:14]  
!     -----------------------------------------------------------------
!$OMP THREADPRIVATE(apad,bpad,d,nexpo3,nmpsrtm,ntyps,radjust,raer,raswca)
!$OMP THREADPRIVATE(raswcb,raswcc,raswcd,raswce,raswcf,rcga,rebcua,rebcub)
!$OMP THREADPRIVATE(rebcuc,rebcud,rebcue,rebcuf,rebcug,rebcuh,rebcui,rebcuj)
!$OMP THREADPRIVATE(reffia,reffib,rexpo3,rflaa0,rflaa1,rflbb0,rflbb1,rflbb2)
!$OMP THREADPRIVATE(rflbb3,rflcc0,rflcc1,rflcc2,rflcc3,rfldd0,rfldd1,rfldd2)
!$OMP THREADPRIVATE(rfldd3,rfuaa0,rfuaa1,rfubb0,rfubb1,rfubb2,rfubb3,rfucc0)
!$OMP THREADPRIVATE(rfucc1,rfucc2,rfucc3,rfueta,rfuetb,rfuetc,rfulio,rhsavi)
!$OMP THREADPRIVATE(rhsra,rhsrb,rhsrc,rhsrd,rhsre,rhsrf,rhsrta,rhsrtb,rlilia)
!$OMP THREADPRIVATE(rlilib,rpdh1,rpdu1,rpiza,rpnh,rpnu,rrasy,rray,rriw,rroma)
!$OMP THREADPRIVATE(rromb,rsusha,rsushc,rsushd,rsushe,rsushf,rsushfa,rsushg)
!$OMP THREADPRIVATE(rsushh,rsushk,rswce,rswcp,rtaua,rtdh2o,rtdumg,rth2o,rtiw)
!$OMP THREADPRIVATE(rtumg,rtweight,rweight,rweigs,rweigv,ryfwca,ryfwcb,ryfwcc,ryfwcd,ryfwce,ryfwcf)
END MODULE yoesw