<div dir="ltr"><div><div><div><div>Hi Adam, <br><br></div>I am guessing I am using your code for stability tools in NEK. I have downloaded the code from <a href="https://nek5000.mcs.anl.gov/toolboxes/" target="_blank">https://nek5000.mcs.anl.gov/toolboxes/</a>.<br><br></div>I am trying to reproduce the ext_cyl example, as is, to get some familiarity. <br><br></div>I have been able to compile the code (after some modifications to compile on the platform I intent to run my analysis).  I have the base flow file "ext_cyl0.f00005" in the simulation folder. When I run the ext_cyl, I get the following error:<br><br><br>******************************************************************************************************************<br>call userchk<br> ARNOLDI:<br> Reading checkpoint<br>       0       0 OPEN: ARPext_cyl0.f00005<br>ARPext_cyl0.f00005<br>byte_read() :: fopen failure2!<br> ERROR: mfi_arnp, wrong header; abort.<br><br>call exitt: dying ...<br>******************************************************************************************************************<br></div><br><div><br><div>I have also pasted below  the parameters printed in the log file for your reference.  <br></div><div><br>******************************************************************************************************************<br></div>122   Parameters from file:/global/u2/r/rtumkur/eig_Setup/ARNOLDI/ext_cyl.rea<br>   1      1.00000     p001 DENSITY<br>   2     -50.0000     p002 VISCOS<br>   7      1.00000     p007 RHOCP<br>   8      1.00000     p008 CONDUCT<br>  11      50003.0     p011 NSTEPS<br>  12    -5.000000E-03 p012 DT<br>  15      2500.00     p015 IOSTEP<br>  17      1.00000     p017<br>  18     0.500000E-01 p018 GRID < 0 --> # cells on screen<br>  19     -1.00000     p019 INTYPE<br>  20      7.00000     p020 NORDER<br>  21     0.100000E-09 p021 DIVERGENCE<br>  22     0.100000E-09 p022 HELMHOLTZ<br>  24     0.100000E-09 p024 TOLREL<br>  25     0.100000E-04 p025 TOLABS<br>  26      1.00000     p026 COURANT/NTAU<br>  27      3.00000     p027 TORDER<br>  28      0.00000     p028 TORDER: mesh velocity (0: p28=p27)<br>  31      -1.00000     p031 NPERT: #perturbation modes<br>  59      0.00000     p059 !=0 --> full Jac. eval. for each el.<br>  65      1.00000     p065 #iofiles (eg, 0 or 64); <0 --> sep. dirs<br>  66      4.00000     p066 output : <0=ascii, else binary<br>  67      4.00000     p067 restart: <0=ascii, else binary<br>  69      1.00000     p069 : PROPER RESTART: 0 - do not use my_full_restart; 1- use it<br>  71      5.00000     p071 : ARNOLDI: checkpoint file number for restart (reading)<br>  72      6.00000     p072 : ARNOLDI: checkpoint file number for saving data<br>  99      3.00000     p099 dealiasing: <0--> off/3--> old/4--> new<br> 102      1.00000     p102 Dump out divergence at each time step<br> 108      1.00000     p108 : ARNOLDI: 0 if no arnoldi<br> 109      40.0000     p109 : ARNOLDI: frequency of calling arn_solve<br> 110      4.00000     p110 : ARNOLDI: maximal number of arnoldi cycles<br> 111      90.0000     p111 : ARNOLDI: size of Krylov space<br> 112      30.0000     p112 : ARNOLDI: number of eigenvectors<br> 113      1.00000E-06 p113 : ARNOLDI: arpack tolerance<br> 119      1.00000     p119 : SPONGE STRENGTH<br> 120      25.0000     p120 : SPONGE TOTAL WIDTH<br> 121      3.00000     p121 : SPONGE DROP WIDTH (INFLOW)<br> 122      5.00000     p122 : SPONGE RISE WIDTH (OUTFLOW)<br>******************************************************************************************************************<br><br></div><div>I'd appreciate if you can let me what am I missing in the set up which might be leading to this issue? Please let me know if you need any other information to understand this issue.<br><br><br></div><div>Regards, <br></div><div>Ravi<br><br></div><div><div><div><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 15, 2016 at 4:06 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi<br>
<br>
To force base flow to be constant you have to set negative param(31) in your ###.rea file. Something like<br>
   -1.00000     p031 NPERT: #perturbation modes<br>
Of course you need the base flow as well. It could be analytical one set in useric, or just a field read during initialisation.<br>
Regards<br>
Adam<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
On 2016-03-15 08:49, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Neks,<br>
<br>
I'm trying to use the linear stability tool of Nek5000.<br>
<br>
To start with, I am trying to run the ext_cyl example. However, I have found the following problem:<br>
<br>
ERROR: arnoldi assumes constant base flow.<br>
<br>
I have found on the README tutorial that the base flow has to be specified on the .rea file. I have included a file in the PRESOLVE/RESTART option, however I still have the same error.<br>
<br>
What should I modify?<br>
<br>
Finally, I have also used a subroutine defined in the .usr file that reads the file containing the base flow. However, the error persist...<br>
<br>
Could you help me, please?<br>
<br>
Cheers<br>
SL<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><span style="font-family:arial narrow,sans-serif"><font size="4">Ravi Kumar R. Tumkur</font>,<b>
</b><font size="4">PhD.</font><br>
Phone: <a value="+12177787538">217 778 7538</a></span></div></div>
</div>