July 2014 - Nek5000-users - mailman3.cels.anl.gov

Re: [Nek5000-users] Vorticity
by nek5000-users＠lists.mcs.anl.gov 11 Oct '16

11 Oct '16

Hi, Thanks for your quick reply. Just I wanna know if I can use the interpolated values of velocity at some arbitrary points in the domain (not necessarily the grid points) as the input arguments to theses subroutines. Regards, Alireza Quoting nek5000-users-request(a)lists.mcs.anl.gov: > Send Nek5000-users mailing list submissions to > nek5000-users(a)lists.mcs.anl.gov > > To subscribe or unsubscribe via the World Wide Web, visit > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > or, via email, send a message with subject or body 'help' to > nek5000-users-request(a)lists.mcs.anl.gov > > You can reach the person managing the list at > nek5000-users-owner(a)lists.mcs.anl.gov > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Nek5000-users digest..." > > > Today's Topics: > > 1. Vorticity (nek5000-users(a)lists.mcs.anl.gov) > 2. Re: Vorticity (nek5000-users(a)lists.mcs.anl.gov) > 3. Re: Vorticity (nek5000-users(a)lists.mcs.anl.gov) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 13 Jul 2011 17:55:35 -0400 > From: nek5000-users(a)lists.mcs.anl.gov > Subject: [Nek5000-users] Vorticity > To: nek5000-users(a)lists.mcs.anl.gov > Message-ID: <20110713175535.1042430oct81vz2v(a)webmail.vt.edu> > Content-Type: text/plain; charset=ISO-8859-1; DelSp="Yes"; > format="flowed" > > Hi, > > I need to calculate the vorticity vector and the rate of strain tensor > at any arbitrary point in the domain. I wonder if these values are > computed somewhere in the code or if I should interpolate them in the > usr file. I appreciate any kind of help. > > Thanks in advance, > > > Regards, > Alireza Karimi > > > > ------------------------------ > > Message: 2 > Date: Wed, 13 Jul 2011 17:09:48 -0500 (CDT) > From: nek5000-users(a)lists.mcs.anl.gov > Subject: Re: [Nek5000-users] Vorticity > To: nek5000-users(a)lists.mcs.anl.gov > Message-ID: <Pine.LNX.4.64.1107131704140.26989(a)jacobi.mcs.anl.gov> > Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed > > Hi Alireza, > > For vorticity computation, you can call > > call comp_vort3(vort,work1,work2,vx,vy,vz) > > from userchk with the appropriate definitions of vorticity and work arrays > like > > common /scren/ vort (lx1,ly1,lz1,lelv,3) ! x y z components > $ , work1(lx1,ly1,lz1,lelv) > $ , work2(lx1,ly1,lz1,lelv) > > Let me check what we have for the calculation of strain tensor. > Best, > Aleks > > > On Wed, 13 Jul 2011, nek5000-users(a)lists.mcs.anl.gov wrote: > >> Hi, >> >> I need to calculate the vorticity vector and the rate of strain >> tensor at any >> arbitrary point in the domain. I wonder if these values are computed >> somewhere in the code or if I should interpolate them in the usr file. I >> appreciate any kind of help. >> >> Thanks in advance, >> >> >> Regards, >> Alireza Karimi >> >> _______________________________________________ >> Nek5000-users mailing list >> Nek5000-users(a)lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > > ------------------------------ > > Message: 3 > Date: Wed, 13 Jul 2011 23:59:08 -0500 (CDT) > From: nek5000-users(a)lists.mcs.anl.gov > Subject: Re: [Nek5000-users] Vorticity > To: nek5000-users(a)lists.mcs.anl.gov > Message-ID: > <600285192.94128.1310619548370.JavaMail.root(a)zimbra.anl.gov> > Content-Type: text/plain; charset=utf-8 > > Hi Alireza, > > > For a strain stress computation you can call subroutine comp_sij > > call comp_sij(sij,nij,vx,vy,vz > $ ,ur,us,ut,vr,vs,vt,wr,ws,wt) > > > (navier5.f:2225) with nij=6 for 3D somewhere from userchk once you > declare the stress component array sij and arrays with local > coordinate r-s-t derivatives of velocity components > > common /csij/ > $ sij(lx1*ly1*lz1,6,lelv) ! 6 components > $ ,ur (lx1*ly1*lz1) > ... > $ ,wt (lx1*ly1*lz1) > > > Best, > Aleks > > > > > > > ----- Original Message ----- > From: nek5000-users(a)lists.mcs.anl.gov > To: nek5000-users(a)lists.mcs.anl.gov > Sent: Wednesday, July 13, 2011 5:09:48 PM > Subject: Re: [Nek5000-users] Vorticity > > Hi Alireza, > > For vorticity computation, you can call > > call comp_vort3(vort,work1,work2,vx,vy,vz) > > from userchk with the appropriate definitions of vorticity and work arrays > like > > common /scren/ vort (lx1,ly1,lz1,lelv,3) ! x y z components > $ , work1(lx1,ly1,lz1,lelv) > $ , work2(lx1,ly1,lz1,lelv) > > Let me check what we have for the calculation of strain tensor. > Best, > Aleks > > > On Wed, 13 Jul 2011, nek5000-users(a)lists.mcs.anl.gov wrote: > >> Hi, >> >> I need to calculate the vorticity vector and the rate of strain >> tensor at any >> arbitrary point in the domain. I wonder if these values are computed >> somewhere in the code or if I should interpolate them in the usr file. I >> appreciate any kind of help. >> >> Thanks in advance, >> >> >> Regards, >> Alireza Karimi >> >> _______________________________________________ >> Nek5000-users mailing list >> Nek5000-users(a)lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > _______________________________________________ > Nek5000-users mailing list > Nek5000-users(a)lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > > ------------------------------ > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users(a)lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > > End of Nek5000-users Digest, Vol 29, Issue 2 > ******************************************** > -- Alireza Karimi PhD Candidate Department of Engineering Science and Mechanics Virginia Polytechnic Institute and State University Blacksburg, VA 24061

