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

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

Hi Neks,

I followed what Fisher posted in the below link and tried to extrude a 2D flow into 3D.

	http://lists.mcs.anl.gov/pipermail/nek5000-users/2013-October/002333.html

85 c-----------------------------------------------------------------------
86       subroutine usrdat2
87       include 'SIZE'
88       include 'TOTAL'
89 
90 
91       nelx=lelx,nely=lely,nelz=lelz
92       parameter (l2d=lx1*ly1*lelx)
93       common /mystuff/ u2d(l2d),v2d(l2d)
94 
95        if (nid.eq.0) then
96           open(33,file='2D.dat')
97           do i=1,l2d
98              read(33,*) u2d(i),v2d(i)
99           enddo
100           close(33)
101        endif
102        nbytes = 8*l2d
103        call bcast(u2d,nbytes)
104        call bcast(v2d,nbytes)
105 
106        call z_average_transpose(vx,u2d)
107        call z_average_transpose(vy,v2d)
108 
109        call outpost(vx,vy,vz,pr,t,'3d.dat')
110        call exitti('Quit in usrdat2.$',nelgt)
111 
112       return
113       end
114 c-----------------------------------------------------------------------
115       subroutine usrdat3


when I compiling the simulation ( with a 2D.rea and new SIZE) it displays errors,

	nek5000/test/test3dimposed/3d/box.f:91.8:
	
      		 nelx=lelx,nely=lely,nelz=lelz                                   
       	1
	Error: Unclassifiable statement at (1)
	make: *** [nek5000] Error 1

Is there something I misunderstand ?

Thanks! 

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