<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>In case the attachments didn't go through, please see below for the .usr file and log file.<br>
</p>
<p>I only used 'userchk' and 'usrdat' in the .usr file for interpolation. The other subroutines are empty.<br>
</p>
<p>(I added several print lines in 'gfldr.f' to see where the interpolation halted.)</p>
<p><br>
</p>
<p>Best regards,</p>
<p>Tony<br>
</p>
<br>
<p>.USR FILE:</p>
<p></p>
<div>c-----------------------------------------------------------------------<br>
subroutine userchk()<br>
include 'SIZE'<br>
include 'TOTAL'<br>
include 'ZPER'<br>
include 'CTIMER'<br>
include 'RESTART'<br>
<br>
character*132 sourcefld<br>
<br>
sourcefld='source.f00001'<br>
call gfldr(sourcefld)<br>
call outpost(vx,vy,vz,pr,t,'rst')<br>
call exitt()<br>
<br>
return<br>
end<br>
c-----------------------------------------------------------------------<br>
subroutine usrdat() ! This routine to modify element vertices<br>
include 'SIZE'<br>
include 'TOTAL'<br>
<br>
parameter (ielx=150)<br>
parameter (iely=22)<br>
parameter (ielz=142)<br>
integer i,e,eg,ex,ey,ez<br>
real cony<br>
<br>
cony = 2.36d0<br>
<br>
do e=1,nelv<br>
eg = lglel(e)<br>
if(eg<=ielx*iely*ielz) then<br>
call get_exyz(ex,ey,ez,eg,ielx,iely,ielz)<br>
<br>
do i=1,8<br>
if(i<=2 .or. i==5 .or. i==6) then<br>
if((yc(i,e)>0.0e0.and.yc(i,e)<1.0e0).OR.<br>
& (yc(i,e)>1.0e0.and.yc(i,e)<2.0e0)) then<br>
yc(i,e) = 2.0e0*Real(ey-1)/Real(iely)<br>
end if<br>
else if(i==3 .or. i==4 .or. i>=7) then<br>
if((yc(i,e)>0.0e0.and.yc(i,e)<1.0e0).OR.<br>
& (yc(i,e)>1.0e0.and.yc(i,e)<2.0e0)) then<br>
yc(i,e) = 2.0e0*Real(ey)/Real(iely)<br>
end if<br>
end if<br>
end do<br>
<br>
end if<br>
end do<br>
<br>
do e=1,nelv<br>
eg = lglel(e)<br>
if(eg==1) print*,yc(1,e),yc(2,e),yc(5,e),yc(6,e)<br>
if(eg==1) print*,yc(3,e),yc(4,e),yc(7,e),yc(8,e)<br>
do i=1,8<br>
if(yc(i,e)>0.0e0 .AND. yc(i,e)<1.0e0) then<br>
yc(i,e) = yc(i,e)-1.0e0<br>
yc(i,e) = tanh(cony*yc(i,e))/tanh(cony)<br>
yc(i,e) = yc(i,e)+1.0e0<br>
else if(yc(i,e)>1.0e0 .AND. yc(i,e)<2.0e0) then<br>
yc(i,e) = yc(i,e)-1.0e0<br>
yc(i,e) = tanh(cony*yc(i,e))/tanh(cony)<br>
yc(i,e) = yc(i,e)+1.0e0<br>
end if<br>
end do<br>
end do<br>
<br>
return<br>
end<br>
c-----------------------------------------------------------------------<br>
</div>
<p></p>
<br>
<br>
<p>LOG FILE:<br>
</p>
<br>
<div>/----------------------------------------------------------\\<br>
| _ __ ______ __ __ ______ ____ ____ ____ |<br>
| / | / // ____// //_/ / ____/ / __ \\/ __ \\/ __ \\ |<br>
| / |/ // __/ / ,< /___ \\ / / / // / / // / / / |<br>
| / /| // /___ / /| | ____/ / / /_/ // /_/ // /_/ / |<br>
| /_/ |_//_____//_/ |_|/_____/ \\___/ \\___/ \\___/ |<br>
| |<br>
|----------------------------------------------------------|<br>
| |<br>
| NEK5000: Open Source Spectral Element Solver |<br>
| COPYRIGHT (c) 2008-2017 UCHICAGO ARGONNE, LLC |<br>
| Version: 17.0.0 |<br>
| Web: http://nek5000.mcs.anl.gov |<br>
| |<br>
\\----------------------------------------------------------/<br>
<br>
<br>
Number of processors: 720<br>
REAL wdsize : 8<br>
INTEGER wdsize : 4<br>
Timer accuracy : 9.54E-08<br>
Reading /work/e01/e01/tony/geometry_nonosci_cavity/interp.rea<br>
<br>
mapping elements to processors<br>
Reading /work/e01/e01/tony/geometry_nonosci_cavity/interp.map<br>
element load imbalance: 1 652 653<br>
done :: mapping elements to processors<br>
<br>
0 objects found<br>
118 Parameters from file:<br>
1 1.00000 p001 DENSITY<br>
2 -3150.00 p002 VISCOS<br>
12 -0.11750E-01 p012 DT<br>
15 500.000 p015 IOSTEP<br>
21 0.100000E-07 p021 DIVERGENCE<br>
22 0.100000E-07 p022 HELMHOLTZ<br>
26 2.00000 p026 COURANT/NTAU<br>
27 3.00000 p027 TORDER<br>
59 0.00000 p059 !=0 --> full Jac. eval. for each el.<br>
63 8.00000 p063 =8 --> force 8-byte output<br>
65 1.00000 p065 #iofiles (eg, 0 or 64); <0 --> sep. dirs<br>
66 6.00000 p066 output : <0=ascii, else binary<br>
67 6.00000 p067 restart: <0=ascii, else binary<br>
93 20.0000 p093 Number of previous pressure solns saved<br>
94 5.00000 p094 start projecting velocity after p94 step<br>
95 5.00000 p095 start projecting pressure after p95 step<br>
99 3.00000 p099 dealiasing: <0--> off/3--> old/4--> new<br>
103 0.10000E-01 p103 weight of stabilizing filter (.01)<br>
<br>
nelgt/nelgv/lelt: 469500 469500 653<br>
lx1 /lx2 /lx3 : 8 8 8<br>
done :: read .rea file 45.222 sec<br>
<br>
setup mesh topology<br>
Right-handed check complete for 469500 elements. OK.<br>
setvert3d: 8 60855000 162267000 60855000 60855000<br>
call usrsetvert<br>
done :: usrsetvert<br>
<br>
gs_setup: 10325502 unique labels shared<br>
pairwise times (avg, min, max): 0.00044986 0.00033021 0.000500798<br>
crystal router : 0.00145524 0.00143301 0.0014781<br>
used all_to_all method: pairwise<br>
handle bytes (avg, min, max): 2.64212e+06 2463020 2796252<br>
buffer bytes (avg, min, max): 491593 352288 632640<br>
setupds time 5.3504E-01 seconds 0 8 60855000 469500<br>
8 max multiplicity<br>
done :: setup mesh topology<br>
<br>
0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000<br>
0.000000000000000E+000<br>
9.090909090909091E-002 9.090909090909091E-002 9.090909090909091E-002<br>
9.090909090909091E-002<br>
call usrdat<br>
done :: usrdat<br>
<br>
generate geometry data<br>
NOTE: All elements deformed , param(59) ^=0<br>
done :: generate geometry data<br>
<br>
call usrdat2<br>
done :: usrdat2<br>
<br>
regenerate geometry data 1<br>
NOTE: All elements deformed , param(59) ^=0<br>
done :: regenerate geometry data 1<br>
<br>
verify mesh topology<br>
0.000000000000000E+000 6.00000000000000 Xrange<br>
-0.120000000000000 2.00000000000000 Yrange<br>
0.000000000000000E+000 49.7000000000000 Zrange<br>
done :: verify mesh topology<br>
<br>
IFTRAN = T<br>
IFFLOW = T<br>
IFHEAT = F<br>
IFSPLIT = T<br>
IFLOMACH = F<br>
IFUSERVP = F<br>
IFUSERMV = F<br>
IFSTRS = F<br>
IFCHAR = T<br>
IFCYCLIC = F<br>
IFAXIS = F<br>
IFMVBD = F<br>
IFMELT = F<br>
IFMODEL = F<br>
IFKEPS = F<br>
IFMOAB = F<br>
IFNEKNEK = F<br>
IFSYNC = T<br>
<br>
IFVCOR = F<br>
IFINTQ = F<br>
IFCWUZ = F<br>
IFSWALL = F<br>
IFGEOM = F<br>
IFSURT = F<br>
IFWCNO = F<br>
<br>
IFTMSH for field 1 = F<br>
IFADVC for field 1 = T<br>
IFNONL for field 1 = F<br>
<br>
Dealiasing enabled, lxd= 12<br>
verify mesh topology<br>
0.000000000000000E+000 6.00000000000000 Xrange<br>
-0.120000000000000 2.00000000000000 Yrange<br>
0.000000000000000E+000 49.7000000000000 Zrange<br>
done :: verify mesh topology<br>
<br>
call usrdat3<br>
done :: usrdat3<br>
</div>
<br>
<div> set initial conditions<br>
nekuic (1) for ifld 1<br>
call nekuic for vel<br>
xyz min 0.0000 -0.12000 0.0000<br>
uvwpt min 0.0000 0.0000 0.0000 0.0000 0.0000<br>
xyz max 6.0000 2.0000 49.700<br>
uvwpt max 0.0000 0.0000 0.0000 0.0000 0.0000<br>
done :: set initial conditions<br>
<br>
gridpoints unique/tot: 162267000 240384000<br>
dofs: 159972750 162104250<br>
<br>
Initialization successfully completed 46.588 sec<br>
<br>
nsteps=0 -> skip time loop<br>
running solver in post processing mode<br>
<br>
call userchk<br>
call gfldr<br>
byte swap: F 6.543210 -2.9312772E+35<br>
sourcefld X 0.000000000000000E+000 6.00000000000000<br>
sourcefld Y -0.120000000000000 2.00000000000000<br>
sourcefld Z 0.000000000000000E+000 49.7000000000000<br>
CHECK interp prepare<br>
CHECK interp read<br>
<br>
</div>
<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 20:51<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Nek5000-users Digest, Vol 97, Issue 33</font> </div>
</div>
<font size="2"><span style="font-size:10pt;"></span></font><font size="2"><span style="font-size:10pt;">
<div class="PlainText"><br>
Message: 3<br>
Date: Tue, 28 Mar 2017 19:51:04 +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] Interpolation between two meshes<br>
Message-ID:<br>
<mailman.3915.1490730701.2967.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi Stefan,<br>
<br>
<br>
Many thanks for your reply. I'm running my simulations on ARCHER (Cray XC30 supercomputer) and I'm using the Intel compiler.<br>
<br>
<br>
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.<br>
<br>
<br>
Please have a look and let me know if you need further information about the case. Thank you very much for you help.<br>
<br>
<br>
Best regards,<br>
<br>
Tony<br>
</div>
</span></font></div>
</div>
</body>
</html>