; IDL Version 4.0.1 (Win32 x86)
; Journal File for user@Microsoft Windows Host
; Working directory: C:\RSI\IDL40
; Date: Mon Apr 22 14:58:25 1996
 
im=lonarr(64,64,17)
rdkimg,'\images\dark\cam0\d1900060.002',dkhb,dkim
dark=lin_smooth(dkim,/down_columns)
az=fltarr(64,64) & za=az
imf=rfindfile('\images\airglow\cam0\filt2\g1902060.*',count=nimf)
img=feb19th(imf(60),dark,az,za)
hlp,az,za
;       1   Float array (64, 64).   Min = 0.000000,  Max = 355.366
;       2   Float array (64, 64).   Min = 0.000000,  Max = 47.7645
im(*,*,0)=img
for i=1,16 do begin img=feb19th(imf(60+i),dark) & im(*,*,i)=img & endfor
window,0,xs=640,ys=512
for i=0,16 do tvscl,rebin(im(*,*,i)<pcentile(im(*,*,i),0.99),128,128,/sample),i
for i=0,16 do tvscl,rebin(im(*,*,i)<500,128,128,/sample),i
for i=0,16 do tvscl,rebin(im(*,*,i)<400,128,128,/sample),i
for i=0,16 do tvscl,rebin(im(*,*,i)<300,128,128,/sample),i
help,im
wdelete,0
help
.RUN C:\IDL\STEELE\FEB19TH.PRO
.RUN C:\IDL\STEELE\FEB19TH.PRO
img=feb19th(imf(60),dark,az,za)
img=feb19th(imf(60),dark,az,za,/debug)
im(*,*,0)=img
for i=1,16 do begin img=feb19th(imf(60+i),dark,/debug) & im(*,*,i)=img & endfor
help
save,im,az,za,filename='feb19th.sav'
doc_library,'rb2ll'
;----- Documentation for RB2LL.PRO
;
doc_library,'rb2ll'
;----- Documentation for RB2LL.PRO
; NAME:
;       RB2LL
; PURPOSE:
;       From range, bearing compute latitude, longitude .
; CATEGORY:
; CALLING SEQUENCE:
;       rb2ll, lng0, lat0, dist, azi, lng1, lat1
; INPUTS:
;       lng0, lat0 = long, lat of starting point (deg).     in
;       dist = range to point of interest in RADIANS.       in
;       azi = azimuth to point of interest (degrees).       in
; KEYWORD PARAMETERS:
; OUTPUTS:
;       lng1, lat1 = long, lat of point of interest (deg).  out
; COMMON BLOCKS:
; NOTES:
;       Notes: A unit sphere is assumed, thus dist is in radians.
;         Useful constants:
;         Radius of Earth (mean) = 6371.23 km = 3958.899 miles.
;         Distance to horizon from height H above surface:
;           For small H: dist = sqrt(2*H/R) in Radians
;           For large H: dist = acos(R/(R+H)) in Radians
;         To plot horizon from lat0, lng0, H:
;           rb2ll,lng0,lat0,dist,makex(0,360,1),plng,plat
;           plots,plng-360.,plat,psym=3
; MODIFICATION HISTORY:
;       R. Sterner, 13 Feb,1991
; Copyright (C) 1994, Johns Hopkins University/Applied Physics Laboratory
; This software may be used, copied, or redistributed as long as it is not
; sold and this copyright notice is reproduced on each copy made.  This
; routine is provided as is without any express or implied warranties
; whatsoever.  Other limitations apply as described in the file disclaimer.txt.
;
doc_library,'ll2rb'
;----- Documentation for LL2RB.PRO
;
doc_library,'ll2rb'
;----- Documentation for LL2RB.PRO
; NAME:
;       LL2RB
; PURPOSE:
;       From latitude, longitude compute range, bearing.
; CATEGORY:
; CALLING SEQUENCE:
;       ll2rb, lng0, lat0, lng1, lat1, dist, azi
; INPUTS:
;       lng0, lat0 = long, lat of reference point (deg).    in
;       lng1, lat1 = long, lat of point of interest (deg).  in
; KEYWORD PARAMETERS:
; OUTPUTS:
;       dist = range to point point of interest (radians).  out
;       azi = azimuth to point of interest (degrees).       out
; COMMON BLOCKS:
; NOTES:
;       Notes: A unit sphere is assumed, thus dist is in radians
;         so to get actual distance multiply dist by radius.
;          Useful constants:
;          Radius of Earth (mean) = 6371.23 km = 3958.899 miles.
; MODIFICATION HISTORY:
;       R. Sterner, 13 Feb,1991
; Copyright (C) 1991, Johns Hopkins University/Applied Physics Laboratory
; This software may be used, copied, or redistributed as long as it is not
; sold and this copyright notice is reproduced on each copy made.  This
; routine is provided as is without any express or implied warranties
; whatsoever.  Other limitations apply as described in the file disclaimer.txt.
;
ll2rb,80.0533,273.4161,74.78,265.12,dist,azi & print,dist,azi
;     0.144646      176.891
ll2rb,273.4161,80.0533,265.12,74.78,dist,azi & print,dist,azi
;    0.0970627      203.009
doc_library,'eqldaz'
;----- Documentation for EQLDAZ.PRO
; Purpose:
; EQLDAZ calculates the distance (degrees and km) and azimuth
; between a position in decimal degrees (ETH latitude, EPH longitude),
; and an array of N decimal lat (TH) and long (PHI) positions.
; All input positions are in geographic coordinates which are
; converted to equidistant latitude coordinates for all internal
; calculations using the method of Brown (1984).
; The first-order great-ellipse correction is applied (equation (36)).
; Author: DRHO (using the results of Brown (1984) GJRAS, 445.)
; input:
; double        eth             initial latitude (degrees, N+)
; double        eph             initial longitude (degrees, E+)
; double        th(n)           array of latitudes (degrees, N+)
; double        phi(n)          array of longitudes (degrees, E+)
; integer       n               number of points array
; output:
; double        xdeg            distance between points (degrees)
; double        dist            distance between points (kilometers)
; double        az              azimuth from initial point (degrees)
;
eqldaz,80.0533,273.4161,74.78,265.12,1,xdeg,dist,az
print,xdeg,dist,az
;       5.5836206
;       620.87092
;      -156.98807