gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Updated comments in the evolution module.

Browse Source
This commit is contained in:
Grzegorz Kowal 2008-12-08 21:07:10 -06:00
parent 1756a144d4
commit 7deada5a90
1 changed files with 14 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
src/evolution.F90
View File

@ -33,12 +33,12 @@ module evolution
contains contains
! !
!====================================================================== !===============================================================================
! !
! evolve: subroutine sweeps over all leaf blocks and performes one step ! evolve: subroutine sweeps over all leaf blocks and performs one step time
! time evolution according to the selected integration scheme ! evolution for each according to the selected integration scheme
! !
!====================================================================== !===============================================================================
! !
subroutine evolve subroutine evolve
@ -50,7 +50,7 @@ module evolution
! !
type(block), pointer :: pblock type(block), pointer :: pblock
! !
!---------------------------------------------------------------------- !-------------------------------------------------------------------------------
! !
! iterate over all blocks and perform one step of time evolution ! iterate over all blocks and perform one step of time evolution
! !
@ -74,16 +74,16 @@ module evolution
! TODO: new time step ! TODO: new time step
! !
!---------------------------------------------------------------------- !-------------------------------------------------------------------------------
! !
end subroutine evolve end subroutine evolve
#ifdef RK2 #ifdef RK2
! !
!====================================================================== !===============================================================================
! !
! evolve_rk2: subroutine evolves the current block using RK2 method ! evolve_rk2: subroutine evolves the current block using RK2 integration
! !
!====================================================================== !===============================================================================
! !
subroutine evolve_rk2(pblock) subroutine evolve_rk2(pblock)
@ -107,7 +107,7 @@ module evolution
! !
real, dimension(nvars,igrids,jgrids,kgrids) :: u1, du real, dimension(nvars,igrids,jgrids,kgrids) :: u1, du
! !
!---------------------------------------------------------------------- !-------------------------------------------------------------------------------
! !
! prepare dxi, dyi, and dzi ! prepare dxi, dyi, and dzi
! !
@ -141,17 +141,18 @@ module evolution
do j = 1, jgrids do j = 1, jgrids
do i = 1, igrids do i = 1, igrids
do q = 1, nvars do q = 1, nvars
pblock%u(q,i,j,k) = 0.5 * (pblock%u(q,i,j,k) + u1(q,i,j,k) + dt*du(q,i,j,k)) pblock%u(q,i,j,k) = 0.5 * (pblock%u(q,i,j,k) + u1(q,i,j,k) &
+ dt*du(q,i,j,k))
end do end do
end do end do
end do end do
end do end do
!---------------------------------------------------------------------- !-------------------------------------------------------------------------------
! !
end subroutine evolve_rk2 end subroutine evolve_rk2
#endif /* RK2 */ #endif /* RK2 */
!====================================================================== !===============================================================================
! !
end module end module