anchor.diff:  !    GLON     Geodetic longitude (degrees)
   !    ALT      Altitude (km)
   !    HR       Modified Apex coordinates reference altitude (km)
  -!    VECFLAG  Set this flag to zero to calculate coordinates only, without base vectors. The base
  -!             vector output arguments will be filled with -9999.
  +!    VECFLAG  Set this flag to zero to calculate coordinates only, without base
  +!             vectors. The base vector output arguments will be filled with
  +!             -9999.
   !
   !  OUTPUT ARGUMENTS
   !    QLAT      Quasi-Dipole latitude (degrees)
   !    QLON      Quasi-Dipole longitude (degrees)
   !    MLAT      Modified Apex latitude (degrees)
   !    MLON      Modified Apex longitude (degrees)
  -!    F1,F2     Quasi-Dipole coordinates base vectors described in Richmond (1995). The two elements
  -!              of each vector contain the geodetic east and north components. F1 points in the QD
  -!              east direction, and F2 points in the QD north direction.
  +!    F1,F2     Quasi-Dipole coordinates base vectors described in
  +!              Richmond (1995). The two elements of each vector contain the
  +!              geodetic east and north components. F1 points in the QD east
  +!              direction, and F2 points in the QD north direction.
   !    F         F = |F1 x F2|, as described in Richmond (1995)
  -!    D1,D2,D3  Modified Apex coordinates base vectors described in Richmond (1995). The three
  -!              elements of each vector contain the geodetic east, north, and up components.
  -!              D1 and D2 are generally perpendicular to the magnetic field and point in the
  -!              magnetic east and downward/equatorward directions. D3 points along the magnetic
  -!              field.
  +!    D1,D2,D3  Modified Apex coordinates base vectors described in
  +!              Richmond (1995). The three elements of each vector contain the
  +!              geodetic east, north, and up components. D1 and D2 are generally
  +!              perpendicular to the magnetic field and point in the magnetic
  +!              east and downward/equatorward directions. D3 points along the
  +!              magnetic field.
   !    D         D = |D1 x D2|, as described in Richmond (1995)
  -!    E1,E2,E3  Modified Apex coordinates base vectors described in Richmond (1995). The three
  -!              elements of each vector contain the geodetic east, north, and up components.
  -!              E1 and E2 are generally perpendicular to the magnetic field and point in the
  -!              magnetic east and downward/equatorward directions. E3 points along the magnetic
  -!              field.
  -!
  -!***************************************************************************************************
  +!    E1,E2,E3  Modified Apex coordinates base vectors described in
  +!              Richmond (1995). The three elements of each vector contain the
  +!              geodetic east, north, and up components. E1 and E2 are generally
  +!              perpendicular to the magnetic field and point in the magnetic
  +!              east and downward/equatorward directions. E3 points along the
  +!              magnetic field.
  +!
  +!*******************************************************************************
   !
   !  APXQ2G
   !  Converts Quasi-Dipole to geodetic coordinates
  
positive.code: 
  end subroutine apxg2all
  
  !***************************************************************************************************
  
  subroutine apxq2g(qlat0,qlon0,alt,prec,glatout,glonout,error)
  
  
negative.code: 
  end subroutine alfbasisinit
  
  !******************************************************************************
  
  subroutine alfbasis(nmax,mmax,theta,P,V,W)
  
  