| how to get the data from the contour line? [message #43632] |
Tue, 19 April 2005 18:36 |
lixiaoyao
Messages: 49
Registered: April 2005
|
Member |
|
|
hello all
I have a file like this to draw the contour using IDL.my probelm is
how to get the data for the contour level is equal to 1.2 case?
who can help me for this question?
set_plot,'ps',/copy
device,filename='isothermal.ps'
device,/landscape
;window,xs=640,ys=480
x=2*findgen(101)/100+0.2
y=2*findgen(101)/100+0.2
xx=rebin(x,n_elements(x),n_elements(y))
yy=rebin(transpose(y),n_elements(x),n_elements(y))
;shade_surf,(exp(2*alog(yy)+1/yy))/(exp(xx^2)-alog(xx)),xx,y y
;contour,(exp(2*alog(yy)+1/yy))/(exp(xx^2)-alog(xx)),yy,xx
shade_surf,exp(xx^2/2-alog(xx))/(exp(2*alog(yy)+1/yy)),yy,xx ,xtitle= "Dimensionless
Radius in log scale",ytitle="Mach number in log scale"
contour,exp(xx^2/2-alog(xx))/(exp(2*alog(yy)+1/yy)),yy,xx,le vel=[0.4,0.6,0.892521,1.2,1.4],/xlog,/ylog,xtitle= "Dimensionless
Radius",ytitle="Mach number",title="Isothermal case"
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
