<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<p>Thank you for your solution, it works great.<br>
</p>
<p>Thank you for your time, <br>
</p>
<p>Best wishes,</p>
<p><br>
</p>
<p>Arnold<br>
</p>
<p><br>
</p>
<p><br>
</p>
<br>
<div class="moz-cite-prefix">Le 30/05/2017 à 15:54, <a
class="moz-txt-link-abbreviated"
href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>
a écrit :<br>
</div>
<blockquote
cite="mid:mailman.7627.1496152467.2967.nek5000-users@lists.mcs.anl.gov"
type="cite">Hi; <br>
The easies solution would be to copy your velocity and pressure to
temporary arrays before reading the file and placing them back
later on. e.g. <br>
<br>
real tmp_vel(lx1*ly1*lz1*lelv,ldim),
tmp_pr(lx2*ly2*lz2*lelv) <br>
common /SCREV/ tmp_vel, tmp_pres <br>
integer n_tot <br>
<br>
<br>
call
opcopy(tmp_vel(1,1),tmp_vel(1,2),tmp_vel(1,ndim),vx,vy,vz) <br>
n_tot = nx2*ny2*nz2*nelv <br>
call copy(tmp_pr,pr) <br>
<br>
! read file and use variables <br>
<br>
call
opcopy(vx,vy,vz,tmp_vel(1,1),tmp_vel(1,2),tmp_vel(1,ndim)) <br>
n_tot = nx2*ny2*nz2*nelv <br>
call copy(pr,tmp_pr) <br>
<br>
<br>
Regards <br>
Adam <br>
<br>
On 30/05/17 15:32, <a class="moz-txt-link-abbreviated"
href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>
wrote: <br>
<blockquote type="cite">Hi Neks, <br>
<br>
I may have a problem when I'm trying to read velocity field
(only velocity field, I don't need pressure.) so I can assign
specific boundary conditions using userbc(). <br>
<br>
My questions are : 1) Is it possible to read velocity fields in
nek that are contained in .f000* files without modifying the
pressure "pr" nor the velocity fields "vx,vy,vz" ? <br>
2) Is there any subroutine
(besides load_fld or mfi) that would allow me to read .f000*
files? <br>
<br>
<br>
When I use the mfi (in ic.f) subroutine, it looks like the
pressure field is modified and hence I'm quickly running into
errors like the following: <br>
<br>
100 **ERROR**: Failed in HMHOLTZ: VELX NaN NaN
1.0000E-12. <br>
<br>
If I use load_fld() I seem to me going back to the initial
conditions. <br>
<br>
Thank you for your help, <br>
<br>
Arnold <br>
_______________________________________________ <br>
Nek5000-users mailing list <br>
<a class="moz-txt-link-abbreviated"
href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a>
<br>
<a class="moz-txt-link-freetext"
href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>
<br>
</blockquote>
<br>
<br>
_______________________________________________ <br>
Nek5000-users mailing list <br>
<a class="moz-txt-link-abbreviated"
href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a>
<br>
<a class="moz-txt-link-freetext"
href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>
<br>
</blockquote>
<br>
</body>
</html>