[cig-commits] [commit] devel,master: added "blockmedian" to GMT script; added output of time stamp file for ispec == 1 (1cb1024)

[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 1cb1024aea1eeb45b8b583e9cdb8f792c13d23a8
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Fri May 23 01:43:53 2014 +0200

    added "blockmedian" to GMT script;
    added output of time stamp file for ispec == 1


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

1cb1024aea1eeb45b8b583e9cdb8f792c13d23a8
 src/meshfem3D/compute_volumes_and_areas.F90        |  2 +-
 .../convert_all_result_files_to_GRD_format.sh      | 41 ++++++++++++++++------
 2 files changed, 31 insertions(+), 12 deletions(-)

diff --git a/src/meshfem3D/compute_volumes_and_areas.F90 b/src/meshfem3D/compute_volumes_and_areas.F90
index 8a15087..90460f0 100644
--- a/src/meshfem3D/compute_volumes_and_areas.F90
+++ b/src/meshfem3D/compute_volumes_and_areas.F90
@@ -306,7 +306,7 @@
   do ispec = 1,nspec
 
     ! print information about number of elements done so far
-    if(myrank == 0 .and. (mod(ispec,NSPEC_DISPLAY_INTERVAL) == 0 .or. ispec == nspec)) then
+    if(myrank == 0 .and. (mod(ispec,NSPEC_DISPLAY_INTERVAL) == 0 .or. ispec == 1 .or. ispec == nspec)) then
        write(IMAIN,*) 'for Roland_Sylvain integrals ',ispec,' elements computed out of ',nspec
        ! write time stamp file to give information about progression of simulation
        write(outputname,"('/timestamp_reg',i1.1,'_ispec',i7.7,'_out_of_',i7.7)") iregion_code,ispec,nspec
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
old mode 100644
new mode 100755
index 95a3630..9aebe01
--- 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,16 +1,35 @@
 #!/bin/bash
 
- surface results_g_x_for_GMT.txt -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
- surface results_g_z_for_GMT.txt -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
- surface results_G_xx_for_GMT.txt -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
- surface results_G_zz_for_GMT.txt -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
- surface results_G_xz_for_GMT.txt -Gresults_G_xz.grd -Rd -I4m -f0x,1y -V
- surface results_G_yz_for_GMT.txt -Gresults_G_yz.grd -Rd -I4m -f0x,1y -V
+ blockmedian -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
+ surface ________temp________.xyz -Gresults_g_y.grd -Rd -I4m -f0x,1y -V
+
+ blockmedian -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
+ 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
+ surface ________temp________.xyz -Gresults_G_xx.grd -Rd -I4m -f0x,1y -V
+
+ blockmedian -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
+ surface ________temp________.xyz -Gresults_G_zz.grd -Rd -I4m -f0x,1y -V
+
+ blockmedian -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
+ surface ________temp________.xyz -Gresults_G_xz.grd -Rd -I4m -f0x,1y -V
+
+ blockmedian -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