<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>