<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: Times New Roman; font-size: 12pt; color: #000000'><P>Hi Paul and Aleks,</P>
<P> </P>
<P>I have a question regarding this. I have implemented this and it is working nicely, however it is dumping </P>
<P>out a .f file every step, rather than every iotime step I have specified in the rea.</P>
<P> </P>
<P>Do you know how to fix this? </P>
<P> </P>
<P>- Michael</P>
<P><BR>----- Original Message -----<BR>From: nek5000-users@lists.mcs.anl.gov<BR>To: nek5000-users@lists.mcs.anl.gov<BR>Sent: Wednesday, June 2, 2010 10:57:41 AM GMT -06:00 US/Canada Central<BR>Subject: Re: [Nek5000-users] Compute the derivatives<BR><BR>Hi JC,<BR><BR>To dump the computed derivatives into separate field files you can use in <BR>userchk<BR><BR> ifxyo = .true. ! for cooordinates<BR> ifpo = .false. ! turn off pressure output<BR> ifto = .false. ! turn off pressure output<BR> call outpost(vxx,vxy,vxz,pr,t,'dvx')<BR> call outpost(vyx,vyy,vyz,pr,t,'dvy')<BR>...<BR><BR>Best,<BR>Aleks<BR><BR><BR><BR>On Wed, 2 Jun 2010, nek5000-users@lists.mcs.anl.gov wrote:<BR><BR>><BR>> dudxyz is somewhat dated and slow...<BR>><BR>> I recommend somthing like:<BR>><BR>> parameter (lt=lx1*ly1*lz1*lelv)<BR>> common /mygrad/ vxx(lt),vxy(lt),vxz(lt)<BR>> $ , vyx(lt),vyy(lt),vyz(lt)<BR>> $ , vzx(lt),vzy(lt),vzz(lt)<BR>><BR>> call gradm1(vxx,vxy,vxz,vx)<BR>> call gradm1(vyx,vyy,vyz,vy)<BR>> call gradm1(vzx,vzy,vzz,vz)<BR>><BR>> Note that putting these arrays into a common block, as<BR>> shown, is a good idea for several reasons, one of which<BR>> is that it guarantees that the data will be byte-aligned with the cache line. <BR>> On some machines (e.g., BG/P) this<BR>> is crucial for proper functionality - for this reason we<BR>> don't mix variable types in a given common block.<BR>><BR>> If you subsequently wish to access elements of vxx etc.<BR>> on an element-by-element basis you can of course do something<BR>> like<BR>><BR>> parameter (lt=lx1*ly1*lz1)<BR>> common /mygrad/ vxx(lt,lelt),vxy(lt,lelt),vxz(lt,lelt)<BR>> :<BR>> :<BR>><BR>> Paul<BR>><BR>> On Wed, 2 Jun 2010, nek5000-users@lists.mcs.anl.gov wrote:<BR>><BR>>> Hi Nek's,<BR>>> <BR>>> I have a home-made global stability code, but in order to use it I first <BR>>> need to compute the derivatives of my base flow: dUdX, dUdY, dVdX, dVdY.<BR>>> How could one compute these derivatives using Nek 5k and drop them along <BR>>> with X Y U V in the blah.f**** file (or perhaps another .f**** file) ? I <BR>>> assume I have to make use of the subroutine *dudxyz <BR>>> (du,u,rm1,sm1,tm1,jm1,imsh,isd)*, but I'm not sure what to put as inputs.<BR>>> <BR>>> Regards,<BR>>> JC<BR>>> <BR>> _______________________________________________<BR>> Nek5000-users mailing list<BR>> Nek5000-users@lists.mcs.anl.gov<BR>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<BR>><BR>_______________________________________________<BR>Nek5000-users mailing list<BR>Nek5000-users@lists.mcs.anl.gov<BR>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<BR></P></div></body></html>