sciapy.level2.igrf.gmag_igrf()

Complexity

Conditions

2

Size

Total Lines	69
Code Lines	36

Duplication

Lines	69
Ratio	100 %

Importance

Changes

0

# -*- coding: utf-8 -*-
def gmag_igrf(date, lat, lon, alt=0.,

This code seems to be duplicated in your project.

Too many arguments (6/5)

This function exceeds the maximum number of variables (33/15).

		centered_dipole=False,

Wrong continued indentation (add 12 spaces).

		igrf_name="IGRF.tab"):

Wrong continued indentation (add 12 spaces).

	"""Centered or eccentric dipole geomagnetic coordinates

	Parameters
	----------
	date: `datetime.datetime` instance
	lat: float
		Geographic latitude in degrees north
	lon: float
		Geographic longitude in degrees east
	alt: float, optional
		Altitude in km. Default: 0.
	centered_dipole: bool, optional
		Returns the centered dipole geomagnetic coordinates
		if set to True, returns the eccentric dipole
		geomagnetic coordinates if set to False.
		Default: False
	igrf_name: str, optional
		Default: "IGRF.tab"

	Returns
	-------
	geomag_latitude: numpy.ndarray or float
		Geomagnetic latitude in eccentric dipole coordinates,
		centered dipole coordinates if `centered_dipole` is True.
	geomag_longitude: numpy.ndarray or float
		Geomagnetic longitude in eccentric dipole coordinates,
		centered dipole coordinates if `centered_dipole` is True.
	"""
	ellip = _ellipsoid()
	glat, glon, grad = _geod_to_spher(lat, lon, ellip, alt)
	(lat_GMP, lon_GMP), _, (dX, dY, dZ), B_0 = gmpole(date, ellip.re, igrf_name)

The name lat_GMP does not conform to the variable naming conventions ((([a-z][a-z0-9_]{2,30})|(_[a-z0-9_]*))$).

The name lon_GMP does not conform to the variable naming conventions ((([a-z][a-z0-9_]{2,30})|(_[a-z0-9_]*))$).

The name dX does not conform to the variable naming conventions ((([a-z][a-z0-9_]{2,30})|(_[a-z0-9_]*))$).

The name dY does not conform to the variable naming conventions ((([a-z][a-z0-9_]{2,30})|(_[a-z0-9_]*))$).

The name dZ does not conform to the variable naming conventions ((([a-z][a-z0-9_]{2,30})|(_[a-z0-9_]*))$).

The name B_0 does not conform to the variable naming conventions ((([a-z][a-z0-9_]{2,30})|(_[a-z0-9_]*))$).

	latr, lonr = np.radians(glat), np.radians(glon)
	lat_GMPr, lon_GMPr = np.radians(lat_GMP), np.radians(lon_GMP)

The name lat_GMPr does not conform to the variable naming conventions ((([a-z][a-z0-9_]{2,30})|(_[a-z0-9_]*))$).

The name lon_GMPr does not conform to the variable naming conventions ((([a-z][a-z0-9_]{2,30})|(_[a-z0-9_]*))$).

	sin_lat_gmag = (np.sin(latr) * np.sin(lat_GMPr)
				+ np.cos(latr) * np.cos(lat_GMPr) * np.cos(lonr - lon_GMPr))

Wrong continued indentation (add 13 spaces).

	lon_gmag_y = np.cos(latr) * np.sin(lonr - lon_GMPr)
	lon_gmag_x = (np.cos(latr) * np.sin(lat_GMPr) * np.cos(lonr - lon_GMPr)
				- np.sin(latr) * np.cos(lat_GMPr))

Wrong continued indentation (add 11 spaces).

	lat_gmag = np.arcsin(sin_lat_gmag)
	lat_gmag_geod = np.arctan(np.tan(lat_gmag) / (1. - ellip.epssq))

	B_r = -2. * B_0 * sin_lat_gmag / (grad / ellip.re)**3

The name B_r does not conform to the variable naming conventions ((([a-z][a-z0-9_]{2,30})|(_[a-z0-9_]*))$).

	B_th = -B_0 * np.cos(lat_gmag) / (grad / ellip.re)**3

The name B_th does not conform to the variable naming conventions ((([a-z][a-z0-9_]{2,30})|(_[a-z0-9_]*))$).

	logging.debug("B_r: %s, B_th: %s, dip lat: %s",
			-B_r, B_th, np.degrees(np.arctan(0.5 * B_r / B_th)))

Wrong continued indentation (add 12 spaces).

	lon_gmag = np.arctan2(lon_gmag_y, lon_gmag_x)
	logging.debug("centered dipole coordinates: "
			"lat_gmag: %s, lon_gmag: %s",

Wrong continued indentation (add 12 spaces).

			np.degrees(lat_gmag), np.degrees(lon_gmag))

Wrong continued indentation (add 12 spaces).

	logging.debug("lat_gmag_geod: %s, lat_GMPr: %s",
			np.degrees(lat_gmag_geod), np.degrees(lat_GMPr))

Wrong continued indentation (add 12 spaces).

	if centered_dipole:
		return (np.degrees(lat_gmag), np.degrees(lon_gmag))

	# eccentric dipole coordinates (shifted dipole)
	phi_ed = np.arctan2((grad * np.cos(lat_gmag) * np.sin(lon_gmag) - dY),
			(grad * np.cos(lat_gmag) * np.cos(lon_gmag) - dX))

Wrong continued indentation (add 18 spaces).

	theta_ed = np.arctan2((grad * np.cos(lat_gmag) * np.cos(lon_gmag) - dX),
			((grad * np.sin(lat_gmag) - dZ) * np.cos(phi_ed)))

Wrong continued indentation (add 20 spaces).

	lat_ed = 0.5 * np.pi - theta_ed
	lat_ed_geod = np.arctan(np.tan(lat_ed) / (1. - ellip.epssq))
	logging.debug("lats ed: %s", np.degrees(lat_ed))
	logging.debug("lats ed geod: %s", np.degrees(lat_ed_geod))
	logging.debug("phis ed: %s", np.degrees(phi_ed))
	return (np.degrees(lat_ed_geod), np.degrees(phi_ed))