[cig-commits] [commit] devel, master: added utils/Roland_Sylvain_gravity/display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT_absolute.sh utils/Roland_Sylvain_gravity/display_g_vector_components_SPECFEM3D_GLOBE_with_GMT_absolute.sh utils/Roland_Sylvain_gravity/display_G_tensor_components_SPECFEM3D_GLOBE_with_GMT_absolute.sh (e0db75e)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Thu Nov 6 08:20:51 PST 2014
Repository : https://github.com/geodynamics/specfem3d_globe
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f
>---------------------------------------------------------------
commit e0db75ee246522ae2b0d3b235a7140fc6d46373d
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Mon Jun 2 16:20:20 2014 +0200
added utils/Roland_Sylvain_gravity/display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT_absolute.sh utils/Roland_Sylvain_gravity/display_g_vector_components_SPECFEM3D_GLOBE_with_GMT_absolute.sh utils/Roland_Sylvain_gravity/display_G_tensor_components_SPECFEM3D_GLOBE_with_GMT_absolute.sh
>---------------------------------------------------------------
e0db75ee246522ae2b0d3b235a7140fc6d46373d
...components_SPECFEM3D_GLOBE_with_GMT_absolute.sh | 120 +++++++++++++++++++++
...onents_SPECFEM3D_GLOBE_with_GMT_differences.sh} | 0
...components_SPECFEM3D_GLOBE_with_GMT_absolute.sh | 69 ++++++++++++
...onents_SPECFEM3D_GLOBE_with_GMT_differences.sh} | 0
...ector_norm_SPECFEM3D_GLOBE_with_GMT_absolute.sh | 33 ++++++
...r_norm_SPECFEM3D_GLOBE_with_GMT_differences.sh} | 0
6 files changed, 222 insertions(+)
diff --git a/utils/Roland_Sylvain_gravity/display_G_tensor_components_SPECFEM3D_GLOBE_with_GMT_absolute.sh b/utils/Roland_Sylvain_gravity/display_G_tensor_components_SPECFEM3D_GLOBE_with_GMT_absolute.sh
new file mode 100755
index 0000000..0b6a501
--- /dev/null
+++ b/utils/Roland_Sylvain_gravity/display_G_tensor_components_SPECFEM3D_GLOBE_with_GMT_absolute.sh
@@ -0,0 +1,120 @@
+#!/bin/bash
+
+# The first plotting command needs - and _only_ needs "-K"
+# All following but the last need both, "-K" and "-O"
+# The last one needs _only_ "-O"
+
+central_meridian=0
+
+gmtset ANNOT_FONT_SIZE_PRIMARY 10p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:mm:ssF
+
+# DK DK for the absolute between the 3D calculation and the 1D reference calculation
+makecpt -T-3000/3000/10 -Z > color2.cpt
+
+################################
+
+file=results_G_xx
+
+## DK DK offset de 4.5cm
+## DK DK -Rd signifie -R-180/180/-90/90
+grdimage new_SPECFEM3D_GLOBE_NEX_160_on_600_cores_whole_Earth_with_topo_and_s40rts/${file}.grd -Ccolor2.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > ${file}_absolute.ps
+
+## DK DK offset de -1.5 inch
+pscoast -Rd -JK$central_meridian/9i -B45g30:."G_xx_absolute": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE gravity calculations by Dimitri Komatitsch" -V -O -K >> ${file}_absolute.ps
+
+psscale -Ccolor2.cpt -D12.5/-1.5/16/0.25h -B2.0:"Absolute value (in Eotvos)": -V -O >> ${file}_absolute.ps
+
+## DK DK convert the final file to PDF
+ps2pdf ${file}_absolute.ps
+rm -f ${file}_absolute.ps
+
+################################
+
+file=results_G_yy
+
+## DK DK offset de 4.5cm
+## DK DK -Rd signifie -R-180/180/-90/90
+grdimage new_SPECFEM3D_GLOBE_NEX_160_on_600_cores_whole_Earth_with_topo_and_s40rts/${file}.grd -Ccolor2.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > ${file}_absolute.ps
+
+## DK DK offset de -1.5 inch
+pscoast -Rd -JK$central_meridian/9i -B45g30:."G_yy_absolute": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE gravity calculations by Dimitri Komatitsch" -V -O -K >> ${file}_absolute.ps
+
+psscale -Ccolor2.cpt -D12.5/-1.5/16/0.25h -B2.0:"Absolute value (in Eotvos)": -V -O >> ${file}_absolute.ps
+
+## DK DK convert the final file to PDF
+ps2pdf ${file}_absolute.ps
+rm -f ${file}_absolute.ps
+
+################################
+
+file=results_G_zz
+
+## DK DK offset de 4.5cm
+## DK DK -Rd signifie -R-180/180/-90/90
+grdimage new_SPECFEM3D_GLOBE_NEX_160_on_600_cores_whole_Earth_with_topo_and_s40rts/${file}.grd -Ccolor2.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > ${file}_absolute.ps
+
+## DK DK offset de -1.5 inch
+pscoast -Rd -JK$central_meridian/9i -B45g30:."G_zz_absolute": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE gravity calculations by Dimitri Komatitsch" -V -O -K >> ${file}_absolute.ps
+
+psscale -Ccolor2.cpt -D12.5/-1.5/16/0.25h -B2.0:"Absolute value (in Eotvos)": -V -O >> ${file}_absolute.ps
+
+## DK DK convert the final file to PDF
+ps2pdf ${file}_absolute.ps
+rm -f ${file}_absolute.ps
+
+################################
+
+file=results_G_xy
+
+## DK DK offset de 4.5cm
+## DK DK -Rd signifie -R-180/180/-90/90
+grdimage new_SPECFEM3D_GLOBE_NEX_160_on_600_cores_whole_Earth_with_topo_and_s40rts/${file}.grd -Ccolor2.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > ${file}_absolute.ps
+
+## DK DK offset de -1.5 inch
+pscoast -Rd -JK$central_meridian/9i -B45g30:."G_xy_absolute": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE gravity calculations by Dimitri Komatitsch" -V -O -K >> ${file}_absolute.ps
+
+psscale -Ccolor2.cpt -D12.5/-1.5/16/0.25h -B2.0:"Absolute value (in Eotvos)": -V -O >> ${file}_absolute.ps
+
+## DK DK convert the final file to PDF
+ps2pdf ${file}_absolute.ps
+rm -f ${file}_absolute.ps
+
+################################
+
+file=results_G_xz
+
+## DK DK offset de 4.5cm
+## DK DK -Rd signifie -R-180/180/-90/90
+grdimage new_SPECFEM3D_GLOBE_NEX_160_on_600_cores_whole_Earth_with_topo_and_s40rts/${file}.grd -Ccolor2.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > ${file}_absolute.ps
+
+## DK DK offset de -1.5 inch
+pscoast -Rd -JK$central_meridian/9i -B45g30:."G_xz_absolute": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE gravity calculations by Dimitri Komatitsch" -V -O -K >> ${file}_absolute.ps
+
+psscale -Ccolor2.cpt -D12.5/-1.5/16/0.25h -B2.0:"Absolute value (in Eotvos)": -V -O >> ${file}_absolute.ps
+
+## DK DK convert the final file to PDF
+ps2pdf ${file}_absolute.ps
+rm -f ${file}_absolute.ps
+
+################################
+
+file=results_G_yz
+
+## DK DK offset de 4.5cm
+## DK DK -Rd signifie -R-180/180/-90/90
+grdimage new_SPECFEM3D_GLOBE_NEX_160_on_600_cores_whole_Earth_with_topo_and_s40rts/${file}.grd -Ccolor2.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > ${file}_absolute.ps
+
+## DK DK offset de -1.5 inch
+pscoast -Rd -JK$central_meridian/9i -B45g30:."G_yz_absolute": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE gravity calculations by Dimitri Komatitsch" -V -O -K >> ${file}_absolute.ps
+
+psscale -Ccolor2.cpt -D12.5/-1.5/16/0.25h -B2.0:"Absolute value (in Eotvos)": -V -O >> ${file}_absolute.ps
+
+## DK DK convert the final file to PDF
+ps2pdf ${file}_absolute.ps
+rm -f ${file}_absolute.ps
+
+################################
+
+# Clean up
+rm -f .gmt* gmt.conf gmt.history color2.cpt
+
diff --git a/utils/Roland_Sylvain_gravity/display_G_tensor_components_SPECFEM3D_GLOBE_with_GMT.sh b/utils/Roland_Sylvain_gravity/display_G_tensor_components_SPECFEM3D_GLOBE_with_GMT_differences.sh
similarity index 100%
rename from utils/Roland_Sylvain_gravity/display_G_tensor_components_SPECFEM3D_GLOBE_with_GMT.sh
rename to utils/Roland_Sylvain_gravity/display_G_tensor_components_SPECFEM3D_GLOBE_with_GMT_differences.sh
diff --git a/utils/Roland_Sylvain_gravity/display_g_vector_components_SPECFEM3D_GLOBE_with_GMT_absolute.sh b/utils/Roland_Sylvain_gravity/display_g_vector_components_SPECFEM3D_GLOBE_with_GMT_absolute.sh
new file mode 100755
index 0000000..c132aab
--- /dev/null
+++ b/utils/Roland_Sylvain_gravity/display_g_vector_components_SPECFEM3D_GLOBE_with_GMT_absolute.sh
@@ -0,0 +1,69 @@
+#!/bin/bash
+
+# The first plotting command needs - and _only_ needs "-K"
+# All following but the last need both, "-K" and "-O"
+# The last one needs _only_ "-O"
+
+central_meridian=0
+
+gmtset ANNOT_FONT_SIZE_PRIMARY 10p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:mm:ssF
+
+# DK DK for absolute norm of the g vector
+makecpt -Chaxby -T-10/10/0.001 > color2.cpt
+
+################################
+
+file=results_g_x
+
+## DK DK offset de 4.5cm
+## DK DK -Rd signifie -R-180/180/-90/90
+grdimage new_SPECFEM3D_GLOBE_NEX_160_on_600_cores_whole_Earth_with_topo_and_s40rts/${file}.grd -Ccolor2.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > ${file}_absolute.ps
+
+## DK DK offset de -1.5 inch
+pscoast -Rd -JK$central_meridian/9i -B45g30:."g_x_absolute": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE gravity calculations by Dimitri Komatitsch" -V -O -K >> ${file}_absolute.ps
+
+psscale -Ccolor2.cpt -D12.5/-1.5/16/0.25h -B0.006:"Absolute value (in m/s2)": -V -O >> ${file}_absolute.ps
+
+## DK DK convert the final file to PDF
+ps2pdf ${file}_absolute.ps
+rm -f ${file}_absolute.ps
+
+################################
+
+file=results_g_y
+
+## DK DK offset de 4.5cm
+## DK DK -Rd signifie -R-180/180/-90/90
+grdimage new_SPECFEM3D_GLOBE_NEX_160_on_600_cores_whole_Earth_with_topo_and_s40rts/${file}.grd -Ccolor2.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > ${file}_absolute.ps
+
+## DK DK offset de -1.5 inch
+pscoast -Rd -JK$central_meridian/9i -B45g30:."g_y_absolute": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE gravity calculations by Dimitri Komatitsch" -V -O -K >> ${file}_absolute.ps
+
+psscale -Ccolor2.cpt -D12.5/-1.5/16/0.25h -B0.006:"Absolute value (in m/s2)": -V -O >> ${file}_absolute.ps
+
+## DK DK convert the final file to PDF
+ps2pdf ${file}_absolute.ps
+rm -f ${file}_absolute.ps
+
+################################
+
+file=results_g_z
+
+## DK DK offset de 4.5cm
+## DK DK -Rd signifie -R-180/180/-90/90
+grdimage new_SPECFEM3D_GLOBE_NEX_160_on_600_cores_whole_Earth_with_topo_and_s40rts/${file}.grd -Ccolor2.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > ${file}_absolute.ps
+
+## DK DK offset de -1.5 inch
+pscoast -Rd -JK$central_meridian/9i -B45g30:."g_z_absolute": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE gravity calculations by Dimitri Komatitsch" -V -O -K >> ${file}_absolute.ps
+
+psscale -Ccolor2.cpt -D12.5/-1.5/16/0.25h -B0.006:"Absolute value (in m/s2)": -V -O >> ${file}_absolute.ps
+
+## DK DK convert the final file to PDF
+ps2pdf ${file}_absolute.ps
+rm -f ${file}_absolute.ps
+
+################################
+
+# Clean up
+rm -f .gmt* gmt.conf gmt.history color2.cpt
+
diff --git a/utils/Roland_Sylvain_gravity/display_g_vector_components_SPECFEM3D_GLOBE_with_GMT.sh b/utils/Roland_Sylvain_gravity/display_g_vector_components_SPECFEM3D_GLOBE_with_GMT_differences.sh
similarity index 100%
rename from utils/Roland_Sylvain_gravity/display_g_vector_components_SPECFEM3D_GLOBE_with_GMT.sh
rename to utils/Roland_Sylvain_gravity/display_g_vector_components_SPECFEM3D_GLOBE_with_GMT_differences.sh
diff --git a/utils/Roland_Sylvain_gravity/display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT_absolute.sh b/utils/Roland_Sylvain_gravity/display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT_absolute.sh
new file mode 100755
index 0000000..5513bf1
--- /dev/null
+++ b/utils/Roland_Sylvain_gravity/display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT_absolute.sh
@@ -0,0 +1,33 @@
+#!/bin/bash
+
+# The first plotting command needs - and _only_ needs "-K"
+# All following but the last need both, "-K" and "-O"
+# The last one needs _only_ "-O"
+
+central_meridian=0
+
+gmtset ANNOT_FONT_SIZE_PRIMARY 10p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:mm:ssF
+
+# DK DK for absolute norm of the g vector
+makecpt -T9.06/9.10/0.001 -Z > color2.cpt
+
+################################
+
+file=results_norm_of_g
+
+## DK DK offset de 4.5cm
+## DK DK -Rd signifie -R-180/180/-90/90
+grdimage new_SPECFEM3D_GLOBE_NEX_160_on_600_cores_whole_Earth_with_topo_and_s40rts/${file}.grd -Ccolor2.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > ${file}_absolute.ps
+
+## DK DK offset de -1.5 inch
+pscoast -Rd -JK$central_meridian/9i -B45g30:."norm_of_g_absolute": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE gravity calculations by Dimitri Komatitsch" -V -O -K >> ${file}_absolute.ps
+
+psscale -Ccolor2.cpt -D12.5/-1.5/16/0.25h -B0.006:"Absolute value (in m/s2)": -V -O >> ${file}_absolute.ps
+
+## DK DK convert the final file to PDF
+ps2pdf ${file}_absolute.ps
+rm -f ${file}_absolute.ps
+
+# Clean up
+rm -f .gmt* gmt.conf gmt.history color2.cpt
+
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_differences.sh
similarity index 100%
rename from utils/Roland_Sylvain_gravity/display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT.sh
rename to utils/Roland_Sylvain_gravity/display_g_vector_norm_SPECFEM3D_GLOBE_with_GMT_differences.sh
More information about the CIG-COMMITS
mailing list