<div dir="ltr"><div><div><div><div><div>Hi Rich,<br><br></div>I have updated the subroutine that dumps out ASCII files to expend the format for larger element counts.<br><br></div>Keeping <br>p66 = -1 ! write <br>p67 = -1 ! read<br>
<br></div>in your .rea file is fine. <br><br></div>Your best option is to re-run the case that produced this restart file so that the header is now correct.<br><br>If you don't want to update the source and re-run the case that produced this restart file, it "should" work if you change the first line of the restart file from: <br>
**** 8 8 8 8.0000000E+01 1500 X Y Z U P NELT,NX,NY,N<br><br></div>to<br><div> 18432 8 8 8 8.800000000E+01 1500 X Y Z U P NELT,NX,NY,N<br><br></div><div>spacing is important. <br>
<br></div><div>Best,<br></div><div>Katie<br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jan 30, 2013 at 5:04 PM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
Hi Rich,<br>
<br>
Try, in usrdat2():<br>
<br>
param(67)=0 ! Set read format<br>
<br>
<br>
(66 is for write format)<br>
<br>
Let me know if this works.<br>
<br>
Cheers,<br>
<br>
Paul<div class="im HOEnZb"><br>
<br>
<br>
<br>
On Wed, 30 Jan 2013, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> wrote:<br>
<br>
</div><div class="HOEnZb"><div class="h5"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Nek heads,<br>
I am having trouble getting Nek5000 to read my restart file. The file is an<br>
ASCII file that was written from a previous run. I set parameter p67 to -1<br>
to read in an ascii restart file. Here is the error message in the .o file:<br>
<br>
set initial conditions<br>
Checking restart options:<br>
sqduct.fld06<br>
<br>
nekuic (1) for ifld 1<br>
Reading checkpoint data<br>
WARNING: No data read in for file<br>
sqduct.fld06<br>
<br>
Failed on element 2, point 0.<br>
Last read dump:<br>
0.000000 0.2641195E-02 0.000000 0.000000 0.000000<br>
Last read dump:<br>
0.000000 0.000000 -0.6859594E-02 0.2343750E-01 -0.2600098E-01<br>
Last read dump:<br>
0.2600098E-01 -0.1855469E-01<br>
0 call exitt 1702a 1<br>
<br>
<br>
Here is the restart call in the .rea file:<br>
<br>
***** NO THERMAL BOUNDARY CONDITIONS *****<br>
1 RESTART OPTIONS FOLLOW<br>
sqduct.fld06<br>
1 INITIAL CONDITIONS *****<br>
C Default<br>
<br>
<br>
Here is the first few lines of the restart file sqduct.fld06:<br>
<br>
<br>
**** 8 8 8 8.0000000E+01 1500 X Y Z U P<br>
NELT,NX,NY,N<br>
0.000 0.000 0.000 0.000 0.000 0.000<br>
0.000 0.000 0.000 0.000 0.000 0.000<br>
0.000 0.000 0.000 0.000 0.000 0.000<br>
0.000 0.000 0.000 0.000 0.000 0.000<br>
0.000 0.000 0.000 0.000 0.000 0.000<br>
0.000 0.000 0.000 0.000 0.000 0.000<br>
<br>
<br>
Note that the problem is a square duct with 24 x 24 x 32 elements, which<br>
equal 18432 total elements. Since this is a 5 digit number, I suppose that<br>
is why it is not shown in the first line in the field file. However, I<br>
tried changing the four *'s to 18432, along with a few other variations,<br>
and no luck, same error. Also, the numbers in the dump after the error<br>
message above don't seem to appear in the restart file. Also, I don't have<br>
any trouble withe a restart file that has only 6 x 6 x 8 = 288 elements in<br>
it. Finally, I am using a version of Nek500 that I downloaded about 11<br>
months ago.<br>
<br>
Help is appreciated.<br>
<br>
Rich Johnson<br>
<br>
<br>
<br>
<br>
<br>
On Tue, Jan 29, 2013 at 10:27 PM, <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a>> wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
On Wed, Jan 30, 2013 at 10:13 AM, <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a>> wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello<br>
I want to identify boundary faces to set their boundary condition type.<br>
Following some of the examples, I created the following function. To find<br>
the location of the face, I need to know the coordinates (xf,yf) of the<br>
face center using which I can figure out where the face lies. Is this<br>
possible to find ?<br>
<br>
Thanks<br>
praveen<br>
<br>
</blockquote>
<br>
Hello<br>
<br>
By looking at this<br>
<br>
<a href="http://nek5000.mcs.anl.gov/index.php/Element_Definition" target="_blank">http://nek5000.mcs.anl.gov/<u></u>index.php/Element_Definition</a><br>
<br>
I made the following code. This is a 2d problem and in blue part below I<br>
am getting face midpoint coordinates. I am not sure I have got the<br>
numbering of the faces correct.<br>
<br>
Thanks<br>
praveen<br>
<br>
C=============================<u></u>==============================<u></u>============<br>
subroutine usrdat2<br>
<br>
include 'SIZE'<br>
include 'TOTAL'<br>
<br>
parameter(XTOL=1e-10)<br>
character*3 cbv<br>
integer e, f<br>
real xmin, xmax, ymax, xf(4), yf(4)<br>
<br>
xmin = -5.0<br>
xmax = +10.0<br>
ymax = 3.0<br>
<br>
c set boundary conditions<br>
do e = 1,nelv<br>
xf(3) = 0.5*(xm1(1, 1, 1,e) + xm1(nx1,1, 1,e))<br>
xf(2) = 0.5*(xm1(nx1,1, 1,e) + xm1(nx1,ny1,1,e))<br>
xf(4) = 0.5*(xm1(1, ny1,1,e) + xm1(nx1,ny1,1,e))<br>
xf(1) = 0.5*(xm1(1, ny1,1,e) + xm1(1, 1, 1,e))<br>
yf(3) = 0.5*(ym1(1, 1, 1,e) + ym1(nx1,1, 1,e))<br>
yf(2) = 0.5*(ym1(nx1,1, 1,e) + ym1(nx1,ny1,1,e))<br>
yf(4) = 0.5*(ym1(1, ny1,1,e) + ym1(nx1,ny1,1,e))<br>
yf(1) = 0.5*(ym1(1, ny1,1,e) + ym1(1, 1, 1,e))<br>
do f = 1,2*ndim<br>
cbv = cbc(f,e,1)<br>
if (cbv.ne.'E ') then ! This is boundary edge<br>
if(abs(xf(f)-xmin).lt.XTOL)<u></u>then<br>
cbc(f,e,1) = 'v ' ! inlet<br>
else if(abs(xf(f)-xmax).lt.XTOL)<u></u>then<br>
cbc(f,e,1) = 'O ' ! outlet<br>
else if(abs(yf(f)-ymax).lt.XTOL)<u></u>then<br>
cbc(f,e,1) = 'SYM' ! top wall<br>
else<br>
cbc(f,e,1) = 'W ' ! bottom wall<br>
endif<br>
endif<br>
enddo<br>
enddo<br>
<br>
return<br>
end<br>
<br>
<br>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
<br>
<br>
</blockquote>
<br>
</blockquote>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</div></div></blockquote></div><br></div>