PRO invcm2wl,wlmin,wlmax,ps=ps

IF N_Params() LT 2 THEN BEGIN
    Print,'Usage: invcm2wl,wlmin,wlmax,ps=ps'
    Print,'(Arguments in Angstroms)'
    Return
    END

wlmin=wlmin>8000<wlmax
wlmax=wlmax<8500>wlmin

data=DDRead('\data\modtran\tape7.txt')
wavnum=Double(Transpose(data(0,*)))
transmission=Transpose(data(1,*))

vacwl=1.0D8/wavnum      ; (1E8 A/cm)/(v cm-1)
VacToAir,vacwl          ; argument changes to wavelength in air
airwl=vacwl

IF N_Elements(ps) GT 0 THEN ps=1
IF Keyword_Set(ps) THEN BEGIN
    psname=String(Round(wlmin/10),Round(wlmax/10) $
                 ,Format='("wl",2I3,".ps")')
    PSOpen,'c:\data\modtran\'+psname
    ENDIF
    
PLOT,airwl,transmission,PSym=-1 $
    ,XRange=[wlmin,wlmax],XStyle=1 $
    ,XTitle='Wavelength in Air ('+Angstrom()+')' $
    ,YRange=[0.7,1.0],YTitle='Atmospheric Transmission ' $
    ,Title='MODTRAN, Vertical Path, 85 km to Sea Level, Sub-Arctic Winter'

IF Keyword_Set(ps) THEN PSClose

Return
END