[Nek5000-users] Forcing nek to read useric routine with PRESOLVE/RESTART OPTIONS set to 1 ??

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

Hi Can,

I would read the blah.f00001 file as you are currently doing, then,
in userchk() add something like:


       n = nx1*ny1*nz1*nelv
       if (istep.eq.0) then
          do i=1,n
             x = xm1(i,1,1,1)
             y = ym1(i,1,1,1)
             z = zm1(i,1,1,1)
             vx(i,1,1,1) = function ux of ( x , y, z )
             vy(i,1,1,1) = function uy of ( x , y, z )
             vz(i,1,1,1) = function uz of ( x , y, z )
          enddo
       endif

for whatever vx,vy,vz functions you want.

Paul


On Fri, 20 Apr 2012, nek5000-users at lists.mcs.anl.gov wrote:

> Hi neks,
> I would like to load a field (blah.f00001) and (before the first time step 
> occur) set the velocity of some region of the domain to 0.
> (I want to introduce a columnar vortex in a turbulent field)
> For this i  need nek to read the useric subroutine after it loads my 
> blah.f00001 but by default if PRESOLVE/RESTART OPTIONS is set to 1 it doesn't 
> read the useric subroutine.
> Is there a way to force it to do this ??
>
> B.regards
> Can
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>