<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style>
<!--
@font-face
{font-family:Calibri}
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New"}
a:link, span.MsoHyperlink
{color:blue;
text-decoration:underline}
a:visited, span.MsoHyperlinkFollowed
{color:purple;
text-decoration:underline}
span.EmailStyle17
{font-family:"Calibri","sans-serif";
color:windowtext}
.MsoChpDefault
{font-family:"Calibri","sans-serif"}
@page WordSection1
{margin:72.0pt 72.0pt 72.0pt 72.0pt}
-->
</style><style type="text/css" id="owaParaStyle"></style>
</head>
<body lang="EN-US" link="blue" vlink="purple" fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div><br>
</div>
Dear Donato,
<div><br>
</div>
<div>I've not looked closely at this yet. However, I suspect it might be related to the projection.</div>
<div>Can you run your example with param 94 = 0 in the .rea file ? This should turn off the</div>
<div>velocity/passive scalar projection.</div>
<div><br>
</div>
<div>Another question ---</div>
<div><br>
</div>
<div>in your .rea file, where it says "IFADVC" there should be two lines like:</div>
<div><br>
</div>
<div>
<div><br>
</div>
<div> T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)</div>
<div> F F T T T T T T T T T T IFTMSH (IF mesh for this field is T mesh)</div>
</div>
<div><br>
</div>
<div><br>
</div>
<div>For your cases, you probably want them:</div>
<div><br>
</div>
<div>
<div>
<div><br class="Apple-interchange-newline">
T T T T T T T T T T T T IFNAV & IFADVC (convection in P.S. fields)</div>
<div> F F T F F F F F F F F F F IFTMSH (IF mesh for this field is T mesh)</div>
</div>
<div><br>
</div>
</div>
<div>This would do the following:</div>
<div><br>
</div>
<div>1) It turns on advection for all fields (fluid, temperature, and passive scalars)</div>
<div><br>
</div>
<div>2) It sets the mesh field to be Fluid</div>
<div> It sets the velocity field to be Fluid</div>
<div> It sets the temperature to be Temperature (i.e., fluid + solid)</div>
<div> It sets the remaining passive scalars to be Fluid</div>
<div><br>
</div>
<div>For conjugate heat transfer, normally only the temperature acts over both</div>
<div>the fluid and solid element set, which is why I set those flags as above.</div>
<div><br>
</div>
<div>Let me know if this resolves the issue for you... </div>
<div><br>
</div>
<div>Regards,</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF122718" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> nek5000-users-bounces@lists.mcs.anl.gov [nek5000-users-bounces@lists.mcs.anl.gov] on behalf of nek5000-users@lists.mcs.anl.gov [nek5000-users@lists.mcs.anl.gov]<br>
<b>Sent:</b> Friday, March 28, 2014 4:34 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Problem with passive scalar fields - subroutine outfld<br>
</font><br>
</div>
<div></div>
<div>
<div class="WordSection1">
<p class="MsoNormal" style="margin-right:35.65pt">Dear Neks, </p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">Sorry for my previous mail.</p>
<p class="MsoNormal" style="margin-right:35.65pt">I hope that the layout will be good now.</p>
<p class="MsoNormal" style="margin-right:35.65pt">_________________________________________________________________</p>
<p class="MsoNormal" style="margin-right:35.65pt">Dear users and developers,</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">I’m simulating a 3D conjugate heat transfer case, considering temperature as a passive scalar.</p>
<p class="MsoNormal" style="margin-right:35.65pt">When I simulate only one temperature field everything is ok.</p>
<p class="MsoNormal" style="margin-right:35.65pt">If I add more passive scalar fields (each one identical to the first one), after the subroutine prepost.f/outfld is called, then some of these fields are not good.</p>
<p class="MsoNormal" style="margin-right:35.65pt">I ran up to 8 additional fields, and results are the followings (“--“ means that the field is not good):</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> ldimt T PS1 PS2 PS3 PS4 PS5 PS6 PS7 PS8</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 1 ok</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 2 ok --</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 3 ok -- ok</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 4 ok -- ok ok</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 5 ok -- -- -- ok</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 6 ok -- -- -- ok ok</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 7 ok -- -- -- ok ok ok</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 8 ok -- -- -- ok ok ok ok</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 9 ok -- -- -- ok ok ok ok ok</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">I’m using PN-PN.</p>
<p class="MsoNormal" style="margin-right:35.65pt">Nothing changes with the number of mesh points or the number of steps or the spatial polynomial order.</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">I saw differences with Visit, with mean temperature profiles, and with average temperature in fluid and solid domain.</p>
<p class="MsoNormal" style="margin-right:35.65pt">Differences are noticeable also in logfile. Consider the following example, in which the first field file is dumped at step 250; at step 251 the changes are noticeable in the number of iterations per field only
for PS1, PS2 and PS3.</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">::</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">Step 250, t= 3.19E+00, DT= 1.09E-02, C= 0.571 7.85E+01 2.6463E-01</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> </span>
Solving for heat</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 250 Hmholtz TEMP: 17 9.2949E-04 4.5939E-03 1.0000E-03</p>
<p class="MsoNormal" style="margin-right:35.65pt"> <span lang="IT">250 Hmholtz PS 1: 17 9.2949E-04 4.5939E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 250 Hmholtz PS 2: 17 9.2949E-04 4.5939E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 250 Hmholtz PS 3: 17 9.2949E-04 4.5939E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 250 Hmholtz PS 4: 17 9.2949E-04 4.5939E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 250 Hmholtz PS 5: 17 9.2949E-04 4.5939E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 250 Hmholtz PS 6: 17 9.2949E-04 4.5939E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 250 Hmholtz PS 7: 17 9.2949E-04 4.5939E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 250 Hmholtz PS 8: 17 9.2949E-04 4.5939E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> </span>250 3.1937E+00 2.0457E-01 Heat done</p>
<p class="MsoNormal" style="margin-right:35.65pt"> Solving for fluid</p>
<p class="MsoNormal" style="margin-right:35.65pt">::</p>
<p class="MsoNormal" style="margin-right:35.65pt">call outfld: ifpsco: T</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 250 3.1937E+00 Write checkpoint:</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 0 250 OPEN: cjh_tc0.f00001</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 250 3.1937E+00 done :: Write checkpoint</p>
<p class="MsoNormal" style="margin-right:35.65pt"> file size = 5.6 MB</p>
<p class="MsoNormal" style="margin-right:35.65pt"> avg data-throughput = 367.7MB/s</p>
<p class="MsoNormal" style="margin-right:35.65pt"> <span lang="IT">io-nodes = 1</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">Step 251, t= 3.20E+00, DT= 1.09E-02, C= 0.566 7.88E+01 2.93E-01</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> </span>
Solving for heat</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 251 Hmholtz TEMP: 17 8.3175E-04 4.6145E-03 1.0000E-03</p>
<p class="MsoNormal" style="margin-right:35.65pt"> <span lang="IT">251 Hmholtz PS 1: 74 9.1166E-04 1.1888E+02 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 251 Hmholtz PS 2: 85 9.4452E-04 9.1977E+02 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 251 Hmholtz PS 3: 91 8.5561E-04 2.3562E+03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 251 Hmholtz PS 4: 17 8.3175E-04 4.6145E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 251 Hmholtz PS 5: 17 8.3175E-04 4.6145E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 251 Hmholtz PS 6: 17 8.3175E-04 4.6145E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 251 Hmholtz PS 7: 17 8.3175E-04 4.6145E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 251 Hmholtz PS 8: 17 8.3175E-04 4.6145E-03 1.0000E-03</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> </span>251 3.2046E+00 3.6832E-01 Heat done</p>
<p class="MsoNormal" style="margin-right:35.65pt">::</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">Looking to the outfld subroutine I wasn’t able to figure out the reason for that.</p>
<p class="MsoNormal" style="margin-right:35.65pt">Does anyone know how to solve this problem?</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">Thanks in advance.</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">Relevant parts of logfile, SIZE, .usr, .box and .rea files are attached at the end of the mail.</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">Best Regards,</p>
<p class="MsoNormal" style="margin-right:35.65pt">Donato</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">______________________________________________________________________</p>
<p class="MsoNormal" style="margin-right:35.65pt">SIZE file</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">parameter (ldimt= 9) ! passive scalars + T</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">______________________________________________________________________</p>
<p class="MsoNormal" style="margin-right:35.65pt">*.box file</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">base.rea</p>
<p class="MsoNormal" style="margin-right:35.65pt">3 spatial dimension ( < 0 --> generate .rea/.re2 pair)</p>
<p class="MsoNormal" style="margin-right:35.65pt">10 number of fields</p>
<p class="MsoNormal" style="margin-right:35.65pt">#======</p>
<p class="MsoNormal" style="margin-right:35.65pt">Box</p>
<p class="MsoNormal" style="margin-right:35.65pt">-10 14 -5 nelx,nely,nelz for Box</p>
<p class="MsoNormal" style="margin-right:35.65pt">0. 6. 1. x0,x1,gain</p>
<p class="MsoNormal" style="margin-right:35.65pt">-1.00 -0.98 -0.94 -0.86 -0.70 -0.50 -0.25 -0.00 0.25</p>
<p class="MsoNormal" style="margin-right:35.65pt">0.50 0.70 0.86 0.94 0.98 1.00 y0,y1,...,yn</p>
<p class="MsoNormal" style="margin-right:35.65pt">0. 3. 1. z0,z1,gain</p>
<p class="MsoNormal" style="margin-right:35.65pt">P ,P ,W ,W ,P ,P vel BCs
</p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,E ,P ,P temp BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,E ,P ,P PS1 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,E ,P ,P PS2 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,E ,P ,P PS3 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,E ,P ,P PS4 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,E ,P ,P PS5 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,E ,P ,P PS6 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,E ,P ,P PS7 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,E ,P ,P PS8 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> </span></p>
<p class="MsoNormal" style="margin-right:35.65pt">Box_up</p>
<p class="MsoNormal" style="margin-right:35.65pt">-10 7 -5 nelx,nely,nelz for Box_up</p>
<p class="MsoNormal" style="margin-right:35.65pt">0. 6. 1. x0,x1,gain</p>
<p class="MsoNormal" style="margin-right:35.65pt">1.00 1.02 1.06 1.14 1.30 1.50 1.75 2.00 y0,y1,...,yn</p>
<p class="MsoNormal" style="margin-right:35.65pt">0. 3. 1. z0,z1,gain</p>
<p class="MsoNormal" style="margin-right:35.65pt"> , , , , , vel BCs
</p>
<p class="MsoNormal" style="margin-right:35.65pt">P ,P ,I ,E ,P ,P temp BCs
</p>
<p class="MsoNormal" style="margin-right:35.65pt">P ,P ,I ,E ,P ,P PS1 BCs
</p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,I ,E ,P ,P PS2 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,I ,E ,P ,P PS3 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,I ,E ,P ,P PS4 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,I ,E ,P ,P PS5 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,I ,E ,P ,P PS6 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,I ,E ,P ,P PS7 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,I ,E ,P ,P PS8 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> </span></p>
<p class="MsoNormal" style="margin-right:35.65pt">Box_down</p>
<p class="MsoNormal" style="margin-right:35.65pt">-10 7 -5 nelx,nely,nelz for Box_down</p>
<p class="MsoNormal" style="margin-right:35.65pt">0. 6. 1. x0,x1,gain</p>
<p class="MsoNormal" style="margin-right:35.65pt">-2.00 -1.75 -1.50 -1.30 -1.14 -1.06 -1.02 -1.00 y0,y1,...,yn</p>
<p class="MsoNormal" style="margin-right:35.65pt">0. 3. 1. z0,z1,gain</p>
<p class="MsoNormal" style="margin-right:35.65pt"> , , , , , vel BCs
</p>
<p class="MsoNormal" style="margin-right:35.65pt">P ,P ,E ,I ,P ,P temp BCs
</p>
<p class="MsoNormal" style="margin-right:35.65pt">P ,P ,E ,I ,P ,P PS1 BCs
</p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,I ,P ,P PS2 BCs </span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,I ,P ,P PS3 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,I ,P ,P PS4 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,I ,P ,P PS5 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,I ,P ,P PS6 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,I ,P ,P PS7 BCs
</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT">P ,P ,E ,I ,P ,P PS8 BCs</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> </span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> </span></p>
<p class="MsoNormal" style="margin-right:35.65pt">______________________________________________________________________</p>
<p class="MsoNormal" style="margin-right:35.65pt">logfile</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">103 Parameters</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 1 1.00000 p1 DENSITY</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 2 -2812.94 p2 Re_b</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 7 1.00000 p7 RHOCP</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 8 1.00000 p8 CONDUCT</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 11 500. p11 NSTEPS</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 12 0.05 p12 DT</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 15 250. p15 IOSTEP</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 21 1.000000E-06 p21 DIVERGENCE</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 22 1.000000E-03 p22 HELMHOLTZ</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 23 8 p23 NPSCAL</p>
<p class="MsoNormal" style="margin-right:35.65pt"> <span lang="IT">24 0.100000E-01 p24 TOLREL</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 25 0.100000E-01 p25 TOLABS</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> </span>26 0.50000 p26 TARGET COURANT NUMBER</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 27 3.00000 p27 TORDER</p>
<p class="MsoNormal" style="margin-right:35.65pt"> <span lang="IT">28 0.000000E+00 p28</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> 30 1.000000E+00 p30 1: user-def props via USERVP subr.</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> </span>52 1. p52 if >1, history point dump freq.(steps)</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 54 1.000000E+00 p54 direction of flow rate -- 1-X, 2-Y, 3-Z</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 55 3.000E+00 p55 volumetric flow rate</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 65 1.00000 p65</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 66 0.00000 p66 write format</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 67 0.00000 p67 read format</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 68 1. p68 iostep for avg_all</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 93 20.000 p93 if>0, # res-proj</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 94 5.00000 p94 if>0, start proj vel and p.s.</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 95 5.00000 p95 if>0, start proj pressure</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 99 4.000000E+00 p99 dealiasing</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 102 1.00000 p102 dump out divergence at each time step</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFTRAN = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFFLOW = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFHEAT = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFSPLIT = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFLOMACH = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFUSERVP = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFUSERMV = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFSTRS = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFCHAR = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFCYCLIC = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFAXIS = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFMVBD = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFMELT = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFMODEL = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFKEPS = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFMOAB = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFNEKNEK = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFSYNC = T</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFVCOR = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFINTQ = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFCWUZ = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFSWALL = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFGEOM = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFSURT = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFWCNO = F</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFTMSH for field 1 = F</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFADVC for field 1 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFNONL for field 1 = F</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFTMSH for field 2 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFADVC for field 2 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFNONL for field 2 = F</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFTMSH for field 3 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFADVC for field 3 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFNONL for field 3 = F</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFTMSH for field 4 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFADVC for field 4 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFNONL for field 4 = F</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFTMSH for field 5 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFADVC for field 5 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFNONL for field 5 = F</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFTMSH for field 6 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFADVC for field 6 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFNONL for field 6 = F</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFTMSH for field 7 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFADVC for field 7 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFNONL for field 7 = F</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFTMSH for field 8 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFADVC for field 8 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFNONL for field 8 = F</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFTMSH for field 9 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFADVC for field 9 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFNONL for field 9 = F</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">IFTMSH for field 10 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFADVC for field 10 = T</p>
<p class="MsoNormal" style="margin-right:35.65pt">IFNONL for field 10 = F</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">Dealiasing enabled, lxd= 6</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">______________________________________________________________________</p>
<p class="MsoNormal" style="margin-right:35.65pt">*.rea file</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> ***** OUTPUT FIELD SPECIFICATION *****</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 5 SPECIFICATIONS FOLLOW</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T COORDINATES</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T VELOCITY</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T PRESSURE</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T TEMPERATURE</p>
<p class="MsoNormal" style="margin-right:35.65pt"> F TEMPERATURE GRADIENT</p>
<p class="MsoNormal" style="margin-right:35.65pt"> 8 PASSIVE SCALARS</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T PS 1</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T PS 2</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T PS 3</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T PS 4</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T PS 5</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T PS 6</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T PS 7</p>
<p class="MsoNormal" style="margin-right:35.65pt"> T PS 8</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt">______________________________________________________________________</p>
<p class="MsoNormal" style="margin-right:35.65pt">*.usr file</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> subroutine uservp (ix,iy,iz,ieg)</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'SIZE'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'TOTAL'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'NEKUSE'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> if (ifield.eq.1) then</p>
<p class="MsoNormal" style="margin-right:35.65pt"> udiff = param(2) ! fluid kinematic viscosity</p>
<p class="MsoNormal" style="margin-right:35.65pt"> utrans = param(1) ! fluid density</p>
<p class="MsoNormal" style="margin-right:35.65pt"> if (ieg .gt. nelgv) then ! properties in the solid</p>
<p class="MsoNormal" style="margin-right:35.65pt"> utrans = param(1) ! solid density</p>
<p class="MsoNormal" style="margin-right:35.65pt"> endif</p>
<p class="MsoNormal" style="margin-right:35.65pt"> else</p>
<p class="MsoNormal" style="margin-right:35.65pt"> utrans = param(7) ! fluid rhocp</p>
<p class="MsoNormal" style="margin-right:35.65pt"> udiff = param(8) ! fluid conductivity</p>
<p class="MsoNormal" style="margin-right:35.65pt"> if (ieg .gt. nelgv) then ! properties in the solid</p>
<p class="MsoNormal" style="margin-right:35.65pt"> udiff = param(8) ! solid conductivity</p>
<p class="MsoNormal" style="margin-right:35.65pt"> utrans = param(7) ! solid rhocp</p>
<p class="MsoNormal" style="margin-right:35.65pt"> endif</p>
<p class="MsoNormal" style="margin-right:35.65pt"> endif</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> return</p>
<p class="MsoNormal" style="margin-right:35.65pt"> end</p>
<p class="MsoNormal" style="margin-right:35.65pt">c------------------------------------------</p>
<p class="MsoNormal" style="margin-right:35.65pt"> subroutine userq (ix,iy,iz,ieg)</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'SIZE'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'TOTAL'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'NEKUSE'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> qvol = 1.</p>
<p class="MsoNormal" style="margin-right:35.65pt"> if (ieg.gt.nelgv) qvol = -1.</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> return</p>
<p class="MsoNormal" style="margin-right:35.65pt"> end</p>
<p class="MsoNormal" style="margin-right:35.65pt">c------------------------------------------</p>
<p class="MsoNormal" style="margin-right:35.65pt"> subroutine userbc (ix,iy,iz,iside,ieg)</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'SIZE'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'TOTAL'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'NEKUSE'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> ux=0.0</p>
<p class="MsoNormal" style="margin-right:35.65pt"> uy=0.0</p>
<p class="MsoNormal" style="margin-right:35.65pt"> uz=0.0</p>
<p class="MsoNormal" style="margin-right:35.65pt"> temp=0.0</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> return</p>
<p class="MsoNormal" style="margin-right:35.65pt"> end</p>
<p class="MsoNormal" style="margin-right:35.65pt">c------------------------------------------</p>
<p class="MsoNormal" style="margin-right:35.65pt"> subroutine useric (ix,iy,iz,ieg)</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'SIZE'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'TOTAL'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> include 'NEKUSE'</p>
<p class="MsoNormal" style="margin-right:35.65pt"> </p>
<p class="MsoNormal" style="margin-right:35.65pt"> one = 1.</p>
<p class="MsoNormal" style="margin-right:35.65pt"> pi = 4.*atan(one)</p>
<p class="MsoNormal" style="margin-right:35.65pt"> ux = (9.*(1.-y**8))/8.-0.6283*Cos(x)*Sin(Pi*y)*Sin(z)</p>
<p class="MsoNormal" style="margin-right:35.65pt"> uy = 0.1*(1. + Cos(Pi*y))*Sin(x)*Sin(z)</p>
<p class="MsoNormal" style="margin-right:35.65pt"> <span lang="IT">uz = 0.3141*Cos(z)*Sin(x)*Sin(Pi*y)</span></p>
<p class="MsoNormal" style="margin-right:35.65pt"><span lang="IT"> </span>temp = 0.d0</p>
<p class="MsoNormal" style="margin-right:35.65pt"> return</p>
<p class="MsoNormal" style="margin-right:35.65pt"> end</p>
<p class="MsoNormal"><span style="font-size:11.0pt; font-family:"Calibri","sans-serif""> </span></p>
</div>
<br>
<br>
</div>
</div>
</div>
</div>
</body>
</html>