* RFE2 7.0 - reference-phase analysis
* Field Precision LLC, Copyright 2010
*
* Value of a complex quantity at the reference phase phi0:
*   Q = Qr*cos(phi0) - Qi*sin(phi0)
*
* Variables
*   &PhiR, &PhiI, &ExzR, &ExzI, &EyrR, &EyrI, &EMag, &PDens
* Stored quantities
*   PhiReal  Real particle of the quasi-static potential (V)
*   PhiImag  Imaginary particle of the quasi-static potential (V)
*   PDens    Time-averaged power density (J/m3)
* Region quantities
*  $RegSigma
*  $RegEpsi
*  $RegPotAmp
*  $RegPotPhase

PROGPARAM
  $DegToRad = 0.0174533
  $Epsi0 = 8.854187E-12
  $TwoPi = 6.283185
END

RUNPARAM
  $Freq = 0.0
END

REALTIMEPARAM
  $Phi0 = 0.0
END

INTERPOLATION
* Potential
  PhiRef = $Phi0 $DegToRad * @ENTER @COS &PhiR * @EXCH @SIN &PhiI * -
* Electric field
  ExRef/EzRef = $Phi0 $DegToRad * @ENTER @COS &ExzR * @EXCH @SIN &ExzI * -
  EyRef/ErRef = $Phi0 $DegToRad * @ENTER @COS &EyrR * @EXCH @SIN &EyrI * -
  EMagRef = $Phi0 $DegToRad * @ENTER @COS &ExzR * @EXCH @SIN &ExzI * - 2 ^ $Phi0 $DegToRad * @ENTER @COS &EyrR * @EXCH @SIN &EyrI * - 2 ^ + @SQRT
* Material properties
  EpsiR = $RegEpsi
  Sigma = $RegSigma
* Conductive current: [Sigma*Er, Sigma*Ei]
  JcxRef/JczRef = $Phi0 $DegToRad * @ENTER @COS &ExzR * @EXCH @SIN &ExzI * - $RegSigma *
  JcyRef/JcrRef = $Phi0 $DegToRad * @ENTER @COS &EyrR * @EXCH @SIN &EyrI * - $RegSigma *

  JcMagRef = $Phi0 $DegToRad * @ENTER @COS &ExzR * @EXCH @SIN &ExzI * - 2 ^ $Phi0 $DegToRad * @ENTER @COS &EyrR * @EXCH @SIN &EyrI * - 2 ^ + @SQRT $RegSigma *
* Displacement current: [-2*pi*Freq*Epsi0*EpsiR*Ei,2*pi*Freq*Epsi0*EpsiR*Er]
  JdxRef/JdzRef = $Phi0 $DegToRad * @ENTER @COS &ExzI * @EXCH @SIN &ExzR * + $TwoPi * $Freq * $Epsi0 * $RegEpsi *
  JdyRef/JdrRef = $Phi0 $DegToRad * @ENTER @COS &EyrI * @EXCH @SIN &EyrR * + $TwoPi * $Freq * $Epsi0 * $RegEpsi *
  JdMagRef = $Phi0 $DegToRad * @ENTER @COS &ExzI * @EXCH @SIN &ExzR * + 2 ^ $Phi0 $DegToRad * @ENTER @COS &EyrI * @EXCH @SIN &EyrR * + 2 ^ + @SQRT $TwoPi * $Freq * $Epsi0 * $RegEpsi *
END

VECTOR
* Vector display of electric field direction
  EVect = $Phi0 $DegToRad * @ENTER @COS &ExzR * @EXCH @SIN &ExzI * -;$Phi0 $DegToRad * @ENTER @COS &EyrR * @EXCH @SIN &EyrI * -
END

VOLUME
END

SURFACE
END

ENDFILE