Hello again
I have two 3D simulations, one control, with material in non-cohesive suspension, grain size of 0.0001m and SETTLING VELOCITY of 0.0001 m/s. The other simulation is the same but with the inclusion of an external file with time-varying SETTLING VELOCITY. In this I include the VITCHU subroutine and another subroutine that I created (VELDEPO) to receive the values of the time-varying SETTLING VELOCITY (WCHU), so that the VITCHU subroutine calls this variable. The problem is that this simulation exceeds the number of interactions for the three tracers (salt, temperature and sediment).
I'm not sure if the VITCHU subroutine is reading the calculated variable (WCHU in the VELDEPO subroutine) and is being entered correctly.
Attached the part of the subroutine that I include in the fortran file and the external file.
Grateful for the attention