open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear Jean-Michel,

Thank you very much for looking into it. I have tried to add modified bord3d.f into my fortran file but it gave me errors all the time I run it. I attached the output with error message. So, silly me, I thought that files might be connected and I have changed it in bord3d.f in Source directory without adding it, and as a result it did not work of course, I have just realized it after your note.

So modified included bord3d does not run. It could be some wrong parameters used in steering file, but error message does not indicate so. I have attached all case files. If you could be so kind to take a look when you get a moment, I would be very grateful. I have to combine this case with tidal run, which works fine without stratification, but first I have to sort the open boundary conditions. Thanks a mil.

Kind Regards!

Violeta

Hello Violeta,

I think I got something. You should add in your parameter file :

PRESCRIBED TRACERS VALUES : 0.D0

The 0.D0 will be replaced by your profile but it will trigger the call to TRA_PROF_Z which was not done. In bord3d.f this call is done only if NTRACER.GT.0 (line 387) and NTRACER is the number of prescribed tracers values given. If you do not give any value in the parameter file it is understood that your values are given in the boundary conditions file, where you have zeroes. Sorry for this, there are many possibilities piling up. Putting print messages in bord3d.f you would have seen that your function tra_prof_z was not called. I hope it is the solution...

With best regards,

Jean-Michel Hervouet

Thank you very much Jean-Michel. I will do so. I was afraid to use that option before as I have two values for one tracer but in that option I have to use only one and I did not think that it was right. But your solution could work, just setting it to 0. Thanks again. I will keep you posted.

Kind Regards!

Violeta

Dear Jean-Michel,

Unfortunately PRESCRIBED TRACERS VALUES : 0.D0 did not solve the problem. It does not run, I still get same error. I thought that maybe it is problems with parallelism and tried to run on serial, but got same result. Is there are anything else I can try? I am ready to experiment until it is sorted, no matter how long it will take. Thank you for your help, it is much appreciated.

Kind Regards!

Violeta

Hello,

I have no other idea. You can put prints in your subroutine TRA_PROF_Z and before the CALL TRA_PROF_Z in bord3d.f to see if we really get there, if it does not clear up things then please send me all the necessary files so that I can run your case myself.

Regards,

JMH

Thank you Jean-Michel. I am probably wrong, but I cannot find CALL TRA_PROF_Z in bord3d.f. I have tried to add it myself but got an error, so I probably did not put it in right place or not enough coding. If it is missing so maybe it is the reason why it does not apply, I hope my guess is right. Please advise.

Thanks.

Kind Regards!

Violeta

Hello,

Excuse me, it is a function, so no CALL, it is just the line :

PROFZ=TRA_PROF_Z(......

JMH

Dear Jean-Michel,

Thank you for your suggestions, the prints did not appear in my output so I am not sure if those subroutines were applied at all. I have attached all files I use for this case. I would be really grateful if you could try to run them, only when you have some time for this. Thank you very much, your help is very appreciated.

Kind Regards!

Violeta

I had to change the cas file extension otherwise it cannot be attached. Here it is. Thanks again.

Kind Regards!

Violeta

Sorry, this time it is attached.

Vi