<div dir="ltr">Dear Paul and Neks,<div><br></div><div>I followed your instructions to extrude a 2D field results to 3D: </div><div>What I did are:</div><div>1. run the 2D simulation to get a 2D ascii field. I actived UPT, <span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">there should be 6 columns. </span>as followed:</div><div><br></div><div><div>1874   6   6   1 6.0000000E+02    0 X Y U P T                      NELT,NX,NY,N</div><div>  0.000      0.000      0.000      0.000      0.000      0.000</div><div>  0.000      0.000      0.000      0.000      0.000      0.000</div><div>  0.000      0.000      0.000      0.000      0.000      0.000</div><div>  0.000      0.000      0.000      0.000      0.000      0.000</div><div>......</div><div>.......     <br></div></div><div><div> 0.000      0.000      0.000      0.000      0.000      0.000</div><div><b><font color="#ff0000"> 0.000      0.000</font></b></div><div>  1.000000E+00  1.550000E+01  9.999999E-01  1.290108E-02  8.919954E-08  0.000000E+00</div><div>  8.825276E-01  1.550000E+01  9.999999E-01  1.281575E-02  1.081013E-07  0.000000E+00</div><div>  6.426157E-01  1.550000E+01  9.999999E-01  1.264053E-02  1.421044E-07  0.000000E+00</div><div>  3.573842E-01  1.550000E+01  1.000000E+00  1.243110E-02  1.706964E-07  0.000000E+00</div></div><div>...........</div><div>..........</div><div><br></div><div>2. <span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"> I get rid the first line,  and the row only contains two-column  </span></div><div><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span><br></span></span></div><div><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span><br></span></span></div><div>3. generated a 3D mesh by n2to3, using 4 levels in 1-spanwise length</div><div><br></div><div>4. setup a 3D simulation by a modified SIZE</div><div><br></div><div><span style="background-color:rgb(238,238,238)">     <font color="#ff0000"> parameter (lelx=1874,lely=1,lelz=4) </font></span><font color="#ff0000">!</font> global tensor mesh dimensions</div><div><br></div><div><br></div><div>5. define usrdat2 as follows</div><div><div>      subroutine usrdat2</div><div>      include 'SIZE'</div><div>      include 'TOTAL'</div><div><br></div><div>       parameter (l2d=lx1*ly1*lelx) ! Number of points in 2D field</div><div>       common /mystuff/ u2d(l2d),v2d(l2d)</div><div><br></div><div><br></div><div><br></div><div>       nelx=lelx</div><div>       nely=lely</div><div>       nelz=lelz</div><div><br></div><div>       if (nid.eq.0) then</div><div>          open(33,file='2D.dat')</div><div>          do i=1,l2d</div><div>                read(33,*) u2d(i),v2d(i)</div><div>          enddo</div><div>          close(33)</div><div>       endif</div><div>       nbytes = 8*l2d</div><div>      call bcast(u2d,nbytes)</div><div>      call bcast(v2d,nbytes)</div><div><br></div><div>       call z_average_transpose(vx,u2d) ! distribute u2d to each plane vx</div><div>       call z_average_transpose(vy,v2d)</div><div><br></div><div>       call outpost(vx,vy,vz,pr,t,'   ')</div><div>       call exitti('Quit in usrdat2.$',nelgt)</div><div><br></div><div><br></div><div>      return</div><div>      end</div></div><div><br></div><div>6 run the 3D simulation, and it appears</div><div><div> call usrdat2</div><div><br></div><div>Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.</div><div><br></div><div>Backtrace for this error:</div><div><br></div><div>Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.</div><div><br></div><div>Backtrace for this error:</div><div><br></div><div>Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.</div><div><br></div><div>Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.</div><div><br></div><div>Backtrace for this error:</div><div><br></div><div>Backtrace for this error:</div><div>#0  0x7FB31EDF4407</div><div>#1  0x7FB31EDF4A1E</div><div>#2  0x7FB31E0F40DF</div><div>#0  0x7F66D6B5B407</div><div>#1  0x7F66D6B5BA1E</div><div>#2  0x7F66D5E5B0DF</div><div>#3  0x4ADA8B in get_exyz_</div><div>#3  0x4ADA8B in get_exyz_</div><div>#4  0x4AE100 in z_average_transpose_</div><div>#4  0x4AE100 in z_average_transpose_</div><div>#5  0x40381B in usrdat2_</div><div>#5  0x40381B in usrdat2_</div><div>#6  0x404017 in nek_init_</div><div>#6  0x404017 in nek_init_</div><div>#0  0x7F9A77280407</div><div>#1  0x7F9A77280A1E</div><div>#2  0x7F9A765800DF</div><div>#3  0x4ADA8B in get_exyz_</div><div>#4  0x4AE100 in z_average_transpose_</div><div>#5  0x40381B in usrdat2_</div><div>#6  0x404017 in nek_init_</div><div>#0  0x7FD15A1FE407</div><div>#1  0x7FD15A1FEA1E</div><div>#2  0x7FD1594FE0DF</div><div>#3  0x4ADA8B in get_exyz_</div><div>#4  0x4AE100 in z_average_transpose_</div><div>#5  0x40381B in usrdat2_</div><div>#6  0x404017 in nek_init_</div><div>-------------------------------------------------------</div><div>Primary job  terminated normally, but 1 process returne</div></div><div><br></div><div>I tried different parameters, unfortunately none of them works, Would anyone like to help me?</div><div><br></div><div>Thanks!</div><div><br></div><div>Wei</div><div><br></div></div>