Welcome, Guest
Username: Password: Remember me

TOPIC: Extracting node data for certain layers

Extracting node data for certain layers 6 years 8 months ago #29078

  • SDAC
  • SDAC's Avatar
Hello,

I want to extract the node data for a specific layers. Is that possible, and is there a feature that readily allows this?

If not, does Telemac identify which nodes belong to which layers? If so it would be simple to extract the needed nodes however I've not seen such a thing.

Many thanks!
The administrator has disabled public write access.

Extracting node data for certain layers 6 years 8 months ago #29106

  • Leballeur
  • Leballeur's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 434
  • Thank you received: 163
Hello,

Several tools can allow you to do that:
- python scripts provided with distribution (v7p3r0/scripts/python27/parsers/parserSELAFIN.py)
- pytel tools: github.com/CNR-Engineering/PyTelTools
- matlab tools: fr.mathworks.com/matlabcentral/fileexchange/25021-telemac-tools
- BlueKenue

To identify the node you want to extract, you can first find its number IPOIN2 in the 2D mesh (geometry file for example of size NPOIN2) and then the number in a given layer IPLAN will be provided by:
IPOIN3 = (IPLAN-1)*NPOIN2 + IPOIN2

Regards,
Laurent
The administrator has disabled public write access.
The following user(s) said Thank You: SDAC

Extracting node data for certain layers 6 years 8 months ago #29157

  • SDAC
  • SDAC's Avatar
Thankyou for your response.

To extract an entire layer, I assume I would have to identify every node in that particular layer? I want to remove certain layers (e.g. the boundary and layers close to the boundary) so I can analyse data only for a certain portion of the channel.

Best regards
The administrator has disabled public write access.

Extracting node data for certain layers 6 years 8 months ago #29162

  • Leballeur
  • Leballeur's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 434
  • Thank you received: 163
Hi,

All the nodes that belong to the layer IPLAN are given by the same formula for IPOIN2 varying between 1 and NPOIN2.
Here, layers are horizontal plans while boundaries are vertical.
The global numbers of your 2D boundary points are provided by the 12th column of your boundary conditions file: NBOR2(IPTFR) (of size NPTFR2). If you replace IPOIN2 by this numbers in the previous formula, you'll have it for each layer IPLAN.

More generally, if you want to select a portion of your mesh, you have first to locate the points enclosed in the suited area, and then apply the formula to these points.

Regards,
Laurent
The administrator has disabled public write access.
The following user(s) said Thank You: SDAC
Moderators: pham

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