Dear all here is my IDL program to make contour plot.
1)The program I am having is that after loading "Loadct,27" my program
is not able make the colorful plot. And also I am not able to save the
colorful plots.
2) Hoe I can assign labeling most of the contour as I am having no. of
Rows=572.
Thanks in advance
readcol, 'spectrum.txt',skipline=3,x,y,z,flag,f='f,f,f'
set_plot,'PS'
sdevice, filename='/home/arvind/Desktop/A.eps'
;*********************************************************** ******************************
tmp = where(flag eq 5.0)
chi2 = reform(x(tmp),11,26)
radius = reform(y(tmp),11,26)
temperature = reform(z(tmp),11,26)
;*********************************************************** ********************************
tmp = where(flag eq 5.5)
chi12 = reform(x(tmp),11,26)
radius1 = reform(y(tmp),11,26)
temperature1 = reform(z(tmp),11,26)
;*********************************************************** *******************************
window,xs=600,ys=900
loadct,27
tvlct,0,0,0,255
!p.charsize=1.5
!p.multi=[0,0,2,0]
;=========================================================== ========================
contour,chi2,temperature,radius,xtit='Teff',ytit='Radius', $
c_label=[1,1,1,1,1,1,1,1],nlevel=100,/xst,/yst,/fill;,Positi on=[0.1,
0.1, 0.9, 0.80]
contour,chi2,temperature,radius,c_label=[1,1,1,1,1,1,1,1],le vel=[0.25,0.5,1,2,5,10]/
100.0,$
c_charsize=1.5,c_color=[255,0,0,0,0,255,255],/over
;=========================================================== ========================
contour,chi2< 0.02,temperature,radius,xtit='Teff',ytit='Radius',xr=[2600,3 100],
$
c_label=[1,1,1,1,1,1],nlevels=100,/xst,/yst,/fill
contour,chi2< 0.02,temperature,radius,c_label=[1,1,1,1,1,1],level=[0.25,0. 45,1,2]/
100.0,$
c_charsize=1.5,c_color=[0,0,255,255],/over
;colorbar, BOTTOM=4,DIVISIONS=5,FORMAT='(f6.4)', POSITION=[0.1, 0.9,
0.9, 0.95], $
;MAX=0.02, MIN=MIN(chi2)
;=========================================================== ========================
!p.charsize=1.0
!p.multi=0
device,/close
set_plot,'x
;=========================================================== =======================
print,'End of Program'
end
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
