#
# set the sum of the strain to zero
#
grdmath becker_cut.grd 0. MUL = sum.grd
grdmath sum.grd 0 MUL = sum2.grd
#   prepare the strain data for a correlation analysis
#
foreach name (`cat stemlist`)
echo $name\_cut.grd
grdmath sum.grd $name\_cut.grd ADD = tmp.grd
mv tmp.grd sum.grd
end
#
# compute average
#
grdmath sum.grd 16 DIV = mean_cut.grd
#
# compute the standard deviation
#
foreach name (`cat stemlist`)
echo $name\_cut.grd
grdmath $name\_cut.grd mean_cut.grd SUB 2 POW = tmp.grd
grdmath sum2.grd tmp.grd ADD = tmp2.grd
mv tmp2.grd sum2.grd
end
#
#   compute the std
#
grdmath sum2.grd 16 DIV SQRT = std_cut.grd