open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,
I noticed some weird behaviour in wind interpolation while working on cases with wind variable in space and time:

• in idwm_t2d.f, if a mesh point is at the same ordinate as a station, or coincident with it, interpolated values are not close to the value at the station
• in meteo.f, wind values are read only if the current time exactly coincides with the times in the ascii data file

I tried to fix these problems. Now wind data is interpolated in time between consecutive lines in the ascii file, and there are no weird spots close to stations.
I also added a column in the ascii file to deal with time-varying, constant in space temperature, to be treated with option 1 (linearized formula) for heat exchange.

This is by no means a general improvement: it is just how I needed my case to work, and I thought it may be useful for someone else. I will attach the modified sources and an example ascii file. Sorry if this is not the right place to post this kind of things!


Have a nice day,
Marco

I found some errors (array allocation/referencing) in the codes I uploaded, they are now fixed.