open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I see "Q for flowrate, SL for water level, U and V for velocities and TR for tracer" in
the user manual(v7p0_t2d page 32-33).

The FUNCTION VIT is for scalar velocity without direction? But I want to prescribe the rotary current at the LIQUID BOUNDARIES.

So how could I prescribe U and V for velocities that vary in time?

Thank you in advance!

Best regards

HJW

Hello,

Actually you can do what you want in subroutine bord.f. Though it is complicated because it deals with many possibilities, its only role is to give the prescribed U, V and H. Then it depends where you get your data and there could be a little programming to read them.

With best regards,

Jean-Michel Hervouet

Hi Jean-Michel Hervouet

Could I do it the same as function sl(q and tr ) ?
Is there an example about prescribe U and V(vary in time)at the liquid boundaries?

Thank you in advance!

Best regards

HJW

HI HJW

Can you say a bit more concrete about where do you get U and V varying in time.

you need to modify bord.f but if the case is for nesting, i have recently done so (modifying bord.f in a way that reads a coarser grid telemac run to use in a finer grid). So mine might help you.

Bests

Homayoon

Hi Homayoon

Sure.
First,the picture is my t2d hydrodynamic computation area.And the result is good.



Second, I want to add the sedimentary compution only in Zone A.
So I make a new geo file for zone A.
I want to take the first computation's velocity(Both U and V component ) result as the second computation's liquid boundary condition.
But I don't know how to prescribe U and V at the liquid boundaries.

Thank you in advance!

Best regards

HJW

Hi

You can start from the attached bord.f. if you run first the bigger grid and connect the results of bigger area to "binary file 2" channel when running zoneA, this bord.f will read it and for each boundary point finds the closest data point in the results file of the bigger grid and use it.

i think you will still need to check and modify the file regarding :

1- i have used this subroutine for spherical coordinates, as your coor. is apparently UTM, some changes might be needed although i do not think so.
2- this routine uses time from beginning 0. i dont know how you deal with time.
3- i use this data on all boundary points but probably you need to change something if you want them only on liquid boundary point.
4- i use free surface and U and V.

having a look at the file will make things clear i guess.

cheers
Homayoon

Hi

Thank you very much for your suggestion.

I have solved my problem by modifying FUNCTION VIT and SUBROUTINE READ_FIC_FRLIQ.
I dont think it is a good idea although my trick works.

I want to add functions(the same as FUNCTION SL ) for U(FUNCTION UHJ) and V(FUNCTION VHJ),
I denfined them in the MODULE INTERFACE_TELEMAC2D and wrote their codes in the FORTRAN FILE.

The error : Return type mismatch of function "uhj"
& function "uhj" has no implicit type

How could I add my funtions ?

Thank you in advance!

Best regards

HJW

Hello,

There is no need to add your functions in the interface. Just put them in your Fortran file and in the subroutine where you use them just add :

DOUBLE PRECISION UHJ,VHJ
EXTERNAL UHJ,VHJ

With best regards,

Jean-Michel Hervouet

Hi

It's OK now!

Thank you very much for your suggestion!

Best regards

HJW