open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everybody,

I am working with the t2d_micropol example to try to undestard how it works.

Currently I am trying to input some discharges in open boundaries to have velocities in the domain and see how it behaves.

I had no problem using open boundaries with free tracer but when I used open boundaries with prescribed tracer, the user defined tracers and the Suspended Matter ones work but the Micropollutant tracers become NaN in the first timestep.

Does anybody know what the problem is?

I'm attaching my files, I hope they upload correctly.

Thanks

Santiago

Hello Santiago,

Thanks for sending your problem. It was a bug in calcs2d_micropol with an unexpected division by 0 in your case.

With the enclosed file, it should fix your problem. It will be integrated in the next release.

The changed lines are:
CALL OVD('X=CXY/Z ' ,T1%R,TN%ADR(IND_CSF)%P%R,TN%ADR(IND_SF)%P%R,
& 1.D0, NPOIN,2,0.D0,EPS)

rather than:
CALL OS( 'X=CXY/Z ' ,X=T1,Y=TN%ADR(IND_CSF)%P,Z=TN%ADR(IND_SF)%P,
& C=1.D0)

Please tell me if you are OK now.

Chi-Tuan

Hello Chi-Tuan

Sorry for the late reply.

Thank you for the fortran file, it seems to be working now. I will let you know if I have another problem.

Santiago