[Nek5000-users] Error in generalev with MHD

[nek5000-users at lists.mcs.anl.gov]

Hello again,
I changed all BC's for the magnetic field to 'v' and used userbc to set 
the value to zero, but the error popped up again. The output reads:


216 Matrix:   0 aa         6       6       6
217   0 aa     7.80579E+01  -6.39449E+01   1.81958E+01 -7.90484E+00   
2.16525E+00   0.00000E+00
218   0 aa    -6.39449E+01   8.67803E+01  -4.74885E+01 1.84117E+01 
-1.07792E+01   5.66872E+00
219   0 aa     1.81958E+01  -4.74885E+01   5.96643E+01 -4.03684E+01   
3.02051E+01  -2.06952E+01
220   0 aa    -7.90484E+00   1.84117E+01  -4.03684E+01 6.64609E+01 
-7.64969E+01   4.76373E+01
221   0 aa     2.16525E+00  -1.07792E+01   3.02051E+01 -7.64969E+01   
1.90127E+02  -1.67410E+02
222   0 aa     0.00000E+00   5.66872E+00  -2.06952E+01 4.76373E+01 
-1.67410E+02      Infinity
223
224 Matrix:   0 bb         6       6       6
225   0 bb     7.19606E-02   9.84283E-03  -2.07087E-03 -1.02002E-03   
7.34397E-04   0.00000E+00
226   0 bb     9.84283E-03   4.07143E-02  -4.22288E-03 7.31377E-03 
-3.21920E-03   0.00000E+00
227   0 bb    -2.07087E-03  -4.22288E-03   1.05630E-01 -1.61560E-02   
6.44826E-03   0.00000E+00
228   0 bb    -1.02002E-03   7.31377E-03  -1.61560E-02 1.05960E-01 
-6.21764E-03   0.00000E+00
229   0 bb     7.34397E-04  -3.21920E-03   6.44826E-03 -6.21764E-03   
5.05023E-02   0.00000E+00
230   0 bb     0.00000E+00   0.00000E+00   0.00000E+00 0.00000E+00 
0.00000E+00           NaN
231
232 Matrix:   0 Aeig       6       6       6
233   0 Aeig           NaN           NaN           NaN NaN NaN           
NaN
234   0 Aeig           NaN           NaN           NaN NaN NaN           
NaN
235   0 Aeig           NaN           NaN           NaN NaN NaN           
NaN
236   0 Aeig           NaN           NaN           NaN NaN NaN           
NaN
237   0 Aeig           NaN           NaN           NaN NaN NaN           
NaN
238   0 Aeig           NaN           NaN           NaN NaN NaN           
NaN
239
240 Matrix:   0 Deig       1       6       6
241   0 Deig           NaN           NaN           NaN NaN NaN           
NaN
242  Error in generalev, info=           5           6           1
243
244 call exitt: dying ...


So the problem is probably the NaN and Inf entries. As a next step, I 
took the gpf example as a starting point, copy pasted my own grid and bc 
descriptions, set the bc's to zero in the usr-file and removed the 
Temperaturefield. This produced the same error. Can this be linked to 
the absence of the temperature field? If I leave out the B-Field the 
case is running just fine, altough the BC's are the same in the *.rea file.

Jan