[Nek5000-users] T and P derivatives

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

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
>