scale_atmos_phy_rd.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_atmos_phy_rd
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  use scale_grid_index
  use scale_tracer
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: atmos_phy_rd_setup

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  abstract interface
     subroutine rd( &
          DENS, RHOT, QTRC,      &
          CZ, FZ,                &
          fact_ocean,            &
          fact_land,             &
          fact_urban,            &
          temp_sfc, albedo_land, &
          solins, cosSZA,        &
          flux_rad,              &
          flux_rad_top,          &
          flux_sfc_dn            )
       use scale_precision
       use scale_grid_index
       use scale_tracer
       implicit none

       real(RP), intent(in)  :: dens        (ka,ia,ja)
       real(RP), intent(in)  :: rhot        (ka,ia,ja)
       real(RP), intent(in)  :: qtrc        (ka,ia,ja,qa)
       real(RP), intent(in)  :: cz          (  ka,ia,ja)    ! UNUSED
       real(RP), intent(in)  :: fz          (0:ka,ia,ja)
       real(RP), intent(in)  :: fact_ocean  (ia,ja)
       real(RP), intent(in)  :: fact_land   (ia,ja)
       real(RP), intent(in)  :: fact_urban  (ia,ja)
       real(RP), intent(in)  :: temp_sfc    (ia,ja)
       real(RP), intent(in)  :: albedo_land (ia,ja,2)
       real(RP), intent(in)  :: solins      (ia,ja)
       real(RP), intent(in)  :: cossza      (ia,ja)
       real(RP), intent(out) :: flux_rad    (ka,ia,ja,2,2,2)
       real(RP), intent(out) :: flux_rad_top(ia,ja,2,2,2)
       real(RP), intent(out) :: flux_sfc_dn (ia,ja,2,2)
     end subroutine rd
  end interface
  procedure(rd), pointer :: atmos_phy_rd => null()
  public :: atmos_phy_rd

  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_rd_setup( RD_TYPE )
    use scale_process, only: &
       prc_mpistop
#define EXTM(pre, name, post) pre ## name ## post
#define NAME(pre, name, post) EXTM(pre, name, post)
#ifdef RD
    use name(scale_atmos_phy_rd_, rd,), only: &
       name(atmos_phy_rd_, rd, _setup), &
       name(atmos_phy_rd_, rd,)
#else
    use scale_atmos_phy_rd_mstrnx, only: &
       atmos_phy_rd_mstrnx_setup, &
       atmos_phy_rd_mstrnx
#endif
    use scale_atmos_phy_rd_mm5sw, only: &
       swinit
    implicit none

    character(len=*), intent(in) :: RD_TYPE
    !---------------------------------------------------------------------------

    select case ( rd_type )
    case ( 'MSTRNX' )
       call atmos_phy_rd_mstrnx_setup( rd_type )
       atmos_phy_rd => atmos_phy_rd_mstrnx
    case ( 'WRF' )
       call atmos_phy_rd_mstrnx_setup( 'MSTRNX' )
       atmos_phy_rd => atmos_phy_rd_mstrnx
       call swinit
    case default
       write(*,*) 'xxx invalid Radiation type(', trim(rd_type), '). CHECK!'
       call prc_mpistop
    end select

    return
  end subroutine atmos_phy_rd_setup

end module scale_atmos_phy_rd