[Nek5000-users] Moab Issues

[nek5000-users at lists.mcs.anl.gov]

Hi All,


Sorry to post so many questions but I am tearing my hair out! I solved my previous issue just a matter of installing the Zoltan tools and using them to partition my mesh instead.


However I have now another obstacle to overcome.


I keep getting the error message:


Parallel Read times:
  0.0771379 PARALLEL READ PART
  0.00908303 PARALLEL RESOLVE_SHARED_ENTS
  0.000237942 PARALLEL RESOLVE_SHARED_SETS
  0.0863409 PARALLEL TOTAL
 ABORT: nelv is invalid in nekmoab_proc_map
 ABORT: nelv is invalid in nekmoab_proc_map
 nelv, lelv =            0          64
 ABORT: nelv is invalid in nekmoab_proc_map
 nelv, lelv =            0          64
 nelv, lelv =            0          64
 ABORT: nelv is invalid in nekmoab_proc_map
 nelv, lelv =            0          64

call exitt: dying ...


This has been asked before and the problem was identified that moab enabled Nek5000 is not compatible with 2D quads. However I have a 3D mesh, proven by:


   type    count   total
------------------------------
   Edge    282      3.7e+02
   Quad   152      2.6e+02
    Hex     64        1.2e+02
1D Side   1376    1.8e+03
 Vertex    765


A very small mesh I know. So the issue seems to be that nelv is less than 64 and in the moab.f routine this brings the above error (I've tried simply commenting this out unsuccessfully). Thus it would seem logical that one defines nelv to be the correct value.


I have tried this by defining it in the geometry section of the .rea file and this brings even more errors. To my knowledge as my mesh is 3D this subroutine shouldn't even be called.


Any help is appreciated!


Thanks for your time, Friedrich.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20141013/9d1b83e9/attachment.html>