Convert from geographic to geomagnetic coordinates

A MATLAB function that converts from geographic to geomagnetic coordinates

This is a MATLAB function that computes the geomagnetic latitudes of a location using the geographic coordinates of the location.

The syntax is:

geom_lat = geog2geom(geog_lat, geog_long, height, decimalYear)

The inputs are:

geog_lat: Geographic latitude, specified as a scalar or vector, in degrees. North latitude is positive and south latitude is negative.

geog_long: Geographic longitude specified as a scalar or vector, in degrees. East longitude is positive and west longitude is negative.

height: Distance from the surface of the Earth, specified as a scalar or vector, in meters.

decimalYear: Year, in decimal format, specified as a scalar or vector. This value can have any fraction of the year that has already passed.

The geog_lat, geog_long, height, and decimalYear arguments must all be the same size (scalars or vectors).

The output is:

geom_lat: Geomagnetic latitude, specified as a scalar or vector, in degrees. North latitude is positive and south latitude is negative.

The size of geom_lat is same as the size of each of the inputs