open TELEMAC-MASCARET The mathematically superior suite of solvers

Recently I was trying to edit USER_CONDIN_TRAC to assign differential values to the drawn grid, but I didn't know the size of the grid node (NPOIN) value. I tried to find the cli file, but I still didn't understand it and couldn't find the answer in the manual.I wish someone could give me some answers.

Hello Xinjie,

Do you mean USER_CONDIN_TRAC for TELEMAC-2D or USER_CONDI3D_TRAC for TELEMAC-3D as you post a question for USER_CONDIN_TRAC in the TELEMAC-3D section?

In 2D, NPOIN is the total number of grid nodes in the subdomain (= total number of the mesh in scalar mode, total number of the partitioned mesh of the subdomain in parallel). You do not have to find the value, it is known in the DECLARATIONS_TELEMAC2D module.

If you want to know more how to implement in TELEMAC, I would recommend you to read the developer guide available here e.g.:
wiki.opentelemac.org/doku.php?id=documentation_v8p4r0
see the Miscellaneous Manuals section.

Hope this helps,

Chi-Tuan

Hello again,

Maybe some examples of use of USER_CONDIN_TRAC may help you. You can find such examples in:
- gaia/mud_conservation-t2d/user_fortran/user_condin_trac.f
- telemac2d/cone/user_fortran/user_condin_trac.f
- telemac2d/cone_diffusion/user_fortran/user_condin_trac.f
- telemac2d/convergence/user_fortran/user_condin_trac.f
- telemac2d/gouttedo/user_fortran/user_condin_trac.f
- telemac2d/flume_tracer/user_fortran/user_condin_trac.f
- telemac2d/siphon/user_fortran/user_condin_trac.f
- telemac2d/tracer_diffusion/user_fortran/user_condin_trac.f
- waqtel/tracer_decay/user_fortran/user_condin_trac.f

Hope this helps,

Chi-Tuan

Hello again!
I have also tried to look at the examples you gave, but they do not seem to solve my problem.
I tried to build a grid myself and read its cli file, but now it makes me even more confused. I know that the two columns on the far right represent the order and the node values, but they don't seem to be related to each other, and the node values don't seem to have any pattern, so how do I get the node values in the circle in Figure 2? When I generalize to an n×n grid, how do I get the values of the nodes in the grid?
Or something else, what else can you do to assign values of different coordinates to the corresponding grid points?
Best regard

Hello!
Thanks for your reply, I'm just giving you an example because I know that 3d and 2d have a lot in common. I may not have described it clearly, but let me try to draw it on a graph, which I put in the following file. I was thinking that if I didn't know the table nodes and their ordering rules, it would be difficult to assign values through fortran.
Regard for you!

Hi Xinjie,

I might not understand your question very well, but for a (unstructured / triangular) mesh the data cannot be simplified to a raster.

For me, the easiest way to now the ID of a node is to open the grid in BlueKenue and double click a node to get the ID (IPOIN2).

Hi jurjendejong,
Thank you very much for your reply! The method you mentioned does know the ID of the node on the grid, but in my previous attempt it only knew the ID of the node on the boundary, and it may not be able to get the ID of the node inside the grid. More importantly, I do not know the order of these nodes, because my current work is to add the detailed temperature to the target region, the temperature distribution is shown in Figure 1. If you have other good ways, please give me some suggestions!

Hi,

If I understand correctly you have temperature on a regular (square) grid and you want those values on the triangular mesh.

You can use the Blue Kenue "MapObject" function to assign interpolated values from the regular grid to a mesh. Please refer to the Blue Kenue reference manual for further information.

I hope this helps... Martin

Hi!
Thank you! I'm sorry for not replying to you in time.Due to some reasons, I have been away from work during this period.
I roughly understand what you mean, because I have also used the map object feature in bluekenue to assign values to a grid, maybe it is feasible in 2d conditions, but I can't achieve it on a 3d grid, I think the temperature is the same.
This means that I still need to customize it through fortran.