[Nek5000-users] Modifying an initial field

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

Hi David,

The best place would probably be in userchk.

You could have a block like

if (istep.eq.0) then

C........do stuff

endif

If you are doing something else in userchk that would depend on your
modification of the field quantities, then you would want  this to come
first.

Just in case, the velocity fields are accessed like

vx(i,j,k,e) ( for u)

and likewise for vy and vz, where i,j,k go from 1 to nx1 (ny1, nz1) and e
goes from 1 to nelv.  Sorry if that was more information than you needed,
just wanted to make sure all my bases were covered :-).

Does this help?

Josh Camp
"All that is necessary for the triumph of evil is that good men do nothing"
-- Edmund Burke
------------------------------
From: nek5000-users at lists.mcs.anl.gov
Sent: 10/4/2011 6:08 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Modifying an initial field

Hi neks,

I am loading an initial condition from a .fld file, but I would like to
modify the velocity field before using it as my initial condition.  (I would
like to superpose a shear onto the existing field.)  Can someone tell me
what part of the .usr file this should be done in?

Best,

David
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20111004/506c308c/attachment.html>