# Data is weighted by sine squared inside a mute zone. # The weight is zero when t < x * slope0 # The weight is one when t > tp + x * slopep # Suggested defaults: slopep = slope0= 1./1.45 sec/km; tp=.150 sec # subroutine mutter1( tp, slope0,slopep, t0,dt, x, data,nt) integer it, nt real t,wt, tp, slope0,slopep, t0,dt, x, data(nt) x = abs( x) do it=1,nt { t = abs( t0+(it-1)*dt) if ( t < x * slope0) wt = 0 else if( t > tp + x * slopep) wt = 1. else wt = sin( 0.5 * 3.14159265 * (t-x*slope0)/(tp+x*(slopep-slope0))) ** 2 data(it) = data(it) * wt } return; end