Merge pull request #1794 from xumi1993/devel

Solve some bugs on file IO

Author: danielpeter

setup/constants.h.in

@@ -285,8 +285,8 @@
 
 ! some old tests (currently unstable; do not remove them though, we might fix this one day)
   integer, parameter :: Ntime_step_dsm = 100
-  integer, parameter :: IIN_veloc_dsm = 51, IIN_tract_dsm = 52
-  integer, parameter :: IIN_displ_axisem = 54
+  integer, parameter :: IIN_veloc_dsm = 59, IIN_tract_dsm = 60
+  integer, parameter :: IIN_displ_axisem = 61
 
 !!-----------------------------------------------------------
 !!

src/specfem3D/couple_with_injection.f90

@@ -1955,6 +1955,7 @@ subroutine ReadFKModelInput(Xmin_box, Xmax_box, Ymin_box, Ymax_box, Zmin_box, Zm
      end select
   !!------------------------------------------------------------------------------------------------------
   enddo
+  close(85)
 
   if (allocated(ilayer_fk_input)) then
 
@@ -4381,4 +4382,4 @@ subroutine read_specfem_file(Veloc_specfem,Tract_specfem,nb,it)
     endif
 
     end subroutine read_specfem_file
-  
+