[Nek5000-users] Coriolis Force

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

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/3312e3ce/attachment.html>