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
    use scale_atmos_phy_rd_mstrnx, only: &
       atmos_phy_rd_mstrnx_setup, &
       atmos_phy_rd_mstrnx
    use scale_atmos_phy_rd_offline, only: &
       atmos_phy_rd_offline_setup, &
       atmos_phy_rd_offline
    use scale_atmos_phy_rd_mm5sw, only: &
       swinit
    implicit none

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

    if( io_l ) write(io_fid_log,*) '*** => ', trim(rd_type), ' is selected.'

    select case( rd_type )
    case('OFF')
       ! do nothing
    case( 'MSTRNX' )
       call atmos_phy_rd_mstrnx_setup( rd_type )
       atmos_phy_rd => atmos_phy_rd_mstrnx
    case( 'OFFLINE' )
       call atmos_phy_rd_offline_setup( rd_type )
       atmos_phy_rd => atmos_phy_rd_offline
    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