MED fichier: Fichier source de MEDmeshNodeCoordinateWr.c

MED fichier

MEDmeshNodeCoordinateWr.c
Aller à la documentation de ce fichier.
1 /* This file is part of MED.
2 *
3 * COPYRIGHT (C) 1999 - 2023 EDF R&D, CEA/DEN
4 * MED is free software: you can redistribute it and/or modify
5 * it under the terms of the GNU Lesser General Public License as published by
6 * the Free Software Foundation, either version 3 of the License, or
7 * (at your option) any later version.
8 *
9 * MED is distributed in the hope that it will be useful,
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 * GNU Lesser General Public License for more details.
13 *
14 * You should have received a copy of the GNU Lesser General Public License
15 * along with MED. If not, see <http://www.gnu.org/licenses/>.
16 */
17
18
19 #include <med.h>
20 #include <med_config.h>
21 #include <med_outils.h>
22 #include <string.h>
23 #include <stdlib.h>
24
45 med_err MEDmeshNodeCoordinateWr(const med_idt fid,
46 const char* const meshname,
47 const med_int numdt,
48 const med_int numit,
49 const med_float dt,
50 const med_switch_mode switchmode,
51 const med_int nentity,
52 const med_float* const coordinate)
53 {
54 med_idt _meshid = 0;
55 med_idt _datagroup1=0,_datagroup2=0,_datagroup3=0,_datagroup_trsf=0;
56 med_idt _dataset_trsf = 0;
57 med_err _ret = -1;
58 med_grid_type _gridtype = MED_UNDEF_GRID_TYPE;
59 med_int _intgridtype = 0;
60 med_int _intmeshtype = 0;
61 med_int _meshdim = 0;
62 char _meshpath [MED_MESH_SUPPORT_GRP_SIZE+MED_NAME_SIZE+1]="";
63 char _geotypename[MED_TAILLE_NOM_ENTITE+1]="";
64 char _datagroup_trsfname[MED_TAILLE_NOM_ENTITE+MED_TAILLE_COOTRF+1]=MED_NOM_NOE MED_COOTRF;
65 med_bool _isasupportmesh=MED_FALSE,_isasoftlink=MED_FALSE;
66 med_int _0 = 0;
67 med_int _medintgeotype = MED_NO_GEOTYPE;
68
69 if (_MEDcheckVersion30(fid) < 0) goto ERROR;
70
71 /* Ouverture du maillage de calcul ou du maillage support */
72 if ((_meshid=_MEDmeshDatagroupOpen(fid,meshname,_meshpath,&_isasupportmesh)) < 0) {
73 MED_ERR_(_ret,MED_ERR_OPEN,MED_ERR_DATAGROUP,MED_ERR_MESH_MSG);
74 SSCRUTE(_meshpath); goto ERROR;
75 }
76
77 /* Lecture du type de maillage (attribut MED_NOM_TYP) */
78 if (_MEDattrEntierLire(_meshid,MED_NOM_TYP,&_intmeshtype) < 0) {
79 MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
80 SSCRUTE(meshname);SSCRUTE(MED_NOM_TYP);ISCRUTE(_intmeshtype);goto ERROR;
81 }
82
83 /*Sortir en erreur si maillage support et numdt,numit != MED_NO_PDT,MED_NO_IT*/
84 /* Accès (création éventuelle) du datagroup numdt,numit */
85
86 if ( _isasupportmesh ) {
87 if (( numdt != MED_NO_DT) && (numit != MED_NO_IT) ) {
88 MED_ERR_(_ret,MED_ERR_INVALID,MED_ERR_PARAMETER,MED_ERR_MESH_MSG);
89 SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);goto ERROR;
90 }
91 } else
92 if ( (_datagroup1 =_MEDmeshAssociatedGroupCr(fid, MED_MESH_GRP,
93 meshname, numdt, numit, dt,MED_FALSE,
94 "." ) ) < 0 ) {
95 MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAssociatedGroupCr");
96 goto ERROR;
97 }
98
99 /* Création de la structure de connectivité assoiciée aux grilles MED_CURVILLINEAR_GRID*/
100 if ( ( (med_mesh_type) _intmeshtype ) != MED_UNSTRUCTURED_MESH ) {
101
102 /* Lecture de l'attribut MED_NOM_GTY */
103 if (_MEDattrEntierLire(_meshid,MED_NOM_GTY,&_intgridtype) < 0) {
104 MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
105 SSCRUTE(meshname);SSCRUTE(MED_NOM_GTY);ISCRUTE_int(_gridtype);goto ERROR;
106 }
107 _gridtype= (med_grid_type) _intgridtype;
108
109 if ( _gridtype == MED_CURVILINEAR_GRID ) {
110
111 if ((_datagroup2 = _MEDdatagroupOuvrir(_datagroup1,MED_NOM_MAI)) < 0) {
112
113 if ((_datagroup2 = _MEDdatagroupCreer(_datagroup1,MED_NOM_MAI)) < 0) {
114 MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,MED_NOM_MAI);
115 SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
116 goto ERROR;
117 }
118
119 if ( _MEDattributeIntWr(_datagroup2,MED_NOM_CGT,&_0) < 0) {
120 MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
121 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGT);
122 goto ERROR;
123 }
124
125 /* Lecture de l'attribut MED_NOM_DIM */
126 if (_MEDattrEntierLire(_meshid,MED_NOM_DIM,&_meshdim) < 0) {
127 MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
128 SSCRUTE(meshname);SSCRUTE(MED_NOM_DIM);ISCRUTE(_meshdim);goto ERROR;
129 }
130
131 switch ( _meshdim ) {
132 case 1 :
133 strcpy(_geotypename,MED_NOM_SE2);
134 _medintgeotype = MED_SEG2;
135 break;
136 case 2 :
137 strcpy(_geotypename,MED_NOM_QU4);
138 _medintgeotype = MED_QUAD4;
139 break;
140 case 3 : strcpy(_geotypename,MED_NOM_HE8);
141 break;
142 _medintgeotype = MED_HEXA8;
143 case 0 : strcpy(_geotypename,MED_NOM_PO1);
144 _medintgeotype = MED_POINT1;
145 break;
146 default :
147 MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_PARAMETER,"");
148 SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
149 goto ERROR;
150 }
151
152 if ((_datagroup3 = _MEDdatagroupCreer(_datagroup2,_geotypename)) < 0) {
153 MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,_geotypename);
154 SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
155 goto ERROR;
156 }
157
158 if ( _MEDattributeStringWr(_datagroup3,MED_NOM_PFL,MED_NAME_SIZE,MED_NO_PROFILE_INTERNAL) < 0) {
159 MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
160 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);
161 goto ERROR;
162 }
163
164 if ( _MEDattributeIntWr(_datagroup3,MED_NOM_CGT,&_0) < 0) {
165 MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
166 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGT);
167 goto ERROR;
168 }
169
170 if ( _MEDattributeIntWr(_datagroup3,MED_NOM_CGS,&_0) < 0) {
171 MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
172 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGS);
173 goto ERROR;
174 }
175
176 if (_MEDattributeIntWr(_datagroup3,MED_NOM_GEO,&_medintgeotype) < 0) {
177 MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_NOM_GEO);
178 ISCRUTE(_medintgeotype);
179 goto ERROR;
180 }
181
182 }
183 } else {
184 MED_ERR_(_ret,MED_ERR_INVALID,MED_ERR_GRIDTYPE,MED_ERR_MESH_MSG);
185 SSCRUTE(meshname);SSCRUTE(MED_NOM_GTY);ISCRUTE_int(_gridtype);goto ERROR;
186
187 }
188
189 }
190
191
192 /*Réinitialisation du lien MED_TRSF_NOM/MED_COORDINATE_TRSF s'il existe */
193 /*Si un dataset de transformation existe une erreur est générée par
194 _MEDsoftlinkDel */
195 if ( !_isasupportmesh)
196 if ((_datagroup_trsf = _MEDdatagroupOuvrir(_datagroup1,_datagroup_trsfname)) > 0) {
197
198 if ( _MEDisasoftlink(_datagroup_trsf, MED_NOM_TRF,MED_FALSE, &_isasoftlink ) < 0 ) {
199 MED_ERR_(_ret,MED_ERR_OPEN,MED_ERR_LINK,MED_NOM_COO);
200 goto ERROR;
201 }
202 if (_isasoftlink ) {
203 if ( _MEDsoftlinkDel(_datagroup_trsf,
204 MED_NOM_TRF,
205 MED_TRUE) < 0 ) {
206 MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDsoftlinkDel");
207 goto ERROR;
208 }
209 } else {
210 if ( (_dataset_trsf = _MEDdatasetOuvrir(_datagroup_trsf,MED_NOM_TRF)) >= 0 ) {
211 _MEDdatasetFermer(_dataset_trsf);
212 MED_ERR_(_ret,MED_ERR_NULL,MED_ERR_DATASET,MED_NOM_TRF);
213 goto ERROR;
214 }
215
216 }
217 }
218
219 /*Ecriture des coordonnées*/
220 if ( _MEDmeshAdvancedWr(fid,
221 meshname,
222 MED_COORDINATE,
223 MED_NO_NAME,
224 MED_INTERNAL_UNDEF,
225 numdt,
226 numit,
227 dt,
228 MED_NODE,
229 MED_NONE,
230 MED_NO_CMODE,
231 MED_UNDEF_STMODE,
232 MED_NO_PROFILE,
233 switchmode,
234 MED_ALL_CONSTITUENT,
235 NULL,
236 nentity,
237 coordinate) < 0 ) {
238 MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAdvancedWr");
239 goto ERROR;
240
241 }
242
243
244 _ret = 0;
245
246 ERROR:
247
248 if (_datagroup_trsf>0) if (_MEDdatagroupFermer(_datagroup_trsf) < 0) {
249 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_datagroup_trsfname);
250 ISCRUTE_id(_datagroup_trsf);
251 }
252
253 if (_datagroup3>0) if (_MEDdatagroupFermer(_datagroup3) < 0) {
254 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_geotypename);
255 ISCRUTE_id(_datagroup3);
256 }
257
258 if (_datagroup2>0) if (_MEDdatagroupFermer(_datagroup2) < 0) {
259 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,MED_NOM_MAI);
260 ISCRUTE_id(_datagroup2);
261 }
262
263 if (_datagroup1>0) if (_MEDdatagroupFermer(_datagroup1) < 0) {
264 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_meshpath);
265 ISCRUTE_id(_datagroup1);ISCRUTE(numdt);ISCRUTE(numit);
266 }
267
268 if (_meshid>0) if (_MEDdatagroupFermer(_meshid) < 0) {
269 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,meshname);
270 ISCRUTE_id(_meshid);
271 }
272
273 return _ret;
274 }
MED_ERR_LINK
#define MED_ERR_LINK
Definition: med_err.h:115
MED_ERR_API
#define MED_ERR_API
Definition: med_err.h:111
MED_ERR_CREATE
#define MED_ERR_CREATE
Definition: med_err.h:27
med_grid_type
med_grid_type
Definition: med.h:139
ISCRUTE_id
#define ISCRUTE_id(entier)
Definition: med_utils.h:319
MED_UNDEF_STMODE
Definition: med.h:110
MED_UNSTRUCTURED_MESH
Definition: med.h:133
_MEDattributeIntWr
#define _MEDattributeIntWr(w, x, y)
Definition: med_hdfi.h:45
MED_TRUE
Definition: med.h:262
MED_NOM_DIM
#define MED_NOM_DIM
Definition: med_outils.h:49
ISCRUTE_int
#define ISCRUTE_int(entier)
Definition: med_utils.h:314
med_idt
hid_t med_idt
Definition: med.h:333
MED_FALSE
Definition: med.h:262
MED_NO_GEOTYPE
#define MED_NO_GEOTYPE
Definition: med.h:234
MED_ERR_
#define MED_ERR_(rt, r1, r2, r3)
Definition: med_utils.h:160
MED_NOM_HE8
#define MED_NOM_HE8
Definition: med_outils.h:98
MED_SEG2
#define MED_SEG2
Definition: med.h:202
MED_ERR_ATTRIBUTE
#define MED_ERR_ATTRIBUTE
Definition: med_err.h:100
_MEDmeshAdvancedWr
MEDC_EXPORT med_err _MEDmeshAdvancedWr(const med_idt fid, const char *const meshname, const med_data_type meddatatype, const char *const datasetname, const med_internal_type datatype, const med_int numdt, const med_int numit, const med_float dt, const med_entity_type entitype, const med_geometry_type geotype, const med_connectivity_mode cmode, const med_storage_mode storagemode, const char *const profilename, const med_switch_mode switchmode, const med_int dimselect, const med_filter *const filter, const med_int nentity, const void *const value)
MED_ERR_INVALID
#define MED_ERR_INVALID
Definition: med_err.h:40
MED_MESH_GRP
#define MED_MESH_GRP
Definition: med_outils.h:160
med_err
herr_t med_err
Definition: med.h:334
med_outils.h
MED_NOM_TRF
#define MED_NOM_TRF
Definition: med_outils.h:57
MED_ERR_CLOSE
#define MED_ERR_CLOSE
Definition: med_err.h:30
MED_NOM_CGS
#define MED_NOM_CGS
Definition: med_outils.h:74
MED_NOM_GEO
#define MED_NOM_GEO
Definition: med_outils.h:133
_MEDmeshAssociatedGroupCr
MEDC_EXPORT med_idt _MEDmeshAssociatedGroupCr(const med_idt id, const char *const rootname, const char *const meshname, const med_int numdt, const med_int numit, const med_float dt, const med_bool justopen, const char *const datagroupname)
MED_NOM_NOE
#define MED_NOM_NOE
Definition: med_outils.h:81
MED_ERR_OPEN
#define MED_ERR_OPEN
Definition: med_err.h:37
med_int
int med_int
Definition: med.h:344
MED_CURVILINEAR_GRID
Definition: med.h:141
_MEDisasoftlink
MEDC_EXPORT med_err _MEDisasoftlink(const med_idt id, const char *const linkname, med_bool linkmustexist, med_bool *const isasoftlink)
ISCRUTE
#define ISCRUTE(entier)
Definition: med_utils.h:313
MED_ERR_WRITE
#define MED_ERR_WRITE
Definition: med_err.h:29
med_config.h
_MEDsoftlinkDel
MEDC_EXPORT med_err _MEDsoftlinkDel(const med_idt id, const char *const softlinkname, med_bool linkmustexist)
_MEDdatasetFermer
MEDC_EXPORT med_err _MEDdatasetFermer(med_idt id)
MED_NO_PROFILE_INTERNAL
#define MED_NO_PROFILE_INTERNAL
Definition: med.h:298
med_bool
med_bool
Definition: med.h:262
MED_ERR_CALL
#define MED_ERR_CALL
Definition: med_err.h:48
MED_HEXA8
#define MED_HEXA8
Definition: med.h:215
MED_NOM_TYP
#define MED_NOM_TYP
Definition: med_outils.h:65
MED_ERR_READ
#define MED_ERR_READ
Definition: med_err.h:28
med_float
double med_float
Definition: med.h:338
MED_COORDINATE
Definition: med.h:151
MED_NO_CMODE
Definition: med.h:257
MED_POINT1
#define MED_POINT1
Definition: med.h:200
MED_NO_DT
#define MED_NO_DT
Definition: med.h:322
MED_NOM_PO1
#define MED_NOM_PO1
Definition: med_outils.h:85
MED_INTERNAL_UNDEF
Definition: med.h:163
MED_NONE
#define MED_NONE
Definition: med.h:233
MED_NOM_MAI
#define MED_NOM_MAI
Definition: med_outils.h:78
_MEDattributeStringWr
MEDC_EXPORT med_err _MEDattributeStringWr(med_idt pid, const char *const attname, const med_size attsize, const char *const val)
MED_TAILLE_NOM_ENTITE
#define MED_TAILLE_NOM_ENTITE
Definition: med_outils.h:38
MED_NOM_GTY
#define MED_NOM_GTY
Definition: med_outils.h:66
MED_ERR_DATAGROUP
#define MED_ERR_DATAGROUP
Definition: med_err.h:99
_MEDdatagroupCreer
MEDC_EXPORT med_idt _MEDdatagroupCreer(med_idt pid, const char *const nom)
MED_COOTRF
#define MED_COOTRF
Definition: med_outils.h:203
MED_ERR_DATASET
#define MED_ERR_DATASET
Definition: med_err.h:98
SSCRUTE
#define SSCRUTE(chaine)
Definition: med_utils.h:323
MED_NOM_CGT
#define MED_NOM_CGT
Definition: med_outils.h:73
med_mesh_type
med_mesh_type
Definition: med.h:133
_MEDdatasetOuvrir
MEDC_EXPORT med_idt _MEDdatasetOuvrir(med_idt pid, char *nom)
MED_NAME_SIZE
#define MED_NAME_SIZE
Definition: med.h:83
_MEDcheckVersion30
MEDC_EXPORT med_err _MEDcheckVersion30(med_idt fid)
_MEDmeshDatagroupOpen
MEDC_EXPORT med_idt _MEDmeshDatagroupOpen(const med_idt fid, const char *const meshname, char *const usedpath, med_bool *const isasupportmesh)
med_switch_mode
med_switch_mode
Definition: med.h:98
MED_TAILLE_COOTRF
#define MED_TAILLE_COOTRF
Definition: med_outils.h:205
_MEDattrEntierLire
#define _MEDattrEntierLire(x, y, z)
Definition: med_hdfi.h:68
MED_ALL_CONSTITUENT
#define MED_ALL_CONSTITUENT
Definition: med.h:301
MED_NODE
Definition: med.h:145
MED_NO_PROFILE
#define MED_NO_PROFILE
Definition: med.h:283
MED_MESH_SUPPORT_GRP_SIZE
#define MED_MESH_SUPPORT_GRP_SIZE
Definition: med_outils.h:165
MED_NOM_PFL
#define MED_NOM_PFL
Definition: med_outils.h:140
med.h
MED_ERR_MESH_MSG
#define MED_ERR_MESH_MSG
Definition: med_err.h:133
MED_NOM_COO
#define MED_NOM_COO
Definition: med_outils.h:56
_MEDdatagroupOuvrir
MEDC_EXPORT med_idt _MEDdatagroupOuvrir(med_idt pid, const char *const nom)
MED_NO_IT
#define MED_NO_IT
Definition: med.h:323
MED_NOM_SE2
#define MED_NOM_SE2
Definition: med_outils.h:86
MED_ERR_RANGE
#define MED_ERR_RANGE
Definition: med_err.h:34
MEDmeshNodeCoordinateWr
med_err MEDmeshNodeCoordinateWr(const med_idt fid, const char *const meshname, const med_int numdt, const med_int numit, const med_float dt, const med_switch_mode switchmode, const med_int nentity, const med_float *const coordinate)
Cette routine permet d'écrire dans un maillage le tableau des coordonnées des noeuds,...
Definition: MEDmeshNodeCoordinateWr.c:45
MED_NOM_QU4
#define MED_NOM_QU4
Definition: med_outils.h:90
MED_ERR_NULL
#define MED_ERR_NULL
Definition: med_err.h:35
MED_ERR_PARAMETER
#define MED_ERR_PARAMETER
Definition: med_err.h:103
MED_UNDEF_GRID_TYPE
Definition: med.h:142
MED_QUAD4
#define MED_QUAD4
Definition: med.h:206
MED_NO_NAME
#define MED_NO_NAME
Definition: med.h:270
_MEDdatagroupFermer
MEDC_EXPORT med_err _MEDdatagroupFermer(med_idt id)
MED_ERR_GRIDTYPE
#define MED_ERR_GRIDTYPE
Definition: med_err.h:119