proway2 / go-igrf Goto Github PK

• Find lat/lon near magnetic poles (both)
• Add unittests to check the correctness of values produced for these points

• Define how the library will be used and create the interface.
• Add unittests for this interface (these may not work at the beginning).

In shval3.c some indexes start with 1. Check if it's possible to switch to the traditional 0-based indexing and fix the issue.

Decide how to return and return warnings (there might be more than one) when e.g. horizontal intensity is weak. Check warn_* variables in C code.

Re-read coeffs. Consider make it similar to those available at C implementation. Add degree and order for h.

Implement analog to subroutine shval3 in C.

If igrf.IGRF called few times it reads SHC again and again every time. Refactor igrf.IGRF and come up with the approach that solves this, so that SHC is read once and used for any number of calculations.

Add How to use section to README.md

When date is higher than end_epoch the factor is calculated the very same way it's done in C implementation. However, this calculation doesn't respect leap years, although when date is lower than end_epoch, factor is calculated with respect to leap years.

Check if it's possible to improve accuracy by calculating factor with leap years when date is higher than end_epoch.

Create an extensive set of fixtures for unittests. Including edge cases.

• Define limits for alt argument.
• Add a check
• Add unittests

Testdata was created with FORTRAN, some values are coarse, especially all SV values. Allowed error (used in test cases) for this values is very high.

• Check if it's possible to create more accurate testdata with C version.
• Generate the new testdata
• Implement a way to digest the new testdata and automatically create tests bases on it.
• Refactor tests - decrease max allowed error for SV fields.

Add tests for latitudes near North and South poles:

• 89,9995
• -89,9995

Calculate:
Declination| (D)
DeclinationSV
Inclination| (I)
InclinationSV
HzIntensity (H)
HorizontalSV
NorthSV (X)
EastSV (Y)
VerticalSV (Z)
TotalIntensity (F)
TotalSV (F)

Switch all values from float32 to float64.

Add meaningful comments for private functions.

Compute geocentric colatitude and radius from geodetic colatitude and height. This is a gg_to_geo(h, gdcolat) from Python.

• Read the model coeffs;
• Find coeffs by interpolation for a certain date.

Add go 1.19 version to the testing flow.

More test cases needed for cases when Horizontal Intensity (H):

• h < 100.0
• h < 1000.0
• h < 5000.0