|
|
|
Re: how to get the data from the contour line? [message #43631 is a reply to message #43630] |
Tue, 19 April 2005 18:58  |
Craig Markwardt
Messages: 1869 Registered: November 1996
|
Senior Member |
|
|
"lixiaoyao" <lixiaoyao5880@yahoo.com> writes:
> 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?
How about reading the documentation for CONTOUR and the keywords
beginning in PATH?
Craig
> 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"
>
--
------------------------------------------------------------ --------------
Craig B. Markwardt, Ph.D. EMAIL: craigmnet@REMOVEcow.physics.wisc.edu
Astrophysics, IDL, Finance, Derivatives | Remove "net" for better response
------------------------------------------------------------ --------------
|
|
|