MED fichier
MEDmeshGridIndexCoordinateWr.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
44 const char* const meshname,
45 const med_int numdt,
46 const med_int numit,
47 const med_float dt,
48 const med_int axis,
49 const med_int indexsize,
50 const med_float * const gridindex)
51{
52 med_access_mode _MED_ACCESS_MODE;
53 med_idt _meshid=0;
54 med_idt _datagroup1=0,_datagroup2=0,_datagroup3=0;
55 med_err _ret = -1;
58 med_int _intgridtype = 0;
59 med_int _intaxistype = 0;
60 med_int _meshdim = 0;
62 char _geotypename[MED_TAILLE_NOM_ENTITE+1]="";
63 med_int _0=0;
64 med_int _medintgeotype = MED_NO_GEOTYPE;
65 char _profilename [MED_NAME_SIZE+1]=MED_SAME_PROFILE_INTERNAL;
66
67
68 /*
69 * On inhibe le gestionnaire d'erreur HDF 5
70 */
72 if (_MEDcheckVersion30(fid) < 0) goto ERROR;
73
74 if ( (_MED_ACCESS_MODE = _MEDmodeAcces(fid) ) == MED_ACC_UNDEF ) {
76 goto ERROR;
77 }
78
79 if ( _MED_ACCESS_MODE == MED_ACC_RDONLY) {
81 ISCRUTE_int(_MED_ACCESS_MODE);
82 goto ERROR;
83 }
84
85 strcat(_meshpath,meshname);
86 if ((_meshid = _MEDdatagroupOuvrir(fid,_meshpath)) < 0) {
88 SSCRUTE(meshname);SSCRUTE(_meshpath); goto ERROR;
89 }
90
91
92 /*
93 * Lecture du type de grille (attribut MED_NOM_GTY)
94 */
95 if (_MEDattrEntierLire(_meshid,MED_NOM_GTY,&_intgridtype) < 0) {
97 SSCRUTE(meshname);SSCRUTE(MED_NOM_GTY);ISCRUTE_int(_gridtype);goto ERROR;
98 }
99 _gridtype=(med_grid_type) _intgridtype;
100
101 /*
102 * Les grilles MED_CURVILINEAR ne sont pas définies en utilisant
103 * cette fonction, elles doivent définir leurs coordonnées des noeuds
104 * comme pour les maillages non structurés.
105 */
106 if ((_gridtype != MED_CARTESIAN_GRID) && (_gridtype != MED_POLAR_GRID)) {
108 SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
109 }
110
111 /* Lecture de l'attribut MED_NOM_REP */
112 if ( _MEDattrEntierLire(_meshid,MED_NOM_REP,&_intaxistype) < 0) {
114 SSCRUTE(meshname);SSCRUTE(MED_NOM_REP);
115 ISCRUTE(_intaxistype);goto ERROR;
116 }
117
118
119 /* Lecture de l'attribut MED_NOM_DIM */
120 if (_MEDattrEntierLire(_meshid,MED_NOM_DIM,&_meshdim) < 0) {
122 SSCRUTE(meshname);SSCRUTE(MED_NOM_DIM);ISCRUTE(_meshdim);goto ERROR;
123 }
124
125 if (axis > _meshdim ) {
127 SSCRUTE(meshname);ISCRUTE(_meshdim);ISCRUTE(axis);goto ERROR;
128 }
129
130 if ( (_datagroup1 =_MEDmeshAssociatedGroupCr(fid, MED_MESH_GRP,
131 meshname, numdt, numit, dt, MED_FALSE,
132 "." ) ) < 0 ) {
133 MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAssociatedGroupCr");
134 goto ERROR;
135 }
136
137 if ((_datagroup2 = _MEDdatagroupOuvrir(_datagroup1,MED_NOM_MAI)) < 0) {
138
139 if ((_datagroup2 = _MEDdatagroupCreer(_datagroup1,MED_NOM_MAI)) < 0) {
141 SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
142 goto ERROR;
143 }
144
145 if ( _MEDattributeIntWr(_datagroup2,MED_NOM_CGT,&_0) < 0) {
147 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGT);
148 goto ERROR;
149 }
150
151 /*
152 Le type géométrique de maille utilisé pour les grilles est fonction
153 de la dimension de la grille == nombre d'axes
154 */
155 switch ( _meshdim ) {
156 case 1 :
157 strcpy(_geotypename,MED_NOM_SE2);
158 _medintgeotype = MED_SEG2;
159 break;
160 case 2 :
161 strcpy(_geotypename,MED_NOM_QU4);
162 _medintgeotype = MED_QUAD4;
163 break;
164 case 3 :
165 strcpy(_geotypename,MED_NOM_HE8);
166 _medintgeotype = MED_HEXA8;
167 break;
168 case 0 : strcpy(_geotypename,MED_NOM_PO1);
169 _medintgeotype = MED_POINT1;
170 break;
171 default :
173 SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
174 goto ERROR;
175 }
176
177 if ((_datagroup3 = _MEDdatagroupCreer(_datagroup2,_geotypename)) < 0) {
178 MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,_geotypename);
179 SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
180 goto ERROR;
181 }
182
185 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);
186 goto ERROR;
187 }
188
189 if (_MEDattributeIntWr(_datagroup3,MED_NOM_GEO,&_medintgeotype) < 0) {
191 ISCRUTE(_medintgeotype);
192 goto ERROR;
193 }
194
195 if ( _MEDattributeIntWr(_datagroup3,MED_NOM_CGT,&_0) < 0) {
197 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGT);
198 goto ERROR;
199 }
200
201 if ( _MEDattributeIntWr(_datagroup3,MED_NOM_CGS,&_0) < 0) {
203 SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGS);
204 goto ERROR;
205 }
206
207 }
208
209
210 /*
211 * On ecrit le tableau d'indice dans un dataset HDF
212 */
213 switch(axis)
214 {
215 case 1 :
216 _datatype = MED_COORDINATE_AXIS1;
217 _profilename[0] = '\0';
218 break;
219 case 2 :
220 _datatype = MED_COORDINATE_AXIS2;
221 break;
222 case 3 :
223 _datatype = MED_COORDINATE_AXIS3;
224 break;
225 default :
227 SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
228 }
229
230 if (_MEDmeshAdvancedWr(fid,
231 meshname,
232 _datatype,
235 numdt,
236 numit,
237 dt,
238 MED_NODE,
239 MED_NONE,
242 _profilename,
245 NULL,
246 indexsize,
247 gridindex) < 0 ) {
248 MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAdvancedWr");
249 goto ERROR;
250 }
251
252 _ret = 0;
253
254 ERROR:
255
256 if (_datagroup3>0) if (_MEDdatagroupFermer(_datagroup3) < 0) {
257 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_geotypename);
258 ISCRUTE_id(_datagroup3);
259 }
260
261 if (_datagroup2>0) if (_MEDdatagroupFermer(_datagroup2) < 0) {
263 ISCRUTE_id(_datagroup2);
264 }
265
266 if (_datagroup1>0) if (_MEDdatagroupFermer(_datagroup1) < 0) {
268 ISCRUTE_id(_datagroup1);
269 }
270
271 if (_meshid>0) if (_MEDdatagroupFermer(_meshid) < 0) {
273 ISCRUTE_id(_meshid);
274 }
275
276 return _ret;
277
278}
#define MED_NAME_SIZE
#define MED_NO_PROFILE_INTERNAL
#define MED_ALL_CONSTITUENT
#define MED_NO_GEOTYPE
#define MED_SAME_PROFILE_INTERNAL
#define MED_NO_NAME
#define MED_ERR_(rt, r1, r2, r3)
#define SSCRUTE(chaine)
#define ISCRUTE(entier)
#define ISCRUTE_id(entier)
#define ISCRUTE_int(entier)
med_err MEDmeshGridIndexCoordinateWr(const med_idt fid, const char *const meshname, const med_int numdt, const med_int numit, const med_float dt, const med_int axis, const med_int indexsize, const med_float *const gridindex)
Cette routine permet l'écriture des coordonnées des noeuds d'un maillage structuré selon un axe du re...
#define MED_ERR_OPEN
Definition med_err.h:37
#define MED_ERR_ACCESSMODE
Definition med_err.h:104
#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_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_FILE_MSG
Definition med_err.h:132
#define MED_ERR_WRITE
Definition med_err.h:29
#define MED_ERR_UNRECOGNIZED
Definition med_err.h:39
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_idt _MEDdatagroupCreer(med_idt pid, const char *const nom)
#define _MEDattributeIntWr(w, x, y)
Definition med_hdfi.h:45
MEDC_EXPORT void _MEDmodeErreurVerrouiller(void)
MEDC_EXPORT med_err _MEDdatagroupFermer(med_idt id)
MEDC_EXPORT med_idt _MEDdatagroupOuvrir(med_idt pid, const char *const nom)
MEDC_EXPORT med_err _MEDcheckVersion30(med_idt fid)
MEDC_EXPORT med_access_mode _MEDmodeAcces(med_idt oid)
#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_NOM_CGT
Definition med_outils.h:73
#define MED_MESH_GRP_SIZE
Definition med_outils.h:161
#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_REP
Definition med_outils.h:61
#define MED_NOM_CGS
Definition med_outils.h:74
#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_GEO
Definition med_outils.h:133
#define MED_NOM_MAI
Definition med_outils.h:78
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)