Three-dimensional subsatellite motion under air drag and oblateness perturbations