pro solareph_sbsp,iy,im,id,utimed,ras,dcs,gsdt,b0ang,peeang,ctrlong,r0 ; routine to compute the solar ephemeris values give the time input ; coordinate in arcseconds from sun center to solar latitude and longitude ; INPUTS: ; ; iy = year (e.g., 2002) ; im = month (1-12) ; id = day of month ; utimed = universal time of day (hours). ; ; ; OUTPUTS: ; ; ras = right ascension disk center(radians). ; dcs = declination disk center(radians). ; gsdt = right ascension from central meridian(radians). ; b0ang = b0-angle(radians). ; peeang = p-angle(radians). ; ctrlong = carrington longitude disk center(radians). ; r0 = solar radius (arcsecs). djd = 0.d iy = long(iy) im = long(im) id = long(id) ; --