susw.F90

Go to the documentation of this file.

SUBROUTINE susw

!**** *SUSW*   - INITIALIZE COMMON YOESW

!     PURPOSE.
!     --------
!           INITIALIZE YOESW, THE COMMON THAT CONTAINS COEFFICIENTS
!           NEEDED TO RUN THE SHORTWAVE RADIATION SUBROUTINES

!**   INTERFACE.
!     ----------
!        *CALL* *SUSW

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

!        IMPLICIT ARGUMENTS :
!        --------------------
!        COMMON YOESW

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

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

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

!     AUTHOR.
!     -------
!        JEAN-JACQUES MORCRETTE *ECMWF*

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

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

#include "tsmbkind.h"

USE yomcst   , ONLY : rg
USE yoesw    , ONLY : apad     ,bpad     ,rray     ,rsun     ,&
            &rpdh1    ,rpdu1    ,rpnh     ,rpnu     ,rswce    ,&
            &rswcp    ,rtdh2o   ,rtdumg   ,rth2o    ,rtumg    ,d

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

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



IMPLICIT NONE
real_b   :: zpdh2o=0.8_jprb
real_b   :: zpdumg=0.75_jprb
real_b   :: zprh2o=30000._jprb
real_b   :: zprumg=30000._jprb

!     LOCAL REAL SCALARS
real_b :: zh2o, zumg

rpdh1=zpdh2o+_one_
rpdu1=zpdumg+_one_
zh2o=_one_/( 10._jprb* rg * rpdh1 )
zumg=_one_/( 10._jprb* rg * rpdu1 )
rpnu = zumg/(zprumg**zpdumg)
rpnh = zh2o/(zprh2o**zpdh2o)

!      RSWCP=0.002*RSWCE


rtdh2o = 0.40_jprb
rtdumg = 0.375_jprb
rth2o = 240._jprb
rtumg = 240._jprb
rswce = 0._jprb
rswcp = 0._jprb

rsun(1) = 0.441676_jprb
d(1,:) = (/ 0.00_jprb, 0.00_jprb, 0.00_jprb /)
!* DERIVED FROM HITRAN APRIL 1991
!       H2O:  Pref=300 hPa, Tref=240K, Pdep=0.8
!       O3 :  unchanged

apad(1, 1, :) = (/&
 &0.912418292e+05_jprb,&
 &0.723613782e+05_jprb,&
 &0.596037057e+04_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb/)
apad(1, 2, :) = (/&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb/)
apad(1, 3, :) = (/&
 &0.925887084e-04_jprb,&
 &0.129353723e-01_jprb,&
 &0.800821928e+00_jprb,&
 &0.242715973e+02_jprb,&
 &0.878331486e+02_jprb,&
 &0.191559725e+02_jprb,&
 &0.000000000e+00_jprb/)

bpad(1, 1, :) = (/&
 &0.912418292e+05_jprb,&
 &0.724555318e+05_jprb,&
 &0.602593328e+04_jprb,&
 &0.100000000e+01_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb/)
bpad(1, 2, :) = (/&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb/)
bpad(1, 3, :) = (/&
 &0.925887084e-04_jprb,&
 &0.131812683e-01_jprb,&
 &0.812706117e+00_jprb,&
 &0.249863591e+02_jprb,&
 &0.931071925e+02_jprb,&
 &0.252233437e+02_jprb,&
 &0.100000000e+01_jprb/)

rray(1,1:6) = (/&
 &.428937e-01_jprb, .890743e+00_jprb,-.288555e+01_jprb,&
 &.522744e+01_jprb,-.469173e+01_jprb, .161645e+01_jprb/)

rsun(2) = 0.558324_jprb

d(2,1:3) = (/ 0.000000000_jprb, 0.000000000_jprb, 0.800000000_jprb /)

!* INTERVAL 2:  0.68 - 4.00 MICRONS
!* DERIVED FROM HITRAN APRIL 1991
!       H2O:  Pref=300 hPa, Tref=240K, Pdep=0.80
!       UMG:  Pref=300 hPa, Tref=240K, Pdep=0.75 (CO2+O2+CH4+N2O+CO)
!       O3 :  unchanged

apad(2, 1, :) = (/&
 &0.376655383e-08_jprb,&
 &0.978576773e-04_jprb,&
 &0.387714006e+00_jprb,&
 &0.118461660e+03_jprb,&
 &0.119079797e+04_jprb,&
 &0.293353397e+03_jprb,&
 &0.000000000e+00_jprb/)
apad(2, 2, :) = (/&
 &0.739646016e-08_jprb,&
 &0.131849595e-03_jprb,&
 &0.437772681e+00_jprb,&
 &0.151345118e+03_jprb,&
 &0.233628890e+04_jprb,&
 &0.797219934e+03_jprb,&
 &0.000000000e+00_jprb/)
apad(2, 3, :) = (/&
 &0.410177786e+03_jprb,&
 &0.672595424e+02_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e-00_jprb,&
 &0.000000000e+00_jprb/)

bpad(2, 1, :) = (/&
 &0.376655383e-08_jprb,&
 &0.979023421e-04_jprb,&
 &0.388611139e+00_jprb,&
 &0.120291383e+03_jprb,&
 &0.130531005e+04_jprb,&
 &0.415049409e+03_jprb,&
 &0.100000000e+01_jprb/)
bpad(2, 2, :) = (/&
 &0.739646016e-08_jprb,&
 &0.131861712e-03_jprb,&
 &0.437949001e+00_jprb,&
 &0.151692730e+03_jprb,&
 &0.237071130e+04_jprb,&
 &0.867914360e+03_jprb,&
 &0.100000000e+01_jprb/)
bpad(2, 3, :) = (/&
 &0.410177786e+03_jprb,&
 &0.731185438e+02_jprb,&
 &0.100000000e+01_jprb,&
 &0.000000000e+00_jprb,&
 &0.000000000e+00_jprb,&
 &0.000000000e+00_jprb,&
 &0.000000000e+00_jprb/)

rray(2,1:6) = (/&
 &.697200e-02_jprb, .173297e-01_jprb,-.850903e-01_jprb,&
 &.248261e+00_jprb,-.302031e+00_jprb, .129662e+00_jprb/)
RETURN
END SUBROUTINE susw