[cig-commits] [commit] devel, master: fixed a typo in convert_all_result_files_to_GRD_format.sh and improved the color scale in display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT.sh (76b1524)

[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 76b152455db13f38c528e1bb996bede0f00353d1
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Fri May 23 17:54:10 2014 +0200

    fixed a typo in convert_all_result_files_to_GRD_format.sh and improved the color scale in display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT.sh


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

76b152455db13f38c528e1bb996bede0f00353d1
 .../convert_all_result_files_to_GRD_format.sh      | 36 +++++++++++++++-------
 ...splay_g_vector_norm_SPECFEM3D_GLOBE_with_GMT.sh |  4 ++-
 2 files changed, 28 insertions(+), 12 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 efa3e6e..b94b103 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,25 +1,39 @@
 #!/bin/bash
 
- surface results_g_x_for_GMT.txt -Gresults_g_x.grd -Rd -I4m -f0x,1y -V
+# 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_y_for_GMT.txt -Gresults_g_y.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_z_for_GMT.txt -Gresults_g_z.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_norm_of_g_for_GMT.txt -Gresults_norm_of_g.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_G_xx_for_GMT.txt -Gresults_G_xx.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_yy_for_GMT.txt -Gresults_G_yy.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_zz_for_GMT.txt -Gresults_G_zz.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_xy_for_GMT.txt -Gresults_G_xy.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_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_xy.grd -Rd -I4m -f0x,1y -V
 
- surface results_G_yz_for_GMT.txt -Gresults_G_yz.grd -Rd -I4m -f0x,1y -V
+ blockmean -Rd -I4m -V results_G_xz_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
 
 # Clean up
-rm -f .gmt*
+rm -f .gmt* ________temp________.xyz
 
diff --git a/utils/Roland_Sylvain_gravity/display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT.sh b/utils/Roland_Sylvain_gravity/display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT.sh
index 5b9a9ae..c03fa7c 100755
--- a/utils/Roland_Sylvain_gravity/display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT.sh
+++ b/utils/Roland_Sylvain_gravity/display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT.sh
@@ -11,7 +11,9 @@ gmtset ANNOT_FONT_SIZE_PRIMARY 10p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:m
 # DK DK for absolute norm of the g vector
 #makecpt  -T9.0/9.2/0.001 -Z > color2.cpt
 # DK DK for the difference between the 3D calculation and the 1D reference calculation
-makecpt  -T-0.03/0.03/0.0001 -Z > color2.cpt
+#makecpt  -T-0.03/0.03/0.0001 -Z > color2.cpt
+#-0.00668716430664 z_max: 0.00217151641846 name: z
+makecpt  -T-0.007/0.007/0.00001 -Z > color2.cpt
 
 ################################