[Nek5000-users] Nek5000+MOAB

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

Hi NEKs,

I am trying to use the MOAB enabled version of NEK.

I looked at the pipe example in /examples/moab and it works.

Than I recreated the .cub mesh using the journal file "pipe.jou" and Trelis.
I converted the mesh using mbconvert blah.cub blah.h5m and it worked, even
if there was a worrysome warning message " Element ids are not contiguous!".
This is the mbsize result of the original input.h5m mesh:

type  count   total                 minimum                 average
            rms                 maximum                std.dev.
------- ------ ------- ----------------------- -----------------------
----------------------- ----------------------- -----------------------
   Edge   4952 5.1e+02                0.034623                 0.10396
            0.11844                 0.19893                0.056762
   Quad   2476      20               0.0017335               0.0082427
          0.0086523                0.036989               0.0026306
    Hex   5496     4.7              0.00014971              0.00085585
          0.0023197                0.018287               0.0021561
1D Side  75856 7.8e+03                0.019193                  0.1028
            0.13066                 0.50001                 0.08065
 Vertex  49005

an this is the same result for the new mesh:

   type  count   total                 minimum                 average
                rms                 maximum                std.dev.
------- ------ ------- ----------------------- -----------------------
----------------------- ----------------------- -----------------------
   Edge   4952 5.1e+02                0.034623                 0.10396
            0.11844                 0.19886                0.056762
   Quad   2476      20               0.0017335               0.0082427
          0.0086523                0.036969               0.0026306
    Hex   5496     4.7              0.00014972              0.00085585
          0.0023244                0.018468               0.0021611
1D Side  75856 8.3e+03                0.019195                    0.11
            0.14682                     0.5                0.097238
 Vertex  49005

If i try to use the new mesh in NEK I observe the following error:

read .rea file
 ABORT: nelv is invalid in nekmoab_proc_map
 nelv, lelv =            0        5496


Probably it is more a MOAB related problem, but if anybody can give me any
hints it will be great.

Bests,

Andrea
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140923/f2e1ef15/attachment.html>