#!/bin/bash
#
R=-R-122/-114/32/37
#
#  make and apply the mask if requested
#
if [ "$#" -eq 2 ] 
then
  echo "mask"
  grdmask CGM_mask.xy -Gmask.grd -Rve_$1.grd -N1/NaN/NaN
  grdmath ve_$1.grd mask.grd MUL = ve.grd
  grdmath vn_$1.grd mask.grd MUL = vn.grd
else
  echo "no_mask"
  cp ve_$1.grd ve.grd
  cp vn_$1.grd vn.grd
fi
#
gmt grdmath ve.grd ve_mean.grd SUB = ve_diff.grd
gmt grdmath vn.grd vn_mean.grd SUB = vn_diff.grd
#
#  gert the standard deviation of each, add them and divided by square root of 2
#
gmt grdinfo -L2 ve_diff.grd | grep stdev | awk '{printf $5}' > tmp.dat
gmt grdinfo -L2 vn_diff.grd | grep stdev | awk '{print(" ",$5)}' >> tmp.data
echo -n $1 " " >> stats.dat
awk '{ print ($1+$2)/1.4142 }' < tmp.dat >> stats.dat