scale_atmos_phy_tb_dummy.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_atmos_phy_tb_dummy
  !-----------------------------------------------------------------------------
  !
  !++ 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_tb_dummy_setup
  public :: atmos_phy_tb_dummy

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_tb_dummy_setup( &
       TYPE_TB,       &
       CDZ, CDX, CDY, &
       CZ             )
    use scale_process, only: &
       prc_mpistop
    implicit none

    character(len=*), intent(in) :: TYPE_TB

    real(RP), intent(in) :: CDZ(ka)
    real(RP), intent(in) :: CDX(ia)
    real(RP), intent(in) :: CDY(ja)
    real(RP), intent(in) :: CZ (ka,ia,ja)
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*) '*** Turbulence Dummy'

    if ( type_tb /= 'DUMMY' ) then
       write(*,*) 'xxx ATMOS_PHY_TB_TYPE is not Dummy. Check!'
       call prc_mpistop
    endif

    return
  end subroutine atmos_phy_tb_dummy_setup

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_tb_dummy( &
       qflx_sgs_momz, qflx_sgs_momx, qflx_sgs_momy, & ! (out)
       qflx_sgs_rhot, qflx_sgs_rhoq,                & ! (out)
       tke,                                         & ! (inout) diagnostic variables
       tke_t, nu_c, ri, pr, n2,                     & ! (out) diagnostic variables
       momz, momx, momy, rhot, dens, qtrc,          & ! (in)
       sflx_mw, sflx_mu, sflx_mv, sflx_sh, sflx_qv, & ! (in)
       gsqrt, j13g, j23g, j33g, mapf, dt            ) ! (in)
    implicit none

    ! SGS flux
    real(RP), intent(out)   :: qflx_sgs_momz(ka,ia,ja,3)
    real(RP), intent(out)   :: qflx_sgs_momx(ka,ia,ja,3)
    real(RP), intent(out)   :: qflx_sgs_momy(ka,ia,ja,3)
    real(RP), intent(out)   :: qflx_sgs_rhot(ka,ia,ja,3)
    real(RP), intent(out)   :: qflx_sgs_rhoq(ka,ia,ja,qa,3)

    real(RP), intent(inout) :: TKE(ka,ia,ja)
    real(RP), intent(out)   :: tke_t(ka,ia,ja) ! tendency TKE

    real(RP), intent(out)   :: nu_C(ka,ia,ja) ! eddy viscosity (center)
    real(RP), intent(out)   :: Ri  (ka,ia,ja) ! Richardson number
    real(RP), intent(out)   :: Pr  (ka,ia,ja) ! Prantle number
    real(RP), intent(out)   :: N2  (ka,ia,ja) ! squared Brunt-Vaisala frequency

    real(RP), intent(in)    :: MOMZ(ka,ia,ja)
    real(RP), intent(in)    :: MOMX(ka,ia,ja)
    real(RP), intent(in)    :: MOMY(ka,ia,ja)
    real(RP), intent(in)    :: RHOT(ka,ia,ja)
    real(RP), intent(in)    :: DENS(ka,ia,ja)
    real(RP), intent(in)    :: QTRC(ka,ia,ja,qa)

    real(RP), intent(in)    :: sflx_mw(ia,ja)
    real(RP), intent(in)    :: sflx_mu(ia,ja)
    real(RP), intent(in)    :: sflx_mv(ia,ja)
    real(RP), intent(in)    :: sflx_sh(ia,ja)
    real(RP), intent(in)    :: sflx_qv(ia,ja)

    real(RP), intent(in)    :: GSQRT(ka,ia,ja,7)
    real(RP), intent(in)    :: J13G (ka,ia,ja,7)
    real(RP), intent(in)    :: J23G (ka,ia,ja,7)
    real(RP), intent(in)    :: J33G
    real(RP), intent(in)    :: MAPF (ia,ja,2,4)
    real(DP), intent(in)    :: dt
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*) '*** Physics step: Turbulence(dummy)'

    ! do nothing
    qflx_sgs_momz(:,:,:,:)   = 0.0_rp
    qflx_sgs_momx(:,:,:,:)   = 0.0_rp
    qflx_sgs_momy(:,:,:,:)   = 0.0_rp
    qflx_sgs_rhot(:,:,:,:)   = 0.0_rp
    qflx_sgs_rhoq(:,:,:,:,:) = 0.0_rp

    tke_t(:,:,:) = 0.0_rp

    tke(:,:,:) = 0.0_rp
    nu_c(:,:,:) = 0.0_rp
    ri(:,:,:) = 0.0_rp
    pr(:,:,:) = 1.0_rp
    n2(:,:,:) = 0.0_rp

    return
  end subroutine atmos_phy_tb_dummy

end module scale_atmos_phy_tb_dummy