[cig-commits] [commit] master: output axis s=0 also in the fluid (d62cbe2)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/axisem/compare/ca4b3002407f17f95bbf7013f7622a43ceb4db35...7a0ed4115e0445ccb0219c894eedf4bf157d8689

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

commit d62cbe2fcdba1d23b7a59c5e9670bb083cff03aa
Author: martinvandriel <martin at vandriel.de>
Date:   Mon Oct 13 17:37:47 2014 +0200

    output axis s=0 also in the fluid


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

d62cbe2fcdba1d23b7a59c5e9670bb083cff03aa
 SOLVER/meshes_io.F90 | 16 ++++++++++++++++
 1 file changed, 16 insertions(+)

diff --git a/SOLVER/meshes_io.F90 b/SOLVER/meshes_io.F90
index ab7e618..69c7917 100644
--- a/SOLVER/meshes_io.F90
+++ b/SOLVER/meshes_io.F90
@@ -609,6 +609,22 @@ subroutine build_kwf_grid()
               endif
               mapping_ijel_ikwf(ipol,jpol,iel + nel_solid) = mapping(idest)
           enddo
+
+          ! add gll points on the axis
+          if (axis_fluid(iel)) then
+              ipol = 0
+
+              ipt = (iel-1)*(npol+1)**2 + jpol*(npol+1) + ipol + 1
+              idest = igloc_fluid(ipt)
+              
+              if (.not. check(idest)) then
+                  ct = ct + 1
+                  check(idest) = .true.
+                  mapping(idest) = ct
+                  kwf_mask(ipol,jpol,iel + nel_solid) = .true.
+              endif
+              mapping_ijel_ikwf(ipol,jpol,iel + nel_solid) = mapping(idest)
+          endif
       enddo
   enddo