<div dir="ltr">Hi Neks,<br><div><br></div><div>I'm running some differentially heated cavity simulations and want to output the average Nusselt over a wall on each output timestep.</div><div><br></div><div>Right now I'm calculating the temperature gradient like this:<br></div><div><br></div>common /grad_t/ dTx(lt), dTy(lt), dTz(lt)<br><div></div><div>parameter (lt=lx1*ly1*lz1*lelv)</div><div></div>if(mod(istep,iostep).eq.0) then<br><div style="margin-left:40px">call gradm1(dTx,dTy,dTz,t)<br>call dsavg(dTx)<br>call dsavg(dTy)<br>call dsavg(dTz)<br>call outpost(dTx,dTy,dTz,pr,t,'grd')<br></div><div>endif</div><div><br></div><div>I would appreciate any suggestion to implement a way to calculate and dump the average Nusselt (-dT/dy) between dumps at each iostep time. Part of the motivation is save disk space by avoiding the  call outpost of each field and then post processing the data.<b><br></b></div><div><br></div><div><br></div><div><div><div class="gmail_signature"><div>Regards,</div><div dir="ltr"><b><br></b></div><div dir="ltr"><b>Nicolás</b></div></div></div>
</div></div>