mod_atmos_phy_bl_vars.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module mod_atmos_phy_bl_vars
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  use scale_atmos_grid_cartesc_index
  use scale_tracer
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: atmos_phy_bl_vars_setup

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  integer,  public              :: qs, qe

  real(RP), public, allocatable :: atmos_phy_bl_rhou_t(:,:,:)   ! tendency RHOU [kg/m2/s2]
  real(RP), public, allocatable :: atmos_phy_bl_rhov_t(:,:,:)   ! tendency RHOV [kg/m2/s2]
  real(RP), public, allocatable :: atmos_phy_bl_rhot_t(:,:,:)   ! tendency RHOT [K*kg/m3/s]

  real(RP), public, allocatable :: atmos_phy_bl_rhoq_t(:,:,:,:) ! tendency rho*QTRC [kg/kg/s]

  real(RP), public, allocatable :: atmos_phy_bl_zi    (:,:)     ! depth of the PBL

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_bl_vars_setup
    use scale_prc, only: &
       prc_abort
    use scale_const, only: &
       undef => const_undef
    implicit none

    !---------------------------------------------------------------------------

    log_newline
    log_info("ATMOS_PHY_BL_vars_setup",*) 'Setup'

    allocate( atmos_phy_bl_rhou_t(ka,ia,ja)    )
    allocate( atmos_phy_bl_rhov_t(ka,ia,ja)    )
    allocate( atmos_phy_bl_rhot_t(ka,ia,ja)    )
    allocate( atmos_phy_bl_rhoq_t(ka,ia,ja,qa) )
    atmos_phy_bl_rhou_t(:,:,:)   = undef
    atmos_phy_bl_rhov_t(:,:,:)   = undef
    atmos_phy_bl_rhot_t(:,:,:)   = undef
    atmos_phy_bl_rhoq_t(:,:,:,:) = undef

    allocate( atmos_phy_bl_zi(ia,ja) )
    atmos_phy_bl_zi(:,:) = 100.0_rp ! tentative

    return
  end subroutine atmos_phy_bl_vars_setup

end module mod_atmos_phy_bl_vars