<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Hello,</p>
<p><br>
</p>
<p>In the userchk subroutine in the .usr file for the turbChannel example, what line of code can I add so rms values can be written in a .dat file such as it is done for umean and reystresses values.</p>
<p><br>
</p>
<p><span id="ms-rterangepaste-start"></span></p>
<div> ! write statistics to file<br>
iostep_avg = param(68)<br>
if (nid.eq.0 .and. istep.gt.0 .and. <br>
& mod(istep,iostep_avg).eq.0) then<br>
write(6,*) 'Dumping statistics ...'<br>
open(unit=56,file='reystresses.dat')<br>
write(56,'(A,1pe14.7)') '#time = ', time<br>
write(56,'(A)')<br>
& '# y y+ R_uu R_vv R_ww R_uv'<br>
open(unit=57,file='means.dat')<br>
write(57,'(A,1pe14.7)') '#time = ', time<br>
write(57,'(A)')<br>
& '# y y+ Umean'<br>
<br>
m = ny1*nely/2<br>
do i=1,m<br>
write(56,3) yy(i)+1<br>
& ,(yy(i)+1)*Re_tau <br>
& , (urms_pl(i)-(uavg_pl(i))**2)/u_tau**2<br>
& , vrms_pl(i)/u_tau**2<br>
& , wrms_pl(i)/u_tau**2<br>
& , uvms_pl(i)/u_tau**2<br>
write(57,3) yy(i) + 1.<br>
& , (yy(i)+1.)*Re_tau <br>
& , uavg_pl(i)/u_tau<br>
3 format(1p15e17.9)<br>
enddo<br>
<br>
close(56)<br>
close(57)<br>
<br>
timel = time<br>
<br>
endif<br>
<br>
return<br>
end</div>
<span id="ms-rterangepaste-end"></span><br>
<p></p>
<p><br>
</p>
<div id="Signature">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div>Thank you,</div>
<div><br>
</div>
<div>Virgilio A. Gomez</div>
<div>
<div><br>
</div>
</div>
<div><br>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>