<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;" dir="ltr">
<p>Hi Stefan,</p>
<p><br>
</p>
<p>Many thanks for your reply. I'm running my simulations on ARCHER <span>(Cray XC30 supercomputer</span>) and I'm using the Intel compiler.</p>
<p><br>
</p>
<p>The current simulation has 469500 elements with polynomial order of 7. The source field I used for interpolation has 383700 elements with the same polynomial order. The size of the source field is about 10GB. I'm not sure if it can be easily transferred.
 But I have attached my .usr and log files herewith (also the gfldr.f with few added checking lines). I ran the interpolation using 720 processors.</p>
<p><br>
</p>
<p>Please have a look and let me know if you need further information about the case. Thank you very much for you help.</p>
<p><br>
</p>
<p>Best regards,</p>
<p>Tony<br>
</p>
<br>
<br>
<div style="color: rgb(19, 19, 18);">
<div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><b>From:</b> nek5000-users-bounces@lists.mcs.anl.gov <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users-request@lists.mcs.anl.gov <nek5000-users-request@lists.mcs.anl.gov><br>
<b>Sent:</b> 28 March 2017 13:19<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Nek5000-users Digest, Vol 97, Issue 31</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Send Nek5000-users mailing list submissions to<br>
        nek5000-users@lists.mcs.anl.gov<br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a previewremoved="true" id="LPlnk826007" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
or, via email, send a message with subject or body 'help' to<br>
        nek5000-users-request@lists.mcs.anl.gov<br>
<br>
You can reach the person managing the list at<br>
        nek5000-users-owner@lists.mcs.anl.gov<br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Nek5000-users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
   1. Re: Interpolation between two meshes<br>
      (nek5000-users@lists.mcs.anl.gov)<br>
   2. Re: Copy velocity slice from one boundary to another<br>
      (nek5000-users@lists.mcs.anl.gov)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Tue, 28 Mar 2017 09:57:33 +0200<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: nek5000-users@lists.mcs.anl.gov      <nek5000-users@lists.mcs.anl.gov><br>
Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
Message-ID:<br>
        <mailman.3869.1490687873.2967.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset=utf-8<br>
<br>
Can you please provide more information about your case (machine, compiler) and include the .usr and log file.<br>
If the case is reasonably small, I can try to run it myself. <br>
 <br>
