[cig-commits] [commit] devel: removed "blockmean", which was smoothing out the results way too much (3027b17)

[cig_noreply at geodynamics.org]

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

On branch  : devel
Link       : https://github.com/geodynamics/specfem3d_globe/compare/c5007a6f6e799a4a534271f14c4518471b727614...0c658eb11718da1e846bf3268964771a4605210d

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

commit 3027b17d5bf59805a67ceaed3a4eaf5e3bec2217
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Fri May 23 17:26:49 2014 +0200

    removed "blockmean", which was smoothing out the results way too much


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

3027b17d5bf59805a67ceaed3a4eaf5e3bec2217
 .../convert_all_result_files_to_GRD_format.sh      | 36 +++++++---------------
 1 file changed, 11 insertions(+), 25 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 d3fa65c..efa3e6e 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,39 +1,25 @@
 #!/bin/bash
 
-# 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).
+ surface results_g_x_for_GMT.txt -Gresults_g_x.grd -Rd -I4m -f0x,1y -V
 
- blockmean -Rd -I4m -V results_g_x_for_GMT.txt > ________temp________.xyz
- surface ________temp________.xyz -Gresults_g_x.grd -Rd -I4m -f0x,1y -V
+ surface results_g_y_for_GMT.txt -Gresults_g_y.grd -Rd -I4m -f0x,1y -V
 
- blockmean -Rd -I4m -V results_g_y_for_GMT.txt > ________temp________.xyz
- surface ________temp________.xyz -Gresults_g_y.grd -Rd -I4m -f0x,1y -V
+ surface results_g_z_for_GMT.txt -Gresults_g_z.grd -Rd -I4m -f0x,1y -V
 
- blockmean -Rd -I4m -V results_g_z_for_GMT.txt > ________temp________.xyz
- surface ________temp________.xyz -Gresults_g_z.grd -Rd -I4m -f0x,1y -V
+ surface results_norm_of_g_for_GMT.txt -Gresults_norm_of_g.grd -Rd -I4m -f0x,1y -V
 
- 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
+ surface results_G_xx_for_GMT.txt -Gresults_G_xx.grd -Rd -I4m -f0x,1y -V
 
- blockmean -Rd -I4m -V results_G_xx_for_GMT.txt > ________temp________.xyz
- surface ________temp________.xyz -Gresults_G_xx.grd -Rd -I4m -f0x,1y -V
+ surface results_G_yy_for_GMT.txt -Gresults_G_yy.grd -Rd -I4m -f0x,1y -V
 
- blockmean -Rd -I4m -V results_G_yy_for_GMT.txt > ________temp________.xyz
- surface ________temp________.xyz -Gresults_G_yy.grd -Rd -I4m -f0x,1y -V
+ surface results_G_zz_for_GMT.txt -Gresults_G_zz.grd -Rd -I4m -f0x,1y -V
 
- blockmean -Rd -I4m -V results_G_zz_for_GMT.txt > ________temp________.xyz
- surface ________temp________.xyz -Gresults_G_zz.grd -Rd -I4m -f0x,1y -V
+ surface results_G_xy_for_GMT.txt -Gresults_G_xy.grd -Rd -I4m -f0x,1y -V
 
- blockmean -Rd -I4m -V results_G_xy_for_GMT.txt > ________temp________.xyz
- surface ________temp________.xyz -Gresults_G_xy.grd -Rd -I4m -f0x,1y -V
+ surface results_G_xz_for_GMT.txt -Gresults_G_xz.grd -Rd -I4m -f0x,1y -V
 
- blockmean -Rd -I4m -V results_G_xy_for_GMT.txt > ________temp________.xyz
- surface ________temp________.xyz -Gresults_G_xz.grd -Rd -I4m -f0x,1y -V
-
- blockmean -Rd -I4m -V results_G_yz_for_GMT.txt > ________temp________.xyz
- surface ________temp________.xyz -Gresults_G_yz.grd -Rd -I4m -f0x,1y -V
+ surface results_G_yz_for_GMT.txt -Gresults_G_yz.grd -Rd -I4m -f0x,1y -V
 
 # Clean up
-rm -f .gmt* ________temp________.xyz
+rm -f .gmt*