<div dir="ltr">Hello,<div><br></div><div style>As a first step, here's a code snippet that would give you the xyz coordinates of the inlet (assuming there is only one inlet):</div><div style><br></div><div style>do iel=1,nelv</div>
<div style> do ifc=1,2*ndim</div><div style> if (cbc(ifc,iel,1).eq.'v ') then</div><div style> call facind(kx1,kx2,ky1,ky2,kz1,kz2,nx1,ny1,nz1,ifc)</div><div style> do iz=kz1,kz2</div><div style> do iy=ky1,ky2</div>
<div style> do ix=kx1,kx2</div><div style> x = xm1(ix,iy,iz,iel)</div><div style> y = ym1(ix,iy,iz,iel)</div><div style> z = zm1(ix,iy,iz,iel)</div><div style> ! Do whatever printing you want to do here</div>
<div style> enddo</div><div style> enddo</div><div style> enddo</div><div style> endif</div><div style> enddo</div><div style>enddo</div><div style><br></div><div style>I don't think nek has any built-in routines for interpolating face data, but one could certainly build one using some of the functions from speclib.f.</div>
<div style><br></div><div style>Josh</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Feb 27, 2013 at 5:52 AM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-SG" link="blue" vlink="purple"><div><p class="MsoNormal"><span lang="EN-GB" style="font-size:12.0pt;font-family:"Times New Roman","serif"">Hello nek5000-users,<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:12.0pt;font-family:"Times New Roman","serif""><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-GB" style="font-size:12.0pt;font-family:"Times New Roman","serif"">I have a data file obtained through a RANS simulation that I would like to use as the inflow condition for my subsequent nek5000 simulation. <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:12.0pt;font-family:"Times New Roman","serif""><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-GB" style="font-size:12.0pt;font-family:"Times New Roman","serif"">May I know whether there is an efficient subroutine which allows me to interpolate the RANS data onto the spectral-element mesh or, at least, allows me to dump out the nodal points at the inlet so that I can perform the interpolation and assignment? <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:12.0pt;font-family:"Times New Roman","serif""><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-GB" style="font-size:12.0pt;font-family:"Times New Roman","serif"">I will be performing a parallel simulation involving 300000 elements. Hence I will be grateful if someone could provide me with an efficient method.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:12.0pt;font-family:"Times New Roman","serif""><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-GB" style="font-size:12.0pt;font-family:"Times New Roman","serif"">Thanks in advance.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman","serif""><u></u> <u></u></span></p><p class="MsoNormal" style="margin-bottom:6.0pt"><b><span style="font-size:18.0pt;font-family:"Times New Roman","serif"">Ng Jee Hann (Mr.)</span></b><span style="font-family:"Times New Roman","serif""><u></u><u></u></span></p>
<p class="MsoNormal" style="margin-right:0cm;margin-bottom:6.0pt;margin-left:17.85pt;line-height:14.0pt"><i><span style="font-size:12.0pt;font-family:"Times New Roman","serif"">PhD Student.<u></u><u></u></span></i></p>
<p class="MsoNormal" style="margin-left:18.0pt;line-height:14.0pt"><span style="font-size:12.0pt;font-family:"Times New Roman","serif"">Department of Mechanical Engineering,<u></u><u></u></span></p><p class="MsoNormal" style="margin-left:18.0pt;line-height:14.0pt">
<span style="font-size:12.0pt;font-family:"Times New Roman","serif"">National University of Singapore,<u></u><u></u></span></p><p class="MsoNormal" style="margin-left:18.0pt;line-height:14.0pt"><span style="font-size:12.0pt;font-family:"Times New Roman","serif"">9 Engineering Drive 1, <u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:18.0pt;line-height:14.0pt"><span style="font-size:12.0pt;font-family:"Times New Roman","serif"">Singapore 117576.<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-GB" style="font-size:12.0pt;font-family:"Times New Roman","serif""><u></u> <u></u></span></p>
</div></div><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" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>Josh Camp<br><br>"All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke
</div>