<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div><br>
</div>
Hi Rif,
<div><br>
</div>
<div>You should use "ux" and "uy" ... which are the point wise values for the velocity that</div>
<div>are passed in via NEKUSE.</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF17996" 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> Tuesday, December 29, 2015 11:28 PM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Coriolis Force<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">Hi Neks,
<div>I tried to include a rotation in my simulation for a flat plate by introducing a Coriolis term in the subroutine
<b>userf</b> like so :</div>
<div><br>
</div>
<div>---------------------------------------------------------------------------------------</div>
<div>subroutine userf (ix,iy,iz,eg)</div>
<div>include 'SIZE'</div>
<div>include 'TOTAL'</div>
<div>include 'NEKUSE'</div>
<div><br>
</div>
<div>integer e,f,eg</div>
<div>e = gllel(eg)</div>
<div><br>
</div>
<div>Omega= 2</div>
<div>ffx = +2.0*Omega*vy</div>
<div>
<div>ffy = -2.0*Omega*vx</div>
</div>
<div>ffz = 0.0</div>
<div><br>
</div>
<div>return</div>
<div>end</div>
<div>-----------------------------------------------------------------------------------------------</div>
<div>PROBLEM: When i tried compiling, it gives me an error</div>
<div>ERROR: Incompatible ranks 0 and 4 in assignment at (1)</div>
<div><br>
</div>
<div>Is there a way out of this? Does this have something to do with assigning a scalar to an array? if so how do I get past this?</div>
<div><br>
</div>
<div><br>
</div>
<div>Rif</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>