Use Print[{}, Format, File] in PostOperation

PotentialFlow/magnus.pro

@@ -257,9 +257,6 @@ Resolution{
 	// The run-time variable $newDmdt is used in Generate[A]
         // whereas $circ, $dmdt, $argV and $phiTrailing are evaluated 
 	// by the PostOperation[Trailing].
-	// A file containing "PhiTrailing = $phiTrailing;" is 
-	// after convergence written on disk to be used 
-	// for formatting purposes in the PostOperation section.  
 
 	DeleteFile["KJiter.txt"];
 
@@ -290,7 +287,6 @@ Resolution{
 	  Print[{$syscount, $circ, $dmdt, $jac, $argV}, 
 		Format "iter = %3g Circ = %5.2f Dmdt = %5.2f jac=%5.2f argV = %5.3e"];
 	}
-	Print[{$phiTrailing}, Format "PhiTrailing=%g;", File "PhiTrailing.txt"];
       EndIf
     }
   }
@@ -380,22 +376,25 @@ PostOperation{
       Echo[Str["l=PostProcessing.NbViews-1;",
 	       "View[l].Name = 'stagnation points';",
 	       "View[l].ScaleType = 2; // log scale"],
-	    File "tmp.geo", LastTimeStepOnly] ;
+	   File "tmp.geo"] ;
 
+      If( Flag_Object == 0 )
+	Print[phi, OnElementsOf Vol_rho, File "phi.pos"];
+      Else
 	// Show the isovalue "phi=$PhiTrailing"
-      // which should be perpendicular to the airfoil
-      // if the Kutta condition is fulfilled
+	// which is perpendicular to the airfoil
+	// iff the Kutta condition is fulfilled
 	Print[phi, OnElementsOf Vol_rho, File "KJ.pos"];
-      Echo[Str["Merge 'PhiTrailing.txt';",
-	       "l=PostProcessing.NbViews-1;",
+        Print[{$phiTrailing, 1.001*$phiTrailing}, Format
+	      Str["l=PostProcessing.NbViews-1;",
 		  "View[l].Name = 'isovalue phiTrailing';",
 		  "View[l].IntervalsType = 3;",
 		  "Mesh.SurfaceEdges = 0; // hide mesh",
 		  "View[l].RangeType = 2; // custom range",
-	       "View[l].CustomMin = PhiTrailing;",
-	       "View[l].CustomMax = 1.001*PhiTrailing;",
-	       "View[l].NbIso = 10;"], 
-	   File "tmp.geo", LastTimeStepOnly] ;
+		  "View[l].CustomMin = %g;",
+		  "View[l].CustomMax = %g;"], 
+	      File "tmp.geo"] ;
+      EndIf
     }
   }
 }
@@ -409,6 +408,7 @@ PostOperation{ // for the Airfoil model
        Print[Dmdt, OnRegion Sur_Cut, File > "KJiter.txt", Format Table,
      	     StoreInVariable $dmdt, 
      	     SendToServer "Output/Dmdt" ];
+       // P_edge = {1.0001,0,0} 
 	Print[phi, OnPoint {1.0001,0,0}, File > "KJiter.txt", Color "Ivory",
 	      StoreInVariable $phiTrailing, 
 	      Format Table, SendToServer "Output/PhiTrailing"]; 
@@ -420,7 +420,6 @@ PostOperation{ // for the Airfoil model
   }
 }
 
-
 DefineConstant[
   R_ = {"PotentialFlow", Name "GetDP/1ResolutionChoices", Visible 0},
   C_ = {"-solve -pos", Name "GetDP/9ComputeCommand", Visible 0},