[Nek5000-users] Extruding 2D .fld file to 3D

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Sun Jul 8 01:30:06 CDT 2018


Dear Paul and Neks,

I followed your instructions to extrude a 2D field results to 3D:
What I did are:
1. run the 2D simulation to get a 2D ascii field. I actived UPT, there
should be 6 columns. as followed:

1874   6   6   1 6.0000000E+02    0 X Y U P T
NELT,NX,NY,N
  0.000      0.000      0.000      0.000      0.000      0.000
  0.000      0.000      0.000      0.000      0.000      0.000
  0.000      0.000      0.000      0.000      0.000      0.000
  0.000      0.000      0.000      0.000      0.000      0.000
......
.......
 0.000      0.000      0.000      0.000      0.000      0.000
* 0.000      0.000*
  1.000000E+00  1.550000E+01  9.999999E-01  1.290108E-02  8.919954E-08
0.000000E+00
  8.825276E-01  1.550000E+01  9.999999E-01  1.281575E-02  1.081013E-07
0.000000E+00
  6.426157E-01  1.550000E+01  9.999999E-01  1.264053E-02  1.421044E-07
0.000000E+00
  3.573842E-01  1.550000E+01  1.000000E+00  1.243110E-02  1.706964E-07
0.000000E+00
...........
..........

2.  I get rid the first line,  and the row only contains two-column


3. generated a 3D mesh by n2to3, using 4 levels in 1-spanwise length

4. setup a 3D simulation by a modified SIZE

      parameter (lelx=1874,lely=1,lelz=4) ! global tensor mesh dimensions


5. define usrdat2 as follows
      subroutine usrdat2
      include 'SIZE'
      include 'TOTAL'

       parameter (l2d=lx1*ly1*lelx) ! Number of points in 2D field
       common /mystuff/ u2d(l2d),v2d(l2d)



       nelx=lelx
       nely=lely
       nelz=lelz

       if (nid.eq.0) then
          open(33,file='2D.dat')
          do i=1,l2d
                read(33,*) u2d(i),v2d(i)
          enddo
          close(33)
       endif
       nbytes = 8*l2d
      call bcast(u2d,nbytes)
      call bcast(v2d,nbytes)

       call z_average_transpose(vx,u2d) ! distribute u2d to each plane vx
       call z_average_transpose(vy,v2d)

       call outpost(vx,vy,vz,pr,t,'   ')
       call exitti('Quit in usrdat2.$',nelgt)


      return
      end

6 run the 3D simulation, and it appears
 call usrdat2

Program received signal SIGFPE: Floating-point exception - erroneous
arithmetic operation.

Backtrace for this error:

Program received signal SIGFPE: Floating-point exception - erroneous
arithmetic operation.

Backtrace for this error:

Program received signal SIGFPE: Floating-point exception - erroneous
arithmetic operation.

Program received signal SIGFPE: Floating-point exception - erroneous
arithmetic operation.

Backtrace for this error:

Backtrace for this error:
#0  0x7FB31EDF4407
#1  0x7FB31EDF4A1E
#2  0x7FB31E0F40DF
#0  0x7F66D6B5B407
#1  0x7F66D6B5BA1E
#2  0x7F66D5E5B0DF
#3  0x4ADA8B in get_exyz_
#3  0x4ADA8B in get_exyz_
#4  0x4AE100 in z_average_transpose_
#4  0x4AE100 in z_average_transpose_
#5  0x40381B in usrdat2_
#5  0x40381B in usrdat2_
#6  0x404017 in nek_init_
#6  0x404017 in nek_init_
#0  0x7F9A77280407
#1  0x7F9A77280A1E
#2  0x7F9A765800DF
#3  0x4ADA8B in get_exyz_
#4  0x4AE100 in z_average_transpose_
#5  0x40381B in usrdat2_
#6  0x404017 in nek_init_
#0  0x7FD15A1FE407
#1  0x7FD15A1FEA1E
#2  0x7FD1594FE0DF
#3  0x4ADA8B in get_exyz_
#4  0x4AE100 in z_average_transpose_
#5  0x40381B in usrdat2_
#6  0x404017 in nek_init_
-------------------------------------------------------
Primary job  terminated normally, but 1 process returne

I tried different parameters, unfortunately none of them works, Would
anyone like to help me?

Thanks!

Wei
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180708/4cba684b/attachment.html>


More information about the Nek5000-users mailing list