1 3

Error when creating own mesh
by nek5000-users＠lists.mcs.anl.gov 01 Mar '16

01 Mar '16

Dear Neks, I am have attempted to use nek5000, firstly by following some of the examples etc. Now I am trying to create my own simulation, which is based upon the ocyl examples. My efforts simply include creating my own grid points in the .rea file to create a slightly different dimension... I believe I have understood how to define the problem however I continually receive errors when creating the genmap it reads as follows.... Input (.rea) file name: ocyl Input mesh tolerance (default 0.2): NOTE: smaller is better, but generous is more forgiving for bad meshes. 0.2 reading .rea file data ... ERROR: error reading 1 11 24 aborting 510 in routine rdbdry. 1 quit All I have done again is change the dimensions of the elements in the .rea file so I do not see why it would now not work? Thank you for any help in advance, Friedrich.

1 7

About governing equations
by nek5000-users＠lists.mcs.anl.gov 27 Aug '14

27 Aug '14

Hi, I read in website that NEK5000 can't solve fully compressible flows. How it changes density in low mach number flows (in which we have thermal dilatation effects)? Thanks

1 1

SFD error message from p66
by nek5000-users＠lists.mcs.anl.gov 12 Aug '14

12 Aug '14

Dear Nek's, I'm trying to program a Selective Frequency Damping following the article of Akervik et. al (Physics of Fluids, 2006) and some Nek5000 posts. I'm able to compile the code, however, when I try to run my simulation, I find the following error message: 0 0.0000E+00 Write checkpoint: 0 0.0000E+00 OPEN: At line 335 of file User/nek5_svn/trunk/nek/prepost.f (unit = 24, file = '') Fortran runtime error: File '' does not exist It seems to be something related to param 66. I've tried with p66=6, 0, 1, -1, but I always get the same error. Do you know what happens? Thanks SL

1 5

Re: [Nek5000-users] Specifying boundary condition from external text file
by nek5000-users＠lists.mcs.anl.gov 28 Jul '14

28 Jul '14

