[Nek5000-users] Output of temperature gradient, scalar dissipation rate and convective flux

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

Hi,

thanks for the following suggestion concerning our post processing where we need temperature gradients, scalar dissipation rate and convective heat fluxes:

     subroutine userchk
        :
     parameter (lt=lx1*ly1*lz1*lelt)
     common /mygrad/ tx(lt),ty(lt),tz(lt)
        :
     if (mod(istep,iostep).eq.0) then ! , say,
        call gradm1(tx,ty,tz,t)
        call outpost(tx,ty,tz,pr,t,'gdt') ! write to gdtblah.f000...
     endif
        :
     return
     end

It has however one shortcoming, namely the unnecessary additional output of the pressure field which we don't need 
for this analysis.

As far as I can see when going down from the routine outpost --> outpost2 --> prepost  this is a fixed
structure and the pressure field cannot be simply substituted by another scalar field (e.g. the scalar dissipation rate)
since it might have a different dimension (e.g. in case of spaces Pn-Pn-2 which we are using).

--- Can this problem be solved without editing the routines prepost?  
--- Is prepost_map inside prepost still necessary? 

Thanks, Joerg.