在利用 SAC 的 rms 命令计算波形的均方根时发现了一个问题。

SAC> r SC_HMS.z
SAC> lh b t0
  
  FILE: SC_HMS.z - 1
 --------------
      b = 1.000000e-02
     t0 = 2.277163e+01

SAC> setbb var1 (&1,t0&)
SAC> setbb var2 (&1,t0&+100)
SAC> mtw %var1% %var2%
SAC> rms to user1
SAC> mtw b t0
SAC> rms to user2
SAC> mtw 10 100
SAC> rms to user3
SAC> mtw b 100
SAC> rms to user4
SAC> mtw b %var1%
SAC> rms to user5
SAC> mtw &1,b& %var1%
SAC> rms to user6
SAC> lh user1 user2 user3 user4 user5 user6

 FILE: SC_HMS.z - 1
 --------------

     user1 = 3.900914e-04
     user2 = 6.736934e-06
     user3 = 4.061756e-04
     user4 = -0.000000e+00
     user5 = -0.000000e+00
     user6 = 6.736934e-06

不难发现, mtw 后面的两个值必须为同一类才能算出正确结果,不知道是不是 SAC 的 bug。