open TELEMAC-MASCARET The mathematically superior suite of solvers

I am using the condin subroutine to add initial salinity variable in the space, this works good without continued computations. But when I use continued computations Telemac uses zero salinity as initial condition which is very strange.

The fortran file is attached, can one of you help me in this problem. The condin has part also for initial H and initial velocity, but normally for continued computations the H and V for the last time step of the previous file is taken, so I left this part in the condin file because I know that telemac will ignore it, can this the source of the problem.

I am trying to model the salinity transport in an estuary, I set my boundary conditions as imposed flow from the side of the river and an imposed head with imposed salinity at the side of the sea. the initial conditions was constant head and constant salinity ( zero at the upstream at the river and 38 at the sea). When I tried to model this case the computations usually stop at certain time as the salinity at the sea increases usually and at the river remains zero.

If someone studied such case, can he or she tell me how he or she did it, and what is the most suitable configurations for such case because I fell that what I am doing is wrong and there is something missing.

Hello,

Normally, the salinity should be in the results file and it is read to continue the computation. So in case of computation continued condin is not called (actually it is called to get some initial values as the depth that is necessary for Boussinesq equations) but then velocities, depth and everything is erased by what is read in the previous computation file.

Could that be your problem ?

With best regards,

Jean-Michel Hervouet

yeah this is my problem. But is you know normally we make the hydrodynamic simulation until reaching a steady state, then we use this file as initial condition for the transport simulation. But is in the hydrodynamic simulation there is no tracer, that means that one must usually start with a zero tracer or just fixed number.

Hello,

OK, so take subroutine rescue.f from telemac-2d library into your FORTRAN FILE. You will find in it a sequence of commands that checks if tracers have been found in the results file and that puts their initial value if not. Just erase all this section and your values set in condin.f will be kept.

With best regards,

Jean-Michel Hervouet

ok, that works very good, but only when I merged the condins subroutine with the rescue subroutine in one subroutine.

thank you very much