Welcome, Guest
Username: Password: Remember me

TOPIC: Prescribe the tracer value of all the surface nodes

Prescribe the tracer value of all the surface nodes 1 year 6 months ago #42614

  • lyybert
  • lyybert's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 24
Could any one tell me how to prescribe the tracer value of all the surface node? I think this should be accomplished by modified the fortran code. I want to locate the paragraph of the code dealing it and the variable storing the value of the tracer concentration.
The administrator has disabled public write access.

Prescribe the tracer value of all the surface nodes 1 year 6 months ago #42618

  • zied.amama
  • zied.amama's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 6
  • Thank you received: 2
Greetings,

if I understand well, you want to prescribe the initial value of a specific tracer on the top layer of your TELEMAC-3D calculation.

To do so, I would suggest you to use the USER_CONDI3D_TRAC fortran file to setup your initial conditions.

You could look into the /examples/waqtel/waq3d_micropol/user_fortran_sed/ user_condi3d_trac.f file where an initial concentration of 1.0 was set to the top layer with the following paragraph:


      DO I=1,NPOIN2
        TA%ADR(1)%P%R(NPOIN3-NPOIN2+I)=1.D0
      ENDDO


Zied
The administrator has disabled public write access.

Prescribe the tracer value of all the surface nodes 1 year 6 months ago #42619

  • lyybert
  • lyybert's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 24
Dear Zied,
Thank you for your reply. I think I did not express clearly. I meant that I need to set the value of tracer concentration on the surface every time step. could you give me a hand?
The administrator has disabled public write access.

Prescribe the tracer value of all the surface nodes 1 year 6 months ago #42630

  • zied.amama
  • zied.amama's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 6
  • Thank you received: 2
Unfortunately I am not sure about the way to do it at all time steps.

You could maybe look at the source_trac.f or user_bord3d.f files to impose source terms at the top layer, but I am not sure if this is the best way to go.

Sorry I cannot be of further help,

Zied
The administrator has disabled public write access.

Prescribe the tracer value of all the surface nodes 1 year 5 months ago #42655

  • lyybert
  • lyybert's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 24
Zied,

Thank you for your help. I think Jean-Paul may have given a direction. Thank you for all of you.

Bert
The administrator has disabled public write access.

Prescribe the tracer value of all the surface nodes 1 year 6 months ago #42631

  • jtravert
  • jtravert's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 36
  • Thank you received: 23
Good afternoon,
I guess that it is possible to do it with Fortran subroutines as described by Zied.
Another way, if you find it too complicated with the subroutines might be to use TelApy python wrapper to interact with the code through Application Programming Interface. If you are not familiar with it you can found the user manual and some examples in the notebook.

You can directly change the state variables or retrieve them during the computation.

I would follow these steps:
1) Retrieve all the nodes number from the surface
2) Run the simulation with a loop one time-step at a time
3) At each time step set the tracer value for all nodes at the free-surface (look into set set_array functions or their parallel implementation mpi_set, mpi_set_array)

Good luck,

Jean-Paul
The administrator has disabled public write access.
The following user(s) said Thank You: zied.amama

Prescribe the tracer value of all the surface nodes 1 year 5 months ago #42654

  • lyybert
  • lyybert's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 24
Jean-Paul,
Thank you for your reply. Using the API sounds good. I am struggling in this way.
The administrator has disabled public write access.

Prescribe the tracer value of all the surface nodes 1 year 5 months ago #42642

  • pham
  • pham's Avatar
  • OFFLINE
  • Administrator
  • Posts: 1559
  • Thank you received: 602
Hello Bert,

You can read the following topic:
www.opentelemac.org/index.php/kunena/21-...g-in-time-telemac-3d

Anyway, you may encounter the same issue as Juli. Prescribing a Dirichlet condition for surface temperature may be too stiff.

Hope this helps,

Chi-Tuan
The administrator has disabled public write access.

Prescribe the tracer value of all the surface nodes 1 year 5 months ago #42656

  • lyybert
  • lyybert's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 24
Chi-Tuan,

I think the post is very useful. Thank you very much.

Best wishes,

Bert
The administrator has disabled public write access.
Moderators: pham

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.