One possible solution is to generate a hidden plot in the same window and then create an axis to it. For example:
x=10d^(findgen(41)*0.125)
y=1.0/(1.0+x*1e-3)
x2=10d^((alog10(x)-15.)/1.3)
p = plot(x,y,/xlog,axis_style=1)
p2= plot(x2,y,/xlog,/current,/hide,axis_style=0)
a = axis('x',location='top', target=p2, axis_range=[x2[0],x2[40]], /log)
|