clf
start_time = [172 22 36 0];
window = 250;
das = 'A988';

[seisTINF,hdr] = segycut(start_time,window,das,'1');

tvec=linspace(0,window,window*hdr.sps);

conv_factor_1INCH_DM_V_V2PA=  1*248.843/(0.01*7/12);
vect_signalPA=-seisTINF*conv_factor_1INCH_DM_V_V2PA*hdr.atod/32;

npoles = 2

high = .25

[B2,A2] = butter(npoles,[high]/(hdr.sps/2),'high')

highpassV = filter(B2,A2,vect_signalPA);

plot(tvec,highpassV);
xlabel('Time (s) after GMT 22:36 Julian Day 172');
ylabel('Pressure (Pa)');
axis tight

set(gcf,'Position',[0 0 700 150])