open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi

I have a question about Subroutine BIL3D, in the calculation of the Flux by diffusion of tracers you have the following:

BOTTOM AND FREE SURFACE
00282 !
00283 IF(ATABOF%ADR(ITRAC)%P%TYPR.NE.'0') THEN
00284 ! WITH MASS-LUMPING LIKE IN DIFF3D
00285 DO I=1,NPOIN2
00286 FLUDI(5+ITRAC) = FLUDI(5+ITRAC)
00287 & + ATABOF%ADR(ITRAC)%P%R(I)*VOLU2D%R(I)
00288 & *TA%ADR(ITRAC)%P%R(I)
00289 ENDDO
00290 ENDIF

I don't understand what is the physical meaning of the variable VOLU2D. I tried to look for the subroutine that calculates this variable in order to understand it and get its units, but I was not able to find it

Thanks

Viviana

I can't give u a direct answer but perhaps this helps as well:

I usually like to start looking here:
docs.opentelemac.org/briefs.html
to find out about the subroutine, (I didn't look into this one before).

Then I lokk at the code, and especially at your VOLU2D which seems to be known by the subroutine without beeing explicitly given.

So I assume it will be known by the telemac3D routine.
Looking docs.opentelemac.org/telemac3d_8f_source.html (or in the files on the pc) I find it's used by 3 subroutines, all having mass as a common theme.

MASBAS2D, ASSED and SED3D.

In the first there is some description :
VOLU2D |<->| INTEGRAL OF TEST FUNCTIONS, WITH MASKING
and
VOLU2D : VOLUME WITH POSSIBLE MASKING
second and third have nearly the same explanation
VOLU2D |-->| INTEGRAL OF TEST FUNCTIONS IN 2D

So I would assume it's some volume and that since in MABAS2D a routine called vectors is called the values are set there.

It might be used and changed elsewhere but thats what can be found the way I usually start searching.
Maybe that's allready sufficient for you to understand bil3d better (?).

And usually SI units are used so m³ would usually be the value for a real volume and I assume (but I'am only a beginner to telemac ) it is one.

Kind regards
csee

oops since it is 2D it's m² .

Thank you very much, usually I also check in the source code, to try to find the variables. I found the vector subroutine where this variable is calculated and it seems the units are m2. Now, the problem is that the units of the flux (FLUDI) are very strange. I was hopping to have W/m2 as units. If someone could help clarifying the units of this flux, I would be grateful.

Thanks

Viviana