<p dir="ltr">Oana, <br>
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?<br>
David</p>
<div class="gmail_quote">On Sep 16, 2015 3:53 PM, <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi David,<br>
<br>
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.<br>
<br>
However the check performed there seems to me not too relevant for your problem. You could in fact just compute<br>
nelx*nely*nelz>lelx*lely*lelz<br>
<br>
Oana<br>
________________________________________<br>
From: <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov">nek5000-users-bounces@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov">nek5000-users-bounces@lists.mcs.anl.gov</a>] on behalf of <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>]<br>
Sent: Wednesday, September 16, 2015 2:14 PM<br>
To: Nek Mailing List<br>
Subject: [Nek5000-users] userchk in turbChannel<br>
<br>
Hi all,<br>
<br>
I'm trying to use a round pipe mesh with this LES formulation but this<br>
section in userchk is tripping me up:<br>
<br>
nelx = 8 ! Number of elements in x,y, and z directions.<br>
nely = 8 ! NOTE, this may vary from one mesh to the next.<br>
nelz = 8 !<br>
<br>
ntot = nx1*ny1*nz1*nelv<br>
<br>
! do some checks<br>
if (istep.eq.0) then<br>
if(mod(nely,2).ne.0) then<br>
if(nid.eq.0) write(6,*) 'ABORT: nely has to be even!'<br>
call exitt<br>
endif<br>
if(nelx.gt.lelx .or. nely.gt.lely .or. nelz.gt.lelz) then<br>
if(nid.eq.0) write(6,*) 'ABORT: nel_xyz > lel_xyz!'<br>
call exitt<br>
endif<br>
endif<br>
<br>
How should I determine nelx and nely for the circular area at the end<br>
of the pipe? I have 56 elements in this area (x-y plane) extruded to<br>
10 layers with n2to3 (560 total elements).<br>
<br>
Thanks,<br>
<br>
David<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote></div>