We ran smoke-beis version 4.8 and found NaN from line 1190 EMSIL in tmpbeis361. Also, there wasn’t NaN in output of normbeis361. Any suggestions where the NaN may generated?
We solved the issue after finding a precision bug in zenith angle calculation source code (…/src/biog/czangle.f) under SMOKE-beis version 4.8. We declare all related parameters that involved in zenith angle calculation with double precision/real8 to replace the original single precision/real4. Then, the calculation reflect correctly and NaN disappear.