Fixing Kp=0 bug for Gallagher

2026-01-08 22:58:05 -05:00 · 2025-11-09 20:28:45 -07:00
parent 7f23f40b30
commit af6762b141
1 changed files with 5 additions and 0 deletions
--- a/src/base/earthhelper.F90
+++ b/src/base/earthhelper.F90
@@ -96,6 +96,11 @@ module earthhelper

        D = 0.0

+        xfn = 0.0_rp
+        yfn = 0.0_rp
+        xfd = 0.0_rp
+        yfd = 0.0_rp
+
        if (present(kpO)) then
            kp = nint(kpO) ! Gallagher model takes integer Kp
        else