adds absorbing boundary VTK output

Author: danielpeter

src/generate_databases/save_arrays_solver.F90

@@ -414,7 +414,7 @@ subroutine save_arrays_solver_mesh()
     call save_arrays_solver_injection_boundary()
   endif
 
-  !! setup wavefield discontinuity interface
+  ! setup wavefield discontinuity interface
   if (IS_WAVEFIELD_DISCONTINUITY) then
     call save_arrays_solver_mesh_wavefield_discontinuity()
   endif
@@ -768,6 +768,32 @@ subroutine save_arrays_solver_files()
       deallocate(iglob_tmp)
     endif
 
+    ! saves absorbing surface points
+    if (num_abs_boundary_faces > 0) then
+      ! saves free surface interface points
+      allocate(iglob_tmp(NGLLSQUARE*num_abs_boundary_faces),stat=ier)
+      if (ier /= 0) call exit_MPI_without_rank('error allocating array 652b')
+      if (ier /= 0) stop 'error allocating array iglob_tmp'
+      inum = 0
+      iglob_tmp(:) = 0
+      do i = 1,num_abs_boundary_faces
+        do j = 1,NGLLSQUARE
+          inum = inum+1
+          iglob_tmp(inum) = ibool(abs_boundary_ijk(1,j,i), &
+                                  abs_boundary_ijk(2,j,i), &
+                                  abs_boundary_ijk(3,j,i), &
+                                  abs_boundary_ispec(i) )
+        enddo
+      enddo
+      filename = prname(1:len_trim(prname))//'abs_boundary'
+      call write_VTK_data_points(nglob_unique, &
+                                 xstore_unique,ystore_unique,zstore_unique, &
+                                 iglob_tmp,NGLLSQUARE*num_abs_boundary_faces, &
+                                 filename)
+
+      deallocate(iglob_tmp)
+    endif
+
     ! acoustic-elastic domains
     if (ACOUSTIC_SIMULATION .and. ELASTIC_SIMULATION) then
       ! saves points on acoustic-elastic coupling interface

src/meshfem3D/save_databases.F90

@@ -601,8 +601,9 @@ subroutine save_databases(nspec,nglob, &
   deallocate(material_index)
 
   ! setting up wavefield discontinuity interface
-  if (IS_WAVEFIELD_DISCONTINUITY) &
+  if (IS_WAVEFIELD_DISCONTINUITY) then
     call write_wavefield_discontinuity_database()
+  endif
 
   ! user output
   if (myrank == 0) then