🐛 [BUG] - Conflict of file unit of `IO_STF` and `IIN_veloc_dsm`.

[xumi1993]

Description

Hi all, I find that the file unit of IO_STF and IIN_veloc_dsm in the constant.h are the same. I am now try to run specfem3D using INJECTION_TECHNIQUE_TYPE = 2 by reading injection wavefield from proc******_sol_axisem with opening IIN_veloc_dsm. but the file pointer on rank 0 has been closed at

specfem3d/src/specfem3D/print_stf_file.f90

Line 195 in 7a1c764

close(IO_STF)

before

specfem3d/src/specfem3D/compute_stacey_viscoelastic.F90

Line 801 in 7a1c764

read(IIN_veloc_dsm) Veloc_axisem, Tract_axisem

Affected SPECFEM3D version

v4.1.1

Your software and hardware environment

CPU on workstation

Reproduction steps

1. prepare `proc******_sol_axisem`
2. run `specfem3D`

Screenshots

Logs

OS

Linux

[xumi1993]

Solved by #1794