scale_atmos_phy_cp.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
module scale_atmos_phy_cp
  !------------------------------------------------------------------------------
  !
  !+++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  use scale_grid_index
  use scale_tracer
  private

  abstract interface
     subroutine cp( &
          DENS,           &
          MOMZ,           &
          MOMX,           &
          MOMY,           &
          RHOT,           &
          QTRC,           &
          DENS_t_CP,      &
          MOMZ_t_CP,      &
          MOMX_t_CP,      &
          MOMY_t_CP,      &
          RHOT_t_CP,      &
          RHOQ_t_CP,      &
          MFLX_cloudbase, &
          SFLX_convrain,  &
          cloudtop,       &
          cloudbase,      &
          cldfrac_dp,     &
          cldfrac_sh,     &
          kf_nca,         &
          kf_w0avg        )
       use scale_precision
       use scale_stdio
       use scale_prof
       use scale_grid_index
       use scale_tracer
       real(RP), intent(in)    :: dens(ka,ia,ja)
       real(RP), intent(in)    :: momx(ka,ia,ja)
       real(RP), intent(in)    :: momy(ka,ia,ja)
       real(RP), intent(in)    :: momz(ka,ia,ja)
       real(RP), intent(in)    :: rhot(ka,ia,ja)
       real(RP), intent(in)    :: qtrc(ka,ia,ja,qa)
       real(RP), intent(inout) :: dens_t_cp(ka,ia,ja)
       real(RP), intent(inout) :: momz_t_cp(ka,ia,ja)
       real(RP), intent(inout) :: momx_t_cp(ka,ia,ja)
       real(RP), intent(inout) :: momy_t_cp(ka,ia,ja)
       real(RP), intent(inout) :: rhot_t_cp(ka,ia,ja)
       real(RP), intent(inout) :: rhoq_t_cp(ka,ia,ja,qa)
       real(RP), intent(inout) :: mflx_cloudbase(ia,ja)
       real(RP), intent(inout) :: sflx_convrain(ia,ja)
       real(RP), intent(inout) :: cloudtop(ia,ja)
       real(RP), intent(inout) :: cloudbase(ia,ja)
       real(RP), intent(inout) :: cldfrac_dp(ka,ia,ja)
       real(RP), intent(inout) :: cldfrac_sh(ka,ia,ja)
       real(RP), intent(inout) :: kf_nca(ia,ja)
       real(RP), intent(inout) :: kf_w0avg(ka,ia,ja)
     end subroutine cp
  end interface
  procedure(cp), pointer :: atmos_phy_cp => null()
  public :: atmos_phy_cp
  public :: atmos_phy_cp_setup

  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  !------------------------------------------------------------------------------
  subroutine atmos_phy_cp_setup( CP_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 CP
    use name(scale_atmos_phy_mp_, cp,), only: &
       name(atmos_phy_cp_, cp, _setup), &
       name(atmos_phy_cp_, cp,), &
#else
    use scale_atmos_phy_cp_kf, only: &
         atmos_phy_cp_kf_setup, &
         atmos_phy_cp_kf
#endif
    implicit none

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

#ifdef CP
    call name(atmos_phy_cp_, cp, _setup)( cp_type )
    atmos_phy_cp => name(atmos_phy_cp_, cp,)
#else
    select case ( cp_type )
    case( 'KF' )
       call atmos_phy_cp_kf_setup( cp_type )
       atmos_phy_cp => atmos_phy_cp_kf
    end select
#endif

    return
  end subroutine atmos_phy_cp_setup

end module scale_atmos_phy_cp