function spacecraft_geodetic, spacecraft_id, jdn, debug=debug
;
; Given spacecraft ID and Julian Day Number, returns the
; spacecrafts geodetic position (latitude, longitude, height)
;
; B. Knapp, 2004-02-10
;
; Show usage?
  if n_params() lt 2 then begin
    print,' '
    print,' Given a Norad (US SpaceCom) spacecraft ID and double-'
    print,' precision Julian Day Number and fraction, this function'
    print,' returns the geodetic coordinates of the spacecraft as an'
    print,' array of length 3 (i.e., [longitude, latitude, height];'
    print,' degrees, km). The spacecraft_id  must be an integer'
    print,' scalar (e.g., 21701=UARS, 27651=SORCE), but the Julian'
    print,' Day Number argument may be a scalar or an array.'
    print,' '
    print,' Usage:'
    print,' '
    print,'     gc = spacecraft_geodetic(spacecraft_id, jdn)'
    return,' '
  endif
;
; Get the ECI coordinates of the spacecraft (rectangular coords)
  spacecraft_pv, spacecraft_id, jdn, eci, $
    tle_path='/mizar/proj/sorce/science_data_system/data/'
;
; Convert the ECI coordinates to geodetic
  return, eci_to_geodetic(jdn, eci, debug=debug)
;
end