From nek5000-users at lists.mcs.anl.gov Tue Jan 8 09:03:47 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 8 Jan 2013 16:03:47 +0100 Subject: [Nek5000-users] prex output - curved side data Message-ID: Dear Neks, the rea files written out from prex mesh generation module contain a block entitled ******curved side data****** which contains edge and element indices. My question: which information contains the array curve (i), i=1,5 which is listed in the columns afterwards? Thanks in advance and best wishes, Joerg. From nek5000-users at lists.mcs.anl.gov Thu Jan 10 11:27:42 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 10 Jan 2013 09:27:42 -0800 (PST) Subject: [Nek5000-users] (no subject) Message-ID: Hi Neks, This question is in regards to boundary conditions. Could anyone tell me both where in the source the boundary conditions are?implemented (I have been looking through nek5_svn/trunk/nek, in many of the .f files). And also where I might be able to find some information about the Stress Boundary condition, "S". I'm looking to test a wall function on a channel and need to apply a Neumann BC in both the?stream wise?and?span wise?directions as well as a Dirichlet BC in the wall-normal direction. Thank you ? Sincerely, Josh Rabau Mechanical Engineering MS Texas A&M University 2013 Mechanical Engineering BS Rochester Institute of Technology 2011 603-381-6685 josh.rabau at yahoo.com jir2830 at .tamu.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Jan 10 12:30:51 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 10 Jan 2013 10:30:51 -0800 (PST) Subject: [Nek5000-users] Stress Boundary Condition Message-ID: Hi Neks, This question is in regards to?boundary conditions. Could anyone tell me both where in the source the boundary conditions are?implemented (I have been looking through nek5_svn/trunk/nek, in many of the .f files). And also where I might be able to find some information about the Stress Boundary condition, "S". I'm looking to test a wall function on a channel and need to apply a Neumann BC in both the?stream wise?and?span wise?directions as well as a Dirichlet BC in the wall-normal direction. Thank you ? Sincerely, Josh Rabau Mechanical Engineering MS Texas A&M University 2013 Mechanical Engineering BS Rochester Institute of Technology 2011 603-381-6685 josh.rabau at yahoo.com jir2830 at .tamu.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Jan 14 10:52:01 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 14 Jan 2013 10:52:01 -0600 (CST) Subject: [Nek5000-users] mesh problems In-Reply-To: References: Message-ID: > > I faced the same trouble with 2.17 million elements. Let me follow your > suggestion and set up a mesh in cell with z between 0 and 3 while keeping the > x-y mesh exactly the same, just to check if the same errors occur. > Hi Joerg, I had no problems with the 2.17 million case. I ran it with P=32768 and the XXt solver. When I ran genmap, I used a tolerance of .01 --- this can be important for certain classes of meshes. If the genmap tolerance is too loose then it's conceivable that you can have spurious connectivities in the graph that represents the sparse matrices, which might be the source of the problem you encountered. Paul From nek5000-users at lists.mcs.anl.gov Mon Jan 14 12:34:28 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 14 Jan 2013 19:34:28 +0100 Subject: [Nek5000-users] mesh problems In-Reply-To: References: Message-ID: Hi Paul, thanks a lot. 1/ we used the AMG preconditioner since it should scale better for bigger processor numbers. I remember of having long waiting times at the beginning with the XXt. I will check it out. 2/ I got about 80 disconnected elements while running genmap. How many did you get? I will repeat the run with tol=0.1. 3/ did you use my file y5.rea and my file zpts_compare? 4/ did you have to enforce the outer elements to radius r=0.5 in the usrdat? Janet and Joerg. Am 14.01.2013 um 17:52 schrieb nek5000-users at lists.mcs.anl.gov: >> >> I faced the same trouble with 2.17 million elements. Let me follow your >> suggestion and set up a mesh in cell with z between 0 and 3 while keeping the >> x-y mesh exactly the same, just to check if the same errors occur. >> > > Hi Joerg, > > I had no problems with the 2.17 million case. > > I ran it with P=32768 and the XXt solver. > > When I ran genmap, I used a tolerance of .01 --- this can > be important for certain classes of meshes. If the genmap tolerance > is too loose then it's conceivable that you can have spurious > connectivities in the graph that represents the sparse matrices, > which might be the source of the problem you encountered. > > Paul > > > > > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Wed Jan 23 16:32:19 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 23 Jan 2013 16:32:19 -0600 (CST) Subject: [Nek5000-users] Nek5000 maketools compile error In-Reply-To: <1564969178.1259587.1358980129436.JavaMail.root@mcs.anl.gov> Message-ID: Hi Neks, I just upgraded to a MacBook Pro running Mountain Lion on an i7 (with gcc, gfortran), and I am trying to move my Nek installation over. When I run the maketools all script, I get: gcc -mcmodel=medium -c -DUNDERSCORE ../../nek/byte.c gfortran -mcmodel=medium -fdefault-real-8 -o /Users/felker/bin/genmap genmap.o byte.o ld: in section __TEXT,__text reloc 61: length != 2 and X86_64_RELOC_GOT not supported for architecture x86_64 collect2: ld returned 1 exit status make[1]: *** [genmap] Error 1 make: *** [all] Error 1 I found that someone else had this error: http://lists.mcs.anl.gov/pipermail/nek5000-users/2012-June/001843.html but there did not appear to be any solution reported. Is there a workaround? Best, Kyle From nek5000-users at lists.mcs.anl.gov Wed Jan 23 21:15:56 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 23 Jan 2013 21:15:56 -0600 Subject: [Nek5000-users] Nek5000 maketools compile error In-Reply-To: References: <1564969178.1259587.1358980129436.JavaMail.root@mcs.anl.gov> Message-ID: Hi Kyle, I have a Mac as well, and I think what I had to do to get the tools to compile correctly was to disable big memory support. In the maketools script, simply comment out the BIGMEM="true" line and see if that works. Josh On Wed, Jan 23, 2013 at 4:32 PM, wrote: > Hi Neks, > > I just upgraded to a MacBook Pro running Mountain Lion on an i7 (with gcc, > gfortran), and I am trying to move my Nek installation over. When I run the > maketools all script, I get: > > gcc -mcmodel=medium -c -DUNDERSCORE ../../nek/byte.c > gfortran -mcmodel=medium -fdefault-real-8 -o /Users/felker/bin/genmap > genmap.o byte.o > ld: in section __TEXT,__text reloc 61: length != 2 and X86_64_RELOC_GOT > not supported for architecture x86_64 > collect2: ld returned 1 exit status > make[1]: *** [genmap] Error 1 > make: *** [all] Error 1 > > I found that someone else had this error: > > http://lists.mcs.anl.gov/pipermail/nek5000-users/2012-June/001843.html > > but there did not appear to be any solution reported. Is there a > workaround? > > Best, > > Kyle > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > -- Josh Camp "All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Jan 24 01:22:58 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 24 Jan 2013 01:22:58 -0600 Subject: [Nek5000-users] Nek5000 maketools compile error In-Reply-To: References: <1564969178.1259587.1358980129436.JavaMail.root@mcs.anl.gov> Message-ID: Hi Josh, That worked for me (after some additional workarounds with the lack of X11 on Mountain Lion). I would be curious to hear what the reason for the issue was on the new OS-- I did not have any similar problems on 10.6. Thanks! Kyle On Jan 23, 2013, at 9:15 PM, nek5000-users at lists.mcs.anl.gov wrote: > Hi Kyle, > > I have a Mac as well, and I think what I had to do to get the tools to compile correctly was to disable big memory support. In the maketools script, simply comment out the BIGMEM="true" line and see if that works. > > Josh > > > On Wed, Jan 23, 2013 at 4:32 PM, wrote: > Hi Neks, > > I just upgraded to a MacBook Pro running Mountain Lion on an i7 (with gcc, gfortran), and I am trying to move my Nek installation over. When I run the maketools all script, I get: > > gcc -mcmodel=medium -c -DUNDERSCORE ../../nek/byte.c > gfortran -mcmodel=medium -fdefault-real-8 -o /Users/felker/bin/genmap genmap.o byte.o > ld: in section __TEXT,__text reloc 61: length != 2 and X86_64_RELOC_GOT not supported for architecture x86_64 > collect2: ld returned 1 exit status > make[1]: *** [genmap] Error 1 > make: *** [all] Error 1 > > I found that someone else had this error: > > http://lists.mcs.anl.gov/pipermail/nek5000-users/2012-June/001843.html > > but there did not appear to be any solution reported. Is there a workaround? > > Best, > > Kyle > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > > > -- > Josh Camp > > "All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Jan 24 03:13:29 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 24 Jan 2013 14:43:29 +0530 Subject: [Nek5000-users] over-integration Message-ID: Hello nek5000 can use over-integration for the quadratic nonlinear terms by setting lxd, lyd, lzd. Is it possible to use this for integration to approximate some integral like applying Biot-Savart law ? I have seen the use of the mass matrices bm1 and bm2 on the velocity and pressure meshes. Is there a mass matrix on the mesh (lxd,lyd,lzd) also which can be used for integration ? I have read about the use of over-integration as a stabilization method. Can you suggest some papers/books which discuss the theory behind this ? Thanks praveen -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Jan 24 05:43:30 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 24 Jan 2013 05:43:30 -0600 Subject: [Nek5000-users] over-integration In-Reply-To: References: Message-ID: The developers can correct me if I'm wrong here, but nek doesn't actually use over-integration, although it does something similar (dealiasing). If we are computing, for example, u * u (u is vx in nek), the steps are : 1) Interpolate u to the higher space (given by lxd, lyd, lzd) (Nek has arrays for this already...for velocity, they are vxd, vyd, vzd) 2) Compute the product ud*ud on the higher order grid 3) Interpolate the solution back to the original space (lx1,ly1,lz1) I'm not aware of a mass matrix on the "dealiasing" mesh. Josh On Thu, Jan 24, 2013 at 3:13 AM, wrote: > Hello > nek5000 can use over-integration for the quadratic nonlinear terms by > setting lxd, lyd, lzd. Is it possible to use this for integration to > approximate some integral like applying Biot-Savart law ? I have seen the > use of the mass matrices bm1 and bm2 on the velocity and pressure meshes. > Is there a mass matrix on the mesh (lxd,lyd,lzd) also which can be used for > integration ? > > I have read about the use of over-integration as a stabilization method. > Can you suggest some papers/books which discuss the theory behind this ? > > Thanks > praveen > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -- Josh Camp "All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Jan 24 06:08:20 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 24 Jan 2013 13:08:20 +0100 Subject: [Nek5000-users] Hpts and IBM bluegene Message-ID: Hi neks, I had some problems with the routine hpts on a cluster IBM bluegene. The simulations hold when the routine is called. I though that it is related to same kind of input output buffer problem. Any hints? I used the gnu compilers mpicc -gcc and mpif77 -gfortran. Best Andrea -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Thu Jan 24 07:25:36 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 24 Jan 2013 07:25:36 -0600 (CST) Subject: [Nek5000-users] over-integration In-Reply-To: References: Message-ID: Hi Praveen, Josh, Overintegration and dealiasing are the same. (Aliasing corresponds to subsampling in the integrand --- overintegration properly samples.) Also, note that to get back to the coarse grid we actually use the transpose of the interpolation operator --- this corresponds to "sampling the test function on the fine grid" ... It is possible to use over integration on any function --- one needs to think about whether it's on the right-hand side or in some implicit substep. Also, we have a forthcoming paper on the topic... Paul On Thu, 24 Jan 2013, nek5000-users at lists.mcs.anl.gov wrote: > The developers can correct me if I'm wrong here, but nek doesn't actually > use over-integration, although it does something similar (dealiasing). If > we are computing, for example, u * u (u is vx in nek), the steps are : > > 1) Interpolate u to the higher space (given by lxd, lyd, lzd) (Nek has > arrays for this already...for velocity, they are vxd, vyd, vzd) > 2) Compute the product ud*ud on the higher order grid > 3) Interpolate the solution back to the original space (lx1,ly1,lz1) > > I'm not aware of a mass matrix on the "dealiasing" mesh. > > Josh > > > On Thu, Jan 24, 2013 at 3:13 AM, wrote: > >> Hello >> nek5000 can use over-integration for the quadratic nonlinear terms by >> setting lxd, lyd, lzd. Is it possible to use this for integration to >> approximate some integral like applying Biot-Savart law ? I have seen the >> use of the mass matrices bm1 and bm2 on the velocity and pressure meshes. >> Is there a mass matrix on the mesh (lxd,lyd,lzd) also which can be used for >> integration ? >> >> I have read about the use of over-integration as a stabilization method. >> Can you suggest some papers/books which discuss the theory behind this ? >> >> Thanks >> praveen >> >> _______________________________________________ >> Nek5000-users mailing list >> Nek5000-users at lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >> >> > > > -- > Josh Camp > > "All that is necessary for the triumph of evil is that good men do nothing" > -- Edmund Burke > From nek5000-users at lists.mcs.anl.gov Fri Jan 25 12:19:16 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 25 Jan 2013 12:19:16 -0600 Subject: [Nek5000-users] Hpts and IBM bluegene In-Reply-To: References: Message-ID: Hi Andrea, What version number of the code are you using? There were some problems with earlier versions. Can you try updating your source code and seeing that the file postpro.f is up to date with the current repository number? --Katie On Thu, Jan 24, 2013 at 6:08 AM, wrote: > Hi neks, > > I had some problems with the routine hpts on a cluster IBM bluegene. > > The simulations hold when the routine is called. I though that it is > related to same kind of input output buffer problem. Any hints? > > I used the gnu compilers mpicc -gcc and mpif77 -gfortran. > > Best > > Andrea > > > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Jan 29 10:48:12 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 29 Jan 2013 17:48:12 +0100 Subject: [Nek5000-users] Average in arbitrary direction Message-ID: Dear all we are trying to compute averages in arbitrary directions, as described in a previous message from May 2011: https://lists.mcs.anl.gov/mailman/htdig/nek5000-users/2011-May/001335.html It seems that we would need a file that was mentioned and linked in that message, however, it is not available any more. Is there a possibility to get that file or a similar example? Thanks in advance Best Regards Alex Gonzalo -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Jan 29 11:09:13 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 29 Jan 2013 11:09:13 -0600 Subject: [Nek5000-users] Average in arbitrary direction In-Reply-To: References: Message-ID: Hi Aleks, That file has reappeared. If you run into any problems or have any questions, please don't hesitate to ask. Hth, Justin On Tue, Jan 29, 2013 at 10:48 AM, wrote: > Dear all > > we are trying to compute averages in arbitrary directions, as described in > a previous message from May 2011: > > https://lists.mcs.anl.gov/mailman/htdig/nek5000-users/2011-May/001335.html > > It seems that we would need a file that was mentioned and linked in that > message, however, it is not available any more. > > Is there a possibility to get that file or a similar example? > > Thanks in advance > > Best Regards > Alex Gonzalo > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Jan 29 22:43:59 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 30 Jan 2013 10:13:59 +0530 Subject: [Nek5000-users] How to identify boundary faces ? Message-ID: Hello I want to identify boundary faces to set their boundary condition type. Following some of the examples, I created the following function. To find the location of the face, I need to know the coordinates (xf,yf) of the face center using which I can figure out where the face lies. Is this possible to find ? Thanks praveen C======================================================================= subroutine usrdat2 include 'SIZE' include 'TOTAL' parameter(XTOL=1e-10) character*3 cbv integer e, f real xmin, xmax, ymax xmin = -5.0 xmax = +10.0 ymax = 3.0 c set boundary conditions do e = 1,nelt do f = 1,2*ndim cbv = cbc(f,e,1) if (cbv.ne.'E ') then ! This is boundary edge c TODO: Get mid-point coordinate of face f = (xf,yf) if(abs(xf-xmin).lt.XTOL)then cbc(f,e,1) = 'v ' ! inlet else if(abs(xf-xmax).lt.XTOL)then cbc(f,e,1) = 'O ' ! outlet else if(abs(yf-ymax).lt.XTOL)then cbc(f,e,1) = 'SYM' ! top wall else cbc(f,e,1) = 'W ' ! bottom wall endif endif enddo enddo return end -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Tue Jan 29 23:27:56 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 30 Jan 2013 10:57:56 +0530 Subject: [Nek5000-users] How to identify boundary faces ? In-Reply-To: References: Message-ID: On Wed, Jan 30, 2013 at 10:13 AM, wrote: > Hello > I want to identify boundary faces to set their boundary condition type. > Following some of the examples, I created the following function. To find > the location of the face, I need to know the coordinates (xf,yf) of the > face center using which I can figure out where the face lies. Is this > possible to find ? > > Thanks > praveen > Hello By looking at this http://nek5000.mcs.anl.gov/index.php/Element_Definition I made the following code. This is a 2d problem and in blue part below I am getting face midpoint coordinates. I am not sure I have got the numbering of the faces correct. Thanks praveen C======================================================================= subroutine usrdat2 include 'SIZE' include 'TOTAL' parameter(XTOL=1e-10) character*3 cbv integer e, f real xmin, xmax, ymax, xf(4), yf(4) xmin = -5.0 xmax = +10.0 ymax = 3.0 c set boundary conditions do e = 1,nelv xf(3) = 0.5*(xm1(1, 1, 1,e) + xm1(nx1,1, 1,e)) xf(2) = 0.5*(xm1(nx1,1, 1,e) + xm1(nx1,ny1,1,e)) xf(4) = 0.5*(xm1(1, ny1,1,e) + xm1(nx1,ny1,1,e)) xf(1) = 0.5*(xm1(1, ny1,1,e) + xm1(1, 1, 1,e)) yf(3) = 0.5*(ym1(1, 1, 1,e) + ym1(nx1,1, 1,e)) yf(2) = 0.5*(ym1(nx1,1, 1,e) + ym1(nx1,ny1,1,e)) yf(4) = 0.5*(ym1(1, ny1,1,e) + ym1(nx1,ny1,1,e)) yf(1) = 0.5*(ym1(1, ny1,1,e) + ym1(1, 1, 1,e)) do f = 1,2*ndim cbv = cbc(f,e,1) if (cbv.ne.'E ') then ! This is boundary edge if(abs(xf(f)-xmin).lt.XTOL)then cbc(f,e,1) = 'v ' ! inlet else if(abs(xf(f)-xmax).lt.XTOL)then cbc(f,e,1) = 'O ' ! outlet else if(abs(yf(f)-ymax).lt.XTOL)then cbc(f,e,1) = 'SYM' ! top wall else cbc(f,e,1) = 'W ' ! bottom wall endif endif enddo enddo return end -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Wed Jan 30 15:40:17 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 30 Jan 2013 14:40:17 -0700 Subject: [Nek5000-users] How to identify boundary faces ? In-Reply-To: References: Message-ID: Dear Nek heads, I am having trouble getting Nek5000 to read my restart file. The file is an ASCII file that was written from a previous run. I set parameter p67 to -1 to read in an ascii restart file. Here is the error message in the .o file: set initial conditions Checking restart options: sqduct.fld06 nekuic (1) for ifld 1 Reading checkpoint data WARNING: No data read in for file sqduct.fld06 Failed on element 2, point 0. Last read dump: 0.000000 0.2641195E-02 0.000000 0.000000 0.000000 Last read dump: 0.000000 0.000000 -0.6859594E-02 0.2343750E-01 -0.2600098E-01 Last read dump: 0.2600098E-01 -0.1855469E-01 0 call exitt 1702a 1 Here is the restart call in the .rea file: ***** NO THERMAL BOUNDARY CONDITIONS ***** 1 RESTART OPTIONS FOLLOW sqduct.fld06 1 INITIAL CONDITIONS ***** C Default Here is the first few lines of the restart file sqduct.fld06: **** 8 8 8 8.0000000E+01 1500 X Y Z U P NELT,NX,NY,N 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 Note that the problem is a square duct with 24 x 24 x 32 elements, which equal 18432 total elements. Since this is a 5 digit number, I suppose that is why it is not shown in the first line in the field file. However, I tried changing the four *'s to 18432, along with a few other variations, and no luck, same error. Also, the numbers in the dump after the error message above don't seem to appear in the restart file. Also, I don't have any trouble withe a restart file that has only 6 x 6 x 8 = 288 elements in it. Finally, I am using a version of Nek500 that I downloaded about 11 months ago. Help is appreciated. Rich Johnson On Tue, Jan 29, 2013 at 10:27 PM, wrote: > > On Wed, Jan 30, 2013 at 10:13 AM, wrote: > >> Hello >> I want to identify boundary faces to set their boundary condition type. >> Following some of the examples, I created the following function. To find >> the location of the face, I need to know the coordinates (xf,yf) of the >> face center using which I can figure out where the face lies. Is this >> possible to find ? >> >> Thanks >> praveen >> > > Hello > > By looking at this > > http://nek5000.mcs.anl.gov/index.php/Element_Definition > > I made the following code. This is a 2d problem and in blue part below I > am getting face midpoint coordinates. I am not sure I have got the > numbering of the faces correct. > > Thanks > praveen > > C======================================================================= > subroutine usrdat2 > > include 'SIZE' > include 'TOTAL' > > parameter(XTOL=1e-10) > character*3 cbv > integer e, f > real xmin, xmax, ymax, xf(4), yf(4) > > xmin = -5.0 > xmax = +10.0 > ymax = 3.0 > > c set boundary conditions > do e = 1,nelv > xf(3) = 0.5*(xm1(1, 1, 1,e) + xm1(nx1,1, 1,e)) > xf(2) = 0.5*(xm1(nx1,1, 1,e) + xm1(nx1,ny1,1,e)) > xf(4) = 0.5*(xm1(1, ny1,1,e) + xm1(nx1,ny1,1,e)) > xf(1) = 0.5*(xm1(1, ny1,1,e) + xm1(1, 1, 1,e)) > yf(3) = 0.5*(ym1(1, 1, 1,e) + ym1(nx1,1, 1,e)) > yf(2) = 0.5*(ym1(nx1,1, 1,e) + ym1(nx1,ny1,1,e)) > yf(4) = 0.5*(ym1(1, ny1,1,e) + ym1(nx1,ny1,1,e)) > yf(1) = 0.5*(ym1(1, ny1,1,e) + ym1(1, 1, 1,e)) > do f = 1,2*ndim > cbv = cbc(f,e,1) > if (cbv.ne.'E ') then ! This is boundary edge > if(abs(xf(f)-xmin).lt.XTOL)then > cbc(f,e,1) = 'v ' ! inlet > else if(abs(xf(f)-xmax).lt.XTOL)then > cbc(f,e,1) = 'O ' ! outlet > else if(abs(yf(f)-ymax).lt.XTOL)then > cbc(f,e,1) = 'SYM' ! top wall > else > cbc(f,e,1) = 'W ' ! bottom wall > endif > endif > enddo > enddo > > return > end > > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Wed Jan 30 17:03:24 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 30 Jan 2013 17:03:24 -0600 (CST) Subject: [Nek5000-users] How to identify boundary faces ? In-Reply-To: References: Message-ID: Hi Praveen, Below is one approach... not tested, but it should work. In general, I would use facint_v() from navier5.f to get the average of (say) xm1 on a particular face. I didn't do that here because we also needed to verify that the face was internal. Note that 'E ' is not a guaranteed way to check because there are some cases where CBC is not defined as 'E ', even when it is in fact an E-E boundary condition. The reason for this is that the E-E information is no longer used in the nek5000 simulation code. This information comes instead from the .map file. Paul c----------------------------------------------------------------------- subroutine usrdat2 include 'SIZE' include 'TOTAL' character*3 cbv integer e,f,fmid(6) real xmin, xmax, ymax common /scrns/ one(lx1,ly1,lz1,lelt) n = nx1*ny1*nz1*nelt call rone (one,n) ifield = 1 call dssum(one,nx1,ny1,nz1) ! Identify shared faces im = (1+nx1)/2 ! Face midpoints jm = (1+ny1)/2 km = (1+nz1)/2 fmid(4) = 1 + nx1*( jm-1) + nx1*ny1*( km-1) ! x = -1 fmid(2) = nx1 + nx1*( jm-1) + nx1*ny1*( km-1) ! x = 1 fmid(1) = im + nx1*( 1-1) + nx1*ny1*( km-1) ! y = -1 fmid(3) = im + nx1*(ny1-1) + nx1*ny1*( km-1) ! y = 1 fmid(5) = im + nx1*( jm-1) + nx1*ny1*( 1-1) ! z = -1 fmid(6) = im + nx1*( jm-1) + nx1*ny1*(nz1-1) ! z = 1 xmin = -5.0 xmax = +10.0 ymax = 3.0 tol = 1.e-5 do e = 1,nelt ! set boundary conditions do f = 1,2*ndim count = one(fmid(f),1,1,e) xf = xm1(fmid(f),1,1,e) yf = ym1(fmid(f),1,1,e) zf = zm1(fmid(f),1,1,e) if (count.lt.1.1) then ! This is boundary edge if(abs(xf-xmin).lt.tol)then cbc(f,e,1) = 'v ' ! inlet else if(abs(xf-xmax).lt.tol)then cbc(f,e,1) = 'O ' ! outlet else if(abs(yf-ymax).lt.tol)then cbc(f,e,1) = 'SYM' ! top wall else cbc(f,e,1) = 'W ' ! bottom wall endif endif enddo enddo return end c----------------------------------------------------------------------- On Wed, 30 Jan 2013, nek5000-users at lists.mcs.anl.gov wrote: > Hello > I want to identify boundary faces to set their boundary condition type. > Following some of the examples, I created the following function. To find > the location of the face, I need to know the coordinates (xf,yf) of the > face center using which I can figure out where the face lies. Is this > possible to find ? > > Thanks > praveen > > C======================================================================= > subroutine usrdat2 > > include 'SIZE' > include 'TOTAL' > > parameter(XTOL=1e-10) > character*3 cbv > integer e, f > real xmin, xmax, ymax > > xmin = -5.0 > xmax = +10.0 > ymax = 3.0 > > c set boundary conditions > do e = 1,nelt > do f = 1,2*ndim > cbv = cbc(f,e,1) > if (cbv.ne.'E ') then ! This is boundary edge > c TODO: Get mid-point coordinate of face f = (xf,yf) > if(abs(xf-xmin).lt.XTOL)then > cbc(f,e,1) = 'v ' ! inlet > else if(abs(xf-xmax).lt.XTOL)then > cbc(f,e,1) = 'O ' ! outlet > else if(abs(yf-ymax).lt.XTOL)then > cbc(f,e,1) = 'SYM' ! top wall > else > cbc(f,e,1) = 'W ' ! bottom wall > endif > endif > enddo > enddo > > return > end > From nek5000-users at lists.mcs.anl.gov Wed Jan 30 17:04:49 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 30 Jan 2013 17:04:49 -0600 (CST) Subject: [Nek5000-users] How to identify boundary faces ? In-Reply-To: References: Message-ID: Hi Rich, Try, in usrdat2(): param(67)=0 ! Set read format (66 is for write format) Let me know if this works. Cheers, Paul On Wed, 30 Jan 2013, nek5000-users at lists.mcs.anl.gov wrote: > Dear Nek heads, > I am having trouble getting Nek5000 to read my restart file. The file is an > ASCII file that was written from a previous run. I set parameter p67 to -1 > to read in an ascii restart file. Here is the error message in the .o file: > > set initial conditions > Checking restart options: > sqduct.fld06 > > nekuic (1) for ifld 1 > Reading checkpoint data > WARNING: No data read in for file > sqduct.fld06 > > Failed on element 2, point 0. > Last read dump: > 0.000000 0.2641195E-02 0.000000 0.000000 0.000000 > Last read dump: > 0.000000 0.000000 -0.6859594E-02 0.2343750E-01 -0.2600098E-01 > Last read dump: > 0.2600098E-01 -0.1855469E-01 > 0 call exitt 1702a 1 > > > Here is the restart call in the .rea file: > > ***** NO THERMAL BOUNDARY CONDITIONS ***** > 1 RESTART OPTIONS FOLLOW > sqduct.fld06 > 1 INITIAL CONDITIONS ***** > C Default > > > Here is the first few lines of the restart file sqduct.fld06: > > > **** 8 8 8 8.0000000E+01 1500 X Y Z U P > NELT,NX,NY,N > 0.000 0.000 0.000 0.000 0.000 0.000 > 0.000 0.000 0.000 0.000 0.000 0.000 > 0.000 0.000 0.000 0.000 0.000 0.000 > 0.000 0.000 0.000 0.000 0.000 0.000 > 0.000 0.000 0.000 0.000 0.000 0.000 > 0.000 0.000 0.000 0.000 0.000 0.000 > > > Note that the problem is a square duct with 24 x 24 x 32 elements, which > equal 18432 total elements. Since this is a 5 digit number, I suppose that > is why it is not shown in the first line in the field file. However, I > tried changing the four *'s to 18432, along with a few other variations, > and no luck, same error. Also, the numbers in the dump after the error > message above don't seem to appear in the restart file. Also, I don't have > any trouble withe a restart file that has only 6 x 6 x 8 = 288 elements in > it. Finally, I am using a version of Nek500 that I downloaded about 11 > months ago. > > Help is appreciated. > > Rich Johnson > > > > > > On Tue, Jan 29, 2013 at 10:27 PM, wrote: > >> >> On Wed, Jan 30, 2013 at 10:13 AM, wrote: >> >>> Hello >>> I want to identify boundary faces to set their boundary condition type. >>> Following some of the examples, I created the following function. To find >>> the location of the face, I need to know the coordinates (xf,yf) of the >>> face center using which I can figure out where the face lies. Is this >>> possible to find ? >>> >>> Thanks >>> praveen >>> >> >> Hello >> >> By looking at this >> >> http://nek5000.mcs.anl.gov/index.php/Element_Definition >> >> I made the following code. This is a 2d problem and in blue part below I >> am getting face midpoint coordinates. I am not sure I have got the >> numbering of the faces correct. >> >> Thanks >> praveen >> >> C======================================================================= >> subroutine usrdat2 >> >> include 'SIZE' >> include 'TOTAL' >> >> parameter(XTOL=1e-10) >> character*3 cbv >> integer e, f >> real xmin, xmax, ymax, xf(4), yf(4) >> >> xmin = -5.0 >> xmax = +10.0 >> ymax = 3.0 >> >> c set boundary conditions >> do e = 1,nelv >> xf(3) = 0.5*(xm1(1, 1, 1,e) + xm1(nx1,1, 1,e)) >> xf(2) = 0.5*(xm1(nx1,1, 1,e) + xm1(nx1,ny1,1,e)) >> xf(4) = 0.5*(xm1(1, ny1,1,e) + xm1(nx1,ny1,1,e)) >> xf(1) = 0.5*(xm1(1, ny1,1,e) + xm1(1, 1, 1,e)) >> yf(3) = 0.5*(ym1(1, 1, 1,e) + ym1(nx1,1, 1,e)) >> yf(2) = 0.5*(ym1(nx1,1, 1,e) + ym1(nx1,ny1,1,e)) >> yf(4) = 0.5*(ym1(1, ny1,1,e) + ym1(nx1,ny1,1,e)) >> yf(1) = 0.5*(ym1(1, ny1,1,e) + ym1(1, 1, 1,e)) >> do f = 1,2*ndim >> cbv = cbc(f,e,1) >> if (cbv.ne.'E ') then ! This is boundary edge >> if(abs(xf(f)-xmin).lt.XTOL)then >> cbc(f,e,1) = 'v ' ! inlet >> else if(abs(xf(f)-xmax).lt.XTOL)then >> cbc(f,e,1) = 'O ' ! outlet >> else if(abs(yf(f)-ymax).lt.XTOL)then >> cbc(f,e,1) = 'SYM' ! top wall >> else >> cbc(f,e,1) = 'W ' ! bottom wall >> endif >> endif >> enddo >> enddo >> >> return >> end >> >> >> _______________________________________________ >> Nek5000-users mailing list >> Nek5000-users at lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >> >> > From nek5000-users at lists.mcs.anl.gov Wed Jan 30 17:23:10 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 30 Jan 2013 17:23:10 -0600 Subject: [Nek5000-users] How to identify boundary faces ? In-Reply-To: References: Message-ID: Hi Rich, I have updated the subroutine that dumps out ASCII files to expend the format for larger element counts. Keeping p66 = -1 ! write p67 = -1 ! read in your .rea file is fine. Your best option is to re-run the case that produced this restart file so that the header is now correct. If you don't want to update the source and re-run the case that produced this restart file, it "should" work if you change the first line of the restart file from: **** 8 8 8 8.0000000E+01 1500 X Y Z U P NELT,NX,NY,N to 18432 8 8 8 8.800000000E+01 1500 X Y Z U P NELT,NX,NY,N spacing is important. Best, Katie On Wed, Jan 30, 2013 at 5:04 PM, wrote: > > > Hi Rich, > > Try, in usrdat2(): > > param(67)=0 ! Set read format > > > (66 is for write format) > > Let me know if this works. > > Cheers, > > Paul > > > > > On Wed, 30 Jan 2013, nek5000-users at lists.mcs.anl.**govwrote: > > Dear Nek heads, >> I am having trouble getting Nek5000 to read my restart file. The file is >> an >> ASCII file that was written from a previous run. I set parameter p67 to -1 >> to read in an ascii restart file. Here is the error message in the .o >> file: >> >> set initial conditions >> Checking restart options: >> sqduct.fld06 >> >> nekuic (1) for ifld 1 >> Reading checkpoint data >> WARNING: No data read in for file >> sqduct.fld06 >> >> Failed on element 2, point 0. >> Last read dump: >> 0.000000 0.2641195E-02 0.000000 0.000000 0.000000 >> Last read dump: >> 0.000000 0.000000 -0.6859594E-02 0.2343750E-01 -0.2600098E-01 >> Last read dump: >> 0.2600098E-01 -0.1855469E-01 >> 0 call exitt 1702a 1 >> >> >> Here is the restart call in the .rea file: >> >> ***** NO THERMAL BOUNDARY CONDITIONS ***** >> 1 RESTART OPTIONS FOLLOW >> sqduct.fld06 >> 1 INITIAL CONDITIONS ***** >> C Default >> >> >> Here is the first few lines of the restart file sqduct.fld06: >> >> >> **** 8 8 8 8.0000000E+01 1500 X Y Z U P >> NELT,NX,NY,N >> 0.000 0.000 0.000 0.000 0.000 0.000 >> 0.000 0.000 0.000 0.000 0.000 0.000 >> 0.000 0.000 0.000 0.000 0.000 0.000 >> 0.000 0.000 0.000 0.000 0.000 0.000 >> 0.000 0.000 0.000 0.000 0.000 0.000 >> 0.000 0.000 0.000 0.000 0.000 0.000 >> >> >> Note that the problem is a square duct with 24 x 24 x 32 elements, which >> equal 18432 total elements. Since this is a 5 digit number, I suppose that >> is why it is not shown in the first line in the field file. However, I >> tried changing the four *'s to 18432, along with a few other variations, >> and no luck, same error. Also, the numbers in the dump after the error >> message above don't seem to appear in the restart file. Also, I don't have >> any trouble withe a restart file that has only 6 x 6 x 8 = 288 elements in >> it. Finally, I am using a version of Nek500 that I downloaded about 11 >> months ago. >> >> Help is appreciated. >> >> Rich Johnson >> >> >> >> >> >> On Tue, Jan 29, 2013 at 10:27 PM, > >> wrote: >> >> >>> On Wed, Jan 30, 2013 at 10:13 AM, > >>> wrote: >>> >>> Hello >>>> I want to identify boundary faces to set their boundary condition type. >>>> Following some of the examples, I created the following function. To >>>> find >>>> the location of the face, I need to know the coordinates (xf,yf) of the >>>> face center using which I can figure out where the face lies. Is this >>>> possible to find ? >>>> >>>> Thanks >>>> praveen >>>> >>>> >>> Hello >>> >>> By looking at this >>> >>> http://nek5000.mcs.anl.gov/**index.php/Element_Definition >>> >>> I made the following code. This is a 2d problem and in blue part below I >>> am getting face midpoint coordinates. I am not sure I have got the >>> numbering of the faces correct. >>> >>> Thanks >>> praveen >>> >>> C=============================**==============================** >>> ============ >>> subroutine usrdat2 >>> >>> include 'SIZE' >>> include 'TOTAL' >>> >>> parameter(XTOL=1e-10) >>> character*3 cbv >>> integer e, f >>> real xmin, xmax, ymax, xf(4), yf(4) >>> >>> xmin = -5.0 >>> xmax = +10.0 >>> ymax = 3.0 >>> >>> c set boundary conditions >>> do e = 1,nelv >>> xf(3) = 0.5*(xm1(1, 1, 1,e) + xm1(nx1,1, 1,e)) >>> xf(2) = 0.5*(xm1(nx1,1, 1,e) + xm1(nx1,ny1,1,e)) >>> xf(4) = 0.5*(xm1(1, ny1,1,e) + xm1(nx1,ny1,1,e)) >>> xf(1) = 0.5*(xm1(1, ny1,1,e) + xm1(1, 1, 1,e)) >>> yf(3) = 0.5*(ym1(1, 1, 1,e) + ym1(nx1,1, 1,e)) >>> yf(2) = 0.5*(ym1(nx1,1, 1,e) + ym1(nx1,ny1,1,e)) >>> yf(4) = 0.5*(ym1(1, ny1,1,e) + ym1(nx1,ny1,1,e)) >>> yf(1) = 0.5*(ym1(1, ny1,1,e) + ym1(1, 1, 1,e)) >>> do f = 1,2*ndim >>> cbv = cbc(f,e,1) >>> if (cbv.ne.'E ') then ! This is boundary edge >>> if(abs(xf(f)-xmin).lt.XTOL)**then >>> cbc(f,e,1) = 'v ' ! inlet >>> else if(abs(xf(f)-xmax).lt.XTOL)**then >>> cbc(f,e,1) = 'O ' ! outlet >>> else if(abs(yf(f)-ymax).lt.XTOL)**then >>> cbc(f,e,1) = 'SYM' ! top wall >>> else >>> cbc(f,e,1) = 'W ' ! bottom wall >>> endif >>> endif >>> enddo >>> enddo >>> >>> return >>> end >>> >>> >>> ______________________________**_________________ >>> 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 > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Wed Jan 30 18:04:01 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 30 Jan 2013 17:04:01 -0700 Subject: [Nek5000-users] How to identify boundary faces ? In-Reply-To: References: Message-ID: Katie Thank you. The alternate heading worked. Which is better than re-running the 2.5 week long DNS calculation that created it. Rich Johnson On Wed, Jan 30, 2013 at 4:23 PM, wrote: > Hi Rich, > > I have updated the subroutine that dumps out ASCII files to expend the > format for larger element counts. > > Keeping > p66 = -1 ! write > p67 = -1 ! read > > in your .rea file is fine. > > Your best option is to re-run the case that produced this restart file so > that the header is now correct. > > If you don't want to update the source and re-run the case that produced > this restart file, it "should" work if you change the first line of the > restart file from: > **** 8 8 8 8.0000000E+01 1500 X Y Z U P > NELT,NX,NY,N > > to > 18432 8 8 8 8.800000000E+01 1500 X Y Z U > P NELT,NX,NY,N > > spacing is important. > > Best, > Katie > > > > On Wed, Jan 30, 2013 at 5:04 PM, wrote: > >> >> >> Hi Rich, >> >> Try, in usrdat2(): >> >> param(67)=0 ! Set read format >> >> >> (66 is for write format) >> >> Let me know if this works. >> >> Cheers, >> >> Paul >> >> >> >> >> On Wed, 30 Jan 2013, nek5000-users at lists.mcs.anl.**govwrote: >> >> Dear Nek heads, >>> I am having trouble getting Nek5000 to read my restart file. The file is >>> an >>> ASCII file that was written from a previous run. I set parameter p67 to >>> -1 >>> to read in an ascii restart file. Here is the error message in the .o >>> file: >>> >>> set initial conditions >>> Checking restart options: >>> sqduct.fld06 >>> >>> nekuic (1) for ifld 1 >>> Reading checkpoint data >>> WARNING: No data read in for file >>> sqduct.fld06 >>> >>> Failed on element 2, point 0. >>> Last read dump: >>> 0.000000 0.2641195E-02 0.000000 0.000000 0.000000 >>> Last read dump: >>> 0.000000 0.000000 -0.6859594E-02 0.2343750E-01 >>> -0.2600098E-01 >>> Last read dump: >>> 0.2600098E-01 -0.1855469E-01 >>> 0 call exitt 1702a 1 >>> >>> >>> Here is the restart call in the .rea file: >>> >>> ***** NO THERMAL BOUNDARY CONDITIONS ***** >>> 1 RESTART OPTIONS FOLLOW >>> sqduct.fld06 >>> 1 INITIAL CONDITIONS ***** >>> C Default >>> >>> >>> Here is the first few lines of the restart file sqduct.fld06: >>> >>> >>> **** 8 8 8 8.0000000E+01 1500 X Y Z U P >>> NELT,NX,NY,N >>> 0.000 0.000 0.000 0.000 0.000 0.000 >>> 0.000 0.000 0.000 0.000 0.000 0.000 >>> 0.000 0.000 0.000 0.000 0.000 0.000 >>> 0.000 0.000 0.000 0.000 0.000 0.000 >>> 0.000 0.000 0.000 0.000 0.000 0.000 >>> 0.000 0.000 0.000 0.000 0.000 0.000 >>> >>> >>> Note that the problem is a square duct with 24 x 24 x 32 elements, which >>> equal 18432 total elements. Since this is a 5 digit number, I suppose >>> that >>> is why it is not shown in the first line in the field file. However, I >>> tried changing the four *'s to 18432, along with a few other variations, >>> and no luck, same error. Also, the numbers in the dump after the error >>> message above don't seem to appear in the restart file. Also, I don't >>> have >>> any trouble withe a restart file that has only 6 x 6 x 8 = 288 elements >>> in >>> it. Finally, I am using a version of Nek500 that I downloaded about 11 >>> months ago. >>> >>> Help is appreciated. >>> >>> Rich Johnson >>> >>> >>> >>> >>> >>> On Tue, Jan 29, 2013 at 10:27 PM, > >>> wrote: >>> >>> >>>> On Wed, Jan 30, 2013 at 10:13 AM, > >>>> wrote: >>>> >>>> Hello >>>>> I want to identify boundary faces to set their boundary condition type. >>>>> Following some of the examples, I created the following function. To >>>>> find >>>>> the location of the face, I need to know the coordinates (xf,yf) of the >>>>> face center using which I can figure out where the face lies. Is this >>>>> possible to find ? >>>>> >>>>> Thanks >>>>> praveen >>>>> >>>>> >>>> Hello >>>> >>>> By looking at this >>>> >>>> http://nek5000.mcs.anl.gov/**index.php/Element_Definition >>>> >>>> I made the following code. This is a 2d problem and in blue part below I >>>> am getting face midpoint coordinates. I am not sure I have got the >>>> numbering of the faces correct. >>>> >>>> Thanks >>>> praveen >>>> >>>> C=============================**==============================** >>>> ============ >>>> subroutine usrdat2 >>>> >>>> include 'SIZE' >>>> include 'TOTAL' >>>> >>>> parameter(XTOL=1e-10) >>>> character*3 cbv >>>> integer e, f >>>> real xmin, xmax, ymax, xf(4), yf(4) >>>> >>>> xmin = -5.0 >>>> xmax = +10.0 >>>> ymax = 3.0 >>>> >>>> c set boundary conditions >>>> do e = 1,nelv >>>> xf(3) = 0.5*(xm1(1, 1, 1,e) + xm1(nx1,1, 1,e)) >>>> xf(2) = 0.5*(xm1(nx1,1, 1,e) + xm1(nx1,ny1,1,e)) >>>> xf(4) = 0.5*(xm1(1, ny1,1,e) + xm1(nx1,ny1,1,e)) >>>> xf(1) = 0.5*(xm1(1, ny1,1,e) + xm1(1, 1, 1,e)) >>>> yf(3) = 0.5*(ym1(1, 1, 1,e) + ym1(nx1,1, 1,e)) >>>> yf(2) = 0.5*(ym1(nx1,1, 1,e) + ym1(nx1,ny1,1,e)) >>>> yf(4) = 0.5*(ym1(1, ny1,1,e) + ym1(nx1,ny1,1,e)) >>>> yf(1) = 0.5*(ym1(1, ny1,1,e) + ym1(1, 1, 1,e)) >>>> do f = 1,2*ndim >>>> cbv = cbc(f,e,1) >>>> if (cbv.ne.'E ') then ! This is boundary edge >>>> if(abs(xf(f)-xmin).lt.XTOL)**then >>>> cbc(f,e,1) = 'v ' ! inlet >>>> else if(abs(xf(f)-xmax).lt.XTOL)**then >>>> cbc(f,e,1) = 'O ' ! outlet >>>> else if(abs(yf(f)-ymax).lt.XTOL)**then >>>> cbc(f,e,1) = 'SYM' ! top wall >>>> else >>>> cbc(f,e,1) = 'W ' ! bottom wall >>>> endif >>>> endif >>>> enddo >>>> enddo >>>> >>>> return >>>> end >>>> >>>> >>>> ______________________________**_________________ >>>> 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 >> > > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: