PRO noxprofile,no2ht,no2dens,no3ht,no3dens,show=show

IF N_Elements(show) EQ 0 THEN show=0

no2ht=[22.2,23.1,24.8,26.3,29.3,32.1,35.7,40.3,42.3,46.8]
no2dens=[2.39,4.12,4.78,4.99,3.77,2.49,1.31,0.495,0.272,0.092]*1.E9

no3ht=Findgen(10)*2+20
no3dens=[0.67,0.71,0.78,0.86,1.02,1.72,3.1,3.4,3.1,2.6]*1E7
no3fht=Findgen(40)*2
r=0.67/0.71
no3fdens=[0.67E7*r^(11-Indgen(10)),no3dens,2.6E7*(0.5^(Indgen(20)+1))]

mtno2=ddread('\usask\893\crossect\dsteele.txt')
f=(101325.*1E-6)/(8.3143*273.15)*1E-5*6.02252E23

IF Keyword_Set(show) THEN BEGIN
    Window,/Free
    win=!D.Window
    Help,win
    Plot_OI,no2dens,no2ht $
           ,XTitle='Concentration (cm!e-3!n)',XRange=[1E6,1E10] $
           ,YTitle='Altitude (km)'
    XYOutS,1E9,40,'[NO!i2!n]',/Data,CharSize=2
    OPlot,no3fdens,no3fht,LineStyle=2
    XYOutS,4E7,34,'[NO!i3!n]',/Data,CharSize=2

    OPlot,f*mtno2(1,*),mtno2(0,*),LineStyle=1
    XYOutS,6E8,5,'[NO!i2!n](Modtran)',/Data
    Print,'Press any key to continue...'
    k=Get_Kbrd(1)
    WDelete,win
    ENDIF

no2ht=Reform(mtno2(0,*))
no2dens=Reform(f*mtno2(1,*))

no3ht=no3fht
no3dens=no3fdens

Return
END