<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal style='margin-right:35.65pt'>Dear Neks, <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>Sorry for my previous mail.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>I hope that the layout will be good now.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>_________________________________________________________________<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>Dear users and developers,<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>I’m simulating a 3D conjugate heat transfer case, considering temperature as a passive scalar.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>When I simulate only one temperature field everything is ok.<o:p></o:p></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.<o:p></o:p></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):<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'> ldimt T PS1 PS2 PS3 PS4 PS5 PS6 PS7 PS8<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 1 ok<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 2 ok --<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 3 ok -- ok<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 4 ok -- ok ok<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 5 ok -- -- -- ok<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 6 ok -- -- -- ok ok<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 7 ok -- -- -- ok ok ok<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 8 ok -- -- -- ok ok ok ok<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 9 ok -- -- -- ok ok ok ok ok<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>I’m using PN-PN.<o:p></o:p></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.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></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.<o:p></o:p></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.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>::<o:p></o:p></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<o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT> </span>Solving for heat<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 250 Hmholtz TEMP: 17 9.2949E-04 4.5939E-03 1.0000E-03<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT> </span>250 3.1937E+00 2.0457E-01 Heat done<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> Solving for fluid<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>::<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>call outfld: ifpsco: T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 250 3.1937E+00 Write checkpoint:<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 0 250 OPEN: cjh_tc0.f00001<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 250 3.1937E+00 done :: Write checkpoint<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> file size = 5.6 MB<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> avg data-throughput = 367.7MB/s<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <span lang=IT>io-nodes = 1<o:p></o:p></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<o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT> </span>Solving for heat<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 251 Hmholtz TEMP: 17 8.3175E-04 4.6145E-03 1.0000E-03<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></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<o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT> </span>251 3.2046E+00 3.6832E-01 Heat done<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>::<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></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.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>Does anyone know how to solve this problem?<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>Thanks in advance.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></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.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>Best Regards,<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>Donato<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>______________________________________________________________________<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>SIZE file<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>parameter (ldimt= 9) ! passive scalars + T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>______________________________________________________________________<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>*.box file<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>base.rea<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>3 spatial dimension ( < 0 --> generate .rea/.re2 pair)<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>10 number of fields<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>#======<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>Box<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>-10 14 -5 nelx,nely,nelz for Box<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>0. 6. 1. x0,x1,gain<o:p></o:p></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<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>0.50 0.70 0.86 0.94 0.98 1.00 y0,y1,...,yn<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>0. 3. 1. z0,z1,gain<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>P ,P ,W ,W ,P ,P vel BCs <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,E ,P ,P temp BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,E ,P ,P PS1 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,E ,P ,P PS2 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,E ,P ,P PS3 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,E ,P ,P PS4 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,E ,P ,P PS5 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,E ,P ,P PS6 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,E ,P ,P PS7 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,E ,P ,P PS8 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT><o:p> </o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'>Box_up<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>-10 7 -5 nelx,nely,nelz for Box_up<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>0. 6. 1. x0,x1,gain<o:p></o:p></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<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>0. 3. 1. z0,z1,gain<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> , , , , , vel BCs <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>P ,P ,I ,E ,P ,P temp BCs <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>P ,P ,I ,E ,P ,P PS1 BCs <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,I ,E ,P ,P PS2 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,I ,E ,P ,P PS3 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,I ,E ,P ,P PS4 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,I ,E ,P ,P PS5 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,I ,E ,P ,P PS6 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,I ,E ,P ,P PS7 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,I ,E ,P ,P PS8 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT><o:p> </o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'>Box_down<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>-10 7 -5 nelx,nely,nelz for Box_down<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>0. 6. 1. x0,x1,gain<o:p></o:p></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<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>0. 3. 1. z0,z1,gain<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> , , , , , vel BCs <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>P ,P ,E ,I ,P ,P temp BCs <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>P ,P ,E ,I ,P ,P PS1 BCs <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,I ,P ,P PS2 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,I ,P ,P PS3 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,I ,P ,P PS4 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,I ,P ,P PS5 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,I ,P ,P PS6 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,I ,P ,P PS7 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT>P ,P ,E ,I ,P ,P PS8 BCs<o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT><o:p> </o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT><o:p> </o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'>______________________________________________________________________<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>logfile<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>103 Parameters<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 1 1.00000 p1 DENSITY<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 2 -2812.94 p2 Re_b<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 7 1.00000 p7 RHOCP<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 8 1.00000 p8 CONDUCT<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 11 500. p11 NSTEPS<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 12 0.05 p12 DT<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 15 250. p15 IOSTEP<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 21 1.000000E-06 p21 DIVERGENCE<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 22 1.000000E-03 p22 HELMHOLTZ<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 23 8 p23 NPSCAL<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <span lang=IT>24 0.100000E-01 p24 TOLREL<o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT> 25 0.100000E-01 p25 TOLABS<o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT> </span>26 0.50000 p26 TARGET COURANT NUMBER<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 27 3.00000 p27 TORDER<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <span lang=IT>28 0.000000E+00 p28<o:p></o:p></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.<o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT> </span>52 1. p52 if >1, history point dump freq.(steps)<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 54 1.000000E+00 p54 direction of flow rate -- 1-X, 2-Y, 3-Z<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 55 3.000E+00 p55 volumetric flow rate<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 65 1.00000 p65<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 66 0.00000 p66 write format<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 67 0.00000 p67 read format<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 68 1. p68 iostep for avg_all<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 93 20.000 p93 if>0, # res-proj<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 94 5.00000 p94 if>0, start proj vel and p.s.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 95 5.00000 p95 if>0, start proj pressure<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 99 4.000000E+00 p99 dealiasing<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 102 1.00000 p102 dump out divergence at each time step<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFTRAN = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFFLOW = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFHEAT = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFSPLIT = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFLOMACH = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFUSERVP = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFUSERMV = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFSTRS = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFCHAR = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFCYCLIC = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFAXIS = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFMVBD = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFMELT = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFMODEL = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFKEPS = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFMOAB = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFNEKNEK = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFSYNC = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFVCOR = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFINTQ = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFCWUZ = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFSWALL = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFGEOM = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFSURT = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFWCNO = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFTMSH for field 1 = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFADVC for field 1 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFNONL for field 1 = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFTMSH for field 2 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFADVC for field 2 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFNONL for field 2 = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFTMSH for field 3 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFADVC for field 3 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFNONL for field 3 = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFTMSH for field 4 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFADVC for field 4 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFNONL for field 4 = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFTMSH for field 5 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFADVC for field 5 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFNONL for field 5 = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFTMSH for field 6 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFADVC for field 6 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFNONL for field 6 = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFTMSH for field 7 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFADVC for field 7 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFNONL for field 7 = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFTMSH for field 8 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFADVC for field 8 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFNONL for field 8 = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFTMSH for field 9 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFADVC for field 9 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFNONL for field 9 = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFTMSH for field 10 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFADVC for field 10 = T<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>IFNONL for field 10 = F<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>Dealiasing enabled, lxd= 6<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>______________________________________________________________________<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>*.rea file<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'> ***** OUTPUT FIELD SPECIFICATION *****<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 5 SPECIFICATIONS FOLLOW<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T COORDINATES<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T VELOCITY<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T PRESSURE<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T TEMPERATURE<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> F TEMPERATURE GRADIENT<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> 8 PASSIVE SCALARS<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T PS 1<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T PS 2<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T PS 3<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T PS 4<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T PS 5<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T PS 6<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T PS 7<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> T PS 8<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'>______________________________________________________________________<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>*.usr file<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'> subroutine uservp (ix,iy,iz,ieg)<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'SIZE'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'TOTAL'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'NEKUSE'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'> if (ifield.eq.1) then<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> udiff = param(2) ! fluid kinematic viscosity<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> utrans = param(1) ! fluid density<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> if (ieg .gt. nelgv) then ! properties in the solid<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> utrans = param(1) ! solid density<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> endif<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> else<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> utrans = param(7) ! fluid rhocp<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> udiff = param(8) ! fluid conductivity<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> if (ieg .gt. nelgv) then ! properties in the solid<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> udiff = param(8) ! solid conductivity<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> utrans = param(7) ! solid rhocp<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> endif<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> endif<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'> return<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> end<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>c------------------------------------------<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> subroutine userq (ix,iy,iz,ieg)<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'SIZE'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'TOTAL'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'NEKUSE'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'> qvol = 1.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> if (ieg.gt.nelgv) qvol = -1.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'> return<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> end<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>c------------------------------------------<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> subroutine userbc (ix,iy,iz,iside,ieg)<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'SIZE'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'TOTAL'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'NEKUSE'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'> ux=0.0<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> uy=0.0<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> uz=0.0<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> temp=0.0<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'> return<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> end<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'>c------------------------------------------<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> subroutine useric (ix,iy,iz,ieg)<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'SIZE'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'TOTAL'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> include 'NEKUSE'<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'><o:p> </o:p></p><p class=MsoNormal style='margin-right:35.65pt'> one = 1.<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> pi = 4.*atan(one)<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> ux = (9.*(1.-y**8))/8.-0.6283*Cos(x)*Sin(Pi*y)*Sin(z)<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> uy = 0.1*(1. + Cos(Pi*y))*Sin(x)*Sin(z)<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> <span lang=IT>uz = 0.3141*Cos(z)*Sin(x)*Sin(Pi*y)<o:p></o:p></span></p><p class=MsoNormal style='margin-right:35.65pt'><span lang=IT> </span>temp = 0.d0<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> return<o:p></o:p></p><p class=MsoNormal style='margin-right:35.65pt'> end<o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'><o:p> </o:p></span></p></div>
<br><br></body></html>