open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello

I have been trying to include the wind and air temperature in a 3D simulation. So I have two text files from where I am reading the data. I used subroutine METEO and I created a new subroutine to read the air temperature using METEO as an example for reading a file. However, the simulation is stopping and in the output file it says:

T_AIR
T_AIR FILE TOO SHORT

Which corresponds to the case when it finishes reading the file. However, in the output file I just see the listing for the wind values but not for the air temperatures. I am attaching my fortran file where I added subroutine Thermal at the end, also in lines 1062 to 1064 I am calling this subroutine from the BORD3D subroutine. I am attaching my case file and the formatted data files and output file too.

Any help is greatly appreciated

Thanks so much

Viviana

Hello,

Difficult to say. You should print AT, AT1 and AT2 after label 10 to see what happens. A possible mistake is the test AT.LT.AT2, in case AT is exactly AT2, maybe AT.LE.AT2 will work, or you could choose a larger final AT2 in the last line of your file.

With best regards,

Jean-Michel Hervouet

Hi

Thanks so much for your quick reply. I was able to fix the problem, for some reason in the subroutine that reads the air temperature the line:
IF(LT.EQ.0) THEN
reads the file

Was not working so it was never really reading the file. So I changed to :

IF(LT.EQ.1) THEN
reads the file
And it works now. In the METEO subroutine the IF(LT.EQ.0) works, but not in the new subroutine. So that was the problem. Also I made sure that my input files has at least 2 lines with AT larger than my total simulation time.

Thanks again

Viviana