!> \file radsw_param.f
!! This file contains SW band parameters setup.
!> \ingroup module_radsw_main
!! @{
!!!!! ============================================================== !!!!!
!!!!! sw-rrtm3 radiation package description !!!!!
!!!!! ============================================================== !!!!!
! !
! this package includes ncep's modifications of the rrtm-sw radiation !
! code from aer inc. !
! !
! the sw-rrtm3 package includes these parts: !
! !
! 'radsw_rrtm3_param.f' !
! 'radsw_rrtm3_datatb.f' !
! 'radsw_rrtm3_main.f' !
! !
! the 'radsw_rrtm3_param.f' contains: !
! !
! 'module_radsw_parameters' -- band parameters set up !
! !
! the 'radsw_rrtm3_datatb.f' contains: !
! !
! 'module_radsw_ref' -- reference temperature and pressure !
! 'module_radsw_cldprtb' -- cloud property coefficients table !
! 'module_radsw_sflux' -- spectral distribution of solar flux !
! 'module_radsw_kgbnn' -- absorption coeffients for 14 !
! bands, where nn = 16-29 !
! !
! the 'radsw_rrtm3_main.f' contains: !
! !
! 'module_radsw_main' -- main sw radiation transfer !
! !
! in the main module 'module_radsw_main' there are only two !
! externally callable subroutines: !
! !
! 'swrad' -- main rrtm3 sw radiation routine !
! 'rswinit' -- initialization routine !
! !
! all the sw radiation subprograms become contained subprograms !
! in module 'module_radsw_main' and many of them are not directly !
! accessable from places outside the module. !
! !
! compilation sequence is: !
! !
! 'radsw_rrtm3_param.f' !
! 'radsw_rrtm3_datatb.f' !
! 'radsw_rrtm3_main.f' !
! !
! and all should be put in front of routines that use sw modules !
! !
! ncep modifications history log: !
! !
! see list in program "radsw_rrtm3_main.f" !
! !
!!!!! ============================================================== !!!!!
!!!!! end descriptions !!!!!
!!!!! ============================================================== !!!!!
!> This module is for specifying the band structures and program parameters
!! used by the RRTMG-SW scheme.
!========================================!
module module_radsw_parameters !
!........................................!
use physparam, only : kind_phys
implicit none
!
public
!
!> derived type for SW fluxes at TOA
type :: topfsw_type
!> total-sky upward flux
real (kind=kind_phys) :: upfxc
!> total-sky downward flux
real (kind=kind_phys) :: dnfxc
!> clear-sky upward flux
real (kind=kind_phys) :: upfx0
end type
!
!> derived type for SW fluxes at surface
type :: sfcfsw_type
!> total-sky upward flux
real (kind=kind_phys) :: upfxc
!> total-sky downward flux
real (kind=kind_phys) :: dnfxc
!> clear-sky upward flux
real (kind=kind_phys) :: upfx0
!> clear-sky downward flux
real (kind=kind_phys) :: dnfx0
end type
!
!> derived type for SW fluxes' column profiles (at layer interfaces)
type :: profsw_type
!> total-sky upward flux
real (kind=kind_phys) :: upfxc
!> total-sky downward flux
real (kind=kind_phys) :: dnfxc
!> clear-sky upward flux
real (kind=kind_phys) :: upfx0
!> clear-sky downward flux
real (kind=kind_phys) :: dnfx0
end type
!
!> derived type for special components of surface SW fluxes
type :: cmpfsw_type
!> total-sky downward flux cover UV-B spectrum
real (kind=kind_phys) :: uvbfc
!> clear-sky downward flux cover UV-B spectrum
real (kind=kind_phys) :: uvbf0
!> total-sky downward flux for near-IR direct beam
real (kind=kind_phys) :: nirbm
!> total-sky downward flux for near-IR diffused part
real (kind=kind_phys) :: nirdf
!> total-sky downward flux for UV+Visible direct
real (kind=kind_phys) :: visbm
!> total-sky downward flux for UV+Visible diffused
real (kind=kind_phys) :: visdf
end type
!
!! \name Parameter constants for SW band structures
!> band range lower index
integer, parameter :: NBLOW = 16
!> band range upper index
integer, parameter :: NBHGH = 29
!> total number of SW bands (14)
integer, parameter :: NBANDS = NBHGH-NBLOW+1
!> total number of g-point in all bands
integer, parameter :: NGPTSW = 112
!> maximum number of g-point in one band
integer, parameter :: NGMAX = 16
!> maximum number of absorbing gases
integer, parameter :: MAXGAS = 7
!> index upper limit of optical depth and transmittance tables
integer, parameter :: NTBMX = 10000
!> SW bands counter starting index (for compatibility with previous
!! SW radiation schemes)
integer, parameter :: NSWSTR = 1
! integer, parameter :: NSWEND = NBANDS
integer, parameter :: NBDSW = NBANDS
!> \name The actual number of g-point for bands 16-29
integer :: NG16, NG17, NG18, NG19, NG20, NG21, NG22,
& NG23, NG24, NG25, NG26, NG27, NG28, NG29
parameter ( NG16=06, NG17=12, NG18=08, NG19=08, NG20=10,
& NG21=10, NG22=02, NG23=10, NG24=08, NG25=06,
& NG26=06, NG27=08, NG28=06, NG29=12)
integer, dimension(NBLOW:NBHGH) :: NG
data NG / NG16, NG17, NG18, NG19, NG20, NG21, NG22,
& NG23, NG24, NG25, NG26, NG27, NG28, NG29 /
!> \name Accumulative starting index for bands 16-29
integer :: NS16, NS17, NS18, NS19, NS20, NS21, NS22,
& NS23, NS24, NS25, NS26, NS27, NS28, NS29
parameter ( NS16=00, NS17=NS16+NG16, NS18=NS17+NG17,
& NS19=NS18+NG18, NS20=NS19+NG19, NS21=NS20+NG20,
& NS22=NS21+NG21, NS23=NS22+NG22, NS24=NS23+NG23,
& NS25=NS24+NG24, NS26=NS25+NG25, NS27=NS26+NG26,
& NS28=NS27+NG27, NS29=NS28+NG28 )
!> array contains values of NS16-NS29
integer, dimension(NBLOW:NBHGH) :: NGS
data NGS / NS16, NS17, NS18, NS19, NS20, NS21, NS22, &
& NS23, NS24, NS25, NS26, NS27, NS28, NS29 /
!> reverse checking of band index for each g-point
integer, dimension(NGPTSW) :: NGB
data NGB(:) / 16,16,16,16,16,16, & ! band 16
& 17,17,17,17,17,17,17,17,17,17,17,17, & ! band 17
& 18,18,18,18,18,18,18,18, & ! band 18
& 19,19,19,19,19,19,19,19, & ! band 19
& 20,20,20,20,20,20,20,20,20,20, & ! band 20
& 21,21,21,21,21,21,21,21,21,21, & ! band 21
& 22,22, & ! band 22
& 23,23,23,23,23,23,23,23,23,23, & ! band 23
& 24,24,24,24,24,24,24,24, & ! band 24
& 25,25,25,25,25,25, & ! band 25
& 26,26,26,26,26,26, & ! band 26
& 27,27,27,27,27,27,27,27, & ! band 27
& 28,28,28,28,28,28, & ! band 28
& 29,29,29,29,29,29,29,29,29,29,29,29 / ! band 29
!> \name Starting/ending wavenumber for each of the SW bands
real (kind=kind_phys), dimension(NBANDS):: wvnum1, wvnum2
data wvnum1(:) / &
& 2600.0, 3251.0, 4001.0, 4651.0, 5151.0, 6151.0, 7701.0, &
& 8051.0,12851.0,16001.0,22651.0,29001.0,38001.0, 820.0 /
data wvnum2(:) / &
& 3250.0, 4000.0, 4650.0, 5150.0, 6150.0, 7700.0, 8050.0, &
& 12850.0,16000.0,22650.0,29000.0,38000.0,50000.0, 2600.0 /
!
!........................................!
end module module_radsw_parameters !
!========================================!
!! @}