Re: how to use Pulverized coal combustion module by using GU
Posted: Tue May 08, 2012 12:48 am
Hello,
Running under a debugger, I can see that the error occurs in cpflux.f90, line 315 due to a zero value of xdfchi = ahetch(ichcor(icla))*exp(-ehetch(ichcor(icla)*4185/(rr*propce(iel,ipctem)))
for cell 1 (iel = 1).
In turn, this is due to ehetch(ichcor(icla)) being very high (13.10e+12), so exp(-....) is close to zero.
This variable is set in your dp_FCP file, so you probably need to use a more lower value (check that your units are correct: this activation energy is in kcal/mol, so 1.3e+13 seems a bit high...)
In case you need more debugging, I recommend building a debug build of the code (see installation documentation), so you can run that version. It is slower, but the backtrace in case of error gives you the file and line number of the error, even when not running under a debugger.
Best regards,
Yvan
Running under a debugger, I can see that the error occurs in cpflux.f90, line 315 due to a zero value of xdfchi = ahetch(ichcor(icla))*exp(-ehetch(ichcor(icla)*4185/(rr*propce(iel,ipctem)))
for cell 1 (iel = 1).
In turn, this is due to ehetch(ichcor(icla)) being very high (13.10e+12), so exp(-....) is close to zero.
This variable is set in your dp_FCP file, so you probably need to use a more lower value (check that your units are correct: this activation energy is in kcal/mol, so 1.3e+13 seems a bit high...)
In case you need more debugging, I recommend building a debug build of the code (see installation documentation), so you can run that version. It is slower, but the backtrace in case of error gives you the file and line number of the error, even when not running under a debugger.
Best regards,
Yvan