MED fichier: Fichier source de MEDmeshGridStructWr.c

MED fichier

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