<pre><span style>Dear Gmsh Team.</span><br style><div style><br>Here is as an attached file a small med 3.0.4 file that cannot be opened in Gmsh 2.5.1-svn. FYI, it can be opened in Salome.<br>The file was created from an initial med file that Gmsh can read by applying some algebraic operations on the fields it contains using MEDLoader (a package in charge of loading from a file or write to a file in MED format a med data structure).</div>

<div style><br></div><div style><div>dimrestriction =0</div><div>mesh=MEDLoader.ReadUMeshFromFile(medfilename,meshname,dimrestriction)</div><div>p1 = MEDLoader.ReadField(ON_NODES,medfilename,meshname,dimrestriction,fieldname,iteration1,order1)</div>

<div>p2 = MEDLoader.ReadField(ON_NODES,medfilename,meshname,dimrestriction,fieldname,iteration2,order2)</div><div>p1.setMesh(mesh)</div><div>p2.setMesh(mesh)</div><div>result=p2-p1</div><div>result.setName(fieldnameresu)</div>

<div>MEDLoader.WriteField(outfile,result,True)</div></div><div style><br></div><div style>The error raised by Gmsh is "Could not read MED families", that can be found twice in Gmsh's sources :</div><div style>

<br></div><div style>Geo/GModelIO_MED.cpp:359:      Msg::Error("Could not read MED families");</div><div style>Geo/GModelIO_MED.cpp:387:    Msg::Error("Could not read MED families");</div><div style> </div>

<div style><div>#if (MED_MAJOR_NUM == 3)</div><div>    if(MEDmeshEntityFamilyNumberRd(fid, meshName, MED_NO_DT, MED_NO_IT, MED_CELL,</div><div>                                   type, &fam[0]) < 0){</div><div>#else</div>

<div>    if(MEDfamLire(fid, meshName, &fam[0], numEle, MED_MAILLE, type) < 0){</div><div>#endif</div><div>      Msg::Error("Could not read MED families");</div><div>      return 0;</div><div>    }</div></div>

<div style><br></div><div style>And then :<br></div><div style><br></div><div style><div>  // read family info</div><div>  med_int numFamilies = MEDnFam(fid, meshName);</div><div>  if(numFamilies < 0){</div><div>    Msg::Error("Could not read MED families");</div>

<div>    return 0;</div><div>  }</div></div></pre><div><br></div>Thank you very much in advance for any help/advice, <div><br></div><div>Best Regards,</div><div><br></div><div>Luca<br><br>
</div>