% Plot results of 2000 May 3 calibration at Palmer Station. See the data file for data, comments.

a=load('PalmerCal2000.txt');

[f,i]=sort(a(:,1));
f=a(i,1);  % Frequency, Hz
V2cal=a(i,5); % Voltage pp at preamp
LR_NS=a(i,6); % Signal back at line receiver, North South
LR_EW=a(i,7); % Signal back at line receiver, East-West
BB_NS=a(i,8); % Signal from I2000 box, broadband monitor, North South
BB_EW=a(i,10); % Signal from I2000 box, broadband monitor, East-West
NB_NS=a(i,9); % Signal from I2000 box, narrowband monitor, North South
NB_EW=a(i,11); % Signal from I2000 box, narrowband monitor, East-West


figure;
loglog(f,LR_NS,'r.-',f,LR_EW,'r.:',f,BB_NS,'g.-',f,BB_EW,'g.:', f,NB_NS,'m.-',f,NB_EW,'m.:');
xlabel('Frequency (Hz)');
ylabel('Voltage, RMS (mV)');
title('Calibration of Palmer VLF receiver');
xlim([min(f) max(f)]);

legend('Line receiver, NS','Line receiver, EW','Broadband, NS','Broadband, EW','Narrowband, NS','Narrowband, EW');