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
  use scale_atmos_phy_mp
  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
       use scale_atmos_phy_mp
       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_MP)
       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
    use scale_atmos_phy_cp_kf, only: &
         atmos_phy_cp_kf_setup, &
         atmos_phy_cp_kf
    implicit none

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

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

    select case( cp_type )
    case('OFF')
       ! do nothing
    case('KF')
       call atmos_phy_cp_kf_setup( cp_type )
       atmos_phy_cp => atmos_phy_cp_kf
    case default
       write(*,*) 'xxx invalid Cumulus parameterization type(', trim(cp_type), '). CHECK!'
       call prc_mpistop
    end select

    return
  end subroutine atmos_phy_cp_setup

end module scale_atmos_phy_cp