[Nek5000-users] Coriolis Force

[nek5000-users at lists.mcs.anl.gov]

Hi Rif,

You should use "ux" and "uy" ... which are the point wise values for the velocity that
are passed in via NEKUSE.

Paul

________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
Sent: Tuesday, December 29, 2015 11:28 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Coriolis Force

Hi Neks,
I tried to include a rotation in my simulation for a flat plate by introducing a Coriolis term in the subroutine userf  like so :

---------------------------------------------------------------------------------------
subroutine userf (ix,iy,iz,eg)
include 'SIZE'
include 'TOTAL'
include 'NEKUSE'

integer e,f,eg
e = gllel(eg)

Omega= 2
ffx = +2.0*Omega*vy
ffy = -2.0*Omega*vx
ffz = 0.0

return
end
-----------------------------------------------------------------------------------------------
PROBLEM: When i tried compiling, it gives me an error
ERROR: Incompatible ranks 0 and 4 in assignment at (1)

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?


Rif
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20151230/1d1bb33d/attachment.html>