<div dir="ltr"><div class="gmail_default" style="font-size:small">Dear Neks,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I am using the MHD example for kinematic Dynamo.</div><div class="gmail_default" style="font-size:small">In the example gpf, the force is added as </div><div class="gmail_default" style="font-size:small">f=-Nabla^2 u - partial u/partial t, where u is</div><div class="gmail_default" style="font-size:small">         ux   = -sys</div><div class="gmail_default">         uy   = -cxc</div><div class="gmail_default">         uz   = sxc + cys</div><div class="gmail_default"><br></div><div class="gmail_default">I want to do it the straight forward approach i.e. directly imposing the velocity in userchk and switching off the iteration Navier-Stokes equation. </div><div class="gmail_default">For this I set the IFNAV to False</div><div class="gmail_default"><div class="gmail_default">  F T F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)</div><div class="gmail_default">  F F F T T T T T T T T T IFTMSH (IF mesh for this field is T mesh)</div><div class="gmail_default"><br></div><div class="gmail_default">IC being the same as imposed velocity</div><div><br></div><div>and in userchk added the following</div><div><br></div><div><div>      if (ifield.eq.1) then<span class="gmail-Apple-tab-span" style="white-space:pre">  </span>! velocity</div><div>         omega = 1.</div><div>      do i=1,nx1*ny1*nz1*nelv   </div><div><span class="gmail-Apple-tab-span" style="white-space:pre">          </span> x = xm1(i,1,1,1); y = ym1(i,1,1,1); z = zm1(i,1,1,1)</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">          </span> </div><div>         st  = sin(omega*time)</div><div>         ct  = cos(omega*time)</div><div>c</div><div>         argxc = x + ct</div><div>         argyc = y + ct</div><div>         argxs = x + st</div><div>         argys = y + st</div><div>c</div><div>         visc = param(2)</div><div>         s    = 3./2.</div><div>         s    = sqrt(s)</div><div>c        </div><div>         cxc = s*cos(argxc)</div><div>         cys = s*cos(argys)</div><div>         sxc = s*sin(argxc)</div><div>         sys = s*sin(argys)</div><div>c        </div><div>         vx(i,1,1,1) = -sys</div><div>         vy(i,1,1,1) = -cxc </div><div>         vz(i,1,1,1) = sxc + cys</div><div>      enddo    </div><div>      endif</div><div><br></div><div>      ibstep = iostep</div><div>      if (istep.gt.0.and.mod(istep,ibstep).eq.0)</div><div>     $   call outpost(bx,by,bz,pm,t,'   ')<span class="gmail-Apple-tab-span" style="white-space:pre">      </span>! odd blah.fld/.f for B-field</div><div>c    $   call outpost(bx,by,bz,pm,t,'mag')<span class="gmail-Apple-tab-span" style="white-space:pre">     </span>!  magblah.fld/.f for B-field</div><div>      if (ifxyo.and.istep.gt.istep) ifxyo = .false.</div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span>  kinerg2 = 0.5*(glsc3(vx,bm1,vx,m)</div><div>     $        +glsc3(vy,bm1,vy,m)</div><div>     $        +glsc3(vz,bm1,vz,m))/volvm1</div><div> </div><div><br></div><div>      n = nx1*ny1*nz1*nelfld(ifldmhd)<span class="gmail-Apple-tab-span" style="white-space:pre">           </span>! Magnetic energy !?nelb</div><div>      energy0 = energy</div><div>      energy  = glsc3(bx,bm1,bx,n)</div><div>     $        + glsc3(by,bm1,by,n) </div><div>     $        + glsc3(bz,bm1,bz,n) </div><div>      energy  = 0.5*energy/volfld(ifldmhd)</div><div> </div><div>      if (istep.gt.1) then</div><div>         ratio  = energy/energy0</div><div>         growth = 0.5*log(ratio)/dt<span class="gmail-Apple-tab-span" style="white-space:pre">               </span>! 1/2 of magnetic energy growth rate</div><div>      if (nid.eq.0) write(6,1) istep,time,growth,energy,</div><div>     $ kinerg1,kinerg2 </div><div>    1    format(i8,1p5e15.6,' growth')</div><div> </div><div>         nstepag = param(120)<span class="gmail-Apple-tab-span" style="white-space:pre">                   </span>! p121 < t_avg < p120</div><div>         nstepsa = param(121)</div><div>         if (istep.lt.nstepag) call</div><div>     $      runtimeavg(Emga,growth,1,nstepsa,1,'gr_Em')<span class="gmail-Apple-tab-span" style="white-space:pre">       </span>! avg. growth</div><div>      endif</div><div><br></div><div>When I calculate the energy of the system, it is decaying. </div><div>Is there anyway I can tamper with / change the velocity field array in nek5000 from userchk?</div><div><br></div><div>Thanking you in advance.</div><div><br></div><div>Cheers,</div><div>Sandeep</div><div><br></div><div>PS: Attached are the .usr and .rea files</div><div> </div></div></div><div><br></div></div>