[cig-commits] [commit] devel: Coupling with DSM (610cc8a)

[cig_noreply at geodynamics.org]

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

On branch  : devel
Link       : https://github.com/geodynamics/specfem3d/compare/8a3f14d7d473f70feb7f073639045daa35c587bc...d759e09dd946c593868753fbb4253d77378fb276

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

commit 610cc8aaf8083135db3c164e37ebf2b08982186c
Author: Clément Durochat <c.durochat at gmail.com>
Date:   Wed Jun 18 18:03:55 2014 +0200

    Coupling with DSM


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

610cc8aaf8083135db3c164e37ebf2b08982186c
 setup/constants.h.in                          |  2 +-
 src/generate_databases/save_arrays_solver.f90 | 32 +++++++++++++++++++++++++++
 2 files changed, 33 insertions(+), 1 deletion(-)

diff --git a/setup/constants.h.in b/setup/constants.h.in
index 3c74111..b4a7be7 100644
--- a/setup/constants.h.in
+++ b/setup/constants.h.in
@@ -157,7 +157,7 @@
 
 !!-----------------------------------------------------------
 !!
-!! Roland_Sylvain gravity calculations
+!! For coupling with DSM by VM
 !!
 !!-----------------------------------------------------------
 
diff --git a/src/generate_databases/save_arrays_solver.f90 b/src/generate_databases/save_arrays_solver.f90
index 95a69f3..587aebe 100644
--- a/src/generate_databases/save_arrays_solver.f90
+++ b/src/generate_databases/save_arrays_solver.f90
@@ -571,6 +571,38 @@
       deallocate(iglob_tmp,v_tmp_i)
     endif
 
+    !! C. DUROCHAT modification : begin !! for coupling with DSM by VM
+    if (COUPLE_WITH_DSM) then
+      !if (num_abs_boundary_faces > 0) then
+      filename = prname(1:len_trim(prname))//'absorb_dsm'
+      open(27,file=filename(1:len_trim(filename)),status='unknown',form='unformatted',iostat=ier)
+      if( ier /= 0 ) stop 'error opening file absorb_dsm'
+      write(27) num_abs_boundary_faces
+      write(27) abs_boundary_ispec
+      write(27) abs_boundary_ijk
+      write(27) abs_boundary_jacobian2Dw
+      write(27) abs_boundary_normal
+      close(27)
+
+      filename = prname(1:len_trim(prname))//'inner'
+      open(27,file=filename(1:len_trim(filename)),status='unknown',form='unformatted',iostat=ier)
+      write(27) ispec_is_inner
+      write(27) ispec_is_elastic
+      close(27)
+        
+      !end if
+
+      !! saves 1. MPI interface
+      if( num_interfaces_ext_mesh >= 1 ) then
+        filename = prname(1:len_trim(prname))//'MPI_1_points'
+        call write_VTK_data_points(nglob, xstore_dummy,ystore_dummy,zstore_dummy, &
+                                   ibool_interfaces_ext_mesh_dummy(1:nibool_interfaces_ext_mesh(1),1), &
+                                   nibool_interfaces_ext_mesh(1), filename)
+      endif
+
+    endif
+    !! C. DUROCHAT modification : end
+
     ! acoustic-poroelastic domains
     if( ACOUSTIC_SIMULATION .and. POROELASTIC_SIMULATION ) then
       ! saves points on acoustic-poroelastic coupling interface