-----Original message-----<br>
> From:nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
> Sent: Tuesday 28th March 2017 4:03<br>
> To: nek5000-users@lists.mcs.anl.gov<br>
> Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
> <br>
> Hi Stefan,<br>
> <br>
> Many thanks for your reply. I have tried the new 'gfldr' subroutine in two ways.<br>
> <br>
> First, I copied the whole 'gfldr.f' into my usr file and ran the interpolation with my legacy nek master code<br>
> (tuned few lines to match the old byte_mpi.f format). The source field was read in successfully and I also confirmed the coordinates from the source field. However, the interpolation halted again at 'findpts' in 'gfldr_intp' subroutine.<br>
> <br>
> Second, I ran the interpolation directly with the latest master nek code. The same thing happened and the interpolation halted at 'findpts' and never moved on.<br>
> <br>
> In addition, I also tried reducing polynomial order (lx1) as Mattias suggested, but with no success.<br>
> <br>
> I really don't know what causes the problem when interpolating between two large meshes. Do you have any workaround?<br>
> <br>
> Please help me out. Thank you very much.<br>
> <br>
> Best regards,<br>
> Tony<br>
> <br>
> <br>
> -----------<br>
> From: nek5000-users-bounces@lists.mcs.anl.gov <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users-request@lists.mcs.anl.gov <nek5000-users-request@lists.mcs.anl.gov><br>
> Sent: 27 March 2017 18:00<br>
> To: nek5000-users@lists.mcs.anl.gov<br>
> Subject: Nek5000-users Digest, Vol 97, Issue 29<br>
> Send Nek5000-users mailing list submissions to<br>
> ??????? nek5000-users@lists.mcs.anl.gov<br>
> <br>
> To subscribe or unsubscribe via the World Wide Web, visit<br>
> <a previewremoved="true" id="LPlnk881824" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a> <<a previewremoved="true" id="LPlnk97559" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>><br>
> or, via email, send a message with subject or body 'help' to<br>
> ??????? nek5000-users-request@lists.mcs.anl.gov<br>
> <br>
> You can reach the person managing the list at<br>
> ??????? nek5000-users-owner@lists.mcs.anl.gov<br>
> <br>
> When replying, please edit your Subject line so it is more specific<br>
> than "Re: Contents of Nek5000-users digest..."<br>
> <br>
> <br>
> Today's Topics:<br>
> <br>
> ?? 1. Re: Interpolation between two meshes<br>
> ????? (nek5000-users@lists.mcs.anl.gov)<br>
> <br>
> <br>
> ----------------------------------------------------------------------<br>
> <br>
> Message: 1<br>
> Date: Mon, 27 Mar 2017 18:21:52 0200<br>
> From: nek5000-users@lists.mcs.anl.gov<br>
> To: nek5000-users@lists.mcs.anl.gov????? <nek5000-users@lists.mcs.anl.gov><br>
> Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
> Message-ID:<br>
> ??????? <mailman.3273.1490631729.2967.nek5000-users@lists.mcs.anl.gov><br>
> Content-Type: text/plain; charset=utf-8<br>
> <br>
> Just give it a whirl.<br>
> It's important for us to see if the latest version fixes your problem. I have some confidence that it will ;)<br>
> ?<br>
> -----Original message-----<br>
> > From:nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
> > Sent: Monday 27th March 2017 18:19<br>
> > To: nek5000-users@lists.mcs.anl.gov<br>
> > Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
> > <br>
> > Hi Stefan and Mattias,<br>
> > <br>
> > Thank you both for your kind advice.<br>
> > <br>
> > Stefan, I had a look at the new 'gfldr' subroutine. It seems that the spectral interpolation procedure is the same as before, however, the latest subroutine breaks down the number of fields to interpolate in one go. I still don't see why the old g2gi hangs<br>
> >? at 'findpts'. Is it mainly because of memory size? (I have been using my own copy of the legacy nek, so I'm quite keen on making it work with g2gi.)<br>
> > <br>
> > Mattias, I will give it a go with your suggestion. By reducing the polynomial order, did you mean the order of the old field or the order of the new geometry field, or both?<br>
> > <br>
> > Thank you again for both your help.<br>
> > <br>
> > Best regards,<br>
> > Tony<br>
> > <br>
> > <br>
> > -----------<br>
> > From: nek5000-users-bounces@lists.mcs.anl.gov <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users-request@lists.mcs.anl.gov <nek5000-users-request@lists.mcs.anl.gov><br>
> > Sent: 27 March 2017 12:57<br>
> > To: nek5000-users@lists.mcs.anl.gov<br>
> > Subject: Nek5000-users Digest, Vol 97, Issue 27<br>
> > Send Nek5000-users mailing list submissions to<br>
> > ??????? nek5000-users@lists.mcs.anl.gov<br>
> > <br>
> > To subscribe or unsubscribe via the World Wide Web, visit<br>
> > <br>
> <a previewremoved="true" id="LPlnk642573" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a> <<a previewremoved="true" id="LPlnk978981" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>> <<a href=""></a>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
 <<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>>><br>
> > or, via email, send a message with subject or body 'help' to<br>
> > ??????? nek5000-users-request@lists.mcs.anl.gov<br>
> > <br>
> > You can reach the person managing the list at<br>
> > ??????? nek5000-users-owner@lists.mcs.anl.gov<br>
> > <br>
> > When replying, please edit your Subject line so it is more specific<br>
> > than "Re: Contents of Nek5000-users digest..."<br>
> > <br>
> > <br>
> > Today's Topics:<br>
> > <br>
> > ?? 1. Interpolation between two meshes (nek5000-users@lists.mcs.anl.gov)<br>
> > ?? 2. Re: Interpolation between two meshes<br>
> > ????? (nek5000-users@lists.mcs.anl.gov)<br>
> > ?? 3. Re: Copy velocity slice from one boundary to another<br>
> > ????? (nek5000-users@lists.mcs.anl.gov)<br>
> > ?? 4. Robin Boundary Condition -- Navier's Slip Length<br>
> > ????? (nek5000-users@lists.mcs.anl.gov)<br>
> > ?? 5. Re: Interpolation between two meshes<br>
> > ????? (nek5000-users@lists.mcs.anl.gov)<br>
> > <br>
> > <br>
> > ----------------------------------------------------------------------<br>
> > <br>
> > Message: 1<br>
> > Date: Mon, 27 Mar 2017 04:37:46 0000<br>
> > From: nek5000-users@lists.mcs.anl.gov<br>
> > To: Nek5000 <nek5000-users@lists.mcs.anl.gov><br>
> > Subject: [Nek5000-users] Interpolation between two meshes<br>
> > Message-ID:<br>
> > ??????? <mailman.3211.1490589471.2967.nek5000-users@lists.mcs.anl.gov><br>
> > Content-Type: text/plain; charset="iso-8859-1"<br>
> > <br>
> > Dear Neks,<br>
> > <br>
> > I have been using 'g2gi' to do interpolation between different meshes. Recently, I have noticed that when I interpolate between two big meshes (say 300000 elements in each mesh), the interpolation hangs at 'call findpts' and never moves on. I tried increasing<br>
> >? number of processors and changing 'lpart', but none of these worked.<br>
> > <br>
> > I'm wondering what could cause the interpolation to hang at 'findpts'. Has anyone encountered similar issue?<br>
> > <br>
> > Hope anyone could give me some advice. Thank you very much in advance.<br>
> > <br>
> > Best regards,<br>
> > Tony<br>
> > -------------- next part --------------<br>
> > An HTML attachment was scrubbed...<br>
> > URL: < <>http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/b238e605/attachment-0001.html <<a href=""></a>http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/b238e605/attachment-0001.html <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/b238e605/attachment-0001.html">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/b238e605/attachment-0001.html</a>>>><br>
> > <br>
> > ------------------------------<br>
> > <br>
> > Message: 2<br>
> > Date: Mon, 27 Mar 2017 09:47:37 0200<br>
> > From: nek5000-users@lists.mcs.anl.gov<br>
> > To: nek5000-users@lists.mcs.anl.gov????? <nek5000-users@lists.mcs.anl.gov><br>
> > Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
> > Message-ID:<br>
> > ??????? <mailman.3212.1490600871.2967.nek5000-users@lists.mcs.anl.gov><br>
> > Content-Type: text/plain; charset=utf-8<br>
> > <br>
> > Please update to the latest master. Note, g2gi was replaced by gfldr. <br>
> > <br>
> > Stefan<br>
> > <br>
> > -----Original message-----<br>
> > > From:nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
> > > Sent: Monday 27th March 2017 6:38<br>
> > > To: Nek5000 <nek5000-users@lists.mcs.anl.gov><br>
> > > Subject: [Nek5000-users] Interpolation between two meshes<br>
> > > <br>
> > > Dear Neks,<br>
> > > <br>
> > > I have been using 'g2gi' to do interpolation between different meshes. Recently, I have?noticed that when I interpolate between two big meshes (say 300000 elements in each mesh), the interpolation hangs at 'call findpts' and never moves on. I tried increasing<br>
> > >? number of processors and changing 'lpart', but none of these worked.<br>
> > > <br>
> > > I'm wondering what could cause the interpolation to hang at 'findpts'. Has anyone encountered similar issue?<br>
> > > <br>
> > > Hope anyone could give me some advice. Thank you very much in advance.<br>
> > > <br>
> > > Best regards,<br>
> > > Tony<br>
> > > _______________________________________________<br>
> > > Nek5000-users mailing list<br>
> > > Nek5000-users@lists.mcs.anl.gov<br>
> > > <br>
> > <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a> <<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>>
 <<a href=""></a>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users <<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>>><br>
> > <br>
> > <br>
> > ------------------------------<br>
> > <br>
> > Message: 3<br>
> > Date: Mon, 27 Mar 2017 16:47:00 0700<br>
> > From: nek5000-users@lists.mcs.anl.gov<br>
> > To: nek5000-users@lists.mcs.anl.gov<br>
> > Subject: Re: [Nek5000-users] Copy velocity slice from one boundary to<br>
> > ??????? another<br>
> > Message-ID:<br>
> > ??????? <mailman.3216.1490608022.2967.nek5000-users@lists.mcs.anl.gov><br>
> > Content-Type: text/plain; charset="utf-8"<br>
> > <br>
> > Thanks a lot for replying, your answer was very useful.<br>
> > <br>
> > I look on multimesh.f file.<br>
> > As i understood, i need set subroutine 'set_inflag' to look for 'P? ' bc.<br>
> > After this i need to call 'get_values(field)' subroutine, which interpolate<br>
> > and put the value in valint(ix,iy,iz,ie,ifld) array.<br>
> > And then this array will be available in userbc in .usr file, to set it on<br>
> > inflow boundary.<br>
> > <br>
> > c-----------------------------------------------------------------------------------------------<br>
> > <br>
> > Vatslav,<br>
> > you will have to use the neknek script that allows you to exchange data<br>
> > between two Nek5000 sessions that are run simultaneously. You should have a<br>
> > look at the neknek examples to understand how to set up your usr file.<br>
> > Some of the subroutines that will be useful to you are located in the<br>
> > following link: <br>
> > <a href="http://www.mcs.anl.gov/research/projects/nek5/docs/">http://www.mcs.anl.gov/research/projects/nek5/docs/</a> <<a href="http://www.mcs.anl.gov/research/projects/nek5/docs/">http://www.mcs.anl.gov/research/projects/nek5/docs/</a>> <<a href=""></a>http://www.mcs.anl.gov/research/projects/nek5/docs/
 <<a href="http://www.mcs.anl.gov/research/projects/nek5/docs/">http://www.mcs.anl.gov/research/projects/nek5/docs/</a>>><br>
> > html/multimesh_8f_source.html<br>
> > You will have to modify some of them to fit your needs since for example<br>
> > the subroutine ?set_inflag? is written to look for ?in? boundary conditions<br>
> > and replace them with ?v? ?. You will need a subroutine that looks for ?in?<br>
> > and substitute it with either ?v? ? or ?P? ? depending on the which pipe<br>
> > you are at.<br>
> > Marco<br>
> > <br>
> > Hi, Neks!<br>
> > I have two pipes in my simulation. They are placed near each other on one<br>
> > axis with empty space between.<br>
> > One of them with periodic boundary conditions on z direction, another one -<br>
> > inflow -> outflow.<br>
> > How can I copy velocity slice from one of boundary on first pipe with<br>
> > periodic bc to the inflow boundary on second pipe?<br>
> > Regards, Vatslav<br>
> > <br>
> > c-----------------------------------------------------------------------------------------------<br>
> > <br>
> > Regards, Vatslav<br>
> > -------------- next part --------------<br>
> > An HTML attachment was scrubbed...<br>
> > URL: < <>http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/92c78de8/attachment-0001.html <<a href=""></a>http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/92c78de8/attachment-0001.html <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/92c78de8/attachment-0001.html">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/92c78de8/attachment-0001.html</a>>>><br>
> > <br>
> > ------------------------------<br>
> > <br>
> > Message: 4<br>
> > Date: Mon, 27 Mar 2017 11:31:53 0200<br>
> > From: nek5000-users@lists.mcs.anl.gov<br>
> > To: nek5000-users@lists.mcs.anl.gov<br>
> > Subject: [Nek5000-users] Robin Boundary Condition -- Navier's Slip<br>
> > ??????? Length<br>
> > Message-ID:<br>
> > ??????? <mailman.3222.1490614366.2967.nek5000-users@lists.mcs.anl.gov><br>
> > Content-Type: text/plain; charset="utf-8"<br>
> > <br>
> > Hi Neks,<br>
> > I wish to simulate a transitional channel flow whose wall respect a<br>
> > "slip-length" boundary condition, such as:<br>
> > <br>
> > $ [u? \delta \frac{\partial u}{\partial n}]_{wall}=0$<br>
> > <br>
> > Where \delta is the slip length<br>
> > <br>
> > ^n<br>
> > |??????? u?????????????? u<br>
> > |????? |-->/???????? |----->/<br>
> > |????? |?? /?????????? |????? /<br>
> > |????? |? /??????????? |???? /<br>
> > |????? | /???????????? |??? /<br>
> > |___|/_______|?? /_____<br>
> > ?????????????????????? | /? ^<br>
> > ?????????????????????? |/?? v \delta<br>
> > <br>
> > No-slip?????? Slip-length<br>
> > <br>
> > (a better sketch can be found here:<br>
> > <a href="http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html">http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html</a> <<a href="http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html">http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html</a>>
 <<a href=""></a>http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html <<a href="http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html">http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html</a>>><br>
>  )<br>
> > <br>
> >? <<a href=""></a>http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html <<a href="http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html">http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html</a>>><br>
> > Figure 1 : Low-friction flows of liquid at nanopatterned ... <<a href=""></a>http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html <<a href="http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html">http://www.nature.com/nmat/journal/v2/n4/fig_tab/nmat857_F1.html</a>>>www.nature.com<br>
> > Nature Materials journal covers a range of topics within materials science, from materials engineering and structural materials (metals, alloys, ceramics, composites ...<br>
> > <br>
> > <br>
> > <br>
> > I would like to know if this Robin boundary condition for velocity field is<br>
> > already implemented in the current code.<br>
> > <br>
> > Thank you<br>
> > <br>
> > Ciccio<br>
> > -------------- next part --------------<br>
> > An HTML attachment was scrubbed...<br>
> > URL: < <>http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/9922e2e7/attachment-0001.html <<a href=""></a>http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/9922e2e7/attachment-0001.html <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/9922e2e7/attachment-0001.html">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/9922e2e7/attachment-0001.html</a>>>><br>
> > <br>
> > ------------------------------<br>
> > <br>
> > Message: 5<br>
> > Date: Mon, 27 Mar 2017 13:56:02 0200<br>
> > From: nek5000-users@lists.mcs.anl.gov<br>
> > To: nek5000-users@lists.mcs.anl.gov<br>
> > Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
> > Message-ID:<br>
> > ??????? <mailman.3230.1490615875.2967.nek5000-users@lists.mcs.anl.gov><br>
> > Content-Type: text/plain; charset="windows-1252"; Format="flowed"<br>
> > <br>
> > Hi Tony,<br>
> > I have encountered the same issue with g2gi when interpolating very <br>
> > large meshes. I don't know what is causing the issue but a workaround is <br>
> > to reduce the polynomial order. You obviously loose accuracy by doing <br>
> > this, but I still find it better than starting from scratch with the new <br>
> > mesh. Hope this helps.<br>
> > <br>
> > Best,<br>
> > Mattias<br>
> > <br>
> > <br>
> > On 03/27/2017 06:37 AM, nek5000-users@lists.mcs.anl.gov wrote:<br>
> > > Dear Neks,<br>
> > ><br>
> > > I have been using 'g2gi' to do interpolation between different meshes. <br>
> > > Recently, I have noticed that when I interpolate between two big <br>
> > > meshes (say 300000 elements in each mesh), the interpolation hangs at <br>
> > > 'call findpts' and never moves on. I tried increasing number of <br>
> > > processors and changing 'lpart', but none of these worked.<br>
> > ><br>
> > > I'm wondering what could cause the interpolation to hang at 'findpts'. <br>
> > > Has anyone encountered similar issue?<br>
> > ><br>
> > > Hope anyone could give me some advice. Thank you very much in advance.<br>
> > ><br>
> > > Best regards,<br>
> > > Tony<br>
> > ><br>
> > ><br>
> > > _______________________________________________<br>
> > > Nek5000-users mailing list<br>
> > > Nek5000-users@lists.mcs.anl.gov<br>
> > > <br>
> > <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a> <<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>>
 <<a href=""></a>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users <<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>>><br>
> > <br>
> > <br>
> > -------------- next part --------------<br>
> > An HTML attachment was scrubbed...<br>
> > URL: < <>http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/031d792f/attachment.html <<a href=""></a>http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/031d792f/attachment.html <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/031d792f/attachment.html">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170327/031d792f/attachment.html</a>>>><br>
> > <br>
> > ------------------------------<br>
> > <br>
> > _______________________________________________<br>
> > Nek5000-users mailing list<br>
> > Nek5000-users@lists.mcs.anl.gov<br>
> > <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a> <<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>>
 <<a href=""></a>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users <<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>>><br>
> > <br>
> > <br>
> > End of Nek5000-users Digest, Vol 97, Issue 27<br>
> > *********************************************<br>
> > _______________________________________________<br>
> > Nek5000-users mailing list<br>
> > Nek5000-users@lists.mcs.anl.gov<br>
> > <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a> <<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>><br>
> <br>
> <br>
> ------------------------------<br>
> <br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> Nek5000-users@lists.mcs.anl.gov<br>
> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a> <<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>><br>
> <br>
> <br>
> End of Nek5000-users Digest, Vol 97, Issue 29<br>
> *********************************************<br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> Nek5000-users@lists.mcs.anl.gov<br>
> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Tue, 28 Mar 2017 12:18:56 +0000<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: "nek5000-users@lists.mcs.anl.gov"<br>
        <nek5000-users@lists.mcs.anl.gov><br>
Subject: Re: [Nek5000-users] Copy velocity slice from one boundary to<br>
        another<br>
Message-ID:<br>
        <mailman.3877.1490703543.2967.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Vatslav,<br>
<br>
what I do and it seems to work for my case is the following. Let assume you have two pipes (pipe 1 with periodic bc and pipe 2 with inlet/outlet), and that the flow is moving from pipe 1 to pipe 2. The solution from pipe 1 is interpolated to pipe 2 and used
 as inlet condition.<br>
<br>
In a separate subroutine, do the following:<br>
<br>
  *   in pipe 1, look for the boundary faces with P boundary condition that will be used for the interpolation and substitute P with inP.<br>
  *   in pipe 2, look for the boundary faces with v boundary condition and substitute v with inV.<br>
<br>
In a subroutine analogous to set_inflag:<br>
<br>
  *   in pipe 1, collect boundary faces with flag inP, set intflag to 1 for this boundary face, and substitute back boundary condition to P.<br>
  *   in pip2, do the same: inV - > v<br>
<br>
After, this step, you have marked all boundary faces that will be involved in the interpolation process in Nek5000. From this point, the steps are similar to what is done in the examples. All of this is done in my usr file.<br>
<br>
Hope this help,<br>
<br>
Marco<br>
<br>
On Mar 27, 2017, at 5:47 AM, nek5000-users@lists.mcs.anl.gov<mailto:nek5000-users@lists.mcs.anl.gov> wrote:<br>
<br>
Thanks a lot for replying, your answer was very useful.<br>
<br>
I look on multimesh.f file.<br>
As i understood, i need set subroutine 'set_inflag' to look for 'P  ' bc.<br>
After this i need to call 'get_values(field)' subroutine, which interpolate and put the value in valint(ix,iy,iz,ie,ifld) array.<br>
And then this array will be available in userbc in .usr file, to set it on inflow boundary.<br>
<br>
c-----------------------------------------------------------------------------------------------<br>
<br>
Vatslav,<br>
you will have to use the neknek script that allows you to exchange data between two Nek5000 sessions that are run simultaneously. You should have a look at the neknek examples to understand how to set up your usr file.<br>
Some of the subroutines that will be useful to you are located in the following link:
<a href="http://www.mcs.anl.gov/research/projects/nek5/docs/html/multimesh_8f_source.html">
http://www.mcs.anl.gov/research/projects/nek5/docs/html/multimesh_8f_source.html</a><br>
You will have to modify some of them to fit your needs since for example the subroutine ?set_inflag? is written to look for ?in? boundary conditions and replace them with ?v  ?. You will need a subroutine that looks for ?in? and substitute it with either ?v 
 ? or ?P  ? depending on the which pipe you are at.<br>
Marco<br>
<br>
Hi, Neks!<br>
I have two pipes in my simulation. They are placed near each other on one axis with empty space between.<br>
One of them with periodic boundary conditions on z direction, another one - inflow -> outflow.<br>
How can I copy velocity slice from one of boundary on first pipe with periodic bc to the inflow boundary on second pipe?<br>
Regards, Vatslav<br>
<br>
c-----------------------------------------------------------------------------------------------<br>
<br>
Regards, Vatslav<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<mailto:Nek5000-users@lists.mcs.anl.gov><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170328/98fafd25/attachment.html">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170328/98fafd25/attachment.html</a>><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
End of Nek5000-users Digest, Vol 97, Issue 31<br>
*********************************************<br>
</div>
</span></font></div>
</div>
</body>
</html>