parallelatsmodule Module Reference

Data Types
type	parallelatstype
	Extends ATS so we can do synchronization. More...

Functions/Subroutines
subroutine	par_ats_submit_delt (this, kstp, kper, dt, sloc, idir)
	@ brief Allow and external caller to submit preferred time step More...
subroutine	par_ats_set_delt (this, kstp, kper, pertim, perlencurrent, delt)
	@ brief Set time step More...

Function/Subroutine Documentation

par_ats_set_delt()

subroutine parallelatsmodule::par_ats_set_delt ( class(parallelatstype)  this, integer(i4b), intent(in)  kstp, integer(i4b), intent(in)  kper, real(dp), intent(inout)  pertim, real(dp), intent(in)  perlencurrent, real(dp), intent(inout)  delt  )

private

Set the time step length (delt) for this time step using the ATS controls, but with synchronization across ranks.

Definition at line 60 of file ParallelAts.f90.

    class(ParallelAtsType) :: this
    integer(I4B), intent(in) :: kstp
    integer(I4B), intent(in) :: kper
    real(DP), intent(inout) :: pertim
    real(DP), intent(in) :: perlencurrent
    real(DP), intent(inout) :: delt
    ! local
    type(MpiWorldType), pointer :: mpi_world
    integer :: ierr
    real(DP) :: global_delt
 
    ! reduce delt over all ranks
    mpi_world => get_mpi_world()
    call mpi_allreduce(delt, global_delt, 1, mpi_double_precision, &
                       mpi_min, mpi_world%comm, ierr)
    call check_mpi(ierr)
 
    delt = global_delt
 
    call this%AtsType%ats_set_delt(kstp, kper, pertim, perlencurrent, delt)
 

Here is the call graph for this function:

par_ats_submit_delt()

subroutine parallelatsmodule::par_ats_submit_delt ( class(parallelatstype)  this, integer(i4b), intent(in)  kstp, integer(i4b), intent(in)  kper, real(dp), intent(in)  dt, character(len=*), intent(in)  sloc, integer(i4b), intent(in), optional  idir  )

private

Submit a preferred time step length. Alternatively, if idir is is passed, then either increase or decrease the submitted time step by the dtadj input variable.

Definition at line 29 of file ParallelAts.f90.

    class(ParallelAtsType) :: this
    integer(I4B), intent(in) :: kstp
    integer(I4B), intent(in) :: kper
    real(DP), intent(in) :: dt
    character(len=*), intent(in) :: sloc
    integer(I4B), intent(in), optional :: idir
    ! local
    type(MpiWorldType), pointer :: mpi_world
    integer :: ierr
    real(DP) :: global_dt
 
    if (.not. this%isAdaptivePeriod(kper)) return
 
    ! reduce dtstable over all ranks
    mpi_world => get_mpi_world()
    call mpi_allreduce(this%dtstable, global_dt, 1, mpi_double_precision, &
                       mpi_min, mpi_world%comm, ierr)
    call check_mpi(ierr)
 
    this%dtstable = global_dt
 
    call this%AtsType%ats_submit_delt(kstp, kper, dt, sloc, idir)
 

Here is the call graph for this function: