!> \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 ! !========================================! !! @}