[Nek5000-users] Moab Issues

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

Hi Friedriech,

It appears that Nek-MOAB does not pickup your elements since nelv=0.

One thing to check is whether your block-ids/material sets are specified correctly in the mesh files .cub & .h5m and in the .rea file -- e.g. the following lines in

nek5_svn/examples/moab_conjht/moab_conjht.rea:

1 2      1 fluid set, 2 other/solid sets
20 10 30      fluid set #20, solid sets #10, 30

correspond to a Cubit/MOAB mesh with 3 block-ids/material sets: 10,20, and 30 which matches the  check with 'mbsize -m' below.


If the problem still persists, please contact me offline at obabko at mcs.nal.gov

Aleks

mbsize -m moab_conjht.h5m

...

Material Set 10:
   type count   total                   minimum                   average                       rms                   maximum                  std.dev.
------- ----- ------- ------------------------- ------------------------- ------------------------- ------------------------- -------------------------
    Hex     1       1                         1                         1                         1                         1                         0
1D Side    12      12                         1                         1                         1                         1                         0
 Vertex     0

Material Set 20:
....
Material Set 30:
...

________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
Sent: Monday, October 13, 2014 11:55 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Moab Issues


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/20141014/db2adea9/attachment.html>