[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