[cig-commits] [commit] devel: added an option (in setup/constants.h.in for now) to display the PMLs in a different color in Postscript vector plots. This option could/should be moved to DATA/Par_file one day. (15a0c3b)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/specfem2d

On branch  : devel
Link       : https://github.com/geodynamics/specfem2d/compare/92314517135e003089b35c0b323b8b78becb0e96...3601a7a6d29d2397e321103ab62c680e8d97805b

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

commit 15a0c3baf69d93d54701ca05af4c3aad9d50ac46
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Wed Feb 26 17:23:33 2014 +0100

    added an option (in setup/constants.h.in for now) to display the PMLs in a different color in Postscript vector plots.
    This option could/should be moved to DATA/Par_file one day.


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

15a0c3baf69d93d54701ca05af4c3aad9d50ac46
 setup/constants.h.in        | 5 +++++
 src/specfem2D/plotpost.F90  | 7 ++++++-
 src/specfem2D/specfem2D.F90 | 6 +++---
 3 files changed, 14 insertions(+), 4 deletions(-)

diff --git a/setup/constants.h.in b/setup/constants.h.in
index 07fe670..fec0cd1 100644
--- a/setup/constants.h.in
+++ b/setup/constants.h.in
@@ -39,6 +39,11 @@
 ! in PostScript files about mesh quality, show all elements above this threshold
   double precision, parameter :: THRESHOLD_POSTSCRIPT = 95.d0 / 100.d0
 
+! display all the PML layers in a different (constant) color in the PostScript vector plots if needed
+!! DK DK these two parameters could/should be moved to DATA/Par_file one day
+  logical, parameter :: DISPLAY_PML_IN_DIFFERENT_COLOR = .true.
+  integer, parameter :: ICOLOR_FOR_PML_DISPLAY = 100
+
 ! option to display only part of the mesh and not the whole mesh,
 ! for instance to analyze Cuthill-McKee mesh partitioning etc.
 ! Possible values are:
diff --git a/src/specfem2D/plotpost.F90 b/src/specfem2D/plotpost.F90
index ea2c39c..44444d0 100644
--- a/src/specfem2D/plotpost.F90
+++ b/src/specfem2D/plotpost.F90
@@ -74,7 +74,7 @@
           coorg_send_ps_free_surface,coorg_recv_ps_free_surface, &
           d1_coorg_send_ps_vector_field,d1_coorg_recv_ps_vector_field, &
           d2_coorg_send_ps_vector_field,d2_coorg_recv_ps_vector_field, &
-          coorg_send_ps_vector_field,coorg_recv_ps_vector_field,US_LETTER)
+          coorg_send_ps_vector_field,coorg_recv_ps_vector_field,US_LETTER,is_PML)
 
 !
 ! PostScript display routine
@@ -99,6 +99,8 @@
   integer ibool(NGLLX,NGLLZ,nspec)
   logical, dimension(nspec) :: poroelastic
 
+  logical, dimension(nspec) :: is_PML
+
   double precision xinterp(pointsdisp,pointsdisp),zinterp(pointsdisp,pointsdisp)
   double precision shapeint(ngnod,pointsdisp,pointsdisp)
   double precision Uxinterp(pointsdisp,pointsdisp)
@@ -2020,6 +2022,9 @@ coorg_recv_ps_vector_field
   imat = kmato(ispec)
   icol = mod(imat - 1,NUM_COLORS) + 1
 
+! display all the PML layers in a different (constant) color if needed
+  if(DISPLAY_PML_IN_DIFFERENT_COLOR .and. is_PML(ispec)) icol = ICOLOR_FOR_PML_DISPLAY
+
   if (  myrank == 0 ) then
     if(meshvect) then
       write(24,680) red(icol),green(icol),blue(icol)
diff --git a/src/specfem2D/specfem2D.F90 b/src/specfem2D/specfem2D.F90
index b2b5035..e3ec6dc 100644
--- a/src/specfem2D/specfem2D.F90
+++ b/src/specfem2D/specfem2D.F90
@@ -8203,7 +8203,7 @@ if(coupled_elastic_poro) then
                       coorg_send_ps_free_surface,coorg_recv_ps_free_surface, &
                       d1_coorg_send_ps_vector_field,d1_coorg_recv_ps_vector_field, &
                       d2_coorg_send_ps_vector_field,d2_coorg_recv_ps_vector_field, &
-                      coorg_send_ps_vector_field,coorg_recv_ps_vector_field,US_LETTER)
+                      coorg_send_ps_vector_field,coorg_recv_ps_vector_field,US_LETTER,is_PML)
 
         else if(imagetype_postscript == 2 .and. p_sv) then
 
@@ -8247,7 +8247,7 @@ if(coupled_elastic_poro) then
                       coorg_send_ps_free_surface,coorg_recv_ps_free_surface, &
                       d1_coorg_send_ps_vector_field,d1_coorg_recv_ps_vector_field, &
                       d2_coorg_send_ps_vector_field,d2_coorg_recv_ps_vector_field, &
-                      coorg_send_ps_vector_field,coorg_recv_ps_vector_field,US_LETTER)
+                      coorg_send_ps_vector_field,coorg_recv_ps_vector_field,US_LETTER,is_PML)
 
         else if(imagetype_postscript == 3 .and. p_sv) then
 
@@ -8291,7 +8291,7 @@ if(coupled_elastic_poro) then
                       coorg_send_ps_free_surface,coorg_recv_ps_free_surface, &
                       d1_coorg_send_ps_vector_field,d1_coorg_recv_ps_vector_field, &
                       d2_coorg_send_ps_vector_field,d2_coorg_recv_ps_vector_field, &
-                      coorg_send_ps_vector_field,coorg_recv_ps_vector_field,US_LETTER)
+                      coorg_send_ps_vector_field,coorg_recv_ps_vector_field,US_LETTER,is_PML)
 
         else if(.not. p_sv) then
           call exit_MPI('cannot draw a SH scalar field as a vector plot, turn PostScript plots off')