MEDfamLire.c
Aller à la documentation de ce fichier.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019 #include <med.h>
00020 #include <med_config.h>
00021 #include <med_outils.h>
00022
00023 #include <stdlib.h>
00024 #include <string.h>
00025
00026 med_err
00027 MEDfamLire(med_idt fid,char *maa, med_int *fam, med_int n,
00028 med_entite_maillage type_ent,med_geometrie_element type_geo)
00029 {
00030 med_idt root;
00031 med_idt maaid = 0;
00032 med_idt entid = 0;
00033 med_idt geoid = 0;
00034 med_err ret = -1;
00035 char chemin[MED_TAILLE_MAA+MED_TAILLE_NOM+1];
00036 char nom_ent[MED_TAILLE_NOM_ENTITE+1];
00037 char nom_geo[MED_TAILLE_NOM_ENTITE+1];
00038 med_int i;
00039 med_entite_maillage _type_ent=type_ent;
00040
00041 if ( type_ent == MED_NOEUD_MAILLE ) _type_ent=MED_NOEUD ;
00042
00043
00044
00045
00046 _MEDmodeErreurVerrouiller();
00047 if (MEDcheckVersion(fid) < 0) return -1;
00048
00049
00050
00051
00052
00053 strcpy(chemin,MED_MAA);
00054 strcat(chemin,maa);
00055 if ((maaid = _MEDdatagroupOuvrir(fid,chemin)) < 0) {
00056 MESSAGE("Impossible d'acceder au maillage");
00057 SSCRUTE(chemin);
00058 goto ERREUR;
00059 }
00060
00061
00062
00063
00064
00065 if (_MEDnomEntite(nom_ent,_type_ent) < 0) {
00066 MESSAGE("L'entité demandée n'est pas une entité <med_entite_maillage> : ");
00067 SSCRUTE(nom_ent);
00068 ISCRUTE(_type_ent);
00069 goto ERREUR;
00070 }
00071
00072
00073
00074
00075 if ((entid = _MEDdatagroupOuvrir(maaid,nom_ent)) < 0) {
00076 MESSAGE("Impossible d'accéder au groupe des entités : ");
00077 SSCRUTE(nom_ent);
00078 ISCRUTE(_type_ent);
00079 goto ERREUR;
00080 }
00081
00082
00083
00084
00085
00086
00087 if ((_type_ent==MED_MAILLE)||(_type_ent==MED_FACE)||(_type_ent==MED_ARETE))
00088 {
00089 if (_MEDnomGeometrie30(nom_geo,type_geo) < 0) {
00090 MESSAGE("L'entité demandée n'est pas un élément <med_geometrie_element> : ");
00091 SSCRUTE(nom_geo);
00092 ISCRUTE(type_geo);
00093 goto ERREUR;
00094 }
00095 if ((geoid = _MEDdatagroupOuvrir(entid,nom_geo)) < 0) {
00096 MESSAGE("Impossible d'accéder au groupe des éléments géométriques : ");
00097 SSCRUTE(nom_geo);
00098 ISCRUTE(type_geo);
00099 goto ERREUR;
00100 }
00101
00102 }
00103 else
00104 geoid = -1;
00105
00106
00107
00108
00109
00110
00111
00112 if (geoid == -1)
00113 root = entid;
00114 else
00115 root = geoid;
00116 #if defined(HAVE_F77INT64)
00117 if (_MEDdatasetNumLire(root,MED_NOM_FAM,MED_INT64,
00118 MED_NO_INTERLACE,1,MED_ALL,
00119 MED_NOPF,MED_NO_PFLMOD,MED_PFL_NON_COMPACT,0,MED_NOPG,0,
00120 (unsigned char *)fam) < 0)
00121 for (i=0; i<n; i++)
00122 *(fam+i) = 0;
00123 #else
00124 if (_MEDdatasetNumLire(root,MED_NOM_FAM,MED_INT32,
00125 MED_NO_INTERLACE,1,MED_ALL,
00126 MED_NOPF,MED_NO_PFLMOD,MED_PFL_NON_COMPACT,0,MED_NOPG,0,
00127 (unsigned char *)fam) < 0)
00128 for (i=0; i<n; i++)
00129 *(fam+i) = 0;
00130 #endif
00131
00132
00133
00134
00135 ret = 0;
00136
00137
00138
00139
00140 ERREUR:
00141 if (geoid > 0)
00142 if (_MEDdatagroupFermer(geoid) < 0)
00143 return -1;
00144 if (entid > 0)
00145 if (_MEDdatagroupFermer(entid) < 0)
00146 return -1;
00147 if (maaid > 0)
00148 if (_MEDdatagroupFermer(maaid) < 0)
00149 return -1;
00150
00151
00152
00153
00154
00155 return (ret);
00156 }