program Fg_16_17;

uses
   Utils;

var
   Lambda,Omega,B: Double;

begin
   Compile_Full_Dataset(325,1917,True,'North');

   Compute_Method(0,180,5,13.50,15.00,0.001,0.0,6.0,0.1,'North',mtSimple,False,True,Lambda,Omega,B);
   if (Lambda-10)<0 then
       Lambda:=Lambda+180;
   Compute_Method(Lambda-10,Lambda+10,1,Omega-0.005,Omega+0.005,0.0001,B-0.2,B+0.2,0.002,'North',mtSimple,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'North','FG16a');

   Compute_Method(0,180,5,13.50,15.00,0.001,-6.0,0.0,0.1,'North',mtSimple,False,True,Lambda,Omega,B);
   if (Lambda-10)<0 then
       Lambda:=Lambda+180;
   Compute_Method(Lambda-10,Lambda+10,1,Omega-0.005,Omega+0.005,0.0001,B-0.2,B+0.2,0.002,'North',mtSimple,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'North','FG16b');

   Compute_Method(0,180,5,13.50,15.00,0.001,0.0,6.0,0.1,'North',mtDouble,False,True,Lambda,Omega,B);
   if (Lambda-10)<0 then
       Lambda:=Lambda+180;
   Compute_Method(Lambda-10,Lambda+10,1,Omega-0.005,Omega+0.005,0.0001,B-0.2,B+0.2,0.002,'North',mtDouble,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'North','FG17a');

   Compute_Method(0,180,5,13.50,15.00,0.001,-6.0,0.0,0.1,'North',mtDouble,False,True,Lambda,Omega,B);
   if (Lambda-10)<0 then
       Lambda:=Lambda+180;
   Compute_Method(Lambda-10,Lambda+10,1,Omega-0.005,Omega+0.005,0.0001,B-0.2,B+0.2,0.002,'North',mtDouble,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'North','FG17b');

   Compile_Full_Dataset(325,1917,False,'South');

   Compute_Method(0,180,5,13.50,15.00,0.001,0.0,6.0,0.1,'South',mtSimple,False,True,Lambda,Omega,B);
   if (Lambda-10)<0 then
       Lambda:=Lambda+180;
   Compute_Method(Lambda-10,Lambda+10,1,Omega-0.005,Omega+0.005,0.0001,B-0.2,B+0.2,0.002,'South',mtSimple,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'South','');

   Compute_Method(0,180,5,13.50,15.00,0.001,-6.0,0.0,0.1,'South',mtSimple,False,True,Lambda,Omega,B);
   if (Lambda-10)<0 then
       Lambda:=Lambda+180;
   Compute_Method(Lambda-10,Lambda+10,1,Omega-0.005,Omega+0.005,0.0001,B-0.2,B+0.2,0.002,'South',mtSimple,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'South','');
end.