MED fichier: Fichier source de MEDmeshGridStructRd.c

MED fichier

MEDmeshGridStructRd.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
39 med_err
40 MEDmeshGridStructRd(const med_idt fid,
41 const char* const meshname,
42 const med_int numdt,
43 const med_int numit,
44 med_int * const gridstruct)
45 {
46 med_idt _meshid=0;
47 med_err _ret = -1;
48 med_data_type _datatype = MED_UNDEF_DATATYPE;
49 med_grid_type _gridtype = MED_UNDEF_GRID_TYPE;
50 med_int _intgridtype = 0;
51 med_int _meshdim = 0;
52 char _meshpath[MED_MESH_GRP_SIZE+MED_NAME_SIZE+1]=MED_MESH_GRP;
53 int _i=0;
54
55 med_int * _ptr;
56
57 /*
58 * On inhibe le gestionnaire d'erreur HDF 5
59 */
60 _MEDmodeErreurVerrouiller();
61
62
63 strcat(_meshpath,meshname);
64 if ((_meshid = _MEDdatagroupOuvrir(fid,_meshpath)) < 0) {
65 MED_ERR_(_ret,MED_ERR_OPEN,MED_ERR_DATAGROUP,MED_ERR_MESH_MSG);
66 SSCRUTE(meshname);SSCRUTE(_meshpath); goto ERROR;
67 }
68
69 /* Lecture de l'attribut MED_NOM_DIM */
70 if (_MEDattrEntierLire(_meshid,MED_NOM_DIM,&_meshdim) < 0) {
71 MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
72 SSCRUTE(meshname);SSCRUTE(MED_NOM_DIM);ISCRUTE(_meshdim);goto ERROR;
73 }
74
75 /*
76 * Lecture du type de grille (attribut MED_NOM_GTY)
77 */
78 if (_MEDattrEntierLire(_meshid,MED_NOM_GTY,&_intgridtype) < 0) {
79 MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
80 SSCRUTE(meshname);SSCRUTE(MED_NOM_GTY);ISCRUTE(_intgridtype);goto ERROR;
81 }
82 _gridtype=(med_grid_type) _intgridtype;
83
84 if (_gridtype != MED_CURVILINEAR_GRID) {
85 MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_GRIDTYPE,MED_ERR_MESH_MSG);
86 SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
87 }
88
89
90 /*
91 * Attribut NBR (taille du tableau d'indices)
92 */
93 for (_i=0;_i<_meshdim;++_i) {
94
95 switch(_i) {
96 case 0 :
97 _datatype = MED_COORDINATE_AXIS1;
98 break;
99 case 1 :
100 _datatype = MED_COORDINATE_AXIS2;
101 break;
102 case 2 :
103 _datatype = MED_COORDINATE_AXIS3;
104 break;
105 default :
106 MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_GRIDTYPE,MED_ERR_MESH_MSG);
107 SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
108 }
109
110 _ptr=gridstruct+_i;
111 if (_MEDmeshAdvancedRd(fid,
112 meshname,
113 _datatype,
114 MED_NO_NAME,
115 MED_INT,
116 numdt,
117 numit,
118 MED_NODE,
119 MED_NONE,
120 MED_NO_CMODE,
121 MED_UNDEF_STMODE,
122 MED_NO_PROFILE,
123 MED_FULL_INTERLACE,
124 MED_ALL_CONSTITUENT,
125 NULL,
126 (unsigned char* const) (_ptr)) < 0 ) {
127 MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAdvancedRd");
128 goto ERROR;
129 }
130
131 }
132
133 _ret = 0;
134
135 ERROR:
136
137 if (_meshid>0) if (_MEDdatagroupFermer(_meshid) < 0) {
138 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_meshpath);
139 ISCRUTE_id(_meshid);
140 }
141
142 return _ret;
143
144 }
MED_ERR_API
#define MED_ERR_API
Definition: med_err.h:111
med_grid_type
med_grid_type
Definition: med.h:139
_MEDmeshAdvancedRd
MEDC_EXPORT med_err _MEDmeshAdvancedRd(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_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, unsigned char *const value)
ISCRUTE_id
#define ISCRUTE_id(entier)
Definition: med_utils.h:319
MED_UNDEF_STMODE
Definition: med.h:110
MED_COORDINATE_AXIS1
Definition: med.h:152
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_UNDEF_DATATYPE
Definition: med.h:154
MED_ERR_
#define MED_ERR_(rt, r1, r2, r3)
Definition: med_utils.h:160
MED_ERR_ATTRIBUTE
#define MED_ERR_ATTRIBUTE
Definition: med_err.h:100
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_ERR_CLOSE
#define MED_ERR_CLOSE
Definition: med_err.h:30
MED_FULL_INTERLACE
Definition: med.h:98
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
ISCRUTE
#define ISCRUTE(entier)
Definition: med_utils.h:313
med_config.h
MED_COORDINATE_AXIS2
Definition: med.h:152
MED_ERR_CALL
#define MED_ERR_CALL
Definition: med_err.h:48
MED_ERR_READ
#define MED_ERR_READ
Definition: med_err.h:28
MED_NO_CMODE
Definition: med.h:257
MED_MESH_GRP_SIZE
#define MED_MESH_GRP_SIZE
Definition: med_outils.h:161
MED_NONE
#define MED_NONE
Definition: med.h:233
MED_NOM_GTY
#define MED_NOM_GTY
Definition: med_outils.h:66
MED_ERR_DATAGROUP
#define MED_ERR_DATAGROUP
Definition: med_err.h:99
MEDmeshGridStructRd
med_err MEDmeshGridStructRd(const med_idt fid, const char *const meshname, const med_int numdt, const med_int numit, med_int *const gridstruct)
Cette routine permet la lecture de la structure (nombre de points sur chaque axe du repère) d'un mail...
Definition: MEDmeshGridStructRd.c:40
MED_INT
Definition: med.h:172
SSCRUTE
#define SSCRUTE(chaine)
Definition: med_utils.h:323
MED_NAME_SIZE
#define MED_NAME_SIZE
Definition: med.h:83
med_data_type
med_data_type
Definition: med.h:151
MED_COORDINATE_AXIS3
Definition: med.h:152
_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.h
MED_ERR_MESH_MSG
#define MED_ERR_MESH_MSG
Definition: med_err.h:133
_MEDdatagroupOuvrir
MEDC_EXPORT med_idt _MEDdatagroupOuvrir(med_idt pid, const char *const nom)
MED_ERR_RANGE
#define MED_ERR_RANGE
Definition: med_err.h:34
_MEDmodeErreurVerrouiller
MEDC_EXPORT void _MEDmodeErreurVerrouiller(void)
MED_UNDEF_GRID_TYPE
Definition: med.h:142
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