Advance all the grids over one time step, including all sources
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | iit |
subroutine advance(iit) use mod_global_parameters use mod_particles, only: handle_particles use mod_source, only: add_split_source integer, intent(in) :: iit integer :: iigrid, igrid, idimsplit ! split source addition call add_split_source(prior=.true.) if (dimsplit) then if ((iit/2)*2==iit .or. typedimsplit=='xy') then ! do the sweeps in order of increasing idim, do idimsplit=1,ndim call advect(idimsplit,idimsplit) end do else ! If the parity of "iit" is odd and typedimsplit=xyyx, ! do sweeps backwards do idimsplit=ndim,1,-1 call advect(idimsplit,idimsplit) end do end if else ! Add fluxes from all directions at once call advect(1,ndim) end if ! split source addition call add_split_source(prior=.false.) if(use_particles) call handle_particles ! do iigrid=1,igridstail; igrid=igrids(iigrid); ! ! end do end subroutine advance