<html><body><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div>Hi Neks,<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div> I am bit stuck on this occasion for past 2 days, the code below shows  a part of the example called ocyl2.usr (from subroutine my_base_meshv(basev,h1,h2,rhs,msk,tmp) ). This code is to generate a base velocity (where basev=1 at wall of object and 0 at the boundary of the domain). I assume that I have a blending coefficient h1.(using Paul's method or volume of elements method). This subroutine working perfectly for the 2D cases. That is I have the base velocity ranges between 0 and 1. But when I tried it on an aero dynamic foil case(its a 3D case where just 1 element in Y direction), the value of base velocity does not even exceed 0. Almsot all values are like -0.0021564... Could you please tell me is there anything that I have to change in this subroutine to make it working. </div><div><br data-mce-bogus="1"></div><div>Tha value of m1 and z0 are just explicity giving the base velocity to the elements right? for example on the W or M.V conditions basev becomes 1 and on boundaries of the domain it becomes 0.0 <br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>      m1 = -1.<br>      z0 =  0.0<br><br>      nface = 2*ndim<br><br><br>      do e=1,nelv<br>      do f=1,nface<br><br>c        Set Dirichlet for mesh velocity on all boundaries<br><br>         if (cbc(f,e,1).ne.'E  ') call facev(msk,e,f,z0,nx1,ny1,nz1)<br><br>c        Inhomogeneous Dirichlet on cylinder only<br><br>         if (cbc(f,e,1).eq.'mv ') then<br>           call fcaver(xavg,xm1,e,f)<br><br><br>           if (xavg.gt.-3.0) call facev(tmp,e,f,m1,nx1,ny1,nz1)<br>         endif<br>      enddo<br>      enddo<br><br>      call axhelm (rhs,tmp,h1,h2,1,1)<br><br>      tol    = 1.e-4<br>      imsh   = 1<br>      ifield = 1<br>      call hmholtz('mshv',basev,rhs,h1,h2,msk,vmult,imsh,tol,200,1)<br>      call sub2(basev,tmp,n)<br><br>c     call outpost(basev,basev,basev,basev,basev,'   ')<br>c     stop<br><br><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Regards<br data-mce-bogus="1"></div><div>Sijo<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div></div></body></html>