#
# set the sum of the velocity to zero
#
grdmath ve_bird.grd 0. MUL = sum.grd
grdmath sum.grd 0 MUL = sum2.grd
#
#   prepare the velocity data for a correlation analysis
#
foreach name (`cat stemlist`)
echo ve_$name.grd
grdmath sum.grd ve_$name.grd ADD = tmp.grd
mv tmp.grd sum.grd
end
#
# compute average
#
grdmath sum.grd 10 DIV = mean.grd
#
# compute the standard deviation
#
foreach name (`cat stemlist`)
echo ve_$name.grd.grd
grdmath ve_bird.grd mean.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 8 DIV SQRT = std.grd
mv mean.grd ve_mean.grd
mv std.grd ve_std.grd
#
#  now do the north compinent
#
# set the sum of the velocity to zero
#
grdmath vn_bird.grd 0. MUL = sum.grd
grdmath sum.grd 0 MUL = sum2.grd
#
#   prepare the vnlocity data for a correlation analysis
#
foreach name (`cat stemlist`)
echo vn_$name.grd
grdmath sum.grd vn_$name.grd ADD = tmp.grd
mv tmp.grd sum.grd
end
#
# compute average
#
grdmath sum.grd 10 DIV = mean.grd
#
# compute the standard deviation
#
foreach name (`cat stemlist`)
echo vn_$name.grd.grd
grdmath vn_bird.grd mean.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 8 DIV SQRT = std.grd
mv mean.grd vn_mean.grd
mv std.grd vn_std.grd
#
#
# set the sum of the velocity to zero
#
grdmath ve_bird.grd 0. MUL = sum.grd
grdmath sum.grd 0 MUL = sum2.grd
#
foreach name (`cat stemlist`)
echo exx_$name.grd
grdmath sum.grd exx_$name.grd ADD = tmp.grd
mv tmp.grd sum.grd
end
#
# compute average
#
grdmath sum.grd 10 DIV = mean.grd
#
# compute the standard deviation
#
foreach name (`cat stemlist`)
echo exx_$name.grd.grd
grdmath exx_bird.grd mean.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 8 DIV SQRT = std.grd
mv mean.grd exx_mean.grd
mv std.grd exx_std.grd
#
# set the sum of the velocity to zero
#
grdmath ve_bird.grd 0. MUL = sum.grd
grdmath sum.grd 0 MUL = sum2.grd
#
foreach name (`cat stemlist`)
echo eyy_$name.grd
grdmath sum.grd eyy_$name.grd ADD = tmp.grd
mv tmp.grd sum.grd
end
#
# compute average
#
grdmath sum.grd 10 DIV = mean.grd
#
# compute the standard deviation
#
foreach name (`cat stemlist`)
echo eyy_$name.grd.grd
grdmath eyy_bird.grd mean.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 8 DIV SQRT = std.grd
mv mean.grd eyy_mean.grd
mv std.grd eyy_std.grd
#
#
# set the sum of the velocity to zero
#
grdmath ve_bird.grd 0. MUL = sum.grd
grdmath sum.grd 0 MUL = sum2.grd
#
foreach name (`cat stemlist`)
echo exy_$name.grd
grdmath sum.grd exy_$name.grd ADD = tmp.grd
mv tmp.grd sum.grd
end
#
# compute average
#
grdmath sum.grd 10 DIV = mean.grd
#
# compute the standard deviation
#
foreach name (`cat stemlist`)
echo exy_$name.grd.grd
grdmath exy_bird.grd mean.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 8 DIV SQRT = std.grd
mv mean.grd exy_mean.grd
mv std.grd exy_std.grd
#
gmt grd2xyz exx_mean.grd > exx.dat
gmt grd2xyz eyy_mean.grd > eyy.dat
gmt grd2xyz exy_mean.grd > exy.dat