python处理地震sac数据_Fortran批量读取SAC地震数据文件，写SAC文件，并转换成SEGY...

[Fortran] 纯文本查看 复制代码Module SACParameter

Integer, parameter :: FileUnit = 100

Character( Len = * ), Parameter :: Directer = '"'//'..\WAV\WAVFormYiliang\yiliang\*.*'//'"'

Type :: EventPara

Integer( kind = 4 ) :: StationNumber

Integer( kind = 4 ) :: FileNumber

Character( Len = 4 ), allocatable :: StationName(:) ! 台站名数组

Character( Len = 4 ), allocatable :: ComponentName(:,:) ! 三分量文件名(台站,分量)

Character( Len = 1 ), allocatable :: stringFileName(:,:) ! 文件名信息数组二维字符数组

Real(kind = 4), allocatable :: StationWaveformLong(:,:) ! 重采样后的波形长度(台站,分量)

Real(kind = 4), allocatable :: StationMaximumAmplitude(:,:)! 归一化后的最大振幅(台站,分量)

Real(kind = 4), allocatable :: StationDELTA(:,:) ! 重采样采样间隔(台站,分量)

End Type EventPara

Type :: SACHeader

Real( Kind = 4 ), Allocatable :: DELTA(:,:), DEPMIN(:,:), DEPMAX(:,:), SCALE(:,:), ODELTA(:,:), B(:,:), E(:,:), &

O(:,:), A(:,:), INTERNAL1(:,:), T0(:,:), T1(:,:), T2(:,:), T3(:,:), T4(:,:), T5(:,:), T6(:,:), T7(:,:), T8(:,:), T9(:,:), F(:,:), &

RESP0(:,:),RESP1(:,:),RESP2(:,:),RESP3(:,:),RESP4(:,:),RESP5(:,:),RESP6(:,:),RESP7(:,:),RESP8(:,:),RESP9(:,:), &

STLA(:,:),STLO(:,:),STEL(:,:),STDP(:,:),EVLA(:,:),EVLO(:,:),EVEL(:,:),EVDP(:,:),MAG(