[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