MESH, GLM-MHD: prolongate proper variables during refinement.

- in subroutine prolong_block() prolongate the magnetic field
   components and scalar potential during the block refinements if the
   GLM-MHD equations are used;
This commit is contained in:
Grzegorz Kowal 2010-12-01 13:54:46 -02:00
parent 9e2f69c745
commit 72d15710c7

View File

@ -819,6 +819,9 @@ module mesh
use variables , only : nfl, nqt
#ifdef MHD
use variables , only : ibx, iby, ibz
#ifdef GLM
use variables , only : iph
#endif /* GLM */
#endif /* MHD */
implicit none
@ -888,6 +891,17 @@ module mesh
call expand_mag(dm, fm, ng, pdata%u(ibx,:,:,:), pdata%u(iby,:,:,:), pdata%u(ibz,:,:,:) &
, u(ibx,:,:,:), u(iby,:,:,:), u(ibz,:,:,:))
#endif /* FLUXCT */
#ifdef GLM
! expand the cell centered magnetic field components
!
do q = ibx, ibz
call expand(dm, fm, ng, pdata%u(q,:,:,:), u(q,:,:,:), 't', 't', 't')
end do
! expand the scalar potential Psi
!
call expand(dm, fm, ng, pdata%u(iph,:,:,:), u(iph,:,:,:), 't', 't', 't')
#endif /* GLM */
#endif /* MHD */
! iterate over all children
!