Hi David, The routine is doing just some sanity checks on your geometry so you shouldn't be concerned. However if you want to import a lot more from that example file it could be that certain routines/operations are hard wired to a geometry generated by genbox (which is the case here). So keep wary and if something doesn't work you need to adapt that to your case. Oana ________________________________ From: [email protected] [[email protected]] on behalf of [email protected] [[email protected]] Sent: Friday, September 18, 2015 1:49 PM To: Nek Mailing List Subject: Re: [Nek5000-users] userchk in turbChannel Oana, Okay, I will keep trying to get around it. Can you tell me what exactly this function is doing? Is it checking for mesh or input errors? David On Sep 16, 2015 3:53 PM, <[email protected]<mailto:[email protected]>> wrote: Hi David, I think that piece of code applies to a perfectly rectangular geometry where both nelx, nely are clearly defined as provided in e.g. genbox. For pipe it is indeed a tricky if not impossible matter. However the check performed there seems to me not too relevant for your problem. You could in fact just compute nelx*nely*nelz>lelx*lely*lelz Oana ________________________________________ From: [email protected]<mailto:[email protected]> [[email protected]<mailto:[email protected]>] on behalf of [email protected]<mailto:[email protected]> [[email protected]<mailto:[email protected]>] Sent: Wednesday, September 16, 2015 2:14 PM To: Nek Mailing List Subject: [Nek5000-users] userchk in turbChannel Hi all, I'm trying to use a round pipe mesh with this LES formulation but this section in userchk is tripping me up: nelx = 8 ! Number of elements in x,y, and z directions. nely = 8 ! NOTE, this may vary from one mesh to the next. nelz = 8 ! ntot = nx1*ny1*nz1*nelv ! do some checks if (istep.eq.0) then if(mod(nely,2).ne.0) then if(nid.eq.0) write(6,*) 'ABORT: nely has to be even!' call exitt endif if(nelx.gt.lelx .or. nely.gt.lely .or. nelz.gt.lelz) then if(nid.eq.0) write(6,*) 'ABORT: nel_xyz > lel_xyz!' call exitt endif endif How should I determine nelx and nely for the circular area at the end of the pipe? I have 56 elements in this area (x-y plane) extruded to 10 layers with n2to3 (560 total elements). Thanks, David _______________________________________________ Nek5000-users mailing list [email protected]<mailto:[email protected]> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users _______________________________________________ Nek5000-users mailing list [email protected]<mailto:[email protected]> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users