<div dir="ltr">Hi all,<div><br></div><div>If I need to add an oscillating source term in the energy equation for a small region in a 2D domain would it be fine to do it this way, with these variables: ?</div><div><br></div><div>c =============================</div><div><div>      subroutine userq  (ix,iy,iz,ieg)</div><div>      include 'SIZE'</div><div>      include 'TOTAL'</div><div>      include 'NEKUSE'</div><div><br></div><div><br></div><div>      real pii,frec,ampli</div><div><br></div><div>      pii=4.0*atan(1.0)</div><div>      frec=0.005</div><div>      ampli=1.0</div><div><br></div><div><br></div><div>      if (x .ge. 0 .and. x .le. 0.02 .and.</div><div>     $ y .ge. 0 .and. y .le. 0.02) then</div><div><br></div><div>          qvol = ampli*sin(2.0*pii*frec*time)</div><div><br></div><div>      else</div><div><br></div><div>          qvol = 0.0</div><div><br></div><div>      end if</div><div><br></div><div><br></div><div>      source = 0.0</div><div><br></div><div>      return</div><div>      end</div></div><div>c ==================================</div><div><br></div><div><br></div><div>Thanks in advance,</div><div><br></div><div>JP.</div></div>