<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:11.0pt;
font-family:"Calibri","sans-serif";}
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>Dear users and developers,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I’m simulating a 3D conjugate heat transfer case, considering temperature as a passive scalar.<o:p></o:p></p><p class=MsoNormal>When I simulate only one temperature field everything is ok. <o:p></o:p></p><p class=MsoNormal>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>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='text-indent:36.0pt'>ldimt T PS1 PS2 PS3 PS4 PS5 PS6 PS7 PS8<o:p></o:p></p><p class=MsoNormal> 1 ok<o:p></o:p></p><p class=MsoNormal> 2 ok --<o:p></o:p></p><p class=MsoNormal> 3 ok -- ok<o:p></o:p></p><p class=MsoNormal> 4 ok -- ok ok<o:p></o:p></p><p class=MsoNormal> 5 ok -- -- -- ok<o:p></o:p></p><p class=MsoNormal> 6 ok -- -- -- ok ok<o:p></o:p></p><p class=MsoNormal> 7 ok -- -- -- ok ok ok<o:p></o:p></p><p class=MsoNormal> 8 ok -- -- -- ok ok ok ok<o:p></o:p></p><p class=MsoNormal> 9 ok -- -- -- ok ok ok ok ok<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I’m using PN-PN.<o:p></o:p></p><p class=MsoNormal>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><o:p> </o:p></p><p class=MsoNormal>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>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><o:p> </o:p></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>::<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>Step 250, t= 3.1936627E+00, DT= 1.0935650E-02, C= 0.571 7.8515E+01 2.6463E-01<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> </span><span style='font-size:8.0pt;font-family:"Courier New"'>Solving for heat<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 250 Hmholtz TEMP: 17 9.2949E-04 4.5939E-03 1.0000E-03<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> </span><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> </span><span style='font-size:8.0pt;font-family:"Courier New"'>250 3.1937E+00 2.0457E-01 Heat done<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> Solving for fluid<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>::</span><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>call outfld: ifpsco: T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 250 3.1937E+00 Write checkpoint:<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 0 250 OPEN: cjh_tc0.f00001<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 250 3.1937E+00 done :: Write checkpoint<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> file size = 5.6 MB<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> avg data-throughput = 367.7MB/s<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> io-nodes = 1<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>Step 251, t= 3.2045983E+00, DT= 1.0935650E-02, C= 0.566 7.8809E+01 2.9378E-01<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> </span><span style='font-size:8.0pt;font-family:"Courier New"'>Solving for heat<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 251 Hmholtz TEMP: 17 8.3175E-04 4.6145E-03 1.0000E-03<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> </span><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 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-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> </span><span style='font-size:8.0pt;font-family:"Courier New"'>251 3.2046E+00 3.6832E-01 Heat done<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>::<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Looking to the outfld subroutine I wasn’t able to figure out the reason for that.<o:p></o:p></p><p class=MsoNormal>Does anyone know how to solve this problem?<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Thanks in advance.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>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><o:p> </o:p></p><p class=MsoNormal>Best Regards,<o:p></o:p></p><p class=MsoNormal>Donato<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>________________________________________________________________________________<o:p></o:p></p><p class=MsoNormal>SIZE file<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><span style='font-size:8.0pt;font-family:"Courier New"'> parameter (ldimt= 9) ! passive scalars + T<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>________________________________________________________________________________<o:p></o:p></p><p class=MsoNormal>*.box file<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>base.rea<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>3 spatial dimension ( < 0 --> generate .rea/.re2 pair)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>10 number of fields<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>#======<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>Box<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>-10 14 -5 nelx,nely,nelz for Box<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>0. 6. 1. x0,x1,gain<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>-1.00 -0.98 -0.94 -0.86 -0.70 -0.50 -0.25 -0.00 0.25<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 0.50 0.70 0.86 0.94 0.98 1.00 y0,y1,...,yn<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>0. 3. 1. z0,z1,gain<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,W ,W ,P ,P vel BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,E ,P ,P temp BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,E ,P ,P PS1 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,E ,P ,P PS2 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,E ,P ,P PS3 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,E ,P ,P PS4 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,E ,P ,P PS5 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,E ,P ,P PS6 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,E ,P ,P PS7 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,E ,P ,P PS8 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>Box_up<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>-10 7 -5 nelx,nely,nelz for Box_up<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>0. 6. 1. x0,x1,gain<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>1.00 1.02 1.06 1.14 1.30 1.50 1.75 2.00 y0,y1,...,yn<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>0. 3. 1. z0,z1,gain<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> , , , , , vel BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,I ,E ,P ,P temp BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,I ,E ,P ,P PS1 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,I ,E ,P ,P PS2 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,I ,E ,P ,P PS3 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,I ,E ,P ,P PS4 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,I ,E ,P ,P PS5 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,I ,E ,P ,P PS6 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,I ,E ,P ,P PS7 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,I ,E ,P ,P PS8 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>Box_down<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>-10 7 -5 nelx,nely,nelz for Box_down<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>0. 6. 1. x0,x1,gain<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>-2.00 -1.75 -1.50 -1.30 -1.14 -1.06 -1.02 -1.00 y0,y1,...,yn<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>0. 3. 1. z0,z1,gain<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> , , , , , vel BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,I ,P ,P temp BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,I ,P ,P PS1 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,I ,P ,P PS2 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,I ,P ,P PS3 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,I ,P ,P PS4 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,I ,P ,P PS5 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,I ,P ,P PS6 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,I ,P ,P PS7 BCs <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>P ,P ,E ,I ,P ,P PS8 BCs<o:p></o:p></span></p><p class=MsoNormal><span lang=IT><o:p> </o:p></span></p><p class=MsoNormal><span lang=IT><o:p> </o:p></span></p><p class=MsoNormal>________________________________________________________________________________<o:p></o:p></p><p class=MsoNormal>logfile<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>103 Parameters<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 1 1.00000 p1 DENSITY<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 2 -2812.94 p2 Re_b<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 7 1.00000 p7 RHOCP<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 8 1.00000 p8 CONDUCT<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 11 500. p11 NSTEPS<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 12 0.05 p12 DT<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> </span><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>15 250. p15 IOSTEP<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 21 1.000000E-06 p21 DIVERGENCE<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 22 1.000000E-03 p22 HELMHOLTZ<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 23 8 p23 NPSCAL<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 24 0.100000E-01 p24 TOLREL<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 25 0.100000E-01 p25 TOLABS<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> </span><span style='font-size:8.0pt;font-family:"Courier New"'>26 0.50000 p26 TARGET COURANT NUMBER<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 27 3.00000 p27 TORDER<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> </span><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>28 0.000000E+00 p28<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> 30 1.000000E+00 p30 1: user-def props via USERVP subr.<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> </span><span style='font-size:8.0pt;font-family:"Courier New"'>52 1. p52 if >1, history point dump frequency (in steps)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 54 1.000000E+00 p54 direction of flow rate -- 1-X, 2-Y, 3-Z<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 55 3.000E+00 p55 volumetric flow rate<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 65 1.00000 p65<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 66 0.00000 p66 write format<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 67 0.00000 p67 read format<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 68 1. p68 iostep for avg_all (if ^= 0, else = iostep)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 93 20.000 p93 if>0, # res-proj (adjust MXPREF in SIZEu accordingly)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 94 5.00000 p94 if>0, start proj vel and p.s. after P094 steps<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 95 5.00000 p95 if>0, start proj pressure after p095 steps<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 99 4.000000E+00 p99 dealiasing<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 102 1.00000 p102 dump out divergence at each time step<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFTRAN = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFFLOW = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFHEAT = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFSPLIT = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFLOMACH = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFUSERVP = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFUSERMV = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFSTRS = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFCHAR = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFCYCLIC = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFAXIS = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFMVBD = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFMELT = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFMODEL = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFKEPS = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFMOAB = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFNEKNEK = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFSYNC = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFVCOR = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFINTQ = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFCWUZ = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFSWALL = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFGEOM = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFSURT = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFWCNO = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFTMSH for field 1 = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFADVC for field 1 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFNONL for field 1 = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFTMSH for field 2 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFADVC for field 2 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFNONL for field 2 = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFTMSH for field 3 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFADVC for field 3 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFNONL for field 3 = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFTMSH for field 4 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFADVC for field 4 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFNONL for field 4 = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFTMSH for field 5 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFADVC for field 5 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFNONL for field 5 = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFTMSH for field 6 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFADVC for field 6 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFNONL for field 6 = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFTMSH for field 7 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFADVC for field 7 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFNONL for field 7 = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFTMSH for field 8 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFADVC for field 8 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFNONL for field 8 = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFTMSH for field 9 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFADVC for field 9 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFNONL for field 9 = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFTMSH for field 10 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFADVC for field 10 = T<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> IFNONL for field 10 = F<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> <o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> Dealiasing enabled, lxd= 6<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>________________________________________________________________________________<o:p></o:p></p><p class=MsoNormal>*.rea file<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> ***** OUTPUT FIELD SPECIFICATION *****<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 5 SPECIFICATIONS FOLLOW<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T COORDINATES<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T VELOCITY<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T PRESSURE<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T TEMPERATURE<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> F TEMPERATURE GRADIENT<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> 8 PASSIVE SCALARS<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T PS 1<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T PS 2<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T PS 3<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T PS 4<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T PS 5<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T PS 6<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T PS 7<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> T PS 8<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>________________________________________________________________________________<o:p></o:p></p><p class=MsoNormal>*.usr file<o:p></o:p></p><p class=MsoNormal><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> subroutine uservp (ix,iy,iz,ieg)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'SIZE'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'TOTAL'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'NEKUSE'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> if (ifield.eq.1) then<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> udiff = param(2) ! fluid kinematic viscosity<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> utrans = param(1) ! fluid density<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> if (ieg .gt. nelgv) then ! properties in the solid<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> utrans = param(1) ! solid density<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> endif<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> else<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> utrans = param(7) ! fluid rhocp<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> udiff = param(8) ! fluid conductivity<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> if (ieg .gt. nelgv) then ! properties in the solid<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> udiff = param(8) ! solid conductivity<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> utrans = param(7) ! solid rhocp<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> endif<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> endif<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> return<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> end<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>c------------------------------------------<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> subroutine userq (ix,iy,iz,ieg)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'SIZE'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'TOTAL'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'NEKUSE'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> qvol = 1.<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> if (ieg.gt.nelgv) qvol = -1.<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> return<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> end<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>c------------------------------------------<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> subroutine userbc (ix,iy,iz,iside,ieg)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'SIZE'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'TOTAL'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'NEKUSE'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> ux=0.0<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> uy=0.0<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> uz=0.0<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> temp=0.0<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> return<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> end<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'>c------------------------------------------<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> subroutine useric (ix,iy,iz,ieg)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'SIZE'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'TOTAL'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> include 'NEKUSE'<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> one = 1.<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> pi = 4.*atan(one)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> ux = (9.*(1.-y**8))/8.-0.6283185307179586*Cos(x)*Sin(Pi*y)*Sin(z)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> uy = 0.1*(1. + Cos(Pi*y))*Sin(x)*Sin(z)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> </span><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'>uz = 0.3141592653589793*Cos(z)*Sin(x)*Sin(Pi*y)<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span lang=IT style='font-size:8.0pt;font-family:"Courier New"'> </span><span style='font-size:8.0pt;font-family:"Courier New"'>temp = 0.d0<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> return<o:p></o:p></span></p><p class=MsoNormal style='margin-left:36.0pt'><span style='font-size:8.0pt;font-family:"Courier New"'> end<o:p></o:p></span></p></div>
<br><br></body></html>