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

PROBLEMS: Add magnetic field to the jet problem.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2015-02-17 10:25:52 -02:00
parent 1930bda0cc
commit e5e0d293cf
2 changed files with 22 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/problems.F90
View File

@ -1273,7 +1273,7 @@ module problems
use coordinates, only : ax, ay, az
use equations , only : prim2cons
use equations , only : nv
use equations , only : idn, ivx, ivy, ivz, ipr
use equations , only : idn, ivx, ivy, ivz, ipr, ibx, iby, ibz, ibp
use parameters , only : get_parameter_real
! local variables are not implicit by default
@ -1289,6 +1289,7 @@ module problems
real(kind=8), save :: djet = 1.00d-01
real(kind=8), save :: damb = 1.00d+01
real(kind=8), save :: pres = 1.00d-02
real(kind=8), save :: bphi = 1.00d-03
real(kind=8), save :: vjet = 0.99d+00
real(kind=8), save :: ljet = 1.00d-00
real(kind=8), save :: rjet = 1.00d+00
@ -1327,6 +1328,7 @@ module problems
call get_parameter_real("djet" , djet)
call get_parameter_real("damb" , damb)
call get_parameter_real("pres" , pres)
call get_parameter_real("bphi" , bphi)
call get_parameter_real("vjet" , vjet)
call get_parameter_real("ljet" , ljet)
call get_parameter_real("rjet" , rjet)
@ -1373,6 +1375,16 @@ module problems
q(ivy,1:im) = 0.0d+00
q(ivz,1:im) = 0.0d+00
! if magnetic field is present, set it to be uniform with the desired strength
! and orientation
!
if (ibx > 0) then
q(ibx,1:im) = 0.0d+00
q(iby,1:im) = 0.0d+00
q(ibz,1:im) = bphi
q(ibp,1:im) = 0.0d+00
end if ! ibx > 0
! set the jet injection
!
if (rr <= max(rm, rjet2)) then

10
src/shapes.F90
View File

@ -568,7 +568,7 @@ module shapes
use coordinates , only : ax, ay, az
use equations , only : prim2cons
use equations , only : nv
use equations , only : idn, ivx, ivy, ivz, ipr
use equations , only : idn, ivx, ivy, ivz, ipr, ibx, iby, ibz, ibp
use parameters , only : get_parameter_real
! local variables are not implicit by default
@ -583,6 +583,7 @@ module shapes
!
real(kind=8), save :: djet = 1.00d-01
real(kind=8), save :: pres = 1.00d-02
real(kind=8), save :: bphi = 1.00d-03
real(kind=8), save :: vjet = 0.99d+00
real(kind=8), save :: ljet = 1.00d+00
real(kind=8), save :: rjet = 1.00d+00
@ -621,6 +622,7 @@ module shapes
!
call get_parameter_real("djet" , djet)
call get_parameter_real("pres" , pres)
call get_parameter_real("bphi" , bphi)
call get_parameter_real("ljet" , ljet)
call get_parameter_real("vjet" , vjet)
call get_parameter_real("rjet" , rjet)
@ -642,6 +644,12 @@ module shapes
qj(ivx) = vjet
qj(ivy) = 0.0d+00
qj(ivz) = 0.0d+00
if (ibx > 0) then
qj(ibx) = 0.0d+00
qj(iby) = 0.0d+00
qj(ibz) = bphi
qj(ibp) = 0.0d+00
end if ! ibx > 0
! prepare block coordinates
!