<div dir="ltr"><div class="gmail_default" style="font-size:small">Hello Neks,</div><div class="gmail_default" style="font-size:small">Does the code go into usrdat2 and modify the geometry first or does it go into userchk first? I initialize the passive scalar field at istep=0</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">i am running simulation with three passive scalars modifying  conj_ht example. To save computational cost I run the fluid simulation and setting the IC of passive scalars to be zero and switching off the advection, so it remains zero and simulation runs quickly. </div><div class="gmail_default" style="font-size:small"><div class="gmail_default">  T     IFHEAT</div><div class="gmail_default">  T     IFTRAN</div><div class="gmail_default">  T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)</div><div class="gmail_default">  F F T T T F F F F F F F IFTMSH (IF mesh for this field is T mesh)</div><div class="gmail_default"><br></div><div class="gmail_default">After the simulation reach steadystate, I use the final fld file to start the passive scalar computations</div></div><div class="gmail_default" style="font-size:small">I have the following in userchk, where the initial conditions are imposed when istep=0</div><div class="gmail_default">      if(istep.eq.0) then<br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default">       x = xm1(i,1,1,1)</div><div class="gmail_default">       y = ym1(i,1,1,1)</div><div class="gmail_default">       R1 = sqrt(x*x+y*y)<br></div><div class="gmail_default">       if((<a href="http://R1.le.Ro">R1.le.Ro</a>).and.(R1.ge.Ri)) then<br></div></div><div class="gmail_default">           t(i,1,1,1,3)  =  fr*y/R1* (sin(2.0*pi*z/4.0)+sin(2.0*2.0*pi*z/4.0))<br></div><div class="gmail_default"><div class="gmail_default"><br></div></div><div class="gmail_default" style="font-size:small">I scale the domain in usrdat2 as following:</div><div class="gmail_default" style="font-size:small"><div class="gmail_default">      call rescale_x(xm1,-8.0,8.0)</div><div class="gmail_default">      call rescale_x(ym1,-8.0,8.0) </div><div class="gmail_default"><br></div><div class="gmail_default">Thank you in advance</div><div class="gmail_default">Sincerely,</div><div class="gmail_default">Sandeep</div></div></div>