Hi all
This time I have read all the manual for contour,I have a question,my
code is
in the following. I will say problem following,please help me.
;DEVICE,RETAIN=2
set_plot,'ps',/copy
device,filename='isothermal.ps'
device,/landscape
MAX=100
x=2*findgen(MAX)/MAX+0.1
y=2*findgen(MAX)/MAX+0.1
F=exp(y^2/2-alog(y))
G=(exp(2*alog(x)+1/x))
R=fltarr(MAX,MAX)
for i=0,MAX-1 do R(i,*)=F
for j=0,MAX-1 do R(*,j)/=G
contour,R,x,y,levels=[0.6,0.7,0.8,0.892521,1.0,1.1,1.2],c_la bels=[1,1,1,1,1,1,1],xtitle='Dimensionless
Radius',ytitle='Mach number',title='Isothermal case'
until now,it give the correct answer.
MAX=100
x=2*findgen(MAX)/MAX+0.1
y=2*findgen(MAX)/MAX+0.1
F=exp(y^2/2-alog(y))
G=(exp(2*alog(x)+1/x))
R=fltarr(MAX,MAX)
for i=0,MAX-1 do R(i,*)=F
for j=0,MAX-1 do R(*,j)/=G
contour,R,x,y,level=[0.892521],path_xy=xy,path_info=info,clo sed=0,/path_double
for I=0,(N_ELEMENTS(info)-1) DO BEGIN
S=[INDGEN(info(I).N),0]
print,xy(*,INFO(I).OFFSET+S)
; plots,xy(*,INFO(I).OFFSET+S),/norm
endfor
this times contour results do not agree with I got from the above code.
I have read all the contour help file,I just can not know to solver
it.IF who knows the question is,please tell me.
Thanks a lot
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