Hello again, I want to provide more information, and more detail of my exact case to hopefully make it easier to get assistance. Using the 2D blasius example as a base, I have created a mesh file "box.box" that allow me to modify the domain. After renaming the necessary files (program is now called 'box') this works fine with the an unmodified blasius.usr (now 'box.usr') code. The next step I took was to change the boundary conditions in userbc() [and useric()] in box.usr. As I understand, the nodes on the inlet plane are called sequentially from bottom to top, by element number ('e') and then by node number ('iy'). The elements are not numbered sequentially from bottom to top, but rather left-to-right, then back to the start, then up, then left-to-right again. Therefore for a domain of 410 elements in the streamwise (x) direction, the element numbers along the inlet are, from bottom to top: 1, 411, 821, etc. There are 12 elements in the wall normal (y) direction. I have an array ('ArrangedInlet') in the subroutine containing the horizontal velocity component, ux, value for each of the nodes along the inlet according to their y coordinate The original data lives in the file called 'InletArranged.out'. This is output by a process in MATLAB so there is freedom in what is included and how it is arranged and formatted. The data points are currently entered manually. This is initially to test, but it is intended that this is read automatically from the external file. However in order to assign the correct velocity value to the correct position along the inlet, I arranged the data in the array such that the rows counted iy from 1 to 8 (lx1) and each is addressing a node: Row 1=Node 1; Row 2 = Node 2, etc. where the nodes go bottom-to-top (y=0 to y=10), The columns contained the data for each element along the inlet - ordered from bottom to top. At the moment each row . Using the same domain as described above (NELX, NELY)=(410, 12): column 1 contains data at each node ('iy') for element 1, column 2 contains data at each node ('iy') for element 411, column 3 contains data at each node ('iy') for element 821, and so on. I got around this by using the following transformation: COL=((e-1)/NELX)+1 where: COL is the column number of the array (1 to 12); e = gllel(eg) is the element number; NELX = 410 is the number of elements in the streamwise (x) direction. This is then read by the program by: ux = ArrangedInlet(iy,COL) And it works! ... but only for a serial case. Also, it currently only works for the single set of data. As soon as I run the case in parallel the addressing changes and this method does not work. The element addressing no longer appears to have a discernible pattern. *So my questions are primarily:* *1) How can I get the values onto the inlet and make it work in parallel? * and *2) Is there a more elegant/efficient method of doing this? IE. a way of calling the inlet nodes only?* I have attached my case "box_2D_Inlet.tar.gz". I have tested it and it does work, however currently it bombs out after 9 iterations, but I'm really only testing the inlet parameters. Looking forward to any assistance. Thanks, Tom. On 17 July 2014 10:37, Thomas Yankos <s3286374(a)student.rmit.edu.au> wrote: > Hi All, > > I'm using NEK5000 for a university project. I've found it to be blissfully > fast and relatively easy to use. > > The problem I'm doing is a blasius boundary layer case with perturbations > (simulating Tollmien–Schlichting waves), initially in 2D as per the blasius > example, then eventually moving to 3D. > I've so far been successful in running the blasius example case, and > modifying the mesh, polynomial order and other input parameters. > > Now however, I want to change the way the boundary conditions are defined. > I have input files that are unique for each time step specifying velocity > and pressures along a line (the inlet) in text format. It would be possible > to re-format these as .csv or another file if necessary. > > *How can I get NEK to read these files as boundary conditions on the inlet > for each time step?* > > I understand it will be a specification in the user userbc subroutine in > blasius.usr however I'm not sure how to actually implement this. > Doing it in 2D will be a precursor to a 3D simulation with the file > specifying the inlet *plane *in a text file, rather than just a *line*. > > Looking forward to a response. > > Warm Regards, > Tom. >

1 1

makenek file
by nek5000-users＠lists.mcs.anl.gov 24 Jul '14

24 Jul '14

hi all, Just beginning with nek. Is it mandatory to have the makenek script in every case folder or is there a way to source this script and place it in the executable path?(along with nek5_svn/trunk/tools/scripts ?) Thanks. Thangam.

1 1

Re: [Nek5000-users] nek installation and running errors
by nek5000-users＠lists.mcs.anl.gov 23 Jul '14

23 Jul '14

