scale_atmos_dyn_tinteg_tracer.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module scale_atmos_dyn_tinteg_tracer
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  use scale_atmos_grid_cartesc_index
  use scale_index
  use scale_tracer
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: atmos_dyn_tinteg_tracer_setup
 
  abstract interface
     subroutine tinteg( &
          QTRC, & ! (out)
          qflx, & ! (out)
          qtrc0, rhoq_t, &! (in)
          dens0, dens, & ! (in)
          mflx_hi, num_diff, & ! (in)
          gsqrt, mapf, & ! (in)
          cdz, rcdz, rcdx, rcdy, & ! (in)
          bnd_w, bnd_e, bnd_s, bnd_n, & ! (in)
          twod, & ! (in)
          dtl, & ! (in)
          flag_fct_tracer, & ! (in)
          flag_fct_along_stream ) ! (in)
       use scale_precision
       use scale_atmos_grid_cartesc_index
       use scale_index
       real(RP), intent(inout) :: QTRC    (KA,IA,JA)
       real(RP), intent(out)   :: qflx    (KA,IA,JA,3)
       real(RP), intent(in)    :: QTRC0   (KA,IA,JA)
       real(RP), intent(in)    :: RHOQ_t  (KA,IA,JA)
       real(RP), intent(in)    :: DENS0   (KA,IA,JA)
       real(RP), intent(in)    :: DENS    (KA,IA,JA)
       real(RP), intent(in)    :: mflx_hi (KA,IA,JA,3)
       real(RP), intent(in)    :: num_diff(KA,IA,JA,3)
       real(RP), intent(in)    :: GSQRT   (KA,IA,JA,7)
       real(RP), intent(in)    :: MAPF    (IA,JA)
       real(RP), intent(in)    :: CDZ(KA)
       real(RP), intent(in)    :: RCDZ(KA)
       real(RP), intent(in)    :: RCDX(IA)
       real(RP), intent(in)    :: RCDY(JA)
       logical,  intent(in)    :: BND_W
       logical,  intent(in)    :: BND_E
       logical,  intent(in)    :: BND_S
       logical,  intent(in)    :: BND_N
       logical,  intent(in)    :: TwoD
       real(RP), intent(in)    :: dtl
       logical,  intent(in)    :: FLAG_FCT_TRACER
       logical,  intent(in)    :: FLAG_FCT_ALONG_STREAM
     end subroutine tinteg
  end interface
  procedure(tinteg), pointer :: atmos_dyn_tinteg_tracer => null()
  public :: atmos_dyn_tinteg_tracer
 
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine atmos_dyn_tinteg_tracer_setup( &
       ATMOS_DYN_Tinteg_tracer_TYPE )
 
    use scale_precision
    use scale_atmos_grid_cartesc_index
    use scale_index
    use scale_prc, only: &
       prc_abort
    use scale_atmos_dyn_tinteg_tracer_euler, only: &
       atmos_dyn_tinteg_tracer_euler_setup, &
       atmos_dyn_tinteg_tracer_euler
    use scale_atmos_dyn_tinteg_tracer_rk3, only: &
       atmos_dyn_tinteg_tracer_rk3_setup, &
       atmos_dyn_tinteg_tracer_rk3
    use scale_atmos_dyn_tinteg_tracer_linrk, only: &
       atmos_dyn_tinteg_tracer_linrk_setup, &
       atmos_dyn_tinteg_tracer_linrk       
    implicit none
    character(len=*), intent(in)  :: atmos_dyn_tinteg_tracer_type
 
    character(len=H_SHORT) :: tinteg_type
    !---------------------------------------------------------------------------
 
    
    if (atmos_dyn_tinteg_tracer_type(1:5) == 'LINRK') then
      tinteg_type = 'LINRK'
    else
      tinteg_type = atmos_dyn_tinteg_tracer_type
   end if
 
    select case( tinteg_type )
    case( 'EULER' )
       call atmos_dyn_tinteg_tracer_euler_setup( &
            atmos_dyn_tinteg_tracer_type )
       atmos_dyn_tinteg_tracer => atmos_dyn_tinteg_tracer_euler
    case( 'RK3WS2002' )
       call atmos_dyn_tinteg_tracer_rk3_setup( &
            atmos_dyn_tinteg_tracer_type )
       atmos_dyn_tinteg_tracer => atmos_dyn_tinteg_tracer_rk3
    case( 'LINRK' )
       call atmos_dyn_tinteg_tracer_linrk_setup( &
            atmos_dyn_tinteg_tracer_type )
       atmos_dyn_tinteg_tracer => atmos_dyn_tinteg_tracer_linrk       
    case( 'OFF', 'NONE' )
       ! do nothing
    case default
       log_error("ATMOS_DYN_Tinteg_tracer_setup",*) 'ATMOS_DYN_TINTEG_TRACER_TYPE is invalid: ', atmos_dyn_tinteg_tracer_type
       call prc_abort
    end select
 
    return
  end subroutine atmos_dyn_tinteg_tracer_setup
 
end module scale_atmos_dyn_tinteg_tracer