Dear Nek heads,<br>I am having trouble getting Nek5000 to read my restart file. The file is an ASCII file that was written from a previous run. I set parameter p67 to -1 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: sqduct.fld06 <br> nekuic (1) for ifld 1<br> Reading checkpoint data<br>
WARNING: No data read in for file sqduct.fld06 <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 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 equal 18432 total elements. Since this is a 5 digit number, I suppose that is why it is not shown in the first line in the field file. However, I tried changing the four *'s to 18432, along with a few other variations, and no luck, same error. Also, the numbers in the dump after the error message above don't seem to appear in the restart file. Also, I don't have any trouble withe a restart file that has only 6 x 6 x 8 = 288 elements in it. Finally, I am using a version of Nek500 that I downloaded about 11 months ago. <br>
<br>Help is appreciated.<br><br>Rich Johnson<br><br><br><br><br><br><div class="gmail_quote">On Tue, Jan 29, 2013 at 10:27 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"><div dir="ltr"><div class="gmail_extra"><div class="im"><br><div class="gmail_quote">On Wed, Jan 30, 2013 at 10:13 AM, <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-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">
Hello<div>I want to identify boundary faces to set their boundary condition type. Following some of the examples, I created the following function. To find the location of the face, I need to know the coordinates (xf,yf) of the face center using which I can figure out where the face lies. Is this possible to find ?</div>
<div><br></div><div>Thanks</div><div>praveen</div></div></blockquote></div><br></div>Hello</div><div class="gmail_extra"><br></div><div class="gmail_extra">By looking at this</div><div class="gmail_extra"><br></div><div class="gmail_extra">
<a href="http://nek5000.mcs.anl.gov/index.php/Element_Definition" target="_blank">http://nek5000.mcs.anl.gov/index.php/Element_Definition</a><br></div><div class="gmail_extra"><br></div><div class="gmail_extra">I made the following code. This is a 2d problem and in blue part below I am getting face midpoint coordinates. I am not sure I have got the numbering of the faces correct.</div>
<div class="gmail_extra"><br></div><div class="gmail_extra">Thanks</div><div class="gmail_extra">praveen</div><div class="gmail_extra"><br></div><div class="gmail_extra"><div class="im"><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace">C=======================================================================</font></div>
<div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> subroutine usrdat2</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"><br>
</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> include 'SIZE'</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> include 'TOTAL'</font></div>
<div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"><br></font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> parameter(XTOL=1e-10)</font></div>
<div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> character*3 cbv</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> integer e, f</font></div>
</div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> real xmin, xmax, ymax, xf(4), yf(4)</font></div><div class="im"><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"><br>
</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> xmin = -5.0</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> xmax = +10.0</font></div>
<div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> ymax = 3.0</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"><br>
</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace">c set boundary conditions</font></div></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> do e = 1,nelv</font></div>
<div class="gmail_extra"><font style color="#0000ff" face="'courier new', monospace"> xf(3) = 0.5*(xm1(1, 1, 1,e) + xm1(nx1,1, 1,e))</font></div><div class="gmail_extra">
<font style color="#0000ff" face="'courier new', monospace"> xf(2) = 0.5*(xm1(nx1,1, 1,e) + xm1(nx1,ny1,1,e))</font></div><div class="gmail_extra"><font style color="#0000ff" face="'courier new', monospace"> xf(4) = 0.5*(xm1(1, ny1,1,e) + xm1(nx1,ny1,1,e))</font></div>
<div class="gmail_extra"><font style color="#0000ff" face="'courier new', monospace"> xf(1) = 0.5*(xm1(1, ny1,1,e) + xm1(1, 1, 1,e))</font></div><div class="gmail_extra">
<font style color="#0000ff" face="'courier new', monospace"> yf(3) = 0.5*(ym1(1, 1, 1,e) + ym1(nx1,1, 1,e))</font></div><div class="gmail_extra"><font style color="#0000ff" face="'courier new', monospace"> yf(2) = 0.5*(ym1(nx1,1, 1,e) + ym1(nx1,ny1,1,e))</font></div>
<div class="gmail_extra"><font style color="#0000ff" face="'courier new', monospace"> yf(4) = 0.5*(ym1(1, ny1,1,e) + ym1(nx1,ny1,1,e))</font></div><div class="gmail_extra">
<font style color="#0000ff" face="'courier new', monospace"> yf(1) = 0.5*(ym1(1, ny1,1,e) + ym1(1, 1, 1,e))</font></div><div class="im"><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> do f = 1,2*ndim</font></div>
<div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> cbv = cbc(f,e,1)</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> if (cbv.ne.'E ') then ! This is boundary edge</font></div>
</div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> if(abs(xf(f)-xmin).lt.XTOL)then</font><br></div><div class="im"><div class="gmail_extra">
<font style color="#ff0000" face="'courier new', monospace"> cbc(f,e,1) = 'v ' ! inlet</font></div></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> else if(abs(xf(f)-xmax).lt.XTOL)then</font></div>
<div class="im">
<div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> cbc(f,e,1) = 'O ' ! outlet</font></div></div><div class="gmail_extra">
<font style color="#ff0000" face="'courier new', monospace"> else if(abs(yf(f)-ymax).lt.XTOL)then</font></div><div class="im"><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> cbc(f,e,1) = 'SYM' ! top wall</font></div>
<div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> else</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> cbc(f,e,1) = 'W ' ! bottom wall</font></div>
<div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> endif</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> endif</font></div>
<div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> enddo</font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> enddo</font></div>
<div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"><br></font></div><div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> return</font></div>
<div class="gmail_extra"><font style color="#ff0000" face="'courier new', monospace"> end</font></div><div><br></div></div></div></div>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br>