diff --git a/src/problems.F90 b/src/problems.F90
index b71de44..28e5ee6 100644
--- a/src/problems.F90
+++ b/src/problems.F90
@@ -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
diff --git a/src/shapes.F90 b/src/shapes.F90
index 9619e4f..5fe4a91 100644
--- a/src/shapes.F90
+++ b/src/shapes.F90
@@ -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
 !