diff options
| author | Dominik Brodowski <linux@dominikbrodowski.net> | 2019-07-09 18:26:01 +0200 |
|---|---|---|
| committer | Dominik Brodowski <linux@dominikbrodowski.net> | 2019-07-09 18:26:01 +0200 |
| commit | bc9a2feba4344e74201281e563332688a4d09cc3 (patch) | |
| tree | 41df72ce67c6c07e3bff10bcb01a3f010caea8dc /nrlmsise-00.c | |
| parent | 4693ab643732f4ee7bf9deda1412452ea805f449 (diff) | |
| download | nrlmsise-00-master.tar.gz | |
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>
Diffstat (limited to 'nrlmsise-00.c')
| -rw-r--r-- | nrlmsise-00.c | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/nrlmsise-00.c b/nrlmsise-00.c index 1abad3b..6b63d29 100644 --- a/nrlmsise-00.c +++ b/nrlmsise-00.c @@ -976,7 +976,7 @@ void gtd7(struct nrlmsise_input *input, struct nrlmsise_flags *flags, struct nrl meso_tgn2[1]=pavgm[8]*pma[9][0]*(1.0+flags->sw[20]*flags->sw[22]*glob7s(pma[9], input, flags))*meso_tn2[3]*meso_tn2[3]/(pow((pma[2][0]*pavgm[2]),2.0)); meso_tn3[0]=meso_tn2[3]; - if (input->alt<zn3[0]) { + if (input->alt<=zn3[0]) { /* LOWER STRATOSPHERE AND TROPOSPHERE (below zn3[0]) * Temperature at nodes and gradients at end nodes * Inverse temperature a linear function of spherical harmonics |