<div dir="ltr">Thanks for the answer.<div><br></div><div>One last question: userdat2 looks like:</div><br><font face="monospace, monospace">do iel=1,nelt<br>   do ifc=1,2*ndim<br>     idss = bc(5,ifc,iel,1)<br>     if (idss.eq.1) cbc(ifc,iel,1)='v '<br>  enddo<br>enddo</font><div><br></div><div><br></div><div>the last value "1" in bc array, and cbc array refers to velocity field, right? So for temperature BC I should do the same but specifying "2"?</div><div><br></div><div>e.g</div><div><br></div><div><span style="font-family:monospace,monospace">do iel=1,nelt</span><br style="font-family:monospace,monospace"><span style="font-family:monospace,monospace">   do ifc=1,2*ndim</span><br style="font-family:monospace,monospace"><span style="font-family:monospace,monospace">     idss = bc(5,ifc,iel,2)</span><br style="font-family:monospace,monospace"><span style="font-family:monospace,monospace">     if (idss.eq.1) cbc(ifc,iel,2)='t '</span><br style="font-family:monospace,monospace"><span style="font-family:monospace,monospace">  enddo</span><br style="font-family:monospace,monospace"><span style="font-family:monospace,monospace">enddo</span><br></div><div><span style="font-family:monospace,monospace"><br></span></div><div><br></div><div><span style="font-family:monospace,monospace"><br></span></div><div><span style="font-family:monospace,monospace">JP.</span></div></div>