open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello!

I am trying to implement a boundary condition with a temperature and salinity profile. I have read to Forum topic by Violeta and JMH on the same issue, and I have tried to use the information there without success.

The only thing I get is that temperature is entering the boundaries with the value 0. I have attached all the files.

Also I want to have different profiles on the two boundaries - how do I implemente this?

Thankful for all help!

Charlotta

Hello Charlotta,

You first need to add the keyword PRESCRIBED TRACERS VALUES in your steering file, giving 4 values (2 boundaries x 2 tracers). Without this keyword the values are taken in the boundary conditions file (tricky I admit...). Then the 4 values will be overwritten by your programmed profiles.

Then to have different profiles on the two boundaries you can use the variable I of function TRA_PROF_Z and nest your two profiles within:

IF(I.EQ.1) THEN
... your profile for boundary 1
ELSEIF(I.EQ.2) THEN
... your profile for boundary 2
ELSE
WRITE(LU,*) 'Hum, boundary ',I,' is not implemented in TRA_PROF_Z'
CALL PLANTE(1)
STOP
ENDIF

With best regards,

Jean-Michel Hervouet

Thanks a lot Jean-Michel!

I finally got it working with two different tracers (temperature and salinity) and different profiles on two different boundaries.

Best regards
Charlotta

And the next step is now to implement time varying profiles on the two different boundaries - how do I do that...????

Thanks,
Charlotta

Hello,

Congratulations! Then it depends on the format of your data, but basically, where you have written your profiles depending on space you can also add a dependency on time. Namely in tra_prof_z.f you have the variable TIME in the list of arguments, so there is no fundamental difference.

With best regards,

Jean-Michel