Welcome to the forums. Please post in English or French.

You are not logged in.

#1 2021-12-28 04:17:23

Yutakachan
Member
Registered: 2018-01-30
Posts: 12

[Q] How I can know number of "points" in an element?

Hello, everyone.

Recently I'm struggling with writing the python code to "slice" the (DKT) shell ELGA field data to output purpose.
In other words, I'd like to do POST_CHAMP(EXTR_COQUE) on the cham_elem by myself.

In document D4.06.05, I found that the data of cham_elem.CELV are ordered as

CELV(ligrel) = suite de CELV(GREL) mis bout à bout
CELV(GREL) = suite de CELV(élément) mis bout à bout
CELV(élément) = suite de CELV(point) mis bout à bout
CELV(point) = suite de CELV(sous-point) mis bout à bout
CELV(sous-point) = suite de CMP (scalaire) mises bout à bout

and I can know
* start address (index) of the data for an element
* total data size of that element
* number of sub-points in that element
from cham_elem.CELD.
So if I can know the number of components in that field data or how many "points" are in that element, I can extract the data for each sub-points, say, INF/MOY/SUP.

But number of components are not set in cham_elem.CELD unless the data type is VARI.
Document says
"Le nombre de points de discrétisation (noeud, point de Gauss, ...) d'un cham_elem sur une
maille est déterminé a priori par le nombre de points défini dans le catalogue du type_elem
associé à la maille (voir .NOLI(1) et .NOLI(2))."
but it looks that cham_elem does not have such objects (NOLI).

Could you please give me some instruction?

Thanks in advance and Best Regards,
Yutaka Nishizawa

Offline