[Nek5000-users] Setics subroutine in load_fld

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Tue Jan 31 11:49:29 CST 2017


Hi Navrose,

I do not know if there is any better method, because I would do it probably
the same way. Another way-- you can probably define another set of routine
say load_fld_pert(), with modified setics and RESTART file, but probably it
is as good as using opcopy().

On Tue, Jan 31, 2017 at 12:25 AM, <nek5000-users at lists.mcs.anl.gov> wrote:

> Many thanks Tanmoy.
>
> So each time load_fld is called it will load the base
> flow variables (I am carrying out computations for incompressible N-S
> equations). What if I want to load the perturbations also from a different
> file? One way that I see it is to use load_fld for the restart file for
> perturbations, copy it in the perturbation variables (vxp, vyp, vzp, prp),
> and then again use load_fld for the base flow file. Is there a direct or
> better method for this?
>
> Thank you.
>
> Best wishes,
> Navrose.
> On Monday 30 January 2017 08:15 PM, nek5000-users at lists.mcs.anl.gov wrote:
>
> Hi Navrose,
>
> You are right. Two things are important in the load_fld() subroutine. The
> include file 'RESTART' and the routine setics. So everytime a load_fld
> subroutine is being called, the variable values read from the current file
> would be stored in vx(i,j,k,e), vy(i,j,k,e), vz(i,j,k,e), pr, t and so on...
>
> Best Regards,
> Tanmoy
>
> On Mon, Jan 30, 2017 at 6:04 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
>
>> Dear Neks,
>>
>>   I am relatively new to NEK5000. I am going through some
>> the .usr files that have been used previously in my research
>> group. I notice that in the *userchk* subroutine, *load_fld *has been
>> called a few times. Upon checking the *load_fld* function in
>> in the source code, I found that *setics *subroutine is called at the
>> end of *load_fld* function.
>>
>> Does this mean that each time *load_fld *function is called, the *useric*
>> subroutine in the .usr file is executed? I would like to mention that
>> no restart file is given in the .rea file.
>>
>> Many thanks.
>>
>> Best regards,
>>
>> Navrose.
>>
>>
>>
>>
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>
>>
>
>
> _______________________________________________
> Nek5000-users mailing listNek5000-users at lists.mcs.anl.govhttps://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170131/b1d09f8e/attachment-0001.html>


More information about the Nek5000-users mailing list