MED fichier
MEDmeshNodeCoordinateWr.c
Aller à la documentation de ce fichier.
1/* This file is part of MED.
2 *
3 * COPYRIGHT (C) 1999 - 2020 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
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;
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]="";
65 int _i=0;
66 med_bool _isasupportmesh=MED_FALSE,_isasoftlink=MED_FALSE;
67 med_int _0 = 0;
68 med_int _medintgeotype = MED_NO_GEOTYPE;
69
70if (_MEDcheckVersion30(fid) < 0) goto ERROR;
71
72 /* Ouverture du maillage de calcul ou du maillage support */
73 if ((_meshid=_MEDmeshDatagroupOpen(fid,meshname,_meshpath,&_isasupportmesh)) < 0) {
75 SSCRUTE(_meshpath); goto ERROR;
76 }
77
78 /* Lecture du type de maillage (attribut MED_NOM_TYP) */
79 if (_MEDattrEntierLire(_meshid,MED_NOM_TYP,&_intmeshtype) < 0) {
81 SSCRUTE(meshname);SSCRUTE(MED_NOM_TYP);ISCRUTE(_intmeshtype);goto ERROR;
82 }
83
84 /*Sortir en erreur si maillage support et numdt,numit != MED_NO_PDT,MED_NO_IT*/
85 /* Accès (création éventuelle) du datagroup numdt,numit */
86
87 if ( _isasupportmesh ) {
88 if (( numdt != MED_NO_DT) && (numit != MED_NO_IT) ) {
90 SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);goto ERROR;
91 }
92 } else
93 if ( (_datagroup1 =_MEDmeshAssociatedGroupCr(fid, MED_MESH_GRP,
94 meshname, numdt, numit, dt,MED_FALSE,
95 "." ) ) < 0 ) {
96 MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAssociatedGroupCr");
97 goto ERROR;
98 }
99
100 /* Création de la structure de connectivité assoiciée aux grilles MED_CURVILLINEAR_GRID*/
101 if ( ( (med_mesh_type) _intmeshtype ) != MED_UNSTRUCTURED_MESH ) {
102
103 /* Lecture de l'attribut MED_NOM_GTY */
104 if (_MEDattrEntierLire(_meshid,MED_NOM_GTY,&_intgridtype) < 0) {
106 SSCRUTE(meshname);SSCRUTE(MED_NOM_GTY);ISCRUTE_int(_gridtype);goto ERROR;
107 }
108 _gridtype= (med_grid_type) _intgridtype;
109
110 if ( _gridtype == MED_CURVILINEAR_GRID ) {
111
112 if ((_datagroup2 = _MEDdatagroupOuvrir(_datagroup1,MED_NOM_MAI)) < 0) {
113
114 if ((_datagroup2 = _MEDdatagroupCreer(_datagroup1,MED_NOM_MAI)) < 0) {
116 SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
117 goto ERROR;
118 }
119
120 if ( _MEDattributeIntWr(_datagroup2,MED_NOM_CGT,&_0) < 0) {
122 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGT);
123 goto ERROR;
124 }
125
126 /* Lecture de l'attribut MED_NOM_DIM */
127 if (_MEDattrEntierLire(_meshid,MED_NOM_DIM,&_meshdim) < 0) {
129 SSCRUTE(meshname);SSCRUTE(MED_NOM_DIM);ISCRUTE(_meshdim);goto ERROR;
130 }
131
132 switch ( _meshdim ) {
133 case 1 :
134 strcpy(_geotypename,MED_NOM_SE2);
135 _medintgeotype = MED_SEG2;
136 break;
137 case 2 :
138 strcpy(_geotypename,MED_NOM_QU4);
139 _medintgeotype = MED_QUAD4;
140 break;
141 case 3 : strcpy(_geotypename,MED_NOM_HE8);
142 break;
143 _medintgeotype = MED_HEXA8;
144 case 0 : strcpy(_geotypename,MED_NOM_PO1);
145 _medintgeotype = MED_POINT1;
146 break;
147 default :
149 SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
150 goto ERROR;
151 }
152
153 if ((_datagroup3 = _MEDdatagroupCreer(_datagroup2,_geotypename)) < 0) {
154 MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,_geotypename);
155 SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
156 goto ERROR;
157 }
158
161 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);
162 goto ERROR;
163 }
164
165 if ( _MEDattributeIntWr(_datagroup3,MED_NOM_CGT,&_0) < 0) {
167 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGT);
168 goto ERROR;
169 }
170
171 if ( _MEDattributeIntWr(_datagroup3,MED_NOM_CGS,&_0) < 0) {
173 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGS);
174 goto ERROR;
175 }
176
177 if (_MEDattributeIntWr(_datagroup3,MED_NOM_GEO,&_medintgeotype) < 0) {
179 ISCRUTE(_medintgeotype);
180 goto ERROR;
181 }
182
183 }
184 } else {
186 SSCRUTE(meshname);SSCRUTE(MED_NOM_GTY);ISCRUTE_int(_gridtype);goto ERROR;
187
188 }
189
190 }
191
192
193 /*Réinitialisation du lien MED_TRSF_NOM/MED_COORDINATE_TRSF s'il existe */
194 /*Si un dataset de transformation existe une erreur est générée par
195 _MEDsoftlinkDel */
196 if ( !_isasupportmesh)
197 if ((_datagroup_trsf = _MEDdatagroupOuvrir(_datagroup1,_datagroup_trsfname)) > 0) {
198
199 if ( _MEDisasoftlink(_datagroup_trsf, MED_NOM_TRF,MED_FALSE, &_isasoftlink ) < 0 ) {
201 goto ERROR;
202 }
203 if (_isasoftlink ) {
204 if ( _MEDsoftlinkDel(_datagroup_trsf,
206 MED_TRUE) < 0 ) {
207 MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDsoftlinkDel");
208 goto ERROR;
209 }
210 } else {
211 if ( (_dataset_trsf = _MEDdatasetOuvrir(_datagroup_trsf,MED_NOM_TRF)) >= 0 ) {
212 _MEDdatasetFermer(_dataset_trsf);
214 goto ERROR;
215 }
216
217 }
218 }
219
220 /*Ecriture des coordonnées*/
221 if ( _MEDmeshAdvancedWr(fid,
222 meshname,
226 numdt,
227 numit,
228 dt,
229 MED_NODE,
230 MED_NONE,
234 switchmode,
236 NULL,
237 nentity,
238 coordinate) < 0 ) {
239 MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAdvancedWr");
240 goto ERROR;
241
242 }
243
244
245 _ret = 0;
246
247 ERROR:
248
249 if (_datagroup_trsf>0) if (_MEDdatagroupFermer(_datagroup_trsf) < 0) {
250 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_datagroup_trsfname);
251 ISCRUTE_id(_datagroup_trsf);
252 }
253
254 if (_datagroup3>0) if (_MEDdatagroupFermer(_datagroup3) < 0) {
255 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_geotypename);
256 ISCRUTE_id(_datagroup3);
257 }
258
259 if (_datagroup2>0) if (_MEDdatagroupFermer(_datagroup2) < 0) {
261 ISCRUTE_id(_datagroup2);
262 }
263
264 if (_datagroup1>0) if (_MEDdatagroupFermer(_datagroup1) < 0) {
266 ISCRUTE_id(_datagroup1);ISCRUTE(numdt);ISCRUTE(numit);
267 }
268
269 if (_meshid>0) if (_MEDdatagroupFermer(_meshid) < 0) {
271 ISCRUTE_id(_meshid);
272 }
273
274 return _ret;
275}
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,...
#define MED_NAME_SIZE
Definition med.h:81
med_switch_mode
Definition med.h:96
med_bool
Definition med.h:260
@ MED_FALSE
Definition med.h:260
@ MED_TRUE
Definition med.h:260
#define MED_QUAD4
Definition med.h:204
@ MED_COORDINATE
Definition med.h:149
#define MED_SEG2
Definition med.h:200
#define MED_NO_PROFILE_INTERNAL
Definition med.h:290
@ MED_UNDEF_STMODE
Definition med.h:108
med_grid_type
Definition med.h:137
@ MED_UNDEF_GRID_TYPE
Definition med.h:140
@ MED_CURVILINEAR_GRID
Definition med.h:139
#define MED_ALL_CONSTITUENT
Definition med.h:293
med_mesh_type
Definition med.h:131
@ MED_UNSTRUCTURED_MESH
Definition med.h:131
#define MED_NO_GEOTYPE
Definition med.h:232
int med_int
Definition med.h:333
#define MED_NO_DT
Definition med.h:311
#define MED_NO_IT
Definition med.h:312
#define MED_NONE
Definition med.h:231
@ MED_INTERNAL_UNDEF
Definition med.h:161
#define MED_NO_PROFILE
Definition med.h:279
@ MED_NODE
Definition med.h:143
#define MED_POINT1
Definition med.h:198
#define MED_HEXA8
Definition med.h:213
double med_float
Definition med.h:327
herr_t med_err
Definition med.h:323
hid_t med_idt
Definition med.h:322
@ MED_NO_CMODE
Definition med.h:255
#define MED_NO_NAME
Definition med.h:266
#define MED_ERR_OPEN
Definition med_err.h:37
#define MED_ERR_LINK
Definition med_err.h:115
#define MED_ERR_CALL
Definition med_err.h:48
#define MED_ERR_DATAGROUP
Definition med_err.h:99
#define MED_ERR_CLOSE
Definition med_err.h:30
#define MED_ERR_MESH_MSG
Definition med_err.h:133
#define MED_ERR_INVALID
Definition med_err.h:40
#define MED_ERR_ATTRIBUTE
Definition med_err.h:100
#define MED_ERR_CREATE
Definition med_err.h:27
#define MED_ERR_NULL
Definition med_err.h:35
#define MED_ERR_API
Definition med_err.h:111
#define MED_ERR_GRIDTYPE
Definition med_err.h:119
#define MED_ERR_PARAMETER
Definition med_err.h:103
#define MED_ERR_READ
Definition med_err.h:28
#define MED_ERR_RANGE
Definition med_err.h:34
#define MED_ERR_WRITE
Definition med_err.h:29
#define MED_ERR_DATASET
Definition med_err.h:98
MEDC_EXPORT med_err _MEDisasoftlink(const med_idt id, const char *const linkname, med_bool linkmustexist, med_bool *const isasoftlink)
MEDC_EXPORT med_err _MEDattributeStringWr(med_idt pid, const char *const attname, const med_size attsize, const char *const val)
#define _MEDattrEntierLire(x, y, z)
Definition med_hdfi.h:68
MEDC_EXPORT med_err _MEDdatasetFermer(med_idt id)
MEDC_EXPORT med_idt _MEDdatagroupCreer(med_idt pid, const char *const nom)
#define _MEDattributeIntWr(w, x, y)
Definition med_hdfi.h:45
MEDC_EXPORT med_err _MEDsoftlinkDel(const med_idt id, const char *const softlinkname, med_bool linkmustexist)
MEDC_EXPORT med_err _MEDdatagroupFermer(med_idt id)
MEDC_EXPORT med_idt _MEDdatasetOuvrir(med_idt pid, char *nom)
MEDC_EXPORT med_idt _MEDdatagroupOuvrir(med_idt pid, const char *const nom)
MEDC_EXPORT med_err _MEDcheckVersion30(med_idt fid)
#define MED_MESH_SUPPORT_GRP_SIZE
Definition med_outils.h:165
#define MED_COOTRF
Definition med_outils.h:203
#define MED_NOM_DIM
Definition med_outils.h:49
#define MED_NOM_QU4
Definition med_outils.h:90
#define MED_TAILLE_NOM_ENTITE
Definition med_outils.h:38
#define MED_TAILLE_COOTRF
Definition med_outils.h:205
#define MED_NOM_COO
Definition med_outils.h:56
#define MED_NOM_CGT
Definition med_outils.h:73
#define MED_MESH_GRP
Definition med_outils.h:160
#define MED_NOM_SE2
Definition med_outils.h:86
#define MED_NOM_PFL
Definition med_outils.h:140
#define MED_NOM_CGS
Definition med_outils.h:74
#define MED_NOM_TYP
Definition med_outils.h:65
#define MED_NOM_PO1
Definition med_outils.h:85
#define MED_NOM_GTY
Definition med_outils.h:66
#define MED_NOM_HE8
Definition med_outils.h:98
#define MED_NOM_TRF
Definition med_outils.h:57
#define MED_NOM_GEO
Definition med_outils.h:133
#define MED_NOM_MAI
Definition med_outils.h:78
#define MED_NOM_NOE
Definition med_outils.h:81
MEDC_EXPORT med_idt _MEDmeshDatagroupOpen(const med_idt fid, const char *const meshname, char *const usedpath, med_bool *const isasupportmesh)
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)
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)
#define MED_ERR_(rt, r1, r2, r3)
Definition med_utils.h:160
#define SSCRUTE(chaine)
Definition med_utils.h:323
#define ISCRUTE(entier)
Definition med_utils.h:313
#define ISCRUTE_id(entier)
Definition med_utils.h:319
#define ISCRUTE_int(entier)
Definition med_utils.h:314