-------------------------------------------------------------------------------
   Computes inclination drift for geosynchronous orbits.
   Use this function to plan north/south stationkeeping.
-------------------------------------------------------------------------------
   Form:
   [dV, dT] = DVIDrift( box, year )
-------------------------------------------------------------------------------

   ------
   Inputs
   ------
   box                 Size of stationkeeping box (deg)
   year                Year of maneuvers


   -------
   Outputs
   -------
   dV                  DeltaV per maneuver (m/sec)
   dT                  Time between maneuvers (days)

-------------------------------------------------------------------------------