mu0 = 4e-7 * Pi;
Group {
  CoilP = Region[{101}];
  CoilM = Region[{102}];
  Air   = Region[{103}];
  Dirichlet_a_0 = Region[{200}];
  Vol_Mag    = Region[{CoilP, CoilM, Air}];
  Vol_S0_Mag = Region[{CoilP, CoilM}];
}
Function {
  nu[Vol_Mag] = 1/mu0;
  js[CoilP] = Vector[0., 0., 1e6];
  js[CoilM] = Vector[0., 0., -1e6];
}
Constraint {
  { Name a ;
    Case {
      { Region Dirichlet_a_0 ; Value 0. ; }
    }
  }
}
modelPath = CurrentDirectory;
Include "/usr/share/doc/getdp/examples/templates/Lib_Magnetostatics_a_phi.pro"
eps = 1e-9;
PostOperation {
  { Name B ; NameOfPostProcessing Magnetostatics_a;
    Operation {
      Print[ az, OnElementsOf Vol_Mag, File "az.pos" ] ;
      Print[ b, OnElementsOf Vol_Mag, File "b.pos" ] ;
    }
  }
}