EVOLUTION, GLM-MHD: update Psi after all RK substeps.

- the scalar potential is updated after each Runge-Kutta substep now;
2010-12-01 15:23:27 -02:00 · 2010-12-01 15:23:27 -02:00 · c65688aaf4
commit c65688aaf4
parent c08cf9e98b
1 changed files with 5 additions and 0 deletions
--- a/src/evolution.F90
+++ b/src/evolution.F90
@ -336,6 +336,11 @@ module evolution
 ! update the solution for the scalar potential Psi
 !
    u1(iph,:,:,:) = pblock%u(iph,:,:,:) + ch2 * dt * du(iph,:,:,:)
+
+! evolve Psi due to the source term
+!
+    decay = exp(- 0.5d0 * alpha_p * cmax * dt / dx_min)
+    u1(iph,:,:,:) = decay * u1(iph,:,:,:)
 #endif /* GLM */
 #endif /* MHD */