;+ ; NAME: lstnow ; ; PURPOSE: ; return current lst, for a given observatory ; ; INPUTS: ; KEYWORDS: ; OUTPUTS: ; ; HISTORY: ; Began 2005-07-14 17:30:20 by Marshall Perrin ;- FUNCTION lstnow,obsname=obsname,longitude=longitude if ~(keyword_set(longitude)) and ~(keyword_set(obsname)) then obsname="lick" if keyword_set(obsname) then begin observatory,obsname,obs longitude=obs.longitude endif if ~(keyword_set(longitude)) then message,"Need to set longitude!" tz_offset = gmt_offsec()/3600. ; julday() doesn't know about timezones, so adjust the ; returned value as necessary jdnow = julday()+tz_offset/24. ct2lst,lst,longitude*(-1),tz_offset,jdnow return,lst end