; the following function converts Mars-solar-orbital (MSO) coordinates to Mars solar local
; time.
; The inputs are:
; xmso, ymso, zmso are the Cartesian components in MSO coordinates
; latss is the latitude of the subsolar point in degrees
;
; The outputs are:
; local time, in the usual 0 to 24 hours
;
; Created by: Robert Lillis (rlillis@SSL.Berkeley.edu)
pro mso2lt, xmso, ymso, zmso, latss, local_time
; the first thing to do is rotate the MSO coordinate system around the
; Y. axis by an amount equal to the subsolar latitude. Reason for
; this is that local time is determined with respect to the planetary
; pole
rotate_y, 1.0*latss*!dtor, xmso, ymso, zmso, x1, y1, z1
cart2latlong, x1, y1, z1, r, solar_lat_rot, solar_elon_rot
; now rotate around the y-axis by the subsolar latitude.
;print,'rotated coordinates: ',solar_lat_rot, solar_elon_rot
local_time = ((solar_elon_rot +180.0) mod 360)/15.0
end