Repository: SHTOOLS/SHTOOLS (language: python)

Issue:
updates for v4.10.2




**Reminders**

- [x] Base all changes on the `develop` branch: the `master` branch is used only when releasing new versions.
- [x] Run `make check` to ensure that the python code follows standard formatting conventions.
- [x] If adding new features, update the docstring to provide all information that is required to use the feature.



Discussion / hints:

Fix patch:
diff --git a/AUTHORS.md b/AUTHORS.md
index 6073a1c24..caed08007 100644
--- a/AUTHORS.md
+++ b/AUTHORS.md
@@ -2,9 +2,9 @@
 
 ## Maintainers
 * Mark Wieczorek / [MarkWieczorek](https://github.com/MarkWieczorek)
-* Matthias Meschede / [MMesch](https://github.com/MMesch)
 
 ## Contributors
+* Matthias Meschede / [MMesch](https://github.com/MMesch)
 * Elliott Sales de Andrade / [QuLogic](https://github.com/QuLogic)
 * Ilya Oshchepkov / [ioshchepkov](https://github.com/ioshchepkov)
 * [xoviat](https://github.com/xoviat)
@@ -20,4 +20,8 @@
 * Erik Schnetter / [eschnett](https://github.com/eschnett)
 * Martin Reinecke / [mreineck](https://github.com/mreineck)
 * Juan Sierra / [jdsierral](https://github.com/jdsierral)
-* Pierre-Louis Phan / [pl-phan](https://github.com/pl-phan)
+* Jonathon Kohler / [kohlerjl](https://github.com/kohlerjl)
+* Pedro Lobo / [lobophf](https://github.com/lobophf)
+* Pierre-Louis Phan / [pl-phan](https://github.com/pl-phan), [pphan-imcce](https://github.com/pphan-imcce)
+* Dominik Kutra / [k-dominik](https://github.com/k-dominik)
+* Orion Poplawski / [opoplawski](https://github.com/opoplawski)
diff --git a/binder/environment.yml b/binder/environment.yml
index a2176bac0..1b4c59b2a 100644
--- a/binder/environment.yml
+++ b/binder/environment.yml
@@ -2,6 +2,6 @@ name: pyshtools
 channels:
     - conda-forge
 dependencies:
-    - python=3.10
+    - python=3.11
     - pyshtools>=4.10.0
     - palettable>=3.3
diff --git a/docs/Gemfile.lock b/docs/Gemfile.lock
index 9701c004e..9f2070e8a 100644
--- a/docs/Gemfile.lock
+++ b/docs/Gemfile.lock
@@ -86,7 +86,7 @@ GEM
       activesupport (>= 2)
       nokogiri (>= 1.4)
     http_parser.rb (0.8.0)
-    i18n (1.12.0)
+    i18n (1.13.0)
       concurrent-ruby (~> 1.0)
     jekyll (3.9.3)
       addressable (~> 2.4)
diff --git a/docs/_config.yml b/docs/_config.yml
index f382e180e..d357a860b 100755
--- a/docs/_config.yml
+++ b/docs/_config.yml
@@ -40,7 +40,7 @@ exclude:
 
 feedback_subject_line: SHTOOLS
 
-feedback_email: mark.wieczorek@oca.eu
+feedback_email: mark.wieczorek@ipgp.fr
 # used as a contact email for the Feedback link in the top navigation bar
 
 # feedback_disable: true
diff --git a/docs/css/printstyles.css b/docs/css/printstyles.css
index 9d3e2292f..8d41d7ef1 100755
--- a/docs/css/printstyles.css
+++ b/docs/css/printstyles.css
@@ -121,7 +121,7 @@ body.frontmatter {counter-reset: page 1}
         content: counter(page, lower-roman);
     }
     @bottom-left {
-        content: "mark.wieczorek@oca.eu";   }
+        content: "mark.wieczorek@ipgp.fr";   }
 }
 
 body.first_page {counter-reset: page 1}
diff --git a/docs/pages/fortran/fdoc/shtools.md b/docs/pages/fortran/fdoc/shtools.md
index 43f069c8a..98841614b 100644
--- a/docs/pages/fortran/fdoc/shtools.md
+++ b/docs/pages/fortran/fdoc/shtools.md
@@ -9,14 +9,14 @@ toc: false
 editdoc: fdoc
 ---
 
-SHTOOLS is an archive of Fortran 95 and Python software that can be used to perform spherical harmonic transforms and reconstructions, rotations of data expressed in spherical harmonics, and multitaper spectral analyses on the sphere.
+SHTOOLS is a Fortran-95 library that can be used for spherical harmonic transforms, multitaper spectral analyses, expansions of gridded data into Slepian basis functions, and standard operations on global gravitational and magnetic field data.
 
 ## Features
 
 * Supports all standard normalizations and phase conventions of the spherical harmonic functions.
 * Use of both regularly sampled geographic grids and grids appropriate for Gauss-Legendre quadrature.
 * Spherical harmonic transforms proven to be accurate up to about degree 2800.
-* Perform localized multitaper spectral analyses, or expand functions in terms of localized Slepian bases.
+* Perform localized multitaper spectral analyses, or expand gridded data in terms of localized Slepian basis functions.
 * Perform basic operations on global gravity and magnetic field data.
 * OpenMP compatible and OpenMP thread-safe versions of the Fortran routines.
 
diff --git a/docs/pages/fortran/fortran-contributors.md b/docs/pages/fortran/fortran-contributors.md
index 42ba8e1d2..28de9c415 100644
--- a/docs/pages/fortran/fortran-contributors.md
+++ b/docs/pages/fortran/fortran-contributors.md
@@ -10,9 +10,9 @@ folder: fortran
 
 ## Maintainers
 * Mark Wieczorek / [MarkWieczorek](https://github.com/MarkWieczorek)
-* Matthias Meschede / [MMesch](https://github.com/MMesch)
 
 ## Contributors
+* Matthias Meschede / [MMesch](https://github.com/MMesch)
 * Elliott Sales de Andrade / [QuLogic](https://github.com/QuLogic)
 * Ilya Oshchepkov / [ioshchepkov](https://github.com/ioshchepkov)
 * [xoviat](https://github.com/xoviat)
@@ -27,4 +27,9 @@ folder: fortran
 * Tristan Brugere
 * Erik Schnetter / [eschnett](https://github.com/eschnett)
 * Martin Reinecke / [mreineck](https://github.com/mreineck)
-* Juan Sierra / [jdsierral](https://github.com/jdsierral)
\ No newline at end of file
+* Juan Sierra / [jdsierral](https://github.com/jdsierral)
+* Jonathon Kohler / [kohlerjl](https://github.com/kohlerjl)
+* Pedro Lobo / [lobophf](https://github.com/lobophf)
+* Pierre-Louis Phan / [pl-phan](https://github.com/pl-phan), [pphan-imcce](https://github.com/pphan-imcce)
+* Dominik Kutra / [k-dominik](https://github.com/k-dominik)
+* Orion Poplawski / [opoplawski](https://github.com/opoplawski)
diff --git a/docs/pages/fortran/fortran-release-notes-v4.md b/docs/pages/fortran/fortran-release-notes-v4.md
index 9edf7a80a..ad2bd3192 100644
--- a/docs/pages/fortran/fortran-release-notes-v4.md
+++ b/docs/pages/fortran/fortran-release-notes-v4.md
@@ -7,6 +7,29 @@ summary:
 toc: true
 folder: fortran
 ---
+## Version 4.10.2
+
+**Bug fixes and minor enhancements**
+
+* Add the optional parameter `lmax` to `pysh.SHWindow.to_shgrid()`.
+* Replace `libtool` by `ar` and `ranlib` in the fortran makefile.
+* Replace `PWD` by `CURDIR` in the main makefile.
+* Add method `.change_units` to the `SHMagCoeffs` class.
+* Copy all routines from the `shtools` module to `backends.shtools` (the top-level shtools module will be deprecated in v4.11).
+* Use `pkg_resources` instead of `setuptools.version.pkg_resources` in `setup.py`.
+* Fix a couple bugs where `_np.int_` was mistakenly `_np.int`.
+* Move unix man pages from section 1 to section 3.
+* Modify fortran test programs to accept command line arguments, such as the location of the example data files and program input files.
+* Remove `pypandoc` as a dependency in `setup.py` and don't convert the readme to reST format for pypi.
+* Add a Ganymede gravity model to the `datasets` module.
+* Minor updates to documentation and refactoring of the project `README`.
+
+**Future deprecation**
+
+The module `pyshtools.shtools` will be deprecated in the v4.11 release. This module represents 1 of 2 possible backends for pyshtools, and has been located at `pyshtools.backends.shtools` since version 4.9. Unless explicitly required, the user should avoid using the `backends` modules directly, and should instead call the routines that are located in the top level modules such as `pyshtools.expand` and `pyshtools.rotate`. Setting the backend by use of the routine `pyshtools.backends.selected_preferred_backend()` determines which backend to use when calling the routines in these top level modules.
+
+M. A. Wieczorek, M. Meschede, T. Brugere, A. Corbin, A. Hattori, A. Kalinin, J. Kohler, D. Kutra, K. Leinweber, P. Lobo, I. Oshchepkov, P.-L. Phan, O. Poplawski, M. Reinecke, E. Sales de Andrade, E. Schnetter, S. Schröder, J. Sierra, A. Vasishta, A. Walker, xoviat, B. Xu (2023). SHTOOLS: Version 4.10.2, Zenodo, doi:[10.5281/zenodo.592762](https://doi.org/10.5281/zenodo.592762)
+
 ## Version 4.10.1
 
 **Bug fixes and minor enhancements**
diff --git a/docs/pages/mydoc/contributors.md b/docs/pages/mydoc/contributors.md
index 9b2a128dc..33d1ae030 100644
--- a/docs/pages/mydoc/contributors.md
+++ b/docs/pages/mydoc/contributors.md
@@ -10,16 +10,16 @@ folder: mydoc
 
 ## Maintainers
 * Mark Wieczorek / [MarkWieczorek](https://github.com/MarkWieczorek)
-* Matthias Meschede / [MMesch](https://github.com/MMesch)
 
 ## Contributors
+* Matthias Meschede / [MMesch](https://github.com/MMesch)
 * Elliott Sales de Andrade / [QuLogic](https://github.com/QuLogic)
 * Ilya Oshchepkov / [ioshchepkov](https://github.com/ioshchepkov)
 * [xoviat](https://github.com/xoviat)
 * Benda Xu / [heroxbd](https://github.com/heroxbd)
 * Katrin Leinweber / [katrinleinweber](https://github.com/katrinleinweber)
 * Andrew Walker / [andreww](https://github.com/andreww)
-* Stefan Schroder / [Nasuyo](https://github.com/Nasuyo)
+* Stefan Schröder / [Nasuyo](https://github.com/Nasuyo)
 * Akihisa Hattori / [HattoriAkihisa](https://github.com/HattoriAkihisa)
 * Aaryaman Vasishta / [jammm](https://github.com/jammm)
 * Alex Kalinin / [alkalinin](https://github.com/alkalinin)
@@ -27,4 +27,9 @@ folder: mydoc
 * Tristan Brugere
 * Erik Schnetter / [eschnett](https://github.com/eschnett)
 * Martin Reinecke / [mreineck](https://github.com/mreineck)
-* Juan Sierra / [jdsierral](https://github.com/jdsierral)
\ No newline at end of file
+* Juan Sierra / [jdsierral](https://github.com/jdsierral)
+* Jonathon Kohler / [kohlerjl](https://github.com/kohlerjl)
+* Pedro Lobo / [lobophf](https://github.com/lobophf)
+* Pierre-Louis Phan / [pl-phan](https://github.com/pl-phan), [pphan-imcce](https://github.com/pphan-imcce)
+* Dominik Kutra / [k-dominik](https://github.com/k-dominik)
+* Orion Poplawski / [opoplawski](https://github.com/opoplawski)
diff --git a/docs/pages/mydoc/pydoc/pyshlocalizedadmitcorr.md b/docs/pages/mydoc/pydoc/pyshlocalizedadmitcorr.md
index 9fdc424f0..2856ffdcd 100644
--- a/docs/pages/mydoc/pydoc/pyshlocalizedadmitcorr.md
+++ b/docs/pages/mydoc/pydoc/pyshlocalizedadmitcorr.md
@@ -59,7 +59,7 @@ lmax : optional, integer, default = min(lmaxgin, lmaxtin)
 :   The maximum spherical harmonic degree of the input functions corresponding to gilm and tilm.
 
 taper_wt : optional, float, dimension (k), default = -1 (not used)
-:   The weights to be applied to the spectral estimates when calculating the admittance, correlation, and their associated errors. This must sum to unity. The default value specifies that taper weights will not be used.
+:   The weights to be applied to the spectral estimates when calculating the admittance, correlation, and their associated errors. This must sum to unity, and is only implemented for the case of mtdef=1. The default value specifies that taper weights will not be used.
 
 mtdef : optional, integer, default = 1
 :   1 (default): Calculate the multitaper spectral estimates Sgt, Sgg and Stt first, and then use these to calculate the admittance and correlation functions. 2: Calculate admittance and correlation spectra using each individual taper, and then average these to obtain the multitaper admittance and correlation functions.
@@ -73,7 +73,7 @@ SHLocalizedAdmitCorr will calculate the localized admittance and degree correlat
 
 The admittance is defined as Sgt/Stt, where Sgt is the localized cross-power spectrum of two functions G and T expressed in spherical harmonics. The localized degree-correlation spectrum is defined as Sgt/sqrt(Sgg Stt), which can possess values between -1 and 1. Two methods are available for calculating the multitaper admittance and correlation functions. When mtdef is 1 (default), the multitaper estimates and errors of Sgt, Stt, and Sgg are calculated by calls to SHMultiTaperSE and SHMultiTaperCSE, and these results are then used to calculate the final admittance and correlation functions. When mtdef is 2, the admitance and correlation are calculated invidivually for each individual taper, and these results are then averaged.
 
-If the optional parameter k1linsig is specified, and only a single taper is being used, the uncertainty in the admittance function will be calculated by assuming the two sets of coefficients are related by a linear degree-dependent transfer function and that the lack of correlation is a result of uncorrelated noise. 
+If the optional parameter k1linsig is specified, and only a single taper is being used, the uncertainty in the admittance function will be calculated by assuming the two sets of coefficients are related by a linear degree-dependent transfer function and that the lack of correlation is a result of uncorrelated noise.
 
 When mtdef is 1, by default, the multitaper spectral estimates are calculated as an unweighted average of the individual tapered estimates. However, if the optional argument taper_wt is specified, a weighted average will be employed using the weights in this array. Minimum variance optimal weights can be obtained from the routines SHMTVarOpt if the form of the underlying global power spectrum is known. Taper weights can not be used when mtdef is 2
 
@@ -83,9 +83,9 @@ This routine assumes that the input functions and tapers are expressed using geo
 
 Wieczorek, M. A. and F. J. Simons, Minimum-variance multitaper spectral estimation on the sphere, J. Fourier Anal. Appl., 13, doi:10.1007/s00041-006-6904-1, 665-692, 2007.
 
-Simons, F. J., F. A. Dahlen and M. A. Wieczorek, Spatiospectral concentration on the sphere, SIAM Review, 48, 504-536, doi:10.1137/S0036144504445765, 2006. 
+Simons, F. J., F. A. Dahlen and M. A. Wieczorek, Spatiospectral concentration on the sphere, SIAM Review, 48, 504-536, doi:10.1137/S0036144504445765, 2006.
 
-Wieczorek, M. A. and F. J. Simons, Localized spectral analysis on the sphere, 
+Wieczorek, M. A. and F. J. Simons, Localized spectral analysis on the sphere,
 Geophys. J. Int., 162, 655-675, 2005.
 
 Simons, M., S. C. Solomon and B. H. Hager, Localization of gravity and topography: constrains on the tectonics and mantle dynamics of Venus, Geophys. J. Int., 131, 24-44, 1997.
diff --git a/docs/pages/mydoc/pydoc/pyshmtvaropt.md b/docs/pages/mydoc/pydoc/pyshmtvaropt.md
index d99ce815f..9fe1abdc5 100644
--- a/docs/pages/mydoc/pydoc/pyshmtvaropt.md
+++ b/docs/pages/mydoc/pydoc/pyshmtvaropt.md
@@ -9,7 +9,7 @@ toc: false
 editdoc: pydoc
 ---
 
-Calculate the theoretical minimum variance of a localized multitaper spectral estimate and the corresponding optimal weights to apply to each localized spectrum.
+Calculate the theoretical minimum variance of a localized multitaper spectral estimate and the corresponding optimal weights to apply to each localized spectrum. This routine only works using the tapers of the spherical cap concentration problem.
 
 ## Usage
 
diff --git a/docs/pages/mydoc/pydoc/select_preferred_backend.md b/docs/pages/mydoc/pydoc/select_preferred_backend.md
index ce2e3d775..c7ba0ac37 100644
--- a/docs/pages/mydoc/pydoc/select_preferred_backend.md
+++ b/docs/pages/mydoc/pydoc/select_preferred_backend.md
@@ -14,11 +14,11 @@ transforms in pyshtools.
 
 ## Usage
 
-select_preferred_backend_module([backend, nthreads])
+select_preferred_backend([backend, nthreads])
 
 ## Parameters
 
-backend : str, optional, default = 'shtools'
+backend : str, optional, default = 'ducc'
 :   Name of the preferred backend, either 'shtools' or 'ducc'.
 
 nthreads : int, optional, default = 1
diff --git a/docs/pages/mydoc/release-notes-v4.md b/docs/pages/mydoc/release-notes-v4.md
index 8f3c60496..d9d82b3cd 100644
--- a/docs/pages/mydoc/release-notes-v4.md
+++ b/docs/pages/mydoc/release-notes-v4.md
@@ -7,6 +7,29 @@ summary:
 toc: true
 folder: mydoc
 ---
+## Version 4.10.2
+
+**Bug fixes and minor enhancements**
+
+* Add the optional parameter `lmax` to `pysh.SHWindow.to_shgrid()`.
+* Replace `libtool` by `ar` and `ranlib` in the fortran makefile.
+* Replace `PWD` by `CURDIR` in the main makefile.
+* Add method `.change_units` to the `SHMagCoeffs` class.
+* Copy all routines from the `shtools` module to `backends.shtools` (the top-level shtools module will be deprecated in v4.11).
+* Use `pkg_resources` instead of `setuptools.version.pkg_resources` in `setup.py`.
+* Fix a couple bugs where `_np.int_` was mistakenly `_np.int`.
+* Move unix man pages from section 1 to section 3.
+* Modify fortran test programs to accept command line arguments, such as the location of the example data files and program input files.
+* Remove `pypandoc` as a dependency in `setup.py` and don't convert the readme to reST format for pypi.
+* Add a Ganymede gravity model to the `datasets` module.
+* Minor updates to documentation and refactoring of the project `README`.
+
+**Future deprecation**
+
+The module `pyshtools.shtools` will be deprecated in the v4.11 release. This module represents 1 of 2 possible backends for pyshtools, and has been located at `pyshtools.backends.shtools` since version 4.9. Unless explicitly required, the user should avoid using the `backends` modules directly, and should instead call the routines that are located in the top level modules such as `pyshtools.expand` and `pyshtools.rotate`. Setting the backend by use of the routine `pyshtools.backends.selected_preferred_backend()` determines which backend to use when calling the routines in these top level modules.
+
+M. A. Wieczorek, M. Meschede, T. Brugere, A. Corbin, A. Hattori, A. Kalinin, J. Kohler, D. Kutra, K. Leinweber, P. Lobo, I. Oshchepkov, P.-L. Phan, O. Poplawski, M. Reinecke, E. Sales de Andrade, E. Schnetter, S. Schröder, J. Sierra, A. Vasishta, A. Walker, xoviat, B. Xu (2023). SHTOOLS: Version 4.10.2, Zenodo, doi:[10.5281/zenodo.592762](https://doi.org/10.5281/zenodo.592762)
+
 ## Version 4.10.1
 
 **Bug fixes and minor enhancements**
diff --git a/examples/fortran/MarsCrustalThickness/MarsCrustalThickness.f95 b/examples/fortran/MarsCrustalThickness/MarsCrustalThickness.f95
index ef5d2373b..fe387ba53 100644
--- a/examples/fortran/MarsCrustalThickness/MarsCrustalThickness.f95
+++ b/examples/fortran/MarsCrustalThickness/MarsCrustalThickness.f95
@@ -44,12 +44,67 @@ program MarsCrustalThickness
                       n_out, iter, j, r1, lmaxp, lmaxt, filter_type, half, &
                       degmax, sampling
     character(120) :: grav_file, moho_out, thick_grid_out, topo_file, &
-                      misfit_file
+                      misfit_file, infile
 
-    print*,  "rho_crust (kg/m3) > "
-    read(*,*) rho_crust
-    print*, "rho_mantle (kg/m3) > "
-    read(*,*) rho_mantle
+    ! Path to example data files may be passed as first argument, or use a default.
+    if (command_argument_count() > 0) then
+        call get_command_argument(1, infile)
+    else
+        infile = "../../ExampleDataFiles"
+    end if
+
+    grav_file = trim(infile) // "/gmm3_120_sha.tab"
+    topo_file = trim(infile) // "/MarsTopo719.shape"
+
+    ! A data input file may be passed as second argument, or else prompt for required settings.
+    if (command_argument_count() > 1) then
+        call get_command_argument(2, infile)
+        open(unit=20, file=infile, action="read")
+        read(20,*) rho_crust
+        read(20,*) rho_mantle
+        read(20,*) filter_type
+        if (filter_type /= 0 ) then
+            read(20,*) half
+        end if
+        read(20,*) r1
+        read(20,*) degmax
+        read(20,*) t0
+        read(20,*) moho_out
+        read(20,*) interval
+        read(20,*) thick_grid_out
+        read(20,*) misfit_file
+        close(20)
+    else
+        print*,  "rho_crust (kg/m3) > "
+        read(*,*) rho_crust
+        print*, "rho_mantle (kg/m3) > "
+        read(*,*) rho_mantle
+
+        print*, "Input filter type (1) Minimum amplitude, (2) minimum curvature, (0) no filter "
+        read(*,*) filter_type
+        if (filter_type /= 0 ) then
+            print*, "Degree at which the filter is 1/2 "
+            read(*,*) half
+        end if
+
+        print*, "Remove degree 1 topo coefficients from Bouguer Correction? (0:no, 1:yes) > "
+        read(*,*) r1
+
+        print*, "maximum degree to compute Moho relief to >"
+        read(*,*) degmax
+
+        print*, "Minimum assumed crustal thickness (km) > "
+        read(*,*) t0
+
+        print*, "Moho spherical harmonic coeficient output filename > "
+        read(*,*) moho_out
+        print*, "Grid spacing for output crustal thickness map (degrees) > "
+        read(*,*) interval
+        print*, "gridded crustal thickness output filename >"
+        read(*,*) thick_grid_out
+        print*, "Gravity misfit spherical harmonic filename >"
+        read(*,*) misfit_file
+    end if
 
     pi = acos(-1.0_dp)
     delta_max = 5.0_dp
@@ -66,39 +121,12 @@ program MarsCrustalThickness
 
     nmax = 7    ! nmax of Wieczorek and Phillips (1998)
 
-    print*, "Input filter type (1) Minimum amplitude, (2) minimum curvature, (0) no filter "
-    read(*,*) filter_type
-    if (filter_type /= 0 ) then
-        print*, "Degree at which the filter is 1/2 " 
-        read(*,*) half
-    end if
-
-    grav_file = "../../ExampleDataFiles/gmm3_120_sha.tab"
-    topo_file = "../../ExampleDataFiles/MarsTopo719.shape"
-
-    print*, "Remove degree 1 topo coefficients from Bouguer Correction? (0:no, 1:yes) > "
-    read(*,*) r1
-
-    print*, "maximum degree to compute Moho relief to >"
-    read(*,*) degmax
-
     lmax = 2 * degmax ! this partially takes into account aliasing problems. Technically, it should be nmax*degmax
 
     nlat = 2 * lmax + 2
     nlong = 2 * nlat
-    print*, "Minimum assumed crustal thickness (km) > "
-    read(*,*) t0
     t0 = t0 * 1.0e3_dp
 
-    print*, "Moho spherical harmonic coeficient output filename > "
-    read(*,*) moho_out
-    print*, "Grid spacing for output crustal thickness map (degrees) > "
-    read(*,*) interval
-    print*, "gridded crustal thickness output filename >"
-    read(*,*) thick_grid_out
-    print*, "Gravity misfit spherical harmonic filename >"
-    read(*,*) misfit_file
-
     call cpu_time(timein)
 
     allocate(topo_grid(nlat, nlong), stat = astat(1))
diff --git a/examples/fortran/SHCilmPlus/TestCilmPlus.f95 b/examples/fortran/SHCilmPlus/TestCilmPlus.f95
index 66be0f38d..d83ff0d56 100644
--- a/examples/fortran/SHCilmPlus/TestCilmPlus.f95
+++ b/examples/fortran/SHCilmPlus/TestCilmPlus.f95
@@ -76,7 +76,13 @@ program TestCilmPlus
 
     !print*, "Spherical Harmonic file > "
     !read(*,*) infile
-    infile = "../../ExampleDataFiles/MarsTopo719.shape"
+    ! Path to example data files may be passed as first argument, or use a default.
+    if (command_argument_count() > 0) then
+        call get_command_argument(1, infile)
+    else
+        infile = "../../ExampleDataFiles"
+    end if
+    infile = trim(infile) // "/MarsTopo719.shape"
 
     !print*, "Maximum degree to read > "
     !read(*,*) lmax
diff --git a/examples/fortran/SHExpandDH/TestExpandDH.f95 b/examples/fortran/SHExpandDH/TestExpandDH.f95
index a8b792420..93a5dd235 100644
--- a/examples/fortran/SHExpandDH/TestExpandDH.f95
+++ b/examples/fortran/SHExpandDH/TestExpandDH.f95
@@ -14,13 +14,19 @@ Program TestExpandDH
 
     implicit none
     integer(int32), parameter :: maxdeg = 120, maxgrid = 800
-    character(80) :: infile
+    character(120) :: infile
     real(dp) :: cilm(2, maxdeg+1, maxdeg+1), grid(maxgrid, maxgrid), &
                 griddh(maxgrid, maxgrid), cilm2(2, maxdeg+1, maxdeg+1), &
                 interval, error, maxerror
     integer(int32) :: lmax, n, nlat, nlong, lmax2, l, m, i
 
-    infile = "../../ExampleDataFiles/MarsTopo719.shape"
+    ! Path to example data files may be passed as first argument, or use a default.
+    if (command_argument_count() > 0) then
+        call get_command_argument(1, infile)
+    else
+        infile = "../../ExampleDataFiles"
+    end if
+    infile = trim(infile) // "/MarsTopo719.shape"
 
     call SHRead(infile, cilm, lmax)
     print*, "Lmax of data file = ", lmax
diff --git a/examples/fortran/SHExpandLSQ/TestSHExpandLSQ.f95 b/examples/fortran/SHExpandLSQ/TestSHExpandLSQ.f95
index 66be46601..8db9a5eb8 100644
--- a/examples/fortran/SHExpandLSQ/TestSHExpandLSQ.f95
+++ b/examples/fortran/SHExpandLSQ/TestSHExpandLSQ.f95
@@ -22,7 +22,7 @@ program TestSHExpandLSQ
                 x, y, z, pi, maxerror, cilm1(2,degmax+1, degmax+1), &
                 dd(dmax), misfit
     integer(int32) :: nmax, lmax, l, i, seed, lmaxfile
-    character(80) :: infile
+    character(120) :: infile
 
     d = 0.0_dp
     dd = 0.0_dp
@@ -30,13 +30,28 @@ program TestSHExpandLSQ
     lat = 0.0_dp
     pi = acos(-1.0_dp)
 
-    infile = "../../ExampleDataFiles/MarsTopo719.shape"
+    ! Path to example data files may be passed as first argument, or use a default.
+    if (command_argument_count() > 0) then
+        call get_command_argument(1, infile)
+    else
+        infile = "../../ExampleDataFiles"
+    end if
+    infile = trim(infile) // "/MarsTopo719.shape"
 
     call SHRead(infile, cilm, lmaxfile)
     print*, "Maximum degree of spherical harmonic file = ", lmaxfile
 
-    print*, "Number of random data points to use > "
-    read(*,*) nmax
+    ! A data input file may be passed as second argument, or else prompt for required settings.
+    if (command_argument_count() > 1) then
+        call get_command_argument(2, infile)
+        open(unit=20, file=infile, action="read")
+        read(20,*) nmax
+        close(20)
+    else
+        print*, "Number of random data points to use > "
+        read(*,*) nmax
+    end if
+
     lmax = floor(sqrt(dble(nmax)) - 1)
     print*, "Maximum spherical harmonic degree for overdetermined least-squares inversion = ", lmax
 
@@ -85,4 +100,4 @@ program TestSHExpandLSQ
 
     end do
 
-end program TestSHExpandLSQ
\ No newline at end of file
+end program TestSHExpandLSQ
diff --git a/examples/fortran/SHLocalizedAdmitCorr/LocalizedAdmitCorr.f95 b/examples/fortran/SHLocalizedAdmitCorr/LocalizedAdmitCorr.f95
index f0b1d08c2..fb15d3b56 100644
--- a/examples/fortran/SHLocalizedAdmitCorr/LocalizedAdmitCorr.f95
+++ b/examples/fortran/SHLocalizedAdmitCorr/LocalizedAdmitCorr.f95
@@ -14,7 +14,7 @@ Program LocalizedAdmitCorr
 
     implicit none
 
-    character(120) :: topography_file, potential_file, outfile
+    character(120) :: infile, topography_file, potential_file, outfile
     real(dp) :: header(8), mpr,  r0_pot, gm, mass, lat, lon, pi, theta0, &
                 alpha, sn, r
     real(dp), allocatable :: grav(:,:,:), admit(:), corr(:), admit_error(:), &
@@ -28,8 +28,14 @@ Program LocalizedAdmitCorr
 
     degmax = 360
 
-    topography_file = "../../ExampleDataFiles/MarsTopo719.shape"
-    potential_file = "../../ExampleDataFiles/gmm3_120_sha.tab"
+    ! Path to example data files may be passed as first argument, or use a default.
+    if (command_argument_count() > 0) then
+        call get_command_argument(1, infile)
+    else
+        infile = "../../ExampleDataFiles"
+    end if
+    topography_file = trim(infile) // "/MarsTopo719.shape"
+    potential_file = trim(infile) // "/gmm3_120_sha.tab"
 
     !--------------------------------------------------------------------------
     !
@@ -78,42 +84,59 @@ Program LocalizedAdmitCorr
     !
     !--------------------------------------------------------------------------
 
-    print*, "Number of tapers to use > "
-    read(*,*) K
-
-    print*, "Latitude of feature of interest (degrees) > "
-    read(*,*) lat
-    print*, "Longitude (degrees) > "
-    read(*,*) lon
-    print*, "Angular radius of localization window (degrees) > "
-    read(*,*) theta0
-    theta0 = theta0*pi / 180.0_dp
-
-    print*, "Create localization window using"
-    print*, "(1) Desired concentration factor, alpha"
-    print*, "(2) Desired (approximate) Shannon number; (Lwin+1) Theta0 / pi)"
-    print*, "(3) Desired Spectral bandwidth (Lwin)"
-    read(*,*) option1
-
-    if (option1==1) then
-        print*, "Input desired concentration factor of the Kth taper > "
-        read(*,*) alpha
-        lwin = SHFindLWin(theta0, 0, alpha, K)
-        print*, "Corresponding spherical harmonic bandwidth = ", lwin
-        print*, "Corresponding approximate Shannon number = ", &
-                (lwin+1) * theta0 / pi
-    else if (option1==2) then
-        print*, "Input Shannon number > "
-        read(*,*) sn
-        lwin = nint(sn * pi / theta0) - 1
-        print*, "Corresponding spherical harmonic bandwidth = ", lwin
-    else if (option1==3) then
-        print*, "Input Lwin > "
-        read(*,*) lwin
-        print*, "Corresponding approximate Shannon number = ", &
-                (lwin+1) * theta0 / pi
+    ! A data input file may be passed as second argument, or else prompt for required settings.
+    if (command_argument_count() > 1) then
+        call get_command_argument(2, infile)
+        open(unit=20, file=infile, action="read")
+        read(20,*) K
+        read(20,*) lat
+        read(20,*) lon
+        read(20,*) theta0
+        read(20,*) option1
+        if (option1==1) then
+            read(20,*) alpha
+        else if (option1==2) then
+            read(20,*) sn
+        else if (option1==3) then
+            read(20,*) lwin
+        else
+            stop
+        end if
+        read(20,*) outfile
+        close(20)
     else
-        stop
+        print*, "Number of tapers to use > "
+        read(*,*) K
+
+        print*, "Latitude of feature of interest (degrees) > "
+        read(*,*) lat
+        print*, "Longitude (degrees) > "
+        read(*,*) lon
+        print*, "Angular radius of localization window (degrees) > "
+        read(*,*) theta0
+        theta0 = theta0*pi / 180.0_dp
+
+        print*, "Create localization window using"
+        print*, "(1) Desired concentration factor, alpha"
+        print*, "(2) Desired (approximate) Shannon number; (Lwin+1) Theta0 / pi)"
+        print*, "(3) Desired Spectral bandwidth (Lwin)"
+        read(*,*) option1
+
+        if (option1==1) then
+            print*, "Input desired concentration factor of the Kth taper > "
+            read(*,*) alpha
+        else if (option1==2) then
+            print*, "Input Shannon number > "
+            read(*,*) sn
+        else if (option1==3) then
+            print*, "Input Lwin > "
+            read(*,*) lwin
+        else
+            stop
+        end if
+
+        print*, "Name of output admittance and correlation file > "
+        read(*,*) outfile
     end if
 
     allocate(tapers(lwin+1, (lwin+1)**2), stat = astat(1))
@@ -137,8 +160,21 @@ Program LocalizedAdmitCorr
 
     print*, "Gravity field evaluated at R = (km) ", r0_pot / 1.e3_dp
 
-    print*, "Name of output admittance and correlation file > "
-    read(*,*) outfile
+    if (option1==1) then
+        lwin = SHFindLWin(theta0, 0, alpha, K)
+        print*, "Corresponding spherical harmonic bandwidth = ", lwin
+        print*, "Corresponding approximate Shannon number = ", &
+                (lwin+1) * theta0 / pi
+    else if (option1==2) then
+        lwin = nint(sn * pi / theta0) - 1
+        print*, "Corresponding spherical harmonic bandwidth = ", lwin
+    else if (option1==3) then
+        print*, "Corresponding approximate Shannon number = ", &
+                (lwin+1) * theta0 / pi
+    else
+        stop
+    end if
+
     open(12, file=outfile)
 
     !--------------------------------------------------------------------------
diff --git a/examples/fortran/SHMag/SHMag.f95 b/examples/fortran/SHMag/SHMag.f95
index 08a3222f3..92926a9de 100644
--- a/examples/fortran/SHMag/SHMag.f95
+++ b/examples/fortran/SHMag/SHMag.f95
@@ -25,7 +25,13 @@ program SHMag
                              total(:,:), pot(:,:)
     integer(int32) :: lmax, lmaxp, n, n_out, nlong, nlat, i, j, astat(6), sampling
 
-    infile = "../../ExampleDataFiles/FSU_mars90.sh"
+    ! Path to example data files may be passed as first argument, or use a default.
+    if (command_argument_count() > 0) then
+        call get_command_argument(1, infile)
+    else
+        infile = "../../ExampleDataFiles"
+    end if
+    infile = trim(infile) // "/FSU_mars90.sh"
 
     f = 1.0_dp / 169.864881_dp    ! Mars flattening = (R_eq - R_p)/R_eq
     mpr = 3389.508e3_dp           ! Mean radius of mars
diff --git a/examples/fortran/SHRotate/TestSHRotate.f95 b/examples/fortran/SHRotate/TestSHRotate.f95
index b02d1c115..b5d8cb34e 100644
--- a/examples/fortran/SHRotate/TestSHRotate.f95
+++ b/examples/fortran/SHRotate/TestSHRotate.f95
@@ -27,26 +27,55 @@ program TestSHRotate
         stop
     end if
 
-    infile = "../../ExampleDataFiles/MarsTopo719.shape"
+    ! Path to example data files may be passed as first argument, or use a default.
+    if (command_argument_count() > 0) then
+        call get_command_argument(1, infile)
+    else
+        infile = "../../ExampleDataFiles"
+    end if
+    infile = trim(infile) // "/MarsTopo719.shape"
     print*, "Input file = ", infile
     call SHRead(infile, cilm1, lmax)
     print*, "Lmax of input file = ", lmax
 
-    print*, "Maximum degree to be used with rotations > "
-    read(*,*) lmax
-
-    print*, "Input Euler Angles"
-    print*, "Alpha (deg) = "
-    read(*,*) alpha
-    print*, "Beta (deg) = "
-    read(*,*) beta
-    print*, "Gamma (deg) = "
-    read(*,*) gamma
-
-    print*, "Do these correspond to "
-    print*, "(1) Rotation of the coordinate system without rotation of the physical body, or"
-    print*, "(2) Rotation of the physical body without rotation of the coordinate system."
-    read(*,*) bodycoord
+    ! A data input file may be passed as second argument, or else prompt for required settings.
+    if (command_argument_count() > 1) then
+        call get_command_argument(2, infile)
+        open(unit=20, file=infile, action="read")
+        read(20,*) lmax
+        read(20,*) alpha
+        read(20,*) beta
+        read(20,*) gamma
+        read(20,*) bodycoord
+        read(20,*) outsh
+        read(20,*) interval
+        read(20,*) outgrid
+        close(20)
+    else
+        print*, "Maximum degree to be used with rotations > "
+        read(*,*) lmax
+
+        print*, "Input Euler Angles"
+        print*, "Alpha (deg) = "
+        read(*,*) alpha
+        print*, "Beta (deg) = "
+        read(*,*) beta
+        print*, "Gamma (deg) = "
+        read(*,*) gamma
+
+        print*, "Do these correspond to "
+        print*, "(1) Rotation of the coordinate system without rotation of the physical body, or"
+        print*, "(2) Rotation of the physical body without rotation of the coordinate system."
+        read(*,*) bodycoord
+
+        print*, "Output file name of rotated spherical harmonics > "
+        read(*,*) outsh
+
+        print*, "Grid spacing for output grid (deg) > "
+        read(*,*) interval
+        print*, "File name of gridded rotated field > "
+        read(*,*) outgrid
+    end if
 
     if (bodycoord ==1) then
         x(1) = alpha ; x(2) = beta ; x(3) = gamma
@@ -59,14 +88,6 @@ program TestSHRotate
 
     x = x * pi / 180.0_dp
 
-    print*, "Output file name of rotated spherical harmonics > "
-    read(*,*) outsh
-
-    print*, "Grid spacing for output grid (deg) > "
-    read(*,*) interval
-    print*, "File name of gridded rotated field > "
-    read(*,*) outgrid
-
     allocate(cilm2(2, lmax+1, lmax+1), stat = astat(2))
     allocate(dj(lmax+1, lmax+1, lmax+1), stat = astat(3))
     nlat = 180.0_dp / interval + 1
diff --git a/man/man3/batohilm.3 b/man/man3/batohilm.3
index 74c39cdc6..0d9b6237e 100644
--- a/man/man3/batohilm.3
+++ b/man/man3/batohilm.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "batohilm" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "batohilm" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH BAtoHilm
 .PP
diff --git a/man/man3/batohilmrhoh.3 b/man/man3/batohilmrhoh.3
index 398de579f..4ccb339fe 100644
--- a/man/man3/batohilmrhoh.3
+++ b/man/man3/batohilmrhoh.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "batohilmrhoh" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "batohilmrhoh" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH BAtoHilmRhoH
 .PP
diff --git a/man/man3/cilmminus.3 b/man/man3/cilmminus.3
index 51e2684d3..98a70794b 100644
--- a/man/man3/cilmminus.3
+++ b/man/man3/cilmminus.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "cilmminus" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "cilmminus" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH CilmMinus
 .PP
diff --git a/man/man3/cilmminusrhoh.3 b/man/man3/cilmminusrhoh.3
index df820b9f4..1d58aa53d 100644
--- a/man/man3/cilmminusrhoh.3
+++ b/man/man3/cilmminusrhoh.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "cilmminusrhoh" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "cilmminusrhoh" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH CilmMinusRhoH
 .PP
diff --git a/man/man3/cilmplus.3 b/man/man3/cilmplus.3
index b3ff195b7..80d30c205 100644
--- a/man/man3/cilmplus.3
+++ b/man/man3/cilmplus.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "cilmplus" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "cilmplus" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH CilmPlus
 .PP
diff --git a/man/man3/cilmplusrhoh.3 b/man/man3/cilmplusrhoh.3
index bba9f1b0d..5afeb6fdd 100644
--- a/man/man3/cilmplusrhoh.3
+++ b/man/man3/cilmplusrhoh.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "cilmplusrhoh" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "cilmplusrhoh" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH CilmPlusRhoH
 .PP
diff --git a/man/man3/computedg82.3 b/man/man3/computedg82.3
index 28592bcde..b6a7411c1 100644
--- a/man/man3/computedg82.3
+++ b/man/man3/computedg82.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "computedg82" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "computedg82" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH ComputeDG82
 .PP
diff --git a/man/man3/computedm.3 b/man/man3/computedm.3
index a6f51cb4f..81dd704f7 100644
--- a/man/man3/computedm.3
+++ b/man/man3/computedm.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "computedm" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "computedm" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH ComputeDM
 .PP
diff --git a/man/man3/computedmap.3 b/man/man3/computedmap.3
index 6be64674d..8f4124be1 100644
--- a/man/man3/computedmap.3
+++ b/man/man3/computedmap.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "computedmap" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "computedmap" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH ComputeDMap
 .PP
diff --git a/man/man3/curve2mask.3 b/man/man3/curve2mask.3
index d31cafa3d..ef168d162 100644
--- a/man/man3/curve2mask.3
+++ b/man/man3/curve2mask.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "curve2mask" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "curve2mask" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH Curve2Mask
 .PP
diff --git a/man/man3/dhaj.3 b/man/man3/dhaj.3
index 103adfa04..5d49e2b76 100644
--- a/man/man3/dhaj.3
+++ b/man/man3/dhaj.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "dhaj" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "dhaj" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH DHaj
 .PP
diff --git a/man/man3/djpi2.3 b/man/man3/djpi2.3
index dd8edeaa1..b19b53a87 100644
--- a/man/man3/djpi2.3
+++ b/man/man3/djpi2.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "djpi2" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "djpi2" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH djpi2
 .PP
diff --git a/man/man3/downcontfilterma.3 b/man/man3/downcontfilterma.3
index cf9718fd6..0b0159d7a 100644
--- a/man/man3/downcontfilterma.3
+++ b/man/man3/downcontfilterma.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "downcontfilterma" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "downcontfilterma" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH DownContFilterMA
 .PP
diff --git a/man/man3/downcontfiltermc.3 b/man/man3/downcontfiltermc.3
index 80dafa6ce..ee75e0ca6 100644
--- a/man/man3/downcontfiltermc.3
+++ b/man/man3/downcontfiltermc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "downcontfiltermc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "downcontfiltermc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH DownContFilterMC
 .PP
diff --git a/man/man3/eigvalsym.3 b/man/man3/eigvalsym.3
index 27e23771f..b8d2c7187 100644
--- a/man/man3/eigvalsym.3
+++ b/man/man3/eigvalsym.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "eigvalsym" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "eigvalsym" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH EigValSym
 .PP
diff --git a/man/man3/eigvalvecsym.3 b/man/man3/eigvalvecsym.3
index 443e1f637..42eae61a3 100644
--- a/man/man3/eigvalvecsym.3
+++ b/man/man3/eigvalvecsym.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "eigvalvecsym" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "eigvalvecsym" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH EigValVecSym
 .PP
diff --git a/man/man3/eigvalvecsymtri.3 b/man/man3/eigvalvecsymtri.3
index b5b390774..5a8a60073 100644
--- a/man/man3/eigvalvecsymtri.3
+++ b/man/man3/eigvalvecsymtri.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "eigvalvecsymtri" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "eigvalvecsymtri" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH EigValVecSymTri
 .PP
diff --git a/man/man3/glqgridcoord.3 b/man/man3/glqgridcoord.3
index 49a3c0f3c..64924b879 100644
--- a/man/man3/glqgridcoord.3
+++ b/man/man3/glqgridcoord.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "glqgridcoord" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "glqgridcoord" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH GLQGridCoord
 .PP
diff --git a/man/man3/makecirclecoord.3 b/man/man3/makecirclecoord.3
index 7c116443d..2524beea3 100644
--- a/man/man3/makecirclecoord.3
+++ b/man/man3/makecirclecoord.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makecirclecoord" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makecirclecoord" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeCircleCoord
 .PP
diff --git a/man/man3/makeellipsecoord.3 b/man/man3/makeellipsecoord.3
index 2165af742..b45c196df 100644
--- a/man/man3/makeellipsecoord.3
+++ b/man/man3/makeellipsecoord.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makeellipsecoord" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makeellipsecoord" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeEllipseCoord
 .PP
diff --git a/man/man3/makegeoidgrid.3 b/man/man3/makegeoidgrid.3
index db6d4c24e..364a2a8f1 100644
--- a/man/man3/makegeoidgrid.3
+++ b/man/man3/makegeoidgrid.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegeoidgrid" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegeoidgrid" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGeoidGrid
 .PP
diff --git a/man/man3/makegradientdh.3 b/man/man3/makegradientdh.3
index 706b1d83e..02ed5c1c6 100644
--- a/man/man3/makegradientdh.3
+++ b/man/man3/makegradientdh.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegradientdh" "3" "2021-09-03" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegradientdh" "1" "2022-05-02" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGradientDH
 .PP
diff --git a/man/man3/makegravgradgriddh.3 b/man/man3/makegravgradgriddh.3
index 4e8ac6592..4df03e191 100644
--- a/man/man3/makegravgradgriddh.3
+++ b/man/man3/makegravgradgriddh.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegravgradgriddh" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegravgradgriddh" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGravGradGridDH
 .PP
diff --git a/man/man3/makegravgriddh.3 b/man/man3/makegravgriddh.3
index ba9b005e2..c02d34a73 100644
--- a/man/man3/makegravgriddh.3
+++ b/man/man3/makegravgriddh.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegravgriddh" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegravgriddh" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGravGridDH
 .PP
diff --git a/man/man3/makegravgridpoint.3 b/man/man3/makegravgridpoint.3
index a7dc45445..4a00556c2 100644
--- a/man/man3/makegravgridpoint.3
+++ b/man/man3/makegravgridpoint.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegravgridpoint" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegravgridpoint" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGravGridPoint
 .PP
diff --git a/man/man3/makegrid2d.3 b/man/man3/makegrid2d.3
index 1b5a6430d..4c8f7cddd 100644
--- a/man/man3/makegrid2d.3
+++ b/man/man3/makegrid2d.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegrid2d" "3" "2021-10-22" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegrid2d" "1" "2021-10-22" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGrid2D
 .PP
diff --git a/man/man3/makegriddh.3 b/man/man3/makegriddh.3
index 796fea073..c80e24370 100644
--- a/man/man3/makegriddh.3
+++ b/man/man3/makegriddh.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegriddh" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegriddh" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGridDH
 .PP
diff --git a/man/man3/makegriddhc.3 b/man/man3/makegriddhc.3
index 0df8427de..3930faf32 100644
--- a/man/man3/makegriddhc.3
+++ b/man/man3/makegriddhc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegriddhc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegriddhc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGridDHC
 .PP
diff --git a/man/man3/makegridglq.3 b/man/man3/makegridglq.3
index b6a595157..aeabab8e1 100644
--- a/man/man3/makegridglq.3
+++ b/man/man3/makegridglq.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegridglq" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegridglq" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGridGLQ
 .PP
diff --git a/man/man3/makegridglqc.3 b/man/man3/makegridglqc.3
index b0a1b88c6..b4495c4b8 100644
--- a/man/man3/makegridglqc.3
+++ b/man/man3/makegridglqc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegridglqc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegridglqc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGridGLQC
 .PP
diff --git a/man/man3/makegridpoint.3 b/man/man3/makegridpoint.3
index dab854d08..3de4d6262 100644
--- a/man/man3/makegridpoint.3
+++ b/man/man3/makegridpoint.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegridpoint" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegridpoint" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGridPoint
 .PP
diff --git a/man/man3/makegridpointc.3 b/man/man3/makegridpointc.3
index 5f2b5de69..7c9265155 100644
--- a/man/man3/makegridpointc.3
+++ b/man/man3/makegridpointc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makegridpointc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makegridpointc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeGridPointC
 .PP
diff --git a/man/man3/makemaggradgriddh.3 b/man/man3/makemaggradgriddh.3
index bdbe66cbd..66c93a905 100644
--- a/man/man3/makemaggradgriddh.3
+++ b/man/man3/makemaggradgriddh.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makemaggradgriddh" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makemaggradgriddh" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeMagGradGridDH
 .PP
diff --git a/man/man3/makemaggriddh.3 b/man/man3/makemaggriddh.3
index 34a167379..990c8cf8f 100644
--- a/man/man3/makemaggriddh.3
+++ b/man/man3/makemaggriddh.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makemaggriddh" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makemaggriddh" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeMagGridDH
 .PP
diff --git a/man/man3/makemaggridpoint.3 b/man/man3/makemaggridpoint.3
index 2b255f11c..4c06d6f1e 100644
--- a/man/man3/makemaggridpoint.3
+++ b/man/man3/makemaggridpoint.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "makemaggridpoint" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "makemaggridpoint" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH MakeMagGridPoint
 .PP
diff --git a/man/man3/normalgravity.3 b/man/man3/normalgravity.3
index 1a0f29d75..eafcd3ce4 100644
--- a/man/man3/normalgravity.3
+++ b/man/man3/normalgravity.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "normalgravity" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "normalgravity" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH NormalGravity
 .PP
diff --git a/man/man3/planetsconstants.3 b/man/man3/planetsconstants.3
index 6086e10fe..bf7a96684 100644
--- a/man/man3/planetsconstants.3
+++ b/man/man3/planetsconstants.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "planetsconstants" "3" "2020-03-09" "Fortran 95" "SHTOOLS 4.10"
+.TH "planetsconstants" "1" "2020-03-09" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlanetsConstants
 .PP
diff --git a/man/man3/plbar.3 b/man/man3/plbar.3
index 6ced93fe4..6d5c61b6e 100644
--- a/man/man3/plbar.3
+++ b/man/man3/plbar.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plbar" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plbar" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlBar
 .PP
diff --git a/man/man3/plbar_d1.3 b/man/man3/plbar_d1.3
index 326b1c857..6e197590e 100644
--- a/man/man3/plbar_d1.3
+++ b/man/man3/plbar_d1.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plbar_d1" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plbar_d1" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlBar_d1
 .PP
diff --git a/man/man3/plegendre.3 b/man/man3/plegendre.3
index 4da4cb2b7..a8312851e 100644
--- a/man/man3/plegendre.3
+++ b/man/man3/plegendre.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plegendre" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plegendre" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PLegendre
 .PP
diff --git a/man/man3/plegendre_d1.3 b/man/man3/plegendre_d1.3
index 734e35ca9..36f33198a 100644
--- a/man/man3/plegendre_d1.3
+++ b/man/man3/plegendre_d1.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plegendre_d1" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plegendre_d1" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PLegendre_d1
 .PP
diff --git a/man/man3/plegendrea.3 b/man/man3/plegendrea.3
index b7d79df5b..ea84bb06c 100644
--- a/man/man3/plegendrea.3
+++ b/man/man3/plegendrea.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plegendrea" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plegendrea" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PLegendreA
 .PP
diff --git a/man/man3/plegendrea_d1.3 b/man/man3/plegendrea_d1.3
index e6c49b1d5..4f7a298b0 100644
--- a/man/man3/plegendrea_d1.3
+++ b/man/man3/plegendrea_d1.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plegendrea_d1" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plegendrea_d1" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PLegendreA_d1
 .PP
diff --git a/man/man3/plmbar.3 b/man/man3/plmbar.3
index d041ab0db..cae88b53d 100644
--- a/man/man3/plmbar.3
+++ b/man/man3/plmbar.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plmbar" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plmbar" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlmBar
 .PP
diff --git a/man/man3/plmbar_d1.3 b/man/man3/plmbar_d1.3
index 74d13406b..2416048e2 100644
--- a/man/man3/plmbar_d1.3
+++ b/man/man3/plmbar_d1.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plmbar_d1" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plmbar_d1" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlmBar_d1
 .PP
diff --git a/man/man3/plmindex.3 b/man/man3/plmindex.3
index 7f4cc122b..9da5cd1ac 100644
--- a/man/man3/plmindex.3
+++ b/man/man3/plmindex.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plmindex" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plmindex" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlmIndex
 .PP
diff --git a/man/man3/plmon.3 b/man/man3/plmon.3
index dc5831430..6b41d3742 100644
--- a/man/man3/plmon.3
+++ b/man/man3/plmon.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plmon" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plmon" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlmON
 .PP
diff --git a/man/man3/plmon_d1.3 b/man/man3/plmon_d1.3
index 16551ba75..50f2b9f4d 100644
--- a/man/man3/plmon_d1.3
+++ b/man/man3/plmon_d1.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plmon_d1" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plmon_d1" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlmON_d1
 .PP
diff --git a/man/man3/plmschmidt.3 b/man/man3/plmschmidt.3
index 7932304e8..907df4990 100644
--- a/man/man3/plmschmidt.3
+++ b/man/man3/plmschmidt.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plmschmidt" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plmschmidt" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlmSchmidt
 .PP
diff --git a/man/man3/plmschmidt_d1.3 b/man/man3/plmschmidt_d1.3
index 8cfc565c5..4485a0c82 100644
--- a/man/man3/plmschmidt_d1.3
+++ b/man/man3/plmschmidt_d1.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plmschmidt_d1" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plmschmidt_d1" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlmSchmidt_d1
 .PP
diff --git a/man/man3/plon.3 b/man/man3/plon.3
index 0499eac92..b74dc3ee8 100644
--- a/man/man3/plon.3
+++ b/man/man3/plon.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plon" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plon" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlON
 .PP
diff --git a/man/man3/plon_d1.3 b/man/man3/plon_d1.3
index 38ee89921..e33d1975d 100644
--- a/man/man3/plon_d1.3
+++ b/man/man3/plon_d1.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plon_d1" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plon_d1" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlON_d1
 .PP
diff --git a/man/man3/plschmidt.3 b/man/man3/plschmidt.3
index 5b788c782..fe107554a 100644
--- a/man/man3/plschmidt.3
+++ b/man/man3/plschmidt.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plschmidt" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plschmidt" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlSchmidt
 .PP
diff --git a/man/man3/plschmidt_d1.3 b/man/man3/plschmidt_d1.3
index 7ad2438d2..a0435529a 100644
--- a/man/man3/plschmidt_d1.3
+++ b/man/man3/plschmidt_d1.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "plschmidt_d1" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "plschmidt_d1" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PlSchmidt_d1
 .PP
diff --git a/man/man3/preglq.3 b/man/man3/preglq.3
index 7cb9a4849..c94b4f691 100644
--- a/man/man3/preglq.3
+++ b/man/man3/preglq.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "preglq" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "preglq" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH PreGLQ
 .PP
diff --git a/man/man3/randomgaussian.3 b/man/man3/randomgaussian.3
index d4637a34b..5c46e5d1c 100644
--- a/man/man3/randomgaussian.3
+++ b/man/man3/randomgaussian.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "randomgaussian" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "randomgaussian" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH RandomGaussian
 .PP
diff --git a/man/man3/randomn.3 b/man/man3/randomn.3
index 844f080cf..ed5c37c44 100644
--- a/man/man3/randomn.3
+++ b/man/man3/randomn.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "randomn" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "randomn" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH RandomN
 .PP
diff --git a/man/man3/shadmitcorr.3 b/man/man3/shadmitcorr.3
index db3867227..1f360fc3f 100644
--- a/man/man3/shadmitcorr.3
+++ b/man/man3/shadmitcorr.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shadmitcorr" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shadmitcorr" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHAdmitCorr
 .PP
diff --git a/man/man3/shbias.3 b/man/man3/shbias.3
index ac730882d..80855a504 100644
--- a/man/man3/shbias.3
+++ b/man/man3/shbias.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shbias" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shbias" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHBias
 .PP
diff --git a/man/man3/shbiasadmitcorr.3 b/man/man3/shbiasadmitcorr.3
index ab3ef149a..e0cf1b21b 100644
--- a/man/man3/shbiasadmitcorr.3
+++ b/man/man3/shbiasadmitcorr.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shbiasadmitcorr" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shbiasadmitcorr" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHBiasAdmitCorr
 .PP
diff --git a/man/man3/shbiask.3 b/man/man3/shbiask.3
index 1f7c57f5c..414d93d7c 100644
--- a/man/man3/shbiask.3
+++ b/man/man3/shbiask.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shbiask" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shbiask" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHBiasK
 .PP
diff --git a/man/man3/shbiaskmask.3 b/man/man3/shbiaskmask.3
index ab4fc6e44..c5d27d27d 100644
--- a/man/man3/shbiaskmask.3
+++ b/man/man3/shbiaskmask.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shbiaskmask" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shbiaskmask" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHBiasKMask
 .PP
diff --git a/man/man3/shcilmtocindex.3 b/man/man3/shcilmtocindex.3
index 83781cba7..55f343679 100644
--- a/man/man3/shcilmtocindex.3
+++ b/man/man3/shcilmtocindex.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shcilmtocindex" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shcilmtocindex" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHCilmToCindex
 .PP
diff --git a/man/man3/shcilmtovector.3 b/man/man3/shcilmtovector.3
index 6e3b7dcc5..77c816dac 100644
--- a/man/man3/shcilmtovector.3
+++ b/man/man3/shcilmtovector.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shcilmtovector" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shcilmtovector" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHCilmToVector
 .PP
diff --git a/man/man3/shcindextocilm.3 b/man/man3/shcindextocilm.3
index 0be76ed1a..5cbb257c5 100644
--- a/man/man3/shcindextocilm.3
+++ b/man/man3/shcindextocilm.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shcindextocilm" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shcindextocilm" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHCindexToCilm
 .PP
diff --git a/man/man3/shconfidence.3 b/man/man3/shconfidence.3
index 5346eec50..d57d32e3f 100644
--- a/man/man3/shconfidence.3
+++ b/man/man3/shconfidence.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shconfidence" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shconfidence" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHConfidence
 .PP
diff --git a/man/man3/shcrosspowerdensityl.3 b/man/man3/shcrosspowerdensityl.3
index d07df39ef..6db941955 100644
--- a/man/man3/shcrosspowerdensityl.3
+++ b/man/man3/shcrosspowerdensityl.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shcrosspowerdensityl" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shcrosspowerdensityl" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHCrossPowerDensityL
 .PP
diff --git a/man/man3/shcrosspowerdensitylc.3 b/man/man3/shcrosspowerdensitylc.3
index 79fcb9571..04b0ad6ac 100644
--- a/man/man3/shcrosspowerdensitylc.3
+++ b/man/man3/shcrosspowerdensitylc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shcrosspowerdensitylc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shcrosspowerdensitylc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHCrossPowerDensityLC
 .PP
diff --git a/man/man3/shcrosspowerl.3 b/man/man3/shcrosspowerl.3
index 12ef85806..ca4ad74ce 100644
--- a/man/man3/shcrosspowerl.3
+++ b/man/man3/shcrosspowerl.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shcrosspowerl" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shcrosspowerl" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHCrossPowerL
 .PP
diff --git a/man/man3/shcrosspowerlc.3 b/man/man3/shcrosspowerlc.3
index 16606ff5f..1936c8007 100644
--- a/man/man3/shcrosspowerlc.3
+++ b/man/man3/shcrosspowerlc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shcrosspowerlc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shcrosspowerlc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHCrossPowerLC
 .PP
diff --git a/man/man3/shcrosspowerspectrum.3 b/man/man3/shcrosspowerspectrum.3
index aeb66d80b..3b1f86902 100644
--- a/man/man3/shcrosspowerspectrum.3
+++ b/man/man3/shcrosspowerspectrum.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shcrosspowerspectrum" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shcrosspowerspectrum" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHCrossPowerSpectrum
 .PP
diff --git a/man/man3/shcrosspowerspectrumc.3 b/man/man3/shcrosspowerspectrumc.3
index 5ee698c32..eed8f7fa5 100644
--- a/man/man3/shcrosspowerspectrumc.3
+++ b/man/man3/shcrosspowerspectrumc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shcrosspowerspectrumc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shcrosspowerspectrumc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHCrossPowerSpectrumC
 .PP
diff --git a/man/man3/shcrosspowerspectrumdensity.3 b/man/man3/shcrosspowerspectrumdensity.3
index b96a94683..5965c116e 100644
--- a/man/man3/shcrosspowerspectrumdensity.3
+++ b/man/man3/shcrosspowerspectrumdensity.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shcrosspowerspectrumdensity" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shcrosspowerspectrumdensity" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHCrossPowerSpectrumDensity
 .PP
diff --git a/man/man3/shcrosspowerspectrumdensityc.3 b/man/man3/shcrosspowerspectrumdensityc.3
index 40d59f3dc..fda398aed 100644
--- a/man/man3/shcrosspowerspectrumdensityc.3
+++ b/man/man3/shcrosspowerspectrumdensityc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shcrosspowerspectrumdensityc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shcrosspowerspectrumdensityc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHCrossPowerSpectrumDensityC
 .PP
diff --git a/man/man3/shctor.3 b/man/man3/shctor.3
index e8acb088d..205822582 100644
--- a/man/man3/shctor.3
+++ b/man/man3/shctor.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shctor" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shctor" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHctor
 .PP
diff --git a/man/man3/shexpanddh.3 b/man/man3/shexpanddh.3
index 414523fd6..146f0fa13 100644
--- a/man/man3/shexpanddh.3
+++ b/man/man3/shexpanddh.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shexpanddh" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shexpanddh" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHExpandDH
 .PP
diff --git a/man/man3/shexpanddhc.3 b/man/man3/shexpanddhc.3
index c2f87d1b2..8b8e26e81 100644
--- a/man/man3/shexpanddhc.3
+++ b/man/man3/shexpanddhc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shexpanddhc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shexpanddhc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHExpandDHC
 .PP
diff --git a/man/man3/shexpandglq.3 b/man/man3/shexpandglq.3
index 450d32438..1b0b785af 100644
--- a/man/man3/shexpandglq.3
+++ b/man/man3/shexpandglq.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shexpandglq" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shexpandglq" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHExpandGLQ
 .PP
diff --git a/man/man3/shexpandglqc.3 b/man/man3/shexpandglqc.3
index a8f5c542d..035052739 100644
--- a/man/man3/shexpandglqc.3
+++ b/man/man3/shexpandglqc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shexpandglqc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shexpandglqc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHExpandGLQC
 .PP
diff --git a/man/man3/shexpandlsq.3 b/man/man3/shexpandlsq.3
index 4a15908d0..d179d06cd 100644
--- a/man/man3/shexpandlsq.3
+++ b/man/man3/shexpandlsq.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shexpandlsq" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shexpandlsq" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHExpandLSQ
 .PP
diff --git a/man/man3/shfindlwin.3 b/man/man3/shfindlwin.3
index 8cccc89a4..0aa730dfc 100644
--- a/man/man3/shfindlwin.3
+++ b/man/man3/shfindlwin.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shfindlwin" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shfindlwin" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHFindLWin
 .PP
diff --git a/man/man3/shglq.3 b/man/man3/shglq.3
index c3d4dfd5d..32dab9504 100644
--- a/man/man3/shglq.3
+++ b/man/man3/shglq.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shglq" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shglq" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHGLQ
 .PP
diff --git a/man/man3/shlocalizedadmitcorr.3 b/man/man3/shlocalizedadmitcorr.3
index f0cf57bed..e34711fbc 100644
--- a/man/man3/shlocalizedadmitcorr.3
+++ b/man/man3/shlocalizedadmitcorr.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shlocalizedadmitcorr" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shlocalizedadmitcorr" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHLocalizedAdmitCorr
 .PP
diff --git a/man/man3/shmagpowerl.3 b/man/man3/shmagpowerl.3
index cad4a7b4b..0384b8119 100644
--- a/man/man3/shmagpowerl.3
+++ b/man/man3/shmagpowerl.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shmagpowerl" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shmagpowerl" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHMagPowerL
 .PP
diff --git a/man/man3/shmagpowerspectrum.3 b/man/man3/shmagpowerspectrum.3
index 3b6bdd03d..3081ca0db 100644
--- a/man/man3/shmagpowerspectrum.3
+++ b/man/man3/shmagpowerspectrum.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shmagpowerspectrum" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shmagpowerspectrum" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHMagPowerSpectrum
 .PP
diff --git a/man/man3/shmtcouplingmatrix.3 b/man/man3/shmtcouplingmatrix.3
index 9993a2dbd..62c080dda 100644
--- a/man/man3/shmtcouplingmatrix.3
+++ b/man/man3/shmtcouplingmatrix.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shmtcouplingmatrix" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shmtcouplingmatrix" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHMTCouplingMatrix
 .PP
diff --git a/man/man3/shmtdebias.3 b/man/man3/shmtdebias.3
index cd831260a..ea905e224 100644
--- a/man/man3/shmtdebias.3
+++ b/man/man3/shmtdebias.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shmtdebias" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shmtdebias" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHMTDebias
 .PP
diff --git a/man/man3/shmtvaropt.3 b/man/man3/shmtvaropt.3
index 20e963f2b..5f9278510 100644
--- a/man/man3/shmtvaropt.3
+++ b/man/man3/shmtvaropt.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shmtvaropt" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shmtvaropt" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHMTVarOpt
 .PP
diff --git a/man/man3/shmultiply.3 b/man/man3/shmultiply.3
index e5486dee7..2252e2b0b 100644
--- a/man/man3/shmultiply.3
+++ b/man/man3/shmultiply.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shmultiply" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shmultiply" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHMultiply
 .PP
diff --git a/man/man3/shmultitapercse.3 b/man/man3/shmultitapercse.3
index aa2e743cc..04a42c37b 100644
--- a/man/man3/shmultitapercse.3
+++ b/man/man3/shmultitapercse.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shmultitapercse" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shmultitapercse" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHMultiTaperCSE
 .PP
diff --git a/man/man3/shmultitapermaskcse.3 b/man/man3/shmultitapermaskcse.3
index e5fe8f510..db099e56b 100644
--- a/man/man3/shmultitapermaskcse.3
+++ b/man/man3/shmultitapermaskcse.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shmultitapermaskcse" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shmultitapermaskcse" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHMultiTaperMaskCSE
 .PP
diff --git a/man/man3/shmultitapermaskse.3 b/man/man3/shmultitapermaskse.3
index c2ed4aac8..8e75566d7 100644
--- a/man/man3/shmultitapermaskse.3
+++ b/man/man3/shmultitapermaskse.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shmultitapermaskse" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shmultitapermaskse" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHMultiTaperMaskSE
 .PP
diff --git a/man/man3/shmultitaperse.3 b/man/man3/shmultitaperse.3
index cbc91f314..dbeeafd5b 100644
--- a/man/man3/shmultitaperse.3
+++ b/man/man3/shmultitaperse.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shmultitaperse" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shmultitaperse" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHMultiTaperSE
 .PP
diff --git a/man/man3/shpowerdensityl.3 b/man/man3/shpowerdensityl.3
index 2629b479e..a16c74cf0 100644
--- a/man/man3/shpowerdensityl.3
+++ b/man/man3/shpowerdensityl.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shpowerdensityl" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shpowerdensityl" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHPowerDensityL
 .PP
diff --git a/man/man3/shpowerdensitylc.3 b/man/man3/shpowerdensitylc.3
index 65470ac0e..b008c137b 100644
--- a/man/man3/shpowerdensitylc.3
+++ b/man/man3/shpowerdensitylc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shpowerdensitylc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shpowerdensitylc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHPowerDensityLC
 .PP
diff --git a/man/man3/shpowerl.3 b/man/man3/shpowerl.3
index 8d973c5f4..88b0c2ec5 100644
--- a/man/man3/shpowerl.3
+++ b/man/man3/shpowerl.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shpowerl" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shpowerl" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHPowerL
 .PP
diff --git a/man/man3/shpowerlc.3 b/man/man3/shpowerlc.3
index 9c258552e..b7f6aacbb 100644
--- a/man/man3/shpowerlc.3
+++ b/man/man3/shpowerlc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shpowerlc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shpowerlc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHPowerLC
 .PP
diff --git a/man/man3/shpowerspectrum.3 b/man/man3/shpowerspectrum.3
index 50073ab2c..9557cfaf6 100644
--- a/man/man3/shpowerspectrum.3
+++ b/man/man3/shpowerspectrum.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shpowerspectrum" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shpowerspectrum" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHPowerSpectrum
 .PP
diff --git a/man/man3/shpowerspectrumc.3 b/man/man3/shpowerspectrumc.3
index 1aa5254cb..b4b76c304 100644
--- a/man/man3/shpowerspectrumc.3
+++ b/man/man3/shpowerspectrumc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shpowerspectrumc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shpowerspectrumc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHPowerSpectrumC
 .PP
diff --git a/man/man3/shpowerspectrumdensity.3 b/man/man3/shpowerspectrumdensity.3
index 30371099d..907c0d156 100644
--- a/man/man3/shpowerspectrumdensity.3
+++ b/man/man3/shpowerspectrumdensity.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shpowerspectrumdensity" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shpowerspectrumdensity" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHPowerSpectrumDensity
 .PP
diff --git a/man/man3/shpowerspectrumdensityc.3 b/man/man3/shpowerspectrumdensityc.3
index 1fa18186d..c7f48e5d0 100644
--- a/man/man3/shpowerspectrumdensityc.3
+++ b/man/man3/shpowerspectrumdensityc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shpowerspectrumdensityc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shpowerspectrumdensityc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHPowerSpectrumDensityC
 .PP
diff --git a/man/man3/shread.3 b/man/man3/shread.3
index f156e2e93..4ff84881a 100644
--- a/man/man3/shread.3
+++ b/man/man3/shread.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shread" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shread" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHRead
 .PP
diff --git a/man/man3/shread2.3 b/man/man3/shread2.3
index a4122f3bf..39389c0ff 100644
--- a/man/man3/shread2.3
+++ b/man/man3/shread2.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shread2" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shread2" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHRead2
 .PP
diff --git a/man/man3/shreadjpl.3 b/man/man3/shreadjpl.3
index 18319a84b..c730223f7 100644
--- a/man/man3/shreadjpl.3
+++ b/man/man3/shreadjpl.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shreadjpl" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shreadjpl" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHReadJPL
 .PP
diff --git a/man/man3/shreturntapers.3 b/man/man3/shreturntapers.3
index fc0f780b4..c6977eee2 100644
--- a/man/man3/shreturntapers.3
+++ b/man/man3/shreturntapers.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shreturntapers" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shreturntapers" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHReturnTapers
 .PP
diff --git a/man/man3/shreturntapersm.3 b/man/man3/shreturntapersm.3
index 159bd6b5a..f40b8269d 100644
--- a/man/man3/shreturntapersm.3
+++ b/man/man3/shreturntapersm.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shreturntapersm" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shreturntapersm" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHReturnTapersM
 .PP
diff --git a/man/man3/shreturntapersmap.3 b/man/man3/shreturntapersmap.3
index 26e2eae34..fae915a2e 100644
--- a/man/man3/shreturntapersmap.3
+++ b/man/man3/shreturntapersmap.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shreturntapersmap" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shreturntapersmap" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHReturnTapersMap
 .PP
diff --git a/man/man3/shrotatecoef.3 b/man/man3/shrotatecoef.3
index 323621aad..d5587d88f 100644
--- a/man/man3/shrotatecoef.3
+++ b/man/man3/shrotatecoef.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shrotatecoef" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shrotatecoef" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHRotateCoef
 .PP
diff --git a/man/man3/shrotaterealcoef.3 b/man/man3/shrotaterealcoef.3
index 95a33755d..c8612ecaa 100644
--- a/man/man3/shrotaterealcoef.3
+++ b/man/man3/shrotaterealcoef.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shrotaterealcoef" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shrotaterealcoef" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHRotateRealCoef
 .PP
diff --git a/man/man3/shrotatetapers.3 b/man/man3/shrotatetapers.3
index 9ae331d71..216c482c3 100644
--- a/man/man3/shrotatetapers.3
+++ b/man/man3/shrotatetapers.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shrotatetapers" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shrotatetapers" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHRotateTapers
 .PP
diff --git a/man/man3/shrtoc.3 b/man/man3/shrtoc.3
index 987b22074..31a24ce0f 100644
--- a/man/man3/shrtoc.3
+++ b/man/man3/shrtoc.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shrtoc" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shrtoc" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHrtoc
 .PP
diff --git a/man/man3/shscouplingmatrix.3 b/man/man3/shscouplingmatrix.3
index 7fd42501d..c10ce1863 100644
--- a/man/man3/shscouplingmatrix.3
+++ b/man/man3/shscouplingmatrix.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shscouplingmatrix" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shscouplingmatrix" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHSCouplingMatrix
 .PP
diff --git a/man/man3/shscouplingmatrixcap.3 b/man/man3/shscouplingmatrixcap.3
index 8debd36bf..b02ff74c8 100644
--- a/man/man3/shscouplingmatrixcap.3
+++ b/man/man3/shscouplingmatrixcap.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shscouplingmatrixcap" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shscouplingmatrixcap" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHSCouplingMatrixCap
 .PP
diff --git a/man/man3/shsjkpg.3 b/man/man3/shsjkpg.3
index 7fbba53ca..2528228d0 100644
--- a/man/man3/shsjkpg.3
+++ b/man/man3/shsjkpg.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shsjkpg" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shsjkpg" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHSjkPG
 .PP
diff --git a/man/man3/shslepianvar.3 b/man/man3/shslepianvar.3
index 626b65e37..3daa8701d 100644
--- a/man/man3/shslepianvar.3
+++ b/man/man3/shslepianvar.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shslepianvar" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shslepianvar" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHSlepianVar
 .PP
diff --git a/man/man3/shtools.3 b/man/man3/shtools.3
index cc48c18d3..230c6961f 100644
--- a/man/man3/shtools.3
+++ b/man/man3/shtools.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,14 +14,14 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shtools" "3" "2020-09-24" "Fortran 95" "SHTOOLS 4.10"
+.TH "shtools" "1" "2022-09-09" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHTOOLS
 .PP
 SHTOOLS is a Fortran-95 library that can be used for spherical harmonic
-transforms, multitaper spectral analyses, expansions of gridded data into
-Slepian basis functions, and standard operations on global gravitational and
-magnetic field data.
+transforms, multitaper spectral analyses, expansions of gridded data
+into Slepian basis functions, and standard operations on global
+gravitational and magnetic field data.
 .SH Features
 .IP \[bu] 2
 Supports all standard normalizations and phase conventions of the
@@ -33,8 +33,8 @@ Gauss-Legendre quadrature.
 Spherical harmonic transforms proven to be accurate up to about degree
 2800.
 .IP \[bu] 2
-Perform localized multitaper spectral analyses, or expand gridded data in
-terms of localized Slepian basis functions.
+Perform localized multitaper spectral analyses, or expand gridded data
+in terms of localized Slepian basis functions.
 .IP \[bu] 2
 Perform basic operations on global gravity and magnetic field data.
 .IP \[bu] 2
diff --git a/man/man3/shvectortocilm.3 b/man/man3/shvectortocilm.3
index 9a6926204..dd4cf433e 100644
--- a/man/man3/shvectortocilm.3
+++ b/man/man3/shvectortocilm.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "shvectortocilm" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "shvectortocilm" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SHVectorToCilm
 .PP
diff --git a/man/man3/slepiancoeffs.3 b/man/man3/slepiancoeffs.3
index e82a6687c..cc11da4b2 100644
--- a/man/man3/slepiancoeffs.3
+++ b/man/man3/slepiancoeffs.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "slepiancoeffs" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "slepiancoeffs" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SlepianCoeffs
 .PP
diff --git a/man/man3/slepiancoeffstosh.3 b/man/man3/slepiancoeffstosh.3
index 84c8b0331..61a0befa6 100644
--- a/man/man3/slepiancoeffstosh.3
+++ b/man/man3/slepiancoeffstosh.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "slepiancoeffstosh" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "slepiancoeffstosh" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SlepianCoeffsToSH
 .PP
diff --git a/man/man3/sphericalcapcoef.3 b/man/man3/sphericalcapcoef.3
index 5df5f4e5a..6b9f3739e 100644
--- a/man/man3/sphericalcapcoef.3
+++ b/man/man3/sphericalcapcoef.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "sphericalcapcoef" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "sphericalcapcoef" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH SphericalCapCoef
 .PP
diff --git a/man/man3/wigner3j.3 b/man/man3/wigner3j.3
index f91ce3438..6023a4a5f 100644
--- a/man/man3/wigner3j.3
+++ b/man/man3/wigner3j.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "wigner3j" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "wigner3j" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH Wigner3j
 .PP
diff --git a/man/man3/yilmindexvector.3 b/man/man3/yilmindexvector.3
index 70689dca7..8ef55aa8b 100644
--- a/man/man3/yilmindexvector.3
+++ b/man/man3/yilmindexvector.3
@@ -1,4 +1,4 @@
-.\" Automatically generated by Pandoc 2.17.1.1
+.\" Automatically generated by Pandoc 2.19.2
 .\"
 .\" Define V font for inline verbatim, using C font in formats
 .\" that render this, and otherwise B font.
@@ -14,7 +14,7 @@
 . ftr VB CB
 . ftr VBI CBI
 .\}
-.TH "yilmindexvector" "3" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
+.TH "yilmindexvector" "1" "2021-02-15" "Fortran 95" "SHTOOLS 4.10"
 .hy
 .SH YilmIndexVector
 .PP
diff --git a/pyshtools/__init__.py b/pyshtools/__init__.py
index 89fc7be11..ad6eef51e 100644
--- a/pyshtools/__init__.py
+++ b/pyshtools/__init__.py
@@ -43,7 +43,7 @@
 from ._version import get_versions as _get_versions
 
 # ---- Import shtools subpackages ----
-from . import shtools
+from . import backends
 from . import constants
 from . import datasets
 from . import shclasses
@@ -54,7 +54,7 @@
 from . import rotate
 from . import gravmag
 from . import utils
-from . import backends
+from .backends import shtools
 
 # ---- Import principal classes into pyshtools namespace
 from .shclasses import SHCoeffs
diff --git a/pyshtools/backends/__init__.py b/pyshtools/backends/__init__.py
index 151da17b1..dffc28635 100644
--- a/pyshtools/backends/__init__.py
+++ b/pyshtools/backends/__init__.py
@@ -24,7 +24,7 @@
 from .backends import preferred_backend_module
 from .backends import backend_module
 from . import ducc0_wrapper  # noqa: F401
-from .. import shtools  # noqa: F401
+from . import shtools  # noqa: F401
 
 del backends  # noqa: F821
 
diff --git a/pyshtools/backends/backends.py b/pyshtools/backends/backends.py
index a4066ff7e..56590a33c 100644
--- a/pyshtools/backends/backends.py
+++ b/pyshtools/backends/backends.py
@@ -19,10 +19,8 @@
 backend_module()             Return a reference to the specified backend
                              module.
 """
-
-
 from . import ducc0_wrapper
-from .. import shtools
+from . import shtools
 
 # defaults
 _preferred_backend = "shtools"
diff --git a/pyshtools/shtools/__init__.py b/pyshtools/backends/shtools.py
similarity index 100%
rename from pyshtools/shtools/__init__.py
rename to pyshtools/backends/shtools.py
diff --git a/pyshtools/expand/__init__.py b/pyshtools/expand/__init__.py
index 1aa4bb557..f1057aa9d 100644
--- a/pyshtools/expand/__init__.py
+++ b/pyshtools/expand/__init__.py
@@ -53,14 +53,14 @@
 spharm_lm      Compute the spherical harmonic function for a specific degree l
                and order m.
 """
-from ..shtools import SHGLQ
-from ..shtools import GLQGridCoord
-from ..shtools import SHExpandLSQ
-from ..shtools import SHExpandWLSQ
-from ..shtools import MakeGrid2D
-from ..shtools import MakeGridPoint
-from ..shtools import MakeGridPointC
-from ..shtools import SHMultiply
+from ..backends.shtools import SHGLQ
+from ..backends.shtools import GLQGridCoord
+from ..backends.shtools import SHExpandLSQ
+from ..backends.shtools import SHExpandWLSQ
+from ..backends.shtools import MakeGrid2D
+from ..backends.shtools import MakeGridPoint
+from ..backends.shtools import MakeGridPointC
+from ..backends.shtools import SHMultiply
 
 from .spharm_functions import spharm
 from .spharm_functions import spharm_lm
diff --git a/pyshtools/gravmag/__init__.py b/pyshtools/gravmag/__init__.py
index ad20140ee..f609d5075 100644
--- a/pyshtools/gravmag/__init__.py
+++ b/pyshtools/gravmag/__init__.py
@@ -55,22 +55,22 @@
 mag_spectrum        Compute the spectrum of either the magnetic potential
                     or magnetic field strength.
 """
-from ..shtools import MakeGravGridDH
-from ..shtools import MakeGravGridPoint
-from ..shtools import MakeGravGradGridDH
-from ..shtools import MakeGeoidGridDH
-from ..shtools import CilmPlusDH
-from ..shtools import CilmMinusDH
-from ..shtools import CilmPlusRhoHDH
-from ..shtools import CilmMinusRhoHDH
-from ..shtools import BAtoHilmDH
-from ..shtools import BAtoHilmRhoHDH
-from ..shtools import DownContFilterMA
-from ..shtools import DownContFilterMC
-from ..shtools import NormalGravity
-from ..shtools import MakeMagGridDH
-from ..shtools import MakeMagGridPoint
-from ..shtools import MakeMagGradGridDH
+from ..backends.shtools import MakeGravGridDH
+from ..backends.shtools import MakeGravGridPoint
+from ..backends.shtools import MakeGravGradGridDH
+from ..backends.shtools import MakeGeoidGridDH
+from ..backends.shtools import CilmPlusDH
+from ..backends.shtools import CilmMinusDH
+from ..backends.shtools import CilmPlusRhoHDH
+from ..backends.shtools import CilmMinusRhoHDH
+from ..backends.shtools import BAtoHilmDH
+from ..backends.shtools import BAtoHilmRhoHDH
+from ..backends.shtools import DownContFilterMA
+from ..backends.shtools import DownContFilterMC
+from ..backends.shtools import NormalGravity
+from ..backends.shtools import MakeMagGridDH
+from ..backends.shtools import MakeMagGridPoint
+from ..backends.shtools import MakeMagGradGridDH
 
 from .mag_spectrum import mag_spectrum
 
diff --git a/pyshtools/legendre/__init__.py b/pyshtools/legendre/__init__.py
index 322c563d5..1982ecdfe 100644
--- a/pyshtools/legendre/__init__.py
+++ b/pyshtools/legendre/__init__.py
@@ -51,22 +51,22 @@
 PlmIndex      Compute the index of an array of Legendre function corresponding
               to degree L and angular order M.
 """
-from ..shtools import PlmBar
-from ..shtools import PlmBar_d1
-from ..shtools import PlBar
-from ..shtools import PlBar_d1
-from ..shtools import PlmON
-from ..shtools import PlmON_d1
-from ..shtools import PlON
-from ..shtools import PlON_d1
-from ..shtools import PlmSchmidt
-from ..shtools import PlmSchmidt_d1
-from ..shtools import PlSchmidt
-from ..shtools import PlSchmidt_d1
-from ..shtools import PLegendreA
-from ..shtools import PLegendreA_d1
-from ..shtools import PLegendre
-from ..shtools import PLegendre_d1
+from ..backends.shtools import PlmBar
+from ..backends.shtools import PlmBar_d1
+from ..backends.shtools import PlBar
+from ..backends.shtools import PlBar_d1
+from ..backends.shtools import PlmON
+from ..backends.shtools import PlmON_d1
+from ..backends.shtools import PlON
+from ..backends.shtools import PlON_d1
+from ..backends.shtools import PlmSchmidt
+from ..backends.shtools import PlmSchmidt_d1
+from ..backends.shtools import PlSchmidt
+from ..backends.shtools import PlSchmidt_d1
+from ..backends.shtools import PLegendreA
+from ..backends.shtools import PLegendreA_d1
+from ..backends.shtools import PLegendre
+from ..backends.shtools import PLegendre_d1
 
 from .legendre_functions import legendre
 from .legendre_functions import legendre_lm
diff --git a/pyshtools/rotate/__init__.py b/pyshtools/rotate/__init__.py
index c8c64a152..97eab86d5 100644
--- a/pyshtools/rotate/__init__.py
+++ b/pyshtools/rotate/__init__.py
@@ -10,9 +10,8 @@
 SHRotateRealCoef  Determine the spherical harmonic coefficients of a real
                   function rotated by three Euler angles.
 """
-from ..shtools import djpi2
-from ..shtools import SHRotateCoef
-
+from ..backends.shtools import djpi2
+from ..backends.shtools import SHRotateCoef
 from ..backends import backend_module, select_preferred_backend
 
 
diff --git a/pyshtools/shclasses/shcoeffs.py b/pyshtools/shclasses/shcoeffs.py
index 0d9ad6a71..9432e6265 100644
--- a/pyshtools/shclasses/shcoeffs.py
+++ b/pyshtools/shclasses/shcoeffs.py
@@ -13,7 +13,6 @@
 from scipy.special import factorial as _factorial
 import xarray as _xr
 
-from .. import shtools as _shtools
 from ..spectralanalysis import spectrum as _spectrum
 from ..spectralanalysis import cross_spectrum as _cross_spectrum
 from ..shio import convert as _convert
@@ -25,6 +24,7 @@
 from ..shio import write_bshc as _write_bshc
 from ..backends import backend_module
 from ..backends import preferred_backend
+from ..backends import shtools as _shtools
 
 
 class SHCoeffs(object):
diff --git a/pyshtools/shclasses/shgravcoeffs.py b/pyshtools/shclasses/shgravcoeffs.py
index 593b57310..9fc416696 100644
--- a/pyshtools/shclasses/shgravcoeffs.py
+++ b/pyshtools/shclasses/shgravcoeffs.py
@@ -31,13 +31,13 @@
 from ..shio import write_bshc as _write_bshc
 from ..shio import read_icgem_gfc as _read_icgem_gfc
 from ..shio import write_icgem_gfc as _write_icgem_gfc
-from ..shtools import CilmPlusRhoHDH as _CilmPlusRhoHDH
-from ..shtools import CilmPlusDH as _CilmPlusDH
-from ..shtools import MakeGravGridDH as _MakeGravGridDH
-from ..shtools import MakeGravGradGridDH as _MakeGravGradGridDH
-from ..shtools import MakeGeoidGridDH as _MakeGeoidGridDH
-from ..shtools import djpi2 as _djpi2
-from ..shtools import MakeGravGridPoint as _MakeGravGridPoint
+from ..backends.shtools import CilmPlusRhoHDH as _CilmPlusRhoHDH
+from ..backends.shtools import CilmPlusDH as _CilmPlusDH
+from ..backends.shtools import MakeGravGridDH as _MakeGravGridDH
+from ..backends.shtools import MakeGravGradGridDH as _MakeGravGradGridDH
+from ..backends.shtools import MakeGeoidGridDH as _MakeGeoidGridDH
+from ..backends.shtools import djpi2 as _djpi2
+from ..backends.shtools import MakeGravGridPoint as _MakeGravGridPoint
 from ..backends import backend_module
 from ..backends import preferred_backend
 
diff --git a/pyshtools/shclasses/shgrid.py b/pyshtools/shclasses/shgrid.py
index 8dd68ac8f..7e1bd5715 100644
--- a/pyshtools/shclasses/shgrid.py
+++ b/pyshtools/shclasses/shgrid.py
@@ -10,7 +10,7 @@
 import tempfile as _tempfile
 from ..backends import backend_module
 from ..backends import preferred_backend
-from .. import shtools as _shtools
+from ..backends import shtools as _shtools
 
 try:
     import cartopy.crs as _ccrs
diff --git a/pyshtools/shclasses/shmagcoeffs.py b/pyshtools/shclasses/shmagcoeffs.py
index 105fd726e..d0cf6013e 100644
--- a/pyshtools/shclasses/shmagcoeffs.py
+++ b/pyshtools/shclasses/shmagcoeffs.py
@@ -26,10 +26,10 @@
 from ..shio import read_bshc as _read_bshc
 from ..shio import write_bshc as _write_bshc
 from ..shio import read_igrf as _read_igrf
-from ..shtools import MakeMagGridDH as _MakeMagGridDH
-from ..shtools import MakeMagGradGridDH as _MakeMagGradGridDH
-from ..shtools import djpi2 as _djpi2
-from ..shtools import MakeMagGridPoint as _MakeMagGridPoint
+from ..backends.shtools import MakeMagGridDH as _MakeMagGridDH
+from ..backends.shtools import MakeMagGradGridDH as _MakeMagGradGridDH
+from ..backends.shtools import djpi2 as _djpi2
+from ..backends.shtools import MakeMagGridPoint as _MakeMagGridPoint
 from ..backends import backend_module
 from ..backends import preferred_backend
 
@@ -103,6 +103,8 @@ class SHMagCoeffs(object):
     set_coeffs()          : Set coefficients in-place to specified values.
     change_ref()          : Return a new class instance referenced to a
                             different reference radius.
+    change_units()        : Return a new class instance with the internal
+                            coefficients converted to either nT or T.
     rotate()              : Rotate the coordinate system used to express the
                             spherical harmonic coefficients and return a new
                             class instance.
@@ -1852,6 +1854,50 @@ def change_ref(self, r0=None, lmax=None):
 
         return clm
 
+    def change_units(self, units=None, lmax=None):
+        """
+        Return a new SHMagCoeffs class instance with the internal coefficients
+        converted to either nT or T.
+
+        Usage
+        -----
+        clm = x.change_units([units, lmax])
+
+        Returns
+        -------
+        clm : SHMagCoeffs class instance.
+
+        Parameters
+        ----------
+        units : str, optional, default = self.units
+            The new units of the spherical harmonic coefficients, which can be
+            either 'T' or 'nT'.
+        lmax : int, optional, default = self.lmax
+            Maximum spherical harmonic degree to output.
+        """
+        if lmax is None:
+            lmax = self.lmax
+        if units is None:
+            units = self.units
+
+        clm = self.pad(lmax)
+
+        if units != self.units:
+            if units.lower() == 't' and self.units.lower() == 'nt':
+                factor = 1.e-9
+            elif units.lower() == 'nt' and self.units.lower() == 't':
+                factor = 1.e9
+            else:
+                raise ValueError("units must be either 'nT' or 'T'. "
+                                 "Input value is {:s}".format(repr(units)))
+
+            clm.coeffs *= factor
+            clm.units = units
+            if clm.errors is not None:
+                clm.errors *= factor
+
+        return clm
+
     # ---- Routines that return different magnetic-related class instances ----
     def expand(self, a=None, f=None, lat=None, colat=None, lon=None,
                degrees=True, lmax=None, lmax_calc=None, sampling=2,
diff --git a/pyshtools/shclasses/shwindow.py b/pyshtools/shclasses/shwindow.py
index 9b2ddee77..454e7ecfb 100644
--- a/pyshtools/shclasses/shwindow.py
+++ b/pyshtools/shclasses/shwindow.py
@@ -9,7 +9,7 @@
 from mpl_toolkits.axes_grid1 import make_axes_locatable as _make_axes_locatable
 import copy as _copy
 
-from .. import shtools as _shtools
+from ..backends import shtools as _shtools
 from ..spectralanalysis import spectrum as _spectrum
 
 from .shcoeffs import SHCoeffs
diff --git a/pyshtools/shclasses/slepian.py b/pyshtools/shclasses/slepian.py
index 72c374a06..6ce82b468 100644
--- a/pyshtools/shclasses/slepian.py
+++ b/pyshtools/shclasses/slepian.py
@@ -9,7 +9,7 @@
 from mpl_toolkits.axes_grid1 import make_axes_locatable as _make_axes_locatable
 import copy as _copy
 
-from .. import shtools as _shtools
+from ..backends import shtools as _shtools
 from ..spectralanalysis import spectrum as _spectrum
 
 from .shcoeffs import SHCoeffs
diff --git a/pyshtools/shclasses/slepiancoeffs.py b/pyshtools/shclasses/slepiancoeffs.py
index f42a9bbf3..f15443d78 100644
--- a/pyshtools/shclasses/slepiancoeffs.py
+++ b/pyshtools/shclasses/slepiancoeffs.py
@@ -6,7 +6,7 @@
 
 from ..backends import backend_module
 from ..backends import preferred_backend
-from .. import shtools as _shtools
+from ..backends import shtools as _shtools
 
 from .shcoeffs import SHCoeffs
 from .shgrid import SHGrid
diff --git a/pyshtools/shio/__init__.py b/pyshtools/shio/__init__.py
index fe708aa5c..9530d1f99 100644
--- a/pyshtools/shio/__init__.py
+++ b/pyshtools/shio/__init__.py
@@ -54,16 +54,16 @@
 SHrtoc           Convert real spherical harmonics to complex form.
 SHctor           Convert complex spherical harmonics to real form.
 """
-from ..shtools import SHRead2
-from ..shtools import SHRead2Error
-from ..shtools import SHReadJPL
-from ..shtools import SHReadJPLError
-from ..shtools import SHCilmToCindex
-from ..shtools import SHCindexToCilm
-from ..shtools import SHCilmToVector
-from ..shtools import SHVectorToCilm
-from ..shtools import SHrtoc
-from ..shtools import SHctor
+from ..backends.shtools import SHRead2
+from ..backends.shtools import SHRead2Error
+from ..backends.shtools import SHReadJPL
+from ..backends.shtools import SHReadJPLError
+from ..backends.shtools import SHCilmToCindex
+from ..backends.shtools import SHCindexToCilm
+from ..backends.shtools import SHCilmToVector
+from ..backends.shtools import SHVectorToCilm
+from ..backends.shtools import SHrtoc
+from ..backends.shtools import SHctor
 
 from .convert import convert
 from .shtools import shread
diff --git a/pyshtools/spectralanalysis/__init__.py b/pyshtools/spectralanalysis/__init__.py
index 186ff3b03..8346230b4 100644
--- a/pyshtools/spectralanalysis/__init__.py
+++ b/pyshtools/spectralanalysis/__init__.py
@@ -89,42 +89,41 @@
 SphericalCapCoef       Calculate the spherical harmonic coefficients of a
                        spherical cap.
 """
-from ..shtools import SHAdmitCorr
-from ..shtools import SHConfidence
+from ..backends.shtools import SHAdmitCorr
+from ..backends.shtools import SHConfidence
+from ..backends.shtools import SHMultiTaperSE
+from ..backends.shtools import SHMultiTaperCSE
+from ..backends.shtools import SHLocalizedAdmitCorr
+from ..backends.shtools import SHReturnTapers
+from ..backends.shtools import SHReturnTapersM
+from ..backends.shtools import ComputeDm
+from ..backends.shtools import ComputeDG82
+from ..backends.shtools import SHFindLWin
+from ..backends.shtools import SHBiasK
+from ..backends.shtools import SHMTCouplingMatrix
+from ..backends.shtools import SHBiasAdmitCorr
+from ..backends.shtools import SHMTDebias
+from ..backends.shtools import SHMTVarOpt
+from ..backends.shtools import SHMTVar
+from ..backends.shtools import SHSjkPG
+from ..backends.shtools import SHMultiTaperMaskSE
+from ..backends.shtools import SHMultiTaperMaskCSE
+from ..backends.shtools import SHReturnTapersMap
+from ..backends.shtools import SHBiasKMask
+from ..backends.shtools import ComputeDMap
+from ..backends.shtools import Curve2Mask
+from ..backends.shtools import SHBias
+from ..backends.shtools import SphericalCapCoef
+from ..backends.shtools import SHRotateTapers
+from ..backends.shtools import SlepianCoeffs
+from ..backends.shtools import SlepianCoeffsToSH
+from ..backends.shtools import SHSCouplingMatrix
+from ..backends.shtools import SHSlepianVar
+from ..backends.shtools import SHSCouplingMatrixCap
 
 from .spectrum import spectrum
 from .cross_spectrum import cross_spectrum
 
-from ..shtools import SHMultiTaperSE
-from ..shtools import SHMultiTaperCSE
-from ..shtools import SHLocalizedAdmitCorr
-from ..shtools import SHReturnTapers
-from ..shtools import SHReturnTapersM
-from ..shtools import ComputeDm
-from ..shtools import ComputeDG82
-from ..shtools import SHFindLWin
-from ..shtools import SHBiasK
-from ..shtools import SHMTCouplingMatrix
-from ..shtools import SHBiasAdmitCorr
-from ..shtools import SHMTDebias
-from ..shtools import SHMTVarOpt
-from ..shtools import SHMTVar
-from ..shtools import SHSjkPG
-from ..shtools import SHMultiTaperMaskSE
-from ..shtools import SHMultiTaperMaskCSE
-from ..shtools import SHReturnTapersMap
-from ..shtools import SHBiasKMask
-from ..shtools import ComputeDMap
-from ..shtools import Curve2Mask
-from ..shtools import SHBias
-from ..shtools import SphericalCapCoef
-from ..shtools import SHRotateTapers
-from ..shtools import SlepianCoeffs
-from ..shtools import SlepianCoeffsToSH
-from ..shtools import SHSCouplingMatrix
-from ..shtools import SHSlepianVar
-from ..shtools import SHSCouplingMatrixCap
-
 
 # ---- Define __all__ for use with: from pyshtools import * ----
 __all__ = ['SHAdmitCorr', 'SHConfidence', 'spectrum', 'cross_spectrum',
diff --git a/pyshtools/utils/__init__.py b/pyshtools/utils/__init__.py
index 16d408550..ba184721b 100644
--- a/pyshtools/utils/__init__.py
+++ b/pyshtools/utils/__init__.py
@@ -13,10 +13,10 @@
 figstyle           Set Matplotlib parameters for creating publication quality
                    graphics.
 """
-from ..shtools import MakeCircleCoord
-from ..shtools import MakeEllipseCoord
-from ..shtools import Wigner3j
-from ..shtools import DHaj
+from ..backends.shtools import MakeCircleCoord
+from ..backends.shtools import MakeEllipseCoord
+from ..backends.shtools import Wigner3j
+from ..backends.shtools import DHaj
 from .figstyle import figstyle
 
 
diff --git a/setup.py b/setup.py
index d408e40be..ad56ccbf0 100644
--- a/setup.py
+++ b/setup.py
@@ -70,7 +70,7 @@
 
 EXTRAS_REQUIRE = {
     'cartopy': ['cython', 'pyshp', 'six', 'shapely', 'cartopy>=0.18.0'],
-    'pygmt': ['pygmt==0.3'],
+    'pygmt': ['pygmt>=0.3'],
     'palettable': ['palettable>=3.3'],
     'ducc': ['ducc0>=0.15']
 }
diff --git a/src/pydoc/select_preferred_backend.md b/src/pydoc/select_preferred_backend.md
index e497ab99f..17ee65403 100644
--- a/src/pydoc/select_preferred_backend.md
+++ b/src/pydoc/select_preferred_backend.md
@@ -5,11 +5,11 @@ transforms in pyshtools.
 
 # Usage
 
-select_preferred_backend_module([backend, nthreads])
+select_preferred_backend([backend, nthreads])
 
 # Parameters
 
-backend : str, optional, default = 'shtools'
+backend : str, optional, default = 'ducc'
 :   Name of the preferred backend, either 'shtools' or 'ducc'.
 
 nthreads : int, optional, default = 1


Test patch:
diff --git a/examples/cpp/test.cpp b/examples/cpp/test.cpp
index a5e7d304c..ced7f47e5 100644
--- a/examples/cpp/test.cpp
+++ b/examples/cpp/test.cpp
@@ -26,7 +26,7 @@ int
 main(int argc, char** argv)
 {
 
-  std::string infile = "../ExampleDataFiles/MarsTopo719.shape";
+  std::string infile;
 
   int lmax = 15;
   int cilm_dim = lmax + 1;
@@ -37,6 +37,12 @@ main(int argc, char** argv)
   
   int exitstatus;
 
+  if( argc > 1 ){
+    infile = argv[1];
+  } else {
+    infile = "../ExampleDataFiles";
+  }
+  infile += "/MarsTopo719.shape";
   shtools::SHRead(infile.c_str(),
                    infile.size(),
                    &cilm[0],