<html><head><meta content="text/html; charset=utf-8" http-equiv="Content-Type"></head><body><div><div style="font-family: Calibri,sans-serif; font-size: 11pt;">Hi David,<br><br>The best place would probably be in userchk.<br>
<br>You could have a block like<br><br>if (istep.eq.0) then<br><br>C........do stuff<br><br>endif<br><br>If you are doing something else in userchk that would depend on your modification of the field quantities, then you would want this to come first.<br>
<br>Just in case, the velocity fields are accessed like<br><br>vx(i,j,k,e) ( for u)<br><br>and likewise for vy and vz, where i,j,k go from 1 to nx1 (ny1, nz1) and e goes from 1 to nelv. Sorry if that was more information than you needed, just wanted to make sure all my bases were covered :-).<br>
<br>Does this help?<br><br>Josh Camp<br>"All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke<br></div></div><hr><span style="font-family: Tahoma,sans-serif; font-size: 10pt; font-weight: bold;">From: </span><span style="font-family: Tahoma,sans-serif; font-size: 10pt;"><a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a></span><br>
<span style="font-family: Tahoma,sans-serif; font-size: 10pt; font-weight: bold;">Sent: </span><span style="font-family: Tahoma,sans-serif; font-size: 10pt;">10/4/2011 6:08 PM</span><br><span style="font-family: Tahoma,sans-serif; font-size: 10pt; font-weight: bold;">To: </span><span style="font-family: Tahoma,sans-serif; font-size: 10pt;"><a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a></span><br>
<span style="font-family: Tahoma,sans-serif; font-size: 10pt; font-weight: bold;">Subject: </span><span style="font-family: Tahoma,sans-serif; font-size: 10pt;">[Nek5000-users] Modifying an initial field</span><br><br></body></html>
Hi neks,<br><br>I am loading an initial condition from a .fld file, but I would like to modify the velocity field before using it as my initial condition. (I would like to superpose a shear onto the existing field.) Can someone tell me what part of the .usr file this should be done in?<br>
<br>Best,<br><br>David<br>