Apparently I did not place the /bin in the executable path. now it works fine! PATH=$PATH:~/.bin export PATH Thanks. Thangam. From: thy_horizon(a)hotmail.com To: nek5000-users(a)lists.mcs.anl.gov Subject: nek installation and running errors Date: Wed, 23 Jul 2014 09:17:07 +0530 Thanks Paul. The issue is my system doesnt recognize the commands 'nek' and 'nekb' etc. Do I have to set the environment variables by sourcing the bash file? Thanks. Thangam. From: nek5000-users-request(a)lists.mcs.anl.gov Subject: Nek5000-users Digest, Vol 65, Issue 13 To: nek5000-users(a)lists.mcs.anl.gov Date: Tue, 22 Jul 2014 12:00:08 -0500 Send Nek5000-users mailing list submissions to nek5000-users(a)lists.mcs.anl.gov To subscribe or unsubscribe via the World Wide Web, visit https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users or, via email, send a message with subject or body 'help' to nek5000-users-request(a)lists.mcs.anl.gov You can reach the person managing the list at nek5000-users-owner(a)lists.mcs.anl.gov When replying, please edit your Subject line so it is more specific than "Re: Contents of Nek5000-users digest..." --Forwarded Message Attachment-- From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov Date: Mon, 21 Jul 2014 19:08:13 +0200 Subject: [Nek5000-users] Extruding 2D .fld file to 3D 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 --Forwarded Message Attachment-- From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov Date: Mon, 21 Jul 2014 18:53:10 -0400 Subject: [Nek5000-users] About governing equations Hi, I read in website that NEK5000 can't solve fully compressible flows. How it changes density in low mach number flows (in which we have thermal dilatation effects)? Thanks --Forwarded Message Attachment-- From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov Date: Tue, 22 Jul 2014 15:32:01 +0530 Subject: [Nek5000-users] nek installation and running errors Hi all, I am trying to download and install NEK5000 and trying to get a hang of it. I have downloaded and tried to the first example from the 'quick start' page. The compilation is successful after the makenek command. However, it does not realize further nek commands such as nek, nekb etc. (command not found error)I went on to issue the command ./nek5000 eddy_uv and following was the output. Number of processors: 1 REAL wdsize : 8 INTEGER wdsize : 4 No file SESSION.NAME; using defaults of PATH=. and SESSION=NEK. Beginning session: ./NEK.rea timer accuracy: 0.0000000E+00 sec read .rea file At line 32 of file /home/thangam/nek5_svn/trunk/nek/connect2.f (unit = 9, file = '') Fortran runtime error: File './NEK.rea' does not exist I have put a copy of the scripts under the /bin directory as well. I dont understand what I am missing and where im going wrong! Thanks. Thangam. --Forwarded Message Attachment-- From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov Date: Tue, 22 Jul 2014 10:10:12 -0500 Subject: Re: [Nek5000-users] nek installation and running errors You should use one of the scripts in trunk/tools/scripts The usual command would be nek blah where blah is the name of your .rea file. To run in background (strongly recommended), nekb blah For mpi there are nekmpi and nekbmpi nekmpi blah 4 would run on 4 ranks, etc. There is some description of this in the manual. Paul On Tue, 22 Jul 2014, nek5000-users(a)lists.mcs.anl.gov wrote: > Hi all, > > > I am trying to download and install NEK5000 and trying to get a hang of it. I have downloaded and tried to the first example from the 'quick start' page. The compilation is successful after the makenek command. However, it does not realize further nek commands such as nek, nekb etc. (command not found error)I went on to issue the command ./nek5000 eddy_uv and following was the output. > > Number of processors: 1 > REAL wdsize : 8 > INTEGER wdsize : 4 > No file SESSION.NAME; using defaults of > PATH=. and SESSION=NEK. > > > Beginning session: > ./NEK.rea > > > timer accuracy: 0.0000000E+00 sec > > read .rea file > At line 32 of file /home/thangam/nek5_svn/trunk/nek/connect2.f (unit = 9, file = '') > Fortran runtime error: File './NEK.rea' does not exist > > > > > I have put a copy of the scripts under the /bin directory as well. > > I dont understand what I am missing and where im going wrong! > > Thanks. > Thangam. >

1 0

nek installation and running errors
by nek5000-users＠lists.mcs.anl.gov 22 Jul '14

22 Jul '14

