[Nek5000-users] Uservp and drag_calc routine

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

Shriram,

I was getting drag values printed for my case... you
need to also set the obj definitions.  I did this through
a routine called set_obj, which was in the same .usr file
sent earlier ?

Paul


On Sat, 14 Aug 2010, nek5000-users at lists.mcs.anl.gov wrote:

> Paul,
>
> Yes. I am using a eddy viscosity model, similar to the one in turbchannel.
>
> I just tried using the torq_calc() but couldn't get it to print the drag
> values to logfile.  This is how I tried.
>
>      logical if_drag_out,if_torq_out
>      real x0(ldim)
>      data x0 /ldim*0/
>
>      scale = 1.0
>      if_drag_out = .true.
>      if_torq_out = .false.
>      call torque_calc(scale,x0,if_drag_out,if_torq_out)
>
> The full_restart.usr file served as a reference for me. Am I missing
> something ?
>
>
> Regards
> Shriram
>