[FLASH-USERS] Nonconvergence in subroutine rieman

[Esha Kundu]

Dear all,


I'm running the flash code radflaHD in 1D. The code runs properly for some time and then got aborted because of the following error:


Nonconvergence in subroutine rieman

 Zone index       =           71
 Zone center      =    47374441964285712.
 Iterations tried =           12
 Pressure error   =    1.0471366062501732E-005
 rieman_tol       =    1.0000000000000001E-005

 pL       =    6.1379631463855528E-016  pR       =   6.1379670668751590E-016
 uL       =    298862323.76308846       uR       =   298862461.02809691
 cL       =    2.0001184805161569E-017  cR       =   2.0001182294886163E-017
 gamma_eL =    1.6662744792216422       gamma_eR =   1.6662744792216424
 gamma_cL =    1.6662744792216422       gamma_cR =   1.6662744792216424

 Iteration history:

   n                    p*
   1    0.100000000000E-29
   2    0.783386173068E-21
   3    0.434967960562E-21
   4    0.544861495611E-21
   5    0.579319637957E-21
   6    0.572926460603E-21
   7    0.574080784292E-21
   8    0.574176977933E-21
   9    0.574128881113E-21
  10    0.574152929523E-21
  11    0.574140905318E-21
  12    0.574146917420E-21

 Terminating execution.

 Driver_abort called.


Does anyone know why the rieman subroutine did not converge?


Regards,

Esha


Department of Astronomy
The Oskar Klein Centre
Stockholm University
10691 Stockholm
Sweden

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20171220/b7e2523a/attachment.htm>