anchor.diff:  !    GLAT     Geodetic latitude (degrees)
   !    GLON     Geodetic longitude (degrees)
   !    ALT      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)
  -!    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
  +!    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)
   !
  -!***************************************************************************************************
  +!*******************************************************************************
   !
   !  APXG2ALL
   !  Converts geodetic to Quasi-Dipole and Modified Apex coordinates
  -!  Must call LOADAPXSH first. All arguments are single precision, except for MAFLAG and VECFLAG,
  -!  which are 4-byte integers.
  +!  Must call LOADAPXSH first. All arguments are single precision, except for
  +!  MAFLAG and VECFLAG, which are 4-byte integers.
   !
   !  CALL APXG2ALL (GLAT, GLON, ALT, HR, VECFLAG, &
   !                 QLAT, QLON, MLAT, MLON, &
  
positive.code: 
      implicit none
  
      character(128)              :: datafilenew, datafilelast=''
      real(4)                     :: epochnew, epochlast=-999.0
      real(8)                     :: we0, we1
      integer(4)                  :: iepoch0, iepoch1, iterm, icoord
  
negative.code: 
  end subroutine apxg2all
  
  !*******************************************************************************
  
  subroutine apxq2g(qlat0,qlon0,alt,prec,glatout,glonout,error)
  
  