<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. Please contact me via wzx1989221@gmail.com.</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> 30 March 2017 13:17<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Nek5000-users Digest, Vol 97, Issue 37</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;"></span></font><font size="2"><span style="font-size:10pt;">
<div class="PlainText"><br>
Today's Topics:<br>
<br>
1. Outflow BC with negative velocity<br>
(nek5000-users@lists.mcs.anl.gov)<br>
2. Re: NekNek with PN-PN (nek5000-users@lists.mcs.anl.gov)<br>
3. Re: NekNek with PN-PN (nek5000-users@lists.mcs.anl.gov)<br>
4. Re: Interpolation between two meshes<br>
(nek5000-users@lists.mcs.anl.gov)<br>
5. Accessing the value of drag and lift coefficients in usrchk<br>
(nek5000-users@lists.mcs.anl.gov)<br>
<br>
<br>
<br>
Message: 4<br>
Date: Thu, 30 Mar 2017 11:46:40 +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.4068.1490867217.2967.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Please let me know how I can contact you off list.<br>
<br>
Stefan?<br>
<br>
On 29 Mar 2017, at 14:59, "nek5000-users@lists.mcs.anl.gov <<a href="mailto:nek5000-users@lists.mcs.anl.gov">mailto:nek5000-users@lists.mcs.anl.gov</a>> " <nek5000-users@lists.mcs.anl.gov <<a href="mailto:nek5000-users@lists.mcs.anl.gov">mailto:nek5000-users@lists.mcs.anl.gov</a>>
> wrote:<br>
<br>
Hi Stefan,<br>
<br>
Thank you very much for your quick reply.<br>
<br>
1. Yes, my source field is in double precision since I need to keep the accuracy.<br>
<br>
2. I added several print lines in gfldr.f to check where the interpolation hangs (see below). I also put 'call nekgsync()' before and after findpts as you suggested. I can confirm that findpts is exactly where the interpolation hangs. The last output in the
logfile is 'CHECK interp findpts before' and it never gets updated again.<br>
<br>
3. I've just tested gfldr() with turbChannel example. Everything worked fine and the?source field was interpolated onto the turbChannel mesh (512 elements). So I moved on increasing the number of elements of turbChannel case.<br>
<br>
First, I used 234300 elements?in the turbChannel example (half of my current simulation size) and the interpolation was done successfully. Then, I used 468600 elements for?the turbChannel example (similar to my current simulation size), the interpolation hangs
again at findpts.<br>
<br>
I guess it might be something to do with the memory. What do you think? Is there any solution to this?<br>
<br>
Best regards,<br>
Tony<br>
<br>
??????? if(nid==0) print*,'CHECK interp findpts before'<br>
??????? call nekgsync()<br>
??????? call findpts(inth_gfldr,<br>
???? &?????????????? rcode,1,<br>
???? &?????????????? proc,1,<br>
???? &?????????????? elid,1,<br>
???? &?????????????? rst,ndim,<br>
???? &?????????????? dist,1,<br>
???? &?????????????? xm1,1,<br>
???? &?????????????? ym1,1,<br>
???? &?????????????? zm1,1,ntot)<br>
??????? call nekgsync()<br>
??????? if(nid==0) print*,'CHECK interp findpts after'<br>
</div>
</span></font></div>
</div>
</body>
</html>