[Nek5000-users] Periodic B.C and Genmap Issue

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu Mar 10 08:46:43 CST 2016


Hi Paul,

It is looks like that the n2to3 not generate the 3D  mesh with periodic
boundary conditions correctly by rotating.  First I tried to generate a 3D
mesh, then after applied genmap it appears the similar error.

Following the discussions some years ago, it seems I need to update the
genmap.f,  Have you update the new version, or Would you like to post the
genmap.f again?
thanks!

Wei

n2to3


 This is the code that establishes proper

 element connectivities, as of Oct., 2006.


 Input old (source) file name:

box

 Input new (output) file name:

boxss

 For ASCII output only: 0; For .rea/.re2: 1

0

 input number of levels: (1, 2, 3,...; < 0 for circular sweep.):

-20

 input z min:

0

 input z max:

90

 input gain (0=custom,1=uniform,other=geometric spacing):

1

 This is for CEM: yes or no:

no

 Enter Z(5) FLUID boundary condition (P,v,O,ect):

P

 this is FLUID cbz: P  P   <---

 F F           0           1  ifheat,ifmhd


         360 elements written to boxss.rea


genmap

Input (.rea) file name:

boxss

Input mesh tolerance (default 0.2):

NOTE: smaller is better, but generous is more forgiving for bad meshes.

0.1

 reading .rea file data ...

 start locglob_lexico:           8         360        2880
0.10000000000000001

 locglob:           1           1        2880

 locglob:           2          59        2880

 locglob:           3          84        2880

 locglob:           1         588        2880

 locglob:           2         588        2880

 locglob:           3         588        2880

 locglob:           1         588        2880

 locglob:           2         588        2880

 locglob:           3         588        2880

 done locglob_lexico:         588         588        2880           8

 start periodic vtx:         360         588

     1     4     4  Matrix:    x0

  1   x0    1.2084E-06  1.2084E-06  1.6729E-06  1.6729E-06

  2   x0   -1.2362E-03 -1.2362E-03  1.2363E-03  1.2363E-03

  3   x0   -0.0000E+00 -0.0000E+00 -0.0000E+00 -0.0000E+00

  4   x0   -2.6465E-03  2.6465E-03  2.6465E-03 -2.6465E-03

   343     4     4  Matrix:    x1

1   x1    1.2084E-06  1.6729E-06  1.6729E-06  1.2084E-06

  2   x1    5.4038E-11 -5.4038E-11 -5.4038E-11  5.4038E-11

  3   x1    1.2362E-03 -1.2363E-03 -1.2363E-03  1.2362E-03

  4   x1   -2.6465E-03 -2.6465E-03  2.6465E-03  2.6465E-03

     1     4     4  Matrix:    z0

  1   z0    1.2084E-06  1.2084E-06  1.6729E-06  1.6729E-06

  2   z0    1.4000E-02  1.4000E-02  1.6473E-02  1.6473E-02

  3   z0   -0.0000E+00 -0.0000E+00 -0.0000E+00 -0.0000E+00

  4   z0   -1.1000E-01 -1.0471E-01 -1.0471E-01 -1.1000E-01

   343     4     4  Matrix:    z1

  1   z1    1.2084E-06  1.6729E-06  1.6729E-06  1.2084E-06

  2   z1   -6.1196E-10 -7.2004E-10 -7.2004E-10 -6.1196E-10

  3   z1   -1.4000E-02 -1.6473E-02 -1.6473E-02 -1.4000E-02

  4   z1   -1.1000E-01 -1.1000E-01 -1.0471E-01 -1.0471E-01

           1 5  3  1.00000000E-03  1.10000000E-01 abort: FACE MATCH FAIL

         343 6  5  1.10000000E-04  3.49664311E-03 abort: FACE MATCH FAIL


           0  quit



====================================================

Hi,

Yes, got it now. I don't know why it wasn't working either. I use the
same thing to compress files, either way it's good now. Thanks Paul.

- Michael

----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
<https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>
To: nek5000-users at lists.mcs.anl.gov
<https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>
Sent: Tuesday, June 22, 2010 10:40:57 AM GMT -06:00 Guadalajara /
Mexico City / Monterrey
Subject: Re: [Nek5000-users] Periodic BC & Genmap issue



Hi Michael,

not sure why gunzip didn't work - it worked here... but I've
put the unzipped version in:
www.mcs.anl.gov/~fischer/genmap.f

Let me know if that works for you.

Cheers,

Paul
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20160310/06594a49/attachment-0001.html>


More information about the Nek5000-users mailing list