%
gpstime=753424982;
gpsstring=num2str(gpstime);
warning off 

f=[0:1/64:7000];
%load Unix version of 1/64 (freq)
load freq1over64
L1=L1IFOparams_753424982;

disp('running L1LSCmodel')
pause(.1)
[olg,sensing,actuation,response]=L1LSCmodel(L1,'darm',f);

%olg is dimensionless; sensing in ASQ-counts/(Lx-Ly)meter; 
%actuation in (Lx-Ly)-meters/DARMCTRL-counts; 
%response is (Lx-Ly)meters/ASQ-count, depsite what LSCmodel.m says

%change all units to strain:
sensing=sensing*L1.darm.armlength;
actuation=actuation/L1.darm.armlength;
response=response/L1.darm.armlength;

  olgfile = strcat('S03V1p1-L1-CAL-OLOOP_GAIN-',gpsstring,'.txt');
  mesg = strcat('Writing file: ',olgfile);
  disp(mesg);
  pause(.1)
  op(:,1)=freq;
  op(:,2)=abs(olg(:,2));
  op(:,3)=angle(olg(:,2));
  save(olgfile,'op','-ASCII','-TABS')

  respfile = strcat('S03V1p1-L1-CAL-RESPONSE-',gpsstring,'.txt');
  mesg = strcat('Writing file: ',respfile);
  disp(mesg);
  pause(.1)
  op(:,2)=abs(response(:,2));
  op(:,3)=angle(response(:,2));
  save(respfile,'op','-ASCII','-TABS')

  sensfile = strcat('S03V1p1-L1-CAL-CAV_GAIN-',gpsstring,'.txt');
  mesg = strcat('Writing file: ',sensfile);
  disp(mesg);
  pause(.1)
  op(:,2)=abs(sensing(:,2));
  op(:,3)=angle(sensing(:,2));
  save(sensfile,'op','-ASCII','-TABS')
  
  actfile = strcat('S03V1p1-L1-CAL-ACTUATION-',gpsstring,'.txt');
  mesg = strcat('Writing file: ',actfile);
  disp(mesg);
  pause(.1)
  op(:,2)=abs(actuation(:,2));
  op(:,3)=angle(actuation(:,2));
  save(actfile,'op','-ASCII','-TABS')