[cig-commits] [commit] devel, master: changed blockmedian to blockmean (according to the GMT manual) (ebc9e2c)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem3d_globe

On branches: devel,master
Link       : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f

>---------------------------------------------------------------

commit ebc9e2c6bcd9b48fcc5437d1471b705eca37ce69
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Fri May 23 01:58:08 2014 +0200

    changed blockmedian to blockmean (according to the GMT manual)


>---------------------------------------------------------------

ebc9e2c6bcd9b48fcc5437d1471b705eca37ce69
 .../convert_all_result_files_to_GRD_format.sh      | 24 +++++++++++++---------
 1 file changed, 14 insertions(+), 10 deletions(-)

diff --git a/utils/Roland_Sylvain_gravity/convert_all_result_files_to_GRD_format.sh b/utils/Roland_Sylvain_gravity/convert_all_result_files_to_GRD_format.sh
index 9aebe01..d3fa65c 100755
--- a/utils/Roland_Sylvain_gravity/convert_all_result_files_to_GRD_format.sh
+++ b/utils/Roland_Sylvain_gravity/convert_all_result_files_to_GRD_format.sh
@@ -1,33 +1,37 @@
 #!/bin/bash
 
- blockmedian -Rd -I4m -V results_g_x_for_GMT.txt > ________temp________.xyz
+# from the GMT manual:
+# As a rule of thumb, we use blockmean for most smooth data (such as potential fields)
+# and blockmedian (or modes) for rough, non-Gaussian data (such as topography).
+
+ blockmean -Rd -I4m -V results_g_x_for_GMT.txt > ________temp________.xyz
  surface ________temp________.xyz -Gresults_g_x.grd -Rd -I4m -f0x,1y -V
 
- blockmedian -Rd -I4m -V results_g_y_for_GMT.txt > ________temp________.xyz
+ blockmean -Rd -I4m -V results_g_y_for_GMT.txt > ________temp________.xyz
  surface ________temp________.xyz -Gresults_g_y.grd -Rd -I4m -f0x,1y -V
 
- blockmedian -Rd -I4m -V results_g_z_for_GMT.txt > ________temp________.xyz
+ blockmean -Rd -I4m -V results_g_z_for_GMT.txt > ________temp________.xyz
  surface ________temp________.xyz -Gresults_g_z.grd -Rd -I4m -f0x,1y -V
 
- blockmedian -Rd -I4m -V results_norm_of_g_for_GMT.txt > ________temp________.xyz
+ blockmean -Rd -I4m -V results_norm_of_g_for_GMT.txt > ________temp________.xyz
  surface ________temp________.xyz -Gresults_norm_of_g.grd -Rd -I4m -f0x,1y -V
 
- blockmedian -Rd -I4m -V results_G_xx_for_GMT.txt > ________temp________.xyz
+ blockmean -Rd -I4m -V results_G_xx_for_GMT.txt > ________temp________.xyz
  surface ________temp________.xyz -Gresults_G_xx.grd -Rd -I4m -f0x,1y -V
 
- blockmedian -Rd -I4m -V results_G_yy_for_GMT.txt > ________temp________.xyz
+ blockmean -Rd -I4m -V results_G_yy_for_GMT.txt > ________temp________.xyz
  surface ________temp________.xyz -Gresults_G_yy.grd -Rd -I4m -f0x,1y -V
 
- blockmedian -Rd -I4m -V results_G_zz_for_GMT.txt > ________temp________.xyz
+ blockmean -Rd -I4m -V results_G_zz_for_GMT.txt > ________temp________.xyz
  surface ________temp________.xyz -Gresults_G_zz.grd -Rd -I4m -f0x,1y -V
 
- blockmedian -Rd -I4m -V results_G_xy_for_GMT.txt > ________temp________.xyz
+ blockmean -Rd -I4m -V results_G_xy_for_GMT.txt > ________temp________.xyz
  surface ________temp________.xyz -Gresults_G_xy.grd -Rd -I4m -f0x,1y -V
 
- blockmedian -Rd -I4m -V results_G_xy_for_GMT.txt > ________temp________.xyz
+ blockmean -Rd -I4m -V results_G_xy_for_GMT.txt > ________temp________.xyz
  surface ________temp________.xyz -Gresults_G_xz.grd -Rd -I4m -f0x,1y -V
 
- blockmedian -Rd -I4m -V results_G_yz_for_GMT.txt > ________temp________.xyz
+ blockmean -Rd -I4m -V results_G_yz_for_GMT.txt > ________temp________.xyz
  surface ________temp________.xyz -Gresults_G_yz.grd -Rd -I4m -f0x,1y -V
 
 # Clean up