<br><font size=2 face="sans-serif"><br>
Hello,</font>
<br>
<br><font size=2 face="sans-serif">I would like to report 2 issues with
MED output format (in Gmsh 2.4.2 and all older versions I have tested)</font>
<br>
<br><font size=2 face="sans-serif">- elements are not oriented correctly,
meaning they probably do not follow the MED element numbering conventions.</font>
<br>
<br><font size=2 face="sans-serif">- Families in MED files do not contain
groups or attributes (attributes should be deprecated soon, so groups are
best).</font>
<br>
<br><font size=2 face="sans-serif">In GModel_IO_MED.cpp, we have  of
the following form:</font>
<br>
<br><font size=2 face="Courier New">if(MEDfamCr(fid, meshName, (char*)familyName.c_str(),
</font>
<br><font size=2 face="Courier New">         
  (med_int)num, 0, 0, 0, 0, (char*)groupName.c_str(),</font>
<br><font size=2 face="Courier New">         
  (med_int)entities[i]->physicals.size()) < 0)</font>
<br><font size=2 face="sans-serif"><br>
The last argument of MEDfamCr should be the number of associated groups,
which is always 1</font>
<br><font size=2 face="sans-serif">if families are mapped 1-to-1 with physical
entities, so I beleive the call should be :</font>
<br>
<br><font size=2 face="Courier New">if(MEDfamCr(fid, meshName, (char*)familyName.c_str(),
</font>
<br><font size=2 face="Courier New">         
  (med_int)num, 0, 0, 0, 0, (char*)groupName.c_str(),</font>
<br><font size=2 face="Courier New">         
  (med_int)entities[i]->physicals.size()) > 0)</font>
<br>
<br>
<br><font size=2 face="sans-serif">(i.e. the sign is wrong in the test
for the original version).</font>
<br>
<br><font size=2 face="sans-serif">Normally, even if entities[i]->physicals.size()
== 0, Families should contain a group, so</font>
<br><font size=2 face="sans-serif">maybe another wy of defining the group
should be sued in this case, and the call</font>
<br><font size=2 face="sans-serif">would simply be :</font>
<br>
<br><font size=2 face="Courier New">if(MEDfamCr(fid, meshName, (char*)familyName.c_str(),
</font>
<br><font size=2 face="Courier New">         
  (med_int)num, 0, 0, 0, 0, (char*)groupName.c_str(), 1)</font>
<br>
<br><font size=2 face="sans-serif">I have not tested the patch or built
a complete GMSH compilation environment,</font>
<br><font size=2 face="sans-serif">but I beleive this could fix the second
bug (i.e. missing group definitons) with MED output.</font>
<br>
<br><font size=2 face="sans-serif">I may try to look into element node
numbering for the first bug in the future, but</font>
<br><font size=2 face="sans-serif">this may have to wait a bit.</font>
<br>
<br><font size=2 face="sans-serif">Best regards,</font>
<br>
<br><font size=2 face="sans-serif">        Yvan
Fournier</font>
<br><font size=2 face="sans-serif">        EDF
R&D</font>
<br><font size=2 face="sans-serif">        Fluid
Dynamcis, Power Generation, and Environment</font>
<br><font size=2 face="sans-serif">        6
quai Watier - BP 49</font>
<br><font size=2 face="sans-serif">        78401
Chatou Cedex</font>
<br><font size=2 face="sans-serif">        FRANCE<br>
</font>
<table>
<tr valign=top>
<td rowspan=2>
<td>
<tr valign=top>
<td>
<tr>
<td valign=top>
<td></table>
<br>