<html><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8" /></head><body style='font-size: 10pt; font-family: Verdana,Geneva,sans-serif'>
<p>Hello everyone,</p>
<p> My aim is to obtain a circular mesh with radius equal to 50.0 and centered at (0,0). I used the following in my input.genbox file :</p>
<p>input.rea<br />2<br />2<br />0.0 0.0 xctr yctr<br />-100 -50 nel_r, nel_theta<br />0.0 50.0 1.0 r0 rn ratio<br />0.0 360.0 1.0 @0 @n ratio<br />v ,O Velocity BC's (here, only at r0 and rn)</p>
<p>By this specification, i think we need to provide BC's at the point (0,0) and at the circular ring of radius 50. I want to specify ux=uy=0 at origin and outflow bc at the periphery of the circle.</p>
<p>I get the following in the box.tmp file :</p>
<p>v ,O Velocity BC's (here, only at r0 and rn)</p>
<p>My userbc subroutine is as follows :</p>
<p>subroutine userbc (ix,iy,iz,iside,ieg)<br />include 'SIZE'<br />include 'TOTAL'<br />include 'NEKUSE'</p>
<p>ux=0.0<br />uy=0.0<br />uz=0.0<br />temp=0.</p>
<p>return<br />end</p>
<p>My code runs into an error (screenshot is attachted). Can somebody point out where is the issue.</p>
<p>Thanks.</p>
<p>Rohan </p>
</body></html>