gtd7: fix output at 32.5 km master
authorDominik Brodowski <linux@dominikbrodowski.net>
Tue, 9 Jul 2019 16:26:01 +0000 (18:26 +0200)
committerDominik Brodowski <linux@dominikbrodowski.net>
Tue, 9 Jul 2019 16:26:01 +0000 (18:26 +0200)
commitbc9a2feba4344e74201281e563332688a4d09cc3
tree41df72ce67c6c07e3bff10bcb01a3f010caea8dc
parent4693ab643732f4ee7bf9deda1412452ea805f449
gtd7: fix output at 32.5 km

The original FORTRAN version, in line 221,

       IF(ALT.GE.ZN3(1)) GOTO 6

and the previous C releases contain a bug for calculations at
32.5 km as the upper boundary for what the model calls lower
stratosphere and troposphere; for this exact altitude,
invalid outputs ("-NAN") may be generated. Fix this issue by
not skipping calculations right at this boundary.

Reported-by: Yoshiaki Ando <y-ando@uec.ac.jp>
Signed-off-by: Dominik Brodowski <linux@dominikbrodowski.net>
DOCUMENTATION
nrlmsise-00.c