Thanks Paul. The issue is my system doesnt recognize the commands 'nek' and 'nekb' etc. Do I have to set the environment variables by sourcing the bash file? Thanks. Thangam. From: nek5000-users-request(a)lists.mcs.anl.gov Subject: Nek5000-users Digest, Vol 65, Issue 13 To: nek5000-users(a)lists.mcs.anl.gov Date: Tue, 22 Jul 2014 12:00:08 -0500 Send Nek5000-users mailing list submissions to nek5000-users(a)lists.mcs.anl.gov To subscribe or unsubscribe via the World Wide Web, visit https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users or, via email, send a message with subject or body 'help' to nek5000-users-request(a)lists.mcs.anl.gov You can reach the person managing the list at nek5000-users-owner(a)lists.mcs.anl.gov When replying, please edit your Subject line so it is more specific than "Re: Contents of Nek5000-users digest..." --Forwarded Message Attachment-- From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov Date: Mon, 21 Jul 2014 19:08:13 +0200 Subject: [Nek5000-users] Extruding 2D .fld file to 3D 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 --Forwarded Message Attachment-- From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov Date: Mon, 21 Jul 2014 18:53:10 -0400 Subject: [Nek5000-users] About governing equations Hi, I read in website that NEK5000 can't solve fully compressible flows. How it changes density in low mach number flows (in which we have thermal dilatation effects)? Thanks --Forwarded Message Attachment-- From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov Date: Tue, 22 Jul 2014 15:32:01 +0530 Subject: [Nek5000-users] nek installation and running errors Hi all, I am trying to download and install NEK5000 and trying to get a hang of it. I have downloaded and tried to the first example from the 'quick start' page. The compilation is successful after the makenek command. However, it does not realize further nek commands such as nek, nekb etc. (command not found error)I went on to issue the command ./nek5000 eddy_uv and following was the output. Number of processors: 1 REAL wdsize : 8 INTEGER wdsize : 4 No file SESSION.NAME; using defaults of PATH=. and SESSION=NEK. Beginning session: ./NEK.rea timer accuracy: 0.0000000E+00 sec read .rea file At line 32 of file /home/thangam/nek5_svn/trunk/nek/connect2.f (unit = 9, file = '') Fortran runtime error: File './NEK.rea' does not exist I have put a copy of the scripts under the /bin directory as well. I dont understand what I am missing and where im going wrong! Thanks. Thangam. --Forwarded Message Attachment-- From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov Date: Tue, 22 Jul 2014 10:10:12 -0500 Subject: Re: [Nek5000-users] nek installation and running errors You should use one of the scripts in trunk/tools/scripts The usual command would be nek blah where blah is the name of your .rea file. To run in background (strongly recommended), nekb blah For mpi there are nekmpi and nekbmpi nekmpi blah 4 would run on 4 ranks, etc. There is some description of this in the manual. Paul On Tue, 22 Jul 2014, nek5000-users(a)lists.mcs.anl.gov wrote: > Hi all, > > > I am trying to download and install NEK5000 and trying to get a hang of it. I have downloaded and tried to the first example from the 'quick start' page. The compilation is successful after the makenek command. However, it does not realize further nek commands such as nek, nekb etc. (command not found error)I went on to issue the command ./nek5000 eddy_uv and following was the output. > > Number of processors: 1 > REAL wdsize : 8 > INTEGER wdsize : 4 > No file SESSION.NAME; using defaults of > PATH=. and SESSION=NEK. > > > Beginning session: > ./NEK.rea > > > timer accuracy: 0.0000000E+00 sec > > read .rea file > At line 32 of file /home/thangam/nek5_svn/trunk/nek/connect2.f (unit = 9, file = '') > Fortran runtime error: File './NEK.rea' does not exist > > > > > I have put a copy of the scripts under the /bin directory as well. > > I dont understand what I am missing and where im going wrong! > > Thanks. > Thangam. >

1 0

nek installation and running errors
by nek5000-users＠lists.mcs.anl.gov 22 Jul '14

22 Jul '14

Hi all, I am trying to download and install NEK5000 and trying to get a hang of it. I have downloaded and tried to the first example from the 'quick start' page. The compilation is successful after the makenek command. However, it does not realize further nek commands such as nek, nekb etc. (command not found error)I went on to issue the command ./nek5000 eddy_uv and following was the output. Number of processors: 1 REAL wdsize : 8 INTEGER wdsize : 4 No file SESSION.NAME; using defaults of PATH=. and SESSION=NEK. Beginning session: ./NEK.rea timer accuracy: 0.0000000E+00 sec read .rea file At line 32 of file /home/thangam/nek5_svn/trunk/nek/connect2.f (unit = 9, file = '') Fortran runtime error: File './NEK.rea' does not exist I have put a copy of the scripts under the /bin directory as well. I dont understand what I am missing and where im going wrong! Thanks. Thangam.

1 1

Extruding 2D .fld file to 3D
by nek5000-users＠lists.mcs.anl.gov 21 Jul '14

21 Jul '14

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

1 0