[Nek5000-users] 2D to 3D

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

Dear Paul,

Thank you for your fast response.

In the userdat2 subroutine I've included the same that it was written in 
the post:

       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='file2D.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)


In the SIZE file I've set lelx=number of elements in 2D file.

I've observed that the code crushes due to the z_average_transpose 
subroutine... What's the problem?

I've also tried to read the 2D data file defining the variables in the 
same way that the z_average_transpose subroutine. The code didn't crush, 
but the results were wrong...


       real u2d(nx1,ny1,lelx,lely),v2d(nx1,ny1,lelx,lely)

       if (nid.eq.0) then
          open(33,file='file2D.dat')
        do e=1,nelt
           eg = lglel(e)
           call get_exyz(ex,ey,ez,eg,nelx,nely,nelz)
           do j=1,ny1
           do i=1,nx1
              do k=1,nz1
                 read(33,*) u2d(i,j,ex,ey), v2d(i,j,ex,ey)
              enddo
           enddo
           enddo
        enddo
        close(33)
        end if


Thank you for your help.
SL






El 24-09-2014 15:07, nek5000-users at lists.mcs.anl.gov escribió:
> Dear SL,
> 
> It depends on what is in the user file and specifically in your user
> dat2 subroutine
> (given that that appears to be the point of failure from the
> information in your email).
> 
> Did you upgrade this routine to reflect the fact that your simulation 
> is 3D ?
> 
> Paul
> 
> ________________________________________
> From: nek5000-users-bounces at lists.mcs.anl.gov
> [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of
> nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
> Sent: Wednesday, September 24, 2014 8:01 AM
> To: Nek5000 Users
> Subject: [Nek5000-users] 2D to 3D
> 
> Hi Nek's
> 
> 
> Following this post:
> 
> http://lists.mcs.anl.gov/pipermail/nek5000-users/2013-October/002333.html
> 
> I'm trying to extrude a 2D file to a 3D file.
> 
> First, I've run a 2D case. Then I've used n2to3 to generate the new 3D
> .rea file, called: new3d.rea. I've launch a new simulation with
> new3d.rea and everything works fine. However, when I include the
> information about the 3D extrusion in the usrdat2 subroutine, I get the
> following error message if I try to run the code. It seems something
> related to the z_average_transpose subroutine. What is wrong? Thanks
> 
> 
>   call usrdat
>   done :: usrdat
> 
>   generate geometry data
>   vol_t,vol_v:   1749.7930352474154        1749.7930352474154
>   done :: generate geometry data
> 
>   call usrdat2
> 
> =====================================================================================
> =   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
> =   EXIT CODE: 8
> =   CLEANING UP REMAINING PROCESSES
> =   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
> =====================================================================================
> APPLICATION TERMINATED WITH THE EXIT STRING: Floating point exception
> (signal 8)
> 
> 
> 
> SL
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users