TELEMAC-MASCARET Forum :: Topic: Vertical Tracer Profile (2/2)

Welcome, Guest

TOPIC: Vertical Tracer Profile

Vertical Tracer Profile 1 month 3 weeks ago #44454

Raphaël OFFLINE Fresh Boarder Posts: 20 Thank you received: 1	Thank you, it finally works with both user_vel_prof_z and user_tra_prof_z at the same time. Raphaël
	The administrator has disabled public write access.

Vertical Tracer Profile 1 month 3 weeks ago #44458

Raphaël
OFFLINE
Fresh Boarder
Posts: 20
Thank you received: 1

Hello

I finally have a problem on how I want to impose a flow of tracers at the border of my artificial island.
I have a liquid boundaries file to impose, as a function of time, my flow rates on the upstream boundary (boundary number 2) but also at the level of the island (boundary number 3), with associated concentration of my tracers (I have 25 tracers, those of the waqtel aed2 module, because I work with a telemac3d/aed2 coupling). Above my test file for liquid boundaries:

File Attachment:

File Name: liquid_boundaries_file_ile.txt
File Size: 5 KB

When I use de user_vel_prof_z only (without user_tra_prof_z.f) my liquid boundaries file is read in terms of flow variation and tracer appearance. At this stage, my tracers are injected on each of my layers, which is normal since I don't use user_tra_prof_z.f.
But when I add user_tra_prof_z.f (on its own or coupled with user_vel_prof_z.f), my tracers are indeed injected at the indicated depth, but they are injected at T0 and continuously throughout the simulation (which is no indicated in the liquid boundaries file). Furthermore, their concentrations are equal to those indicated in user_tra_prof_z.f and not to those indicated in my liquid boundaries file. In other words, the liquid boundaries files is not read for the tracers. But in the other hand, current speeds do vary with the variations indicated in the liquid bondaries file...

Thanks again for your help

Raphaël

The administrator has disabled public write access.

Vertical Tracer Profile 1 month 1 week ago #44500

pham
OFFLINE
Administrator
Posts: 1468
Thank you received: 563

Hello Raphaël,

You can define a tracers vertical profile in USER_TRA_PROF_Z, but if you have a look at BORD3D subroutine which calls TRA_PROF_Z function (which finally may call USER_TRA_PROF_Z), you can understand that what you can see is exactly what is currently implemented:

                TABORL%ADR(ITRAC)%P%R(IBORD)=PROFZ

I have to investigate what would be the best fix to deal with other choices for TRACERS VERTICAL PROFILES but currently, if you want to have PROFZ times what is inside LIQUID BOUNDARIES FILE, I think you can change the conditions:

              IF(IPROF.EQ.2.OR.IPROF.EQ.4.OR.IPROF.EQ.0) THEN
!               Rouse concentrations profiles (IPROF=2 or 4) or given by user (IPROF=0)
                TABORL%ADR(ITRAC)%P%R(IBORD)=PROFZ
              ELSEIF(IPROF.EQ.3) THEN
!               Normalised concentrations profiles (IPROF=3)
                TABORL%ADR(ITRAC)%P%R(IBORD)=
     &          TABORL%ADR(ITRAC)%P%R(IBORD)*PROFZ

with

              IF(IPROF.EQ.2.OR.IPROF.EQ.4) THEN
!               Rouse concentrations profiles (IPROF=2 or 4) or given by user (IPROF=0)
                TABORL%ADR(ITRAC)%P%R(IBORD)=PROFZ
              ELSEIF(IPROF.EQ.3.OR.IPROF.EQ.0) THEN
!               Normalised concentrations profiles (IPROF=3)
                TABORL%ADR(ITRAC)%P%R(IBORD)=
     &          TABORL%ADR(ITRAC)%P%R(IBORD)*PROFZ

see the test for IPROF.EQ.0

Hope this helps,

Chi-Tuan

The administrator has disabled public write access.

Vertical Tracer Profile 1 month 5 days ago #44529

Raphaël OFFLINE Fresh Boarder Posts: 20 Thank you received: 1	Hello Thanks a lot, it works ! Raphaël
	The administrator has disabled public write access. The following user(s) said Thank You: pham