function [w] = hrthesh(z)

m=find(z*60 <= 52); % z<=0.667
z(m)=0;
w=z;

if 1
m=find(z*60 >= 150); % z>=3.34
z(m)=0;
w=z;
end