scale_atmos_dyn_tinteg_tracer_euler.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "inc_openmp.h"
module scale_atmos_dyn_tinteg_tracer_euler
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  use scale_grid_index
  use scale_index
  use scale_tracer

#ifdef DEBUG
  use scale_debug, only: &
     check
  use scale_const, only: &
     undef  => const_undef, &
     iundef => const_undef2
#endif
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: atmos_dyn_tinteg_tracer_euler_setup
  public :: atmos_dyn_tinteg_tracer_euler

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains

  !-----------------------------------------------------------------------------
  subroutine atmos_dyn_tinteg_tracer_euler_setup( &
       tinteg_type )
    use scale_process, only: &
       prc_mpistop
    use scale_const, only: &
       undef => const_undef
    use scale_comm, only: &
       comm_vars8_init
    implicit none

    character(len=*) :: tinteg_type

    integer :: iv
    !---------------------------------------------------------------------------

    if ( tinteg_type /= 'EULER' ) then
       write(*,*) 'xxx TINTEG_LARGE_TYPE is not EULER. Check!'
       call prc_mpistop
    end if

    return
  end subroutine atmos_dyn_tinteg_tracer_euler_setup

  !-----------------------------------------------------------------------------
  subroutine atmos_dyn_tinteg_tracer_euler( &
       QTRC, & ! (out)
       qtrc0, rhoq_t, &! (in)
       dens0, dens, & ! (in)
       mflx_hi, num_diff, & ! (in)
       gsqrt, mapf, & ! (in)
       cdz, rcdz, rcdx, rcdy, & ! (in)
       dtl, & ! (in)
       flag_fct_tracer, & ! (in)
       flag_fct_along_stream ) ! (in)
    use scale_comm, only: &
       comm_vars8, &
       comm_wait
    use scale_atmos_dyn_tstep_tracer, only: &
       atmos_dyn_tstep_tracer
    implicit none
    real(RP), intent(inout) :: QTRC    (ka,ia,ja)
    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)
    real(RP), intent(in)    :: dtl
    logical,  intent(in)    :: FLAG_FCT_TRACER
    logical,  intent(in)    :: FLAG_FCT_ALONG_STREAM

    !------------------------------------------------------------------------
    ! Start RK
    !------------------------------------------------------------------------

    call atmos_dyn_tstep_tracer( &
         qtrc, & ! (out)
         qtrc, qtrc0, rhoq_t, &! (in)
         dens0, dens, & ! (in)
         mflx_hi, num_diff, & ! (in)
         gsqrt, mapf, & ! (in)
         cdz, rcdz, rcdx, rcdy, & ! (in)
         dtl, & ! (in)
         flag_fct_tracer, flag_fct_along_stream ) ! (in)

    return
  end subroutine atmos_dyn_tinteg_tracer_euler

end module scale_atmos_dyn_tinteg_tracer_euler