program Fg_13_15;

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,0.0,0.0,'North',mtSimple,False,False,Lambda,Omega,B);
   Compute_Omega_Vector;
   Write_Vector('FG13');
   if (Lambda-10)<0 then
      Lambda:=Lambda+180;
   Compute_Method(Lambda-10,Lambda+10,1,Omega-0.005,Omega+0.005,0.0001,0.0,0.0,0.0,'North',mtSimple,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'North','FG14');

   Compute_Method(0,180,5,13.50,15.00,0.001,3.40,3.40,0.0,'North',mtSimple,False,False,Lambda,Omega,B);
   Compute_Omega_Vector;
   Write_Vector('FG15');
   if (Lambda-10)<0 then
      Lambda:=Lambda+180;
   Compute_Method(Lambda-10,Lambda+10,1,Omega-0.005,Omega+0.005,0.0001,3.40,3.40,0.0,'North',mtSimple,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'North','');

   Compute_Method(0,180,5,13.50,15.00,0.001,-3.40,-3.40,0.0,'North',mtSimple,False,False,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,-3.40,-3.40,0.0,'North',mtSimple,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'North','');

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

   Compute_Method(0,180,5,13.50,15.00,0.001,0.0,0.0,0.0,'South',mtSimple,False,False,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,0.0,0.0,0.0,'South',mtSimple,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'South','');

   Compute_Method(0,180,5,13.50,15.00,0.001,3.39,3.39,0.0,'South',mtSimple,False,False,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,3.39,3.39,0.0,'South',mtSimple,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'South','');

   Compute_Method(0,180,5,13.50,15.00,0.001,-3.39,-3.39,0.0,'South',mtSimple,False,False,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,-3.39,-3.39,0.0,'South',mtSimple,False,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'South','');
end.