program Fg_10_12;

uses
   Utils;

var
   Lambda,Omega,B,ROmega: Double;

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

   Compute_Method(0,180,5,14.15,14.45,0.001,3.36,3.36,0.0,'North',mtSimple,False,False,Lambda,Omega,B);
   Compute_Omega_Vector;
   Write_Vector('FG10a');

   Compute_Method(0,180,5,14.15,14.45,0.001,3.43,3.43,0.0,'North',mtSimple,False,False,Lambda,Omega,B);
   Compute_Omega_Vector;
   Write_Vector('FG10b');

   Compute_Method(0,180,5,14.15,14.45,0.001,3.50,3.50,0.0,'North',mtSimple,False,False,Lambda,Omega,B);
   Compute_Omega_Vector;
   Write_Vector('FG10c');

   Write_Corrections('North','FG11');

   ROmega:=14.272-3.40*SlopeOfCorrections;
   WriteLn('Reparametrized Omega = ',ROmega:12:6);
   Compute_Method(0,180,5,ROmega,ROmega,0.0,3.0,5.0,0.02,'North',mtSimple,False,True,Lambda,Omega,B);
   Compute_B_Vector;
   Write_Vector('FG12');
end.