MEDmeshGridIndexCoordinateWr.c
Aller à la documentation de ce fichier.
/* This file is part of MED.
*
* COPYRIGHT (C) 1999 - 2021 EDF R&D, CEA/DEN
* MED is free software: you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* MED is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with MED. If not, see <http://www.gnu.org/licenses/>.
*/
#include <med.h>
#include <med_config.h>
#include <med_outils.h>
#include <string.h>
#include <stdlib.h>
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)
{
med_access_mode _MED_ACCESS_MODE;
med_idt _meshid=0;
med_idt _datagroup1=0,_datagroup2=0,_datagroup3=0;
med_err _ret = -1;
med_data_type _datatype = MED_UNDEF_DATATYPE;
med_grid_type _gridtype = MED_UNDEF_GRID_TYPE;
med_int _intgridtype = 0;
med_int _intaxistype = 0;
med_int _meshdim = 0;
char _meshpath[MED_MESH_GRP_SIZE+MED_NAME_SIZE+1]=MED_MESH_GRP;
char _geotypename[MED_TAILLE_NOM_ENTITE+1]="";
med_int _0=0;
med_int _medintgeotype = MED_NO_GEOTYPE;
char _profilename [MED_NAME_SIZE+1]=MED_SAME_PROFILE_INTERNAL;
/*
* On inhibe le gestionnaire d'erreur HDF 5
*/
_MEDmodeErreurVerrouiller();
if (_MEDcheckVersion30(fid) < 0) goto ERROR;
if ( (_MED_ACCESS_MODE = _MEDmodeAcces(fid) ) == MED_ACC_UNDEF ) {
MED_ERR_(_ret,MED_ERR_UNRECOGNIZED,MED_ERR_ACCESSMODE,MED_ERR_FILE_MSG);
goto ERROR;
}
if ( _MED_ACCESS_MODE == MED_ACC_RDONLY) {
MED_ERR_(_ret,MED_ERR_INVALID,MED_ERR_ACCESSMODE,MED_ERR_FILE_MSG);
ISCRUTE_int(_MED_ACCESS_MODE);
goto ERROR;
}
strcat(_meshpath,meshname);
if ((_meshid = _MEDdatagroupOuvrir(fid,_meshpath)) < 0) {
MED_ERR_(_ret,MED_ERR_OPEN,MED_ERR_DATAGROUP,MED_ERR_MESH_MSG);
SSCRUTE(meshname);SSCRUTE(_meshpath); goto ERROR;
}
/*
* Lecture du type de grille (attribut MED_NOM_GTY)
*/
if (_MEDattrEntierLire(_meshid,MED_NOM_GTY,&_intgridtype) < 0) {
MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
SSCRUTE(meshname);SSCRUTE(MED_NOM_GTY);ISCRUTE_int(_gridtype);goto ERROR;
}
_gridtype=(med_grid_type) _intgridtype;
/*
* Les grilles MED_CURVILINEAR ne sont pas définies en utilisant
* cette fonction, elles doivent définir leurs coordonnées des noeuds
* comme pour les maillages non structurés.
*/
if ((_gridtype != MED_CARTESIAN_GRID) && (_gridtype != MED_POLAR_GRID)) {
MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_GRIDTYPE,MED_ERR_MESH_MSG);
SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
}
/* Lecture de l'attribut MED_NOM_REP */
if ( _MEDattrEntierLire(_meshid,MED_NOM_REP,&_intaxistype) < 0) {
MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
SSCRUTE(meshname);SSCRUTE(MED_NOM_REP);
ISCRUTE(_intaxistype);goto ERROR;
}
/* Lecture de l'attribut MED_NOM_DIM */
if (_MEDattrEntierLire(_meshid,MED_NOM_DIM,&_meshdim) < 0) {
MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
SSCRUTE(meshname);SSCRUTE(MED_NOM_DIM);ISCRUTE(_meshdim);goto ERROR;
}
if (axis > _meshdim ) {
MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
SSCRUTE(meshname);ISCRUTE(_meshdim);ISCRUTE(axis);goto ERROR;
}
if ( (_datagroup1 =_MEDmeshAssociatedGroupCr(fid, MED_MESH_GRP,
meshname, numdt, numit, dt, MED_FALSE,
"." ) ) < 0 ) {
MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAssociatedGroupCr");
goto ERROR;
}
if ((_datagroup2 = _MEDdatagroupOuvrir(_datagroup1,MED_NOM_MAI)) < 0) {
if ((_datagroup2 = _MEDdatagroupCreer(_datagroup1,MED_NOM_MAI)) < 0) {
MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,MED_NOM_MAI);
SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
goto ERROR;
}
if ( _MEDattributeIntWr(_datagroup2,MED_NOM_CGT,&_0) < 0) {
MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGT);
goto ERROR;
}
/*
Le type géométrique de maille utilisé pour les grilles est fonction
de la dimension de la grille == nombre d'axes
*/
switch ( _meshdim ) {
case 1 :
strcpy(_geotypename,MED_NOM_SE2);
_medintgeotype = MED_SEG2;
break;
case 2 :
strcpy(_geotypename,MED_NOM_QU4);
_medintgeotype = MED_QUAD4;
break;
case 3 :
strcpy(_geotypename,MED_NOM_HE8);
_medintgeotype = MED_HEXA8;
break;
case 0 : strcpy(_geotypename,MED_NOM_PO1);
_medintgeotype = MED_POINT1;
break;
default :
MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_PARAMETER,"");
SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
goto ERROR;
}
if ((_datagroup3 = _MEDdatagroupCreer(_datagroup2,_geotypename)) < 0) {
MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,_geotypename);
SSCRUTE(meshname);ISCRUTE(numdt);ISCRUTE(numit);
goto ERROR;
}
if ( _MEDattributeStringWr(_datagroup3,MED_NOM_PFL,MED_NAME_SIZE,MED_NO_PROFILE_INTERNAL) < 0) {
MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);
goto ERROR;
}
if (_MEDattributeIntWr(_datagroup3,MED_NOM_GEO,&_medintgeotype) < 0) {
MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_NOM_GEO);
ISCRUTE(_medintgeotype);
goto ERROR;
}
if ( _MEDattributeIntWr(_datagroup3,MED_NOM_CGT,&_0) < 0) {
MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGT);
goto ERROR;
}
if ( _MEDattributeIntWr(_datagroup3,MED_NOM_CGS,&_0) < 0) {
MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
SSCRUTE(meshname);ISCRUTE(numit);ISCRUTE(numdt);SSCRUTE(MED_NOM_CGS);
goto ERROR;
}
}
/*
* On ecrit le tableau d'indice dans un dataset HDF
*/
switch(axis)
{
case 1 :
_datatype = MED_COORDINATE_AXIS1;
_profilename[0] = '\0';
break;
case 2 :
_datatype = MED_COORDINATE_AXIS2;
break;
case 3 :
_datatype = MED_COORDINATE_AXIS3;
break;
default :
MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_GRIDTYPE,MED_ERR_MESH_MSG);
SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
}
if (_MEDmeshAdvancedWr(fid,
meshname,
_datatype,
MED_NO_NAME,
MED_INTERNAL_UNDEF,
numdt,
numit,
dt,
MED_NODE,
MED_NONE,
MED_NO_CMODE,
MED_UNDEF_STMODE,
_profilename,
MED_FULL_INTERLACE,
MED_ALL_CONSTITUENT,
NULL,
indexsize,
gridindex) < 0 ) {
MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAdvancedWr");
goto ERROR;
}
_ret = 0;
ERROR:
if (_datagroup3>0) if (_MEDdatagroupFermer(_datagroup3) < 0) {
MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_geotypename);
ISCRUTE_id(_datagroup3);
}
if (_datagroup2>0) if (_MEDdatagroupFermer(_datagroup2) < 0) {
MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,MED_NOM_MAI);
ISCRUTE_id(_datagroup2);
}
if (_datagroup1>0) if (_MEDdatagroupFermer(_datagroup1) < 0) {
MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,meshname);
ISCRUTE_id(_datagroup1);
}
if (_meshid>0) if (_MEDdatagroupFermer(_meshid) < 0) {
MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_meshpath);
ISCRUTE_id(_meshid);
}
return _ret;
}
MED_ERR_CREATE
#define MED_ERR_CREATE
Definition: med_err.h:27
MED_FULL_INTERLACE
Definition: med.h:96
MED_ERR_API
#define MED_ERR_API
Definition: med_err.h:111
MED_ERR_GRIDTYPE
#define MED_ERR_GRIDTYPE
Definition: med_err.h:119
med_grid_type
med_grid_type
Definition: med.h:137
_MEDdatagroupFermer
MEDC_EXPORT med_err _MEDdatagroupFermer(med_idt id)
MED_NOM_DIM
#define MED_NOM_DIM
Definition: med_outils.h:49
ISCRUTE_id
#define ISCRUTE_id(entier)
Definition: med_utils.h:319
MED_COORDINATE_AXIS2
Definition: med.h:150
_MEDmodeAcces
MEDC_EXPORT med_access_mode _MEDmodeAcces(med_idt oid)
ISCRUTE_int
#define ISCRUTE_int(entier)
Definition: med_utils.h:314
MED_NO_GEOTYPE
#define MED_NO_GEOTYPE
Definition: med.h:232
med_idt
hid_t med_idt
Definition: med.h:331
MED_ERR_
#define MED_ERR_(rt, r1, r2, r3)
Definition: med_utils.h:160
_MEDattributeIntWr
#define _MEDattributeIntWr(w, x, y)
Definition: med_hdfi.h:45
MED_CARTESIAN_GRID
Definition: med.h:137
MED_SEG2
#define MED_SEG2
Definition: med.h:200
MED_NOM_REP
#define MED_NOM_REP
Definition: med_outils.h:61
MED_ACC_RDONLY
Definition: med.h:120
MED_NOM_HE8
#define MED_NOM_HE8
Definition: med_outils.h:98
MED_ERR_INVALID
#define MED_ERR_INVALID
Definition: med_err.h:40
MED_UNDEF_STMODE
Definition: med.h:108
MED_ERR_ATTRIBUTE
#define MED_ERR_ATTRIBUTE
Definition: med_err.h:100
med_err
herr_t med_err
Definition: med.h:332
MED_MESH_GRP
#define MED_MESH_GRP
Definition: med_outils.h:160
MED_NO_CMODE
Definition: med.h:255
_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_FALSE
Definition: med.h:260
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_ERR_OPEN
#define MED_ERR_OPEN
Definition: med_err.h:37
MED_NOM_GEO
#define MED_NOM_GEO
Definition: med_outils.h:133
med_int
int med_int
Definition: med.h:342
MED_UNDEF_DATATYPE
Definition: med.h:152
MED_ERR_WRITE
#define MED_ERR_WRITE
Definition: med_err.h:29
_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_ERR_CALL
#define MED_ERR_CALL
Definition: med_err.h:48
MED_NO_PROFILE_INTERNAL
#define MED_NO_PROFILE_INTERNAL
Definition: med.h:296
med.h
ISCRUTE
#define ISCRUTE(entier)
Definition: med_utils.h:313
MED_HEXA8
#define MED_HEXA8
Definition: med.h:213
MED_UNDEF_GRID_TYPE
Definition: med.h:140
med_float
double med_float
Definition: med.h:336
MED_ERR_READ
#define MED_ERR_READ
Definition: med_err.h:28
MED_NOM_PO1
#define MED_NOM_PO1
Definition: med_outils.h:85
MED_POINT1
#define MED_POINT1
Definition: med.h:198
MED_NONE
#define MED_NONE
Definition: med.h:231
MED_MESH_GRP_SIZE
#define MED_MESH_GRP_SIZE
Definition: med_outils.h:161
med_access_mode
med_access_mode
Definition: med.h:120
MED_NOM_GTY
#define MED_NOM_GTY
Definition: med_outils.h:66
MED_NODE
Definition: med.h:143
MED_NOM_MAI
#define MED_NOM_MAI
Definition: med_outils.h:78
MED_ERR_ACCESSMODE
#define MED_ERR_ACCESSMODE
Definition: med_err.h:104
_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
_MEDdatagroupCreer
MEDC_EXPORT med_idt _MEDdatagroupCreer(med_idt pid, const char *const nom)
SSCRUTE
#define SSCRUTE(chaine)
Definition: med_utils.h:323
MED_ERR_DATAGROUP
#define MED_ERR_DATAGROUP
Definition: med_err.h:99
MED_NOM_CGT
#define MED_NOM_CGT
Definition: med_outils.h:73
MED_NAME_SIZE
#define MED_NAME_SIZE
Definition: med.h:81
_MEDcheckVersion30
MEDC_EXPORT med_err _MEDcheckVersion30(med_idt fid)
med_data_type
med_data_type
Definition: med.h:149
med_outils.h
med_config.h
MED_ALL_CONSTITUENT
#define MED_ALL_CONSTITUENT
Definition: med.h:299
MED_NOM_PFL
#define MED_NOM_PFL
Definition: med_outils.h:140
MEDmeshGridIndexCoordinateWr
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...
Definition: MEDmeshGridIndexCoordinateWr.c:43
_MEDattrEntierLire
#define _MEDattrEntierLire(x, y, z)
Definition: med_hdfi.h:68
MED_ERR_RANGE
#define MED_ERR_RANGE
Definition: med_err.h:34
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
MED_COORDINATE_AXIS3
Definition: med.h:150
MED_NOM_SE2
#define MED_NOM_SE2
Definition: med_outils.h:86
MED_NOM_QU4
#define MED_NOM_QU4
Definition: med_outils.h:90
MED_ERR_PARAMETER
#define MED_ERR_PARAMETER
Definition: med_err.h:103
MED_INTERNAL_UNDEF
Definition: med.h:161
MED_SAME_PROFILE_INTERNAL
#define MED_SAME_PROFILE_INTERNAL
Definition: med.h:297
MED_ERR_FILE_MSG
#define MED_ERR_FILE_MSG
Definition: med_err.h:132
_MEDdatagroupOuvrir
MEDC_EXPORT med_idt _MEDdatagroupOuvrir(med_idt pid, const char *const nom)
MED_COORDINATE_AXIS1
Definition: med.h:150
MED_QUAD4
#define MED_QUAD4
Definition: med.h:204
_MEDmodeErreurVerrouiller
MEDC_EXPORT void _MEDmodeErreurVerrouiller(void)
MED_NO_NAME
#define MED_NO_NAME
Definition: med.h:268
MED_POLAR_GRID
Definition: med.h:138
MED_ACC_UNDEF
Definition: med.h:124