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 }