GMT

Dear All,

I run the provided example ftp://topex.ucsd.edu/pub/sandwell/strain/gpsgridder_tests.tar
I would like to make strain-rate map as on Figure 3 of article https://doi.org/10.1002/2016GL070340.

I try to run gpsgridder based on 02_grid_GPS.sh (from gpsgridder_tests.tar) with additional commands:

*% Creating ASCII nodefile using command
awk '{print $1, $2}' tmp_u.txt > nodes.txt

gmt gpsgridder -Nnodes.txt -Gss_new.out blk.lluv -fg -Vl -Fd8 -Cn700+eigen.txt -S0.5

It's not clear what is the output. "the resulting w values are appended to each record and written to the file given in -G" (as notice in Man Pages).

How to get the max shear strain from gpsgridder output to make strain-rate map?

Attached are gpsgridder -N output and ss_new.out mapped.

Thanks,
Oksana