; Determine geographic coordinates of PACE pole at given altitude. IF KEYWORD_SET(rot) THEN BEGIN rdcols,'~cnsr3/pace/pace.pole',n,alt,glatp,glongp ; for 950101 wht=WHERE(alt EQ height,nwht) IF nwht EQ 1 THEN BEGIN gltp=glatp(wht) glnp=glongp(wht) ENDIF ELSE BEGIN polrec3d,relc(90.)+alt,90.-glatp,glongp,xp,yp,zp,/degrees xpole=interpol(xp,alt,height) ypole=interpol(yp,alt,height) zpole=interpol(zp,alt,height) recpol3d,xpole,ypole,zpole,rpole,colatpole,glnp,/degrees gltp=90.-colatpole ; PACE pole at (gltp,glnp) ENDELSE glnp=glnp-(glnp GT 180.)*360. PACEpolecoords=[glnp,gltp,height] ENDIF ELSE PACEpolecoords=[280.73,80.57,300.]