<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"Malgun Gothic";
panose-1:2 11 5 3 2 0 0 2 0 4;}
@font-face
{font-family:"\@Malgun Gothic";}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="blue" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hi Phillip, </p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am using the “Blasius” example test outputting the streamfunction as a scalar.
</p>
<p class="MsoNormal">I included my code below. I am not sure what I am doing wrong.</p>
<p class="MsoNormal">In userchk, I added “compute streamfunction from vorticity” from the psi_omega example, but it does not compute the psi properly… (max_psi = min_psi = 0 in every step)</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">c-----------------------------------------------------------------------</p>
<p class="MsoNormal"> subroutine userchk</p>
<p class="MsoNormal"> include 'SIZE'</p>
<p class="MsoNormal"> include 'TOTAL'</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> common /myblas/ ub(lx1,ly1,lz1,lelt),vb(lx1,ly1,lz1,lelt)</p>
<p class="MsoNormal"> common /blasiusr/ u(lx1,ly1,lz1,lelv),v(lx1,ly1,lz1,lelv)</p>
<p class="MsoNormal"> common /blasiusg/ bin(lx1,ly1,lz1,lelv)</p>
<p class="MsoNormal"> common /scrns/ dum(lx1*ly1*lz1*lelv)</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> common /xtream/ psi(lx1*ly1*lz1*lelt)</p>
<p class="MsoNormal"> $ , rhs(lx1*ly1*lz1*lelt)</p>
<p class="MsoNormal"> $ , h1 (lx1*ly1*lz1*lelt)</p>
<p class="MsoNormal"> $ , h2 (lx1*ly1*lz1*lelt)</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> integer bin</p>
<p class="MsoNormal"> real normu(50),normv(50)</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> real xmax,xmin,dx</p>
<p class="MsoNormal"> save xmax,xmin,dx</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> n = nx1*ny1*nz1*nelv</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">c. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .</p>
<p class="MsoNormal">c Compute streamfunction from vorticity</p>
<p class="MsoNormal"> call col2(rhs,bm1,n)</p>
<p class="MsoNormal"> call rone (h1 ,n)</p>
<p class="MsoNormal"> call rzero(h2 ,n)</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> tol = param(22)</p>
<p class="MsoNormal"> imsh = 1</p>
<p class="MsoNormal"> call hmholtz('psi ',psi,rhs,h1,h2,v1mask,vmult,imsh,tol,200,1)</p>
<p class="MsoNormal">c. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .</p>
<p class="MsoNormal"> if (istep.eq.0) call outpost(vx,vy,vz,pr,psi,' ')</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">c Compute delta*</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> do i=1,n</p>
<p class="MsoNormal"> dum(i) = 1.-vx(i,1,1,1)</p>
<p class="MsoNormal"> enddo</p>
<p class="MsoNormal"> tot_deficit = glsc2(dum,bm1,n)</p>
<p class="MsoNormal"> xmin = glmin(xm1,n)</p>
<p class="MsoNormal"> xmax = glmax(xm1,n)</p>
<p class="MsoNormal"> delta_star = tot_deficit/(xmax-xmin)</p>
<p class="MsoNormal"> psimin = glmin(psi,n)</p>
<p class="MsoNormal"> psimax = glmax(psi,n)</p>
<p class="MsoNormal"> if (nid.eq.0) write(6,1) istep,time,delta_star,psimin,psimax</p>
<p class="MsoNormal"> 1 format(i8,1p4e13.5,' blasius: delta*')</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> return</p>
<p class="MsoNormal"> end</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>