[Nek5000-users] T and P derivatives

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

thank you,

i was just hoping to save same operations.

francesco


On 11/26/2010 03:48 PM, nek5000-users at lists.mcs.anl.gov wrote:
>
> 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
>>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users