open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello All,
I am trying to define salinity and temperature profile by modifying(condim_f) for the specific area of the domain. Without considering the modified condim file the hydrodynamic works well. However, when using the modified fortran file it calculates the initial mass of water in the domain as infinity and shows exceeding maximum iterations. I attached here my CAS file, the modified condim subroutine also screen shots of the error.
your help is very much appreciated!

Hello Sanaz,

I think this is not that you expected to do:
But fix this probably does not solve your issue...
Best regards,
J

Hello Jacques,
I have the temperature and salinity profile in the water column as below?
Surface Temp C Salinity
1m 12.38 13.37
2m 12.56 19.88
3m 11.00 20.57
4m 10.85 20.98
5m 10.79 20.95

I have seen few posts how to implement temperature value in the water column by layer but not depth.
I tried also this one but still giving me error.Sorry I have a limited Fortran knowledge and will be thankful of you could help me where to start.


Thanks,
Sanaz

Dear Sanaz,

Your fortran sounds good except DEPTH is not an INTEGER but a DOUBLE PRECISION. But this does not fix your issue.

As Chi-Tuân suggest you and as for the steering file, you can simplify your condim (only 2 layers, no zones...) in order to debug your case.

Best regards,

Hello,

There are 2 messages in error.jpg that should be investigated first:
- "INITIAL MASS OF WATER IN THE DOMAIN : Infinity"!
- "TPXO : WARNING:CONSTITUENT ID s1 IS NOT ALLOWED".

You seem to use the new TPXO V9 global solution. Contrary to previous versions, it has 2 extra constituents (2N2+S1), I have just read. The problem is that S1 has not been fully expected by TELEMAC yet.

To see if it is your main problem, can you download an older version of this tidal atlas or a regional V8 or V7, or even The regional solution Atlantic Ocean volkov.oce.orst.edu/tides/AO.html as you seem to work on this area.

If it still fails, I would suggest you to strongly simplify your STEERING FILE, e.g. by starting with the less parameters in your STEERING FILE or the tide example and after, write other parameters step by step.

Regards,

Chi-Tuan

Dear Sanaz,

While discussing with Chi-Tuân... there is one issue in your condim!
You erase the T3D default initializations. You need to add your initialization conditions in addition to the default initialisation of T3D (H, U, V, W, TRAC...)

That is why in the latest 8.0 version we split condim in 5 subroutines.

If you use 8.0 version try to replace the name of subroutine CONDIM with CONDI3DTRAC (this subroutine erase only initialization for tracers, basically that is you want to do).

If you use the trunk (development) version, try to replace the name of subroutine CONDIM by USER_CONDI3D_TRAC.

Best regards,
J

Thanks a mil.Problem solved!