
%   Plot the pressure at Bayreuth (50deg4min N 9deg1minE)
%   16-20 May, 1910.

Time  = [ 16+(7/24); 16+(14/24); 16+(21/24); 
          17+(7/24); 17+(14/24); 17+(21/24); 
          18+(7/24); 18+(14/24); 18+(21/24); 
          19+(7/24); 19+(14/24); 19+(21/24); 
          20+(7/24); 20+(14/24); 20+(21/24);
          21+(7/24); 21+(14/24); 21+(21/24);
          22+(7/24); 22+(14/24); 22+(21/24)  ];
       
PressKahl = [46.9; 46.5; 46.5;
         46.5; 46.1; 46.0; 
         46.7; 47.4; 47.2; 
         46.3; 44.9; 45.1;
         45.3; 44.8; 46.0; 
         47.2; 47.3; 49.2;
         50.5; 50.3; 51.4  ];

PressBay = [ 25.5; 24.9; 25.1;
             25.2; 24.9; 25.6;
             26.3; 25.9; 26.2;
             25.8; 25.0; 25.3;
             25.6; 24.7; 25.6; 
             26.8; 27.0; 28.2;
             30.3; 29.7; 31.1 ];

PressmmHg = 700 + PressBay;
PresshPa = (4/3)*PressmmHg;

plot(Time,PressmmHg,'o',Time,PressmmHg)
axis([16,23.5,720,736])
grid on
xlabel('16-23 May 1910')
%  ylabel('Pressure')
title('Surface Pressure (mm Hg) at Bayreuth', ...
      'FontWeight','bold', 'FontSize',12);
print -deps Bayreuth-Hg.eps
pause

plot(Time,PresshPa,'o',Time,PresshPa)
axis([16,23.5,960,980]);
grid on
xlabel('16-23 May 1910')
%  ylabel('Pressure')
title('Surface Pressure (hPa) at Bayreuth', ...
      'FontWeight','bold', 'FontSize',12);

print -deps Bayreuth-hPa.eps

