read_dataset overwrites data in case a variable is not found
In the routine read_dataset, in case a variable is not found, data from the previously read variable is used for this variable. This issue was encountered using a PREVIOUS COMPUTATION FILE, with only 2 of the 23 tracers present in TELEMAC-3D. In that situation, the data from the last read tracer was set as initial condition for tracers 3 to 23, instead of using the values from the .cas file. Attached is a modified version of the file that solves the issue. read_dataset.f