Changes between Version 52 and Version 53 of ModulesOnAstroBear

Timestamp:

04/24/17 16:24:10 (8 years ago)

Author:

Jonathan

Comment:

—

ModulesOnAstroBear

@@ -432 +432 @@
   END SUBROUTINE
 }}}
+
+We can also use the ApplyEMF or ApplyBField routines just as we would use ApplyTemplate but the subroutine should return either the Vector Potential (EMF) or the B field respectively.  You should only use ApplyBField if your B-Field is constant - otherwise use the vector potential with ApplyEMF.
+
+{{{
+  SUBROUTINE ProblemGridInit(Info)
+    CALL ApplyEMF(myEMF, Info)
+  END SUBROUTINE
+
+  SUBROUTINE myEMF(pos, t, dt, q)
+    REAL(KIND=qPREC), DIMENSION(:) :: pos
+    REAL(KIND=qPREC), DIMENSION(:) :: q
+    REAL(KIND=qPREC) :: t, dt
+    q=(/0d0,0d0,cos(pos(1))/)
+  END SUBROUTINE myEMF
+}}}
+
+
+
+
  * Note that if {{{MaintainAuxArrays}}} is false then we don't have to worry about aux fields and we can just initialize the q values with volume averaged (or cell-centered).  Currently this only happens in 1D.
  * If {{{MaintainAuxArrays}}} is true then we must initialize all of Info%aux surrounding the '''core''' region.  So first we set the Bx and Bz aux fields to zero.