fix boundary condition for a-v formulation

Cohomology/indheat.pro

@@ -16,11 +16,11 @@ Group {
   Cut1TO = Region[{Cut1TO}];
   Cut2TO = Region[{Cut2TO}];
 
-  BdOmega = Region[{2004}];
   Term1 = Region[{2002}];
   Term2 = Region[{2003}];
   Terms = Region[{Term1, Term2}];
 
+  BdOmega = Region[{2004, 2002, 2003}];
 
   DefineConstant[ Cut1AV = {2005, Min 1, Step 1, Closed 1,
       Name "Cohomology/H^1(M_c,S_j)/Cut 1 for AV formulation"} ];
@@ -521,6 +521,7 @@ PostProcessing {
   { Name MagDynAV ; NameOfPostProcessing MagDynAV ;
     Operation {
       Print[ e, OnElementsOf Omega_c , File "eAV.pos"] ;
+      Print[ a, OnElementsOf Omega , File "aAV.pos"] ;
       Print[ b, OnElementsOf Omega , File "bAV.pos"] ;
       Print[ v, OnElementsOf Omega_c , File "vAV.pos"] ;
       Print[ j, OnElementsOf Omega_c , File "jAV.pos"] ;