[Nek5000-users] distribute a 1d profile to 3d domain

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

Hi Mirko,

I once did a similar operation but for duplicating a z-profile (resulted from the planar_average_z() routine) with the following code.

Aleks

c-----------------------------------------------------------------------
      subroutine fill3d_1d_t(u,ua)
c
c     Fill 3d array with 1d t-values by replicating them in r-s planes
c
      include 'SIZE'
      include 'PARALLEL'
      include 'ZPER'  ! define nelx,y,z here!
c
      real u(nx1,ny1,nx1,nelv),ua(nz1,nelz)
      integer e,eg,ex,ey,ez
c
      do e=1,nelt
c
         eg = lglel(e)
         call get_exyz(ex,ey,ez,eg,nelx,nely,nelz)
c
         do k=1,nz1
         do j=1,ny1,node
         do i=1,nx1
            u(i,j,k,e) = ua(k,ez)
         enddo
         enddo
         enddo
      enddo
c
      return
      end
c-----------------------------------------------------------------------




________________________________
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: Monday, April 27, 2015 5:52 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] distribute a 1d profile to 3d domain

Dear nek' users

I've a 1d profile u(y) computed with "planar_average_s" subroutine. So now i want to redistribute this profile to the whole domain in order to obtain a 3d field u(x,y,z) How can i do this?

Thanks

--

Mr. Mirko Farano

Ph.D. student at Politecnico di Bari
Tel.: +39-0883529289 / +39-0805963462
Mob.: +39-3202342719
Department of Mechanics Mathematics and Management -
Fluid Machinery and Energy Systems
Via Re David, 200 - 70125 BARI (Italy)
mirko.farano at poliba.it<mailto:mirko.farano at poliba.it>

In joint supervision with the
Ecole Nationale Supérieure d'Arts et Métiers de Paris - Paristech
Tel.: +33-144246436
Mob.: +33-789822134
Dynfluid Laboratory
151 Boulevard de l'Hôpital - 75013 Paris (France)
mirko.farano at ensam.eu<mailto:mirko.farano at ensam.eu>

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