* MagView 3.1 - standard analysis
* Stored quantities
*    Phi: reduced potential (A)
*    Hx: applied magnetic field (A/m)
*    Hy:
*    Hz:
*    Psi: dual potential (A)
*    MuR: relative magnetic permeability
! OPERATOR SET
!   &Bx, &By, &Bz, &BMag, &Hx, &Hy, &Hz, &HMag, &Q[N]

HEADER
  ProgramName = MagView 3.0
  ManualName = magnum.pdf
  WebSiteURL = http://www.fieldp.com
  DataSuffix = GOU
  About = MagView 3.0;Copyright 2009;Field Precision LLC;Albuquerque, NM USA
  NRecord = 3
  SliceDisplay = 4
  ScanDisplay = 4
  VectorDisplay = 1
  PlaneDisplay = 4
  SurfaceDisplay = 4
END

STOREDQUANT
  Phi N
  Hx  N
  Hy  N
  Hz  N
  Psi N
  Mu E
END

REGIONQUANT
  $RegType
  $HCVectX
  $HCVectY
  $HCVectZ
END

PROGPARAM
  $Mu0 = 1.256637E-6
* 1/Mu0
  $IMu0 = 7.957748E5
* Origin for torque calculations, specify in meters
  $xt = 100.00
  $yt = 0.00
  $zt = 0.00
END

RUNPARAM
END

REALTIMEPARAM
END

INTERPOLATION
  Bx = &Bx
  By = &By
  Bz = &Bz
  BMag = &BMag
  Hx = &Hx
  Hy = &Hy
  Hz = &Hz
  HMag = &HMag
  MuR = &Q[6]
  Gamma = 1.0 &Q[6] /
* ===================================================
* Option: Applied magnetic field
*  HxApp = &Q[2]
*  HyApp = &Q[3]
*  HzApp = &Q[4]
*  HMagApp = &Q[2] 2 ^ &Q[3] 2 ^ + &Q[4] 2 ^ + @SQRT
* ===================================================
END

VECTOR
  BVector = &Bx;&By;&Bz
  HVector = &Hx;&Hy;&Hz
* ===================================================
* Option: Applied magnetic field
*  HVectApp = &Q[2];&Q[3];&Q[4]
* ===================================================
* Option: Magnetization direction for permanent magnets
  PMagDirect = $HCVectX;$HCVectY;$HCVectZ
* ===================================================
END

VOLUME
* Fixed-mu materials only
  MagEnergy = 0.5 &BMag * &BMag * &Q[6] $Mu0 * /
END

SURFACE
* Magnetic flux
  MagFlux = &Bx;&By;&Bz
* Force components
  FxSurf = &Bx 2 ^ &BMag 2 ^ 2.0 / - $IMu0 *;&Bx &By * $IMu0 *;&Bx &Bz * $IMu0 *
  FySurf = &By &Bx * $IMu0 *;&By 2 ^ &BMag 2 ^ 2.0 / - $IMu0 *;&By &Bz * $IMu0 *
  FzSurf = &Bz &Bx * $IMu0 *;&Bz &By * $IMu0 *;&Bz 2 ^ &BMag 2 ^ 2.0 / - $IMu0 *
* Torque components
  TxSurf = $y $yt - &Bz &Bx * * $z $zt - &By &Bx * * - $IMu0 *;$y $yt - &Bz &By * * $z $zt - &By 2 ^ &BMag 2 ^ 2.0 / - * - $IMu0 *;$y $yt - &Bz 2 ^ &BMag 2 ^ 2.0 / - * $z $zt - &By &Bz * * - $IMu0 *
  TySurf = $z $zt - &Bx 2 ^ &BMag 2 ^ 2.0 / - * $x $xt - &Bz &Bx * * - $IMu0 *;$z $zt - &Bx &By * * $x $xt - &Bz &By * * - $IMu0 *;$z $zt - &Bx &Bz * * $x $xt - &Bz 2 ^ &BMag 2 ^ 2.0 / - * - $IMu0 *
  TzSurf = $x $xt - &By &Bx * * $y $yt - &Bx 2 ^ &BMag 2 ^ 2.0 / - * - $IMu0 *;$x $xt - &By 2 ^ &BMag 2 ^ 2.0 / - * $y $yt - &Bx &By * * - $IMu0 *;$x $xt - &By &Bz * * $y $yt - &Bx &Bz * * - $IMu0 *
END

ENDFILE

-------------------------------------------------------------------------------
                                  NOTES
-------------------------------------------------------------------------------

s11: &Bx 2 ^ &BMag 2 ^ 2.0 / -
s12: &Bx &By *
s13: &Bx &Bz *

s21: &By &Bx *
s22: &By 2 ^ &BMag 2 ^ 2.0 / -
s23: &By &Bz *

s31: &Bz &Bx *
s32: &Bz &By *
s33: &Bz 2 ^ &BMag 2 ^ 2.0 / -

FxSurf = s11 $IMu0 *;s12 $IMu0 *;s13 $IMu0 *
FySurf = S21 $IMu0 *;S22 $IMu0 *;s23 $IMu0 *
FzSurf = s31 $IMu0 *;s32 $IMu0 *;s33 $IMu0 *

TxSurf =
    $y $yt - S31 * $z $zt - S21 * - $IMu0 *;
    $y $yt - S32 * $z $zt - S22 * - $IMu0 *;
    $y $yt - S33 * $z $zt - S23 * - $IMu0 *;
TySurf =
    $z $zt - S11 * $x $xt - S31 * - $IMu0 *;
    $z $zt - S12 * $x $xt - S32 * - $IMu0 *;
    $z $zt - S13 * $x $xt - S33 * - $IMu0 *;
TzSurf =
    $x $xt - S21 * $y $yt - S11 * - $IMu0 *;
    $x $xt - S22 * $y $yt - S12 * - $IMu0 *;
    $x $xt - S23 * $y $yt - S13 * - $IMu0 *;