import commands
mean=commands.getstatusoutput('gmt grdinfo -L1 -L2 -C topo1_tra.nc')[1].split('\t')[13]
2019-09-17