python-libmed #215425 Could not retrieve the number of node entities in field [open]
Lecture d'un fichier MED qui contient un champ, produit par Code_Aster. Impossible d'en extraire les valeurs:
_MEDfieldnValue30.c [103] : Erreur à l'ouverture du groupe
_MEDfieldnValue30.c [103] : du champ résultat
_MEDfieldnValue30.c [104] : _path = "/CHA/MODES1__DEPL/-0000000000000000001-0000000000000000001/"
MEDfieldnValueWithProfile.c [65] : Erreur d'appel de l'API
MEDfieldnValueWithProfile.c [65] : du champ résultat
MEDfieldnValueWithProfile.c [66] : fieldname = "MODES1__DEPL"
MEDfieldnValueWithProfile.c [66] : "_MEDfieldnValue" = "_MEDfieldnValue"
Traceback (most recent call last):
File "/local/users/dallolio/devel/workspace/testmed/src/test_mednx1.py", line 23, in <module>
print fd.values()
File "/local/users/dallolio/devel/workspace/testmed/src/libmed/field.py", line 151, in values
etypes = self.entities()
File "/local/users/dallolio/devel/workspace/testmed/src/libmed/field.py", line 123, in entities
return filter(filter_func, ['node'] + msh.elements_types())
File "/local/users/dallolio/devel/workspace/testmed/src/libmed/field.py", line 122, in <lambda>
numit)
File "_field.pyx", line 157, in libmed._field.read_number_of_values_with_profile (/home/dallolio/.pyxbld/temp.linux-x86_64-2.7/pyrex/libmed/_field.c:5205)
libmed._utils.MEDError: Could not retrieve the number of node entities in field MODES1__DEPL.
Si utile je peux mettre à disposition le (petit) fichier MED et le code python pour la lecture. | |
| priority | normal |
|---|---|
| type | bug |
| done in | <not specified> |
| closed by | <not specified> |
Comments
-
2014/02/24 08:33, written by dlaxalde
-
2014/02/25 17:08, written by luca.dallolio@gmail.com
add commentOui, les données pour reproduire le problème seraient utiles. Merci.
J'ai examiné le problème avec un accès direct par la bibliothèque MED, ce qui m'a donné des informations (je crois) sur l'origine du problème dans libmed : le champ contient plusieurs étapes temporelles, le message d'erreur (je crois) indique que l'étape n'est pas bien affectée. Voici un code MED (qui fonctionne avec les mêmes données) qui peut peut-être aider :
from med.medfile import * # @UnusedWildImport from med.medmesh import * # @UnusedWildImport from med.medfamily import * # @UnusedWildImport from med.medfield import * # @UnusedWildImport from med.medenum import * # @UnusedWildImport from med.medprofile import * # @UnusedWildImport from med.medlocalization import * # @UnusedWildImport from med.medlink import * # @UnusedWildImport if __name__ == '__main__': fid= MEDfileOpen("bd05barV4.rmed",MED_ACC_RDONLY) print fid mm,m,r = MEDfileNumVersionRd(fid) print "Version Du Fichier (num) :",mm,m,r comment = MEDfileCommentRd(fid) print comment maa, sdim, mdim, meshtype, desc, dtunit, sort, nstep, repere, axisname, axisunit = MEDmeshInfo(fid, 1) print "Maillage de nom : |%s| , de dimension : %d , et de type %s"%(maa,mdim,meshtype) print "\t -Dimension de l'espace : %d"%(sdim) print "\t -Description du maillage : |%s|"%(desc) print "\t -Noms des axes : |%s|"%(axisname) print "\t -Unites des axes : |%s|"%(axisunit) print "\t -Type de repere : %s"%(repere) print "\t -Nombre d'etapes de calcul : %d"%(nstep) print "\t -Unite des dates : |%s|"%(dtunit) ncha = MEDnField(fid) print "Nombre de champs : %d "%(ncha) numchamp = 1 # /* Lecture du nombre de composantes */ ncomp = MEDfieldnComponent(fid,numchamp) # /* Lecture du type du champ, des noms des composantes et du nom de l'unite */ nomcha,_meshname,_local,typcha,comp,unit,_dtunit,_ncstp = MEDfieldInfo(fid,numchamp) print "Nom du champ : |%s| de type %s"%(nomcha,typcha) print "Nombre de composantes : |%d|"%(ncomp) print "Nom des composantes : |%s|"%(comp) print "Unites des composantes : |%s| "%(unit) print "Unites des donnes : |%s| "%(_dtunit) print "Le maillage associe est |%s|"%(_meshname) print "Nombre de sequences de calcul |%d|"%(_ncstp) dt_unit="unknown" type_geo = MED_NO_GEOTYPE numpdt = 5 entite = MED_NODE usermode = MED_COMPACT_STMODE stockage = MED_FULL_INTERLACE numdt, numo, dt, nmesh, meshname,localmesh, meshnumdt, meshnumit = MEDfield23ComputingStepMeshInfo(fid,nomcha,numpdt) nummesh = 1 _nprofile, meshname, pflname, locname = MEDfield23nProfile(fid,nomcha,numdt,numo,entite,type_geo, nummesh ) numprofile = 1 nval,pflname, pflsize, locname, ngauss = MEDfield23nValueWithProfile(fid, nomcha, numdt, numo, entite, type_geo, meshname, numprofile, usermode ) print " +Pas de Temps n.%d (%f) [%s], n. d'ordre %d, avec %d valeur(s) par entite."%(numdt,dt,dt_unit,numo,ngauss) if (typcha.val == MED_FLOAT64): val=MEDFLOAT(ncomp*nval*ngauss) else: val=MEDINT (ncomp*nval*ngauss) MEDfieldValueWithProfileRd(fid, nomcha, numdt,numo, entite,type_geo, usermode, pflname, stockage, MED_ALL_CONSTITUENT, val) print "Values (debut du champ) : %s" % val[1:10] MEDfileClose(fid)