unified api for couple_with_injection

Author: nqdu

src/gpu/transfer_fields_cuda.cu

@@ -324,6 +324,29 @@ void FC_FUNC_(transfer_wavefield_discontinuity_to_device,
 
 }
 
+/**
+ * @author nqdu
+ * @brief copy injection field from host to device
+ * @param size_face no. of points used = NDIM*num_abs_boundary_faces*NGLLSQUARE
+ * @param veloc/tract_inj velocity/traction
+ * @param Mesh_pointer mesh_ struct
+ */
+extern EXTERN_LANG
+void FC_FUNC_(transfer_injection_field_to_device,
+              TRANSFER_INJECTION_FIELD_TO_DEVICE)(
+                         int* size_face,
+                         realw* veloc_inj,realw* traction_inj,
+                         long* Mesh_pointer) {
+
+  TRACE("transfer_injection_field_to_device");
+
+  Mesh* mp = (Mesh*)(*Mesh_pointer); //get mesh pointer out of fortran integer container
+
+  gpuMemcpy_todevice_realw(mp->d_veloc_inj, veloc_inj, (*size_face));
+  gpuMemcpy_todevice_realw(mp->d_tract_inj, traction_inj, (*size_face));
+
+}
+
 /* ----------------------------------------------------------------------------------------------- */
 
 extern EXTERN_LANG

src/specfem3D/iterate_time.F90

@@ -274,6 +274,8 @@ subroutine iterate_time()
           call update_displ_Newmark()
         endif
 
+        if(COUPLE_WITH_INJECTION_TECHNIQUE) call fetch_injection_wavefield()
+
         ! computes acoustic domain (first)
         if (ACOUSTIC_SIMULATION) call compute_forces_acoustic_forward_calling()
         ! computes elastic domain

src/specfem3D/prepare_gpu.f90

@@ -100,6 +100,7 @@ subroutine prepare_GPU()
                                 NB_RUNS_ACOUSTIC_GPU, &
                                 FAULT_SIMULATION, &
                                 IS_WAVEFIELD_DISCONTINUITY, &
+                                COUPLE_WITH_INJECTION_TECHNIQUE,&
                                 UNDO_ATTENUATION_AND_OR_PML, &
                                 PML_CONDITIONS)