[Nek5000-users] T and P derivatives

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Fri Nov 26 08:48:30 CST 2010


PS ... slight bug in what I sent... Should read:

>         call mappr (pm1,pr,px,py) ! interpolate P(m2) --> P(m1)
>         call gradm1(px,py,pz,pm1)


On Fri, 26 Nov 2010, nek5000-users at lists.mcs.anl.gov wrote:

>
> Hi Francesco,
>
> Don't waste too much time looking - we're happy to answer!
>
> As Stefan commented, we don't store the gradients, we just
> compute them.  Typical usage would be
>
>
>      subroutine userchk
>      include 'SIZE'
>      include 'TOTAL'
>      common /mystuff/ uij(lx1*ly1*lz1*lelt,ldim,ldim)
>     $               , px(lx1*ly1,lz1*lelt)
>     $               , py(lx1*ly1,lz1*lelt)
>     $               , pz(lx1*ly1,lz1*lelt)
>
>      call gradm1(uij(1,1,1),uij(1,1,2),uij(1,1,3),vx)
>      call gradm1(uij(1,2,1),uij(1,2,2),uij(1,2,3),vy)
>      call gradm1(uij(1,3,1),uij(1,3,2),uij(1,3,3),vz)
>
>      if (lx2.ne.lx1) then ! use something different for mesh 2
>         call mappr (pm1,pr,px,py) ! interpolate P(m2) --> P(m1)
>         call gradm1(px,py,pz,pr)
>      else
>         call gradm1(px,py,pz,pr)
>      endif
>
>
> --- Paul
>
>
> On Thu, 25 Nov 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
>> Hi all,
>> 
>> Sorry to bug you but i have already lost enough time looking for the 
>> temperature and pressure gradients.
>> Does anybody knows where are they stored?
>> And the velocity time derivatives ?
>> 
>> thank you all
>> 
>> francesco
>> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>



More information about the Nek5000-users mailing list