I want to rotate the 1,2 components of OBS to the R, T components, but the SAC command running results show that the two components are not horizontal. Could anyone can teach me how to correct the two horizontal components?
You must change CMPAZ and CMPINC of BHN, BHE components before rotating as follows:
% sac SAC> r BK.PKD.HHN.SAC # read the north component seismogram SAC> ch CMPAZ 0 CMPINC 90 # change the incident and azimuth of the instrument to be north and horizontal SAC> wh # save the header (the above three lines are not needed if the headers are already in) SAC> r BK.PKD.HHE.SAC # read the east component seismogram SAC> ch CMPAZ 90 CMPINC 90 # change the incident and azimuth of the instrument to be east and horizontal SAC> wh # save the header (the above three lines are not needed if the headers are already in) SAC> r BK.PKD.HHN.SAC BK.PKD.HHE.SAC # read two horizontal component seismograms SAC> rotate to GCP # rotate to the "great circle path" SAC> w BK.PKD.HHR.SAC BK.PKD.HHT.SAC # save the rotated data (first radial, second transverse) SAC> r BK.PKD.HHZ.SAC BK.PKD.HHR.SAC BK.PKD.HHT.SAC # read three-component seismograms (vertical, radial and transverse) SAC> qdp off # turn off "quick and dirty plot" SAC> xlim 300 1500 # zoom in SAC> p1 # plot them out, try to identify Love and Rayleigh waves