28 Jul
2015
28 Jul
'15
9:49 a.m.
Hi again, I was able to solve the problem: I wasn't aware that the forcing in userf is also added to the induction equation. To avoid this, one has to use: ffx=0.0d0 ffy=0.0d0 ffz=0.0d0 if (ifield .eq. 1) then c----velocity forcing here ffx=... ffy=... ffz=... endif as in the mhd example case. My example seems to be running fine (for now). Cheers, Jan