comp.lang.idl-pvwave archive: archive

Home » Public Forums » archive » CONTOUR problem

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

CONTOUR problem [message #62135]

Tue, 26 August 2008 16:21

xiao zhang
Messages: 81
Registered: June 2008

Member

HI~ everyone~ i use the contour pro to display a 2Darray named s11.
And I wrote it like this:

CONTOUR,s11,rangle,vangle,xrange=[0,180],yrange=[0,90],ytitl e='Viewing
zenith angle',xtitle='Relative azimuth angle',/C_LABELS

But the xrange and yrange do not work and it only labeled on line in
the image , Any one see where the problem is ?

Thank you very much~

The whole program is like this if u want to see it:

rangle=findgen(180)
vangle=findgen(90)
y1=rangle*3.1415926*2/360
y2=vangle*3.1415926*2/360
;plot,y
thta1=22.5*3.1415926*2/360
thta2=37.5*3.1415926*2/360
thta3=52.5*3.1415926*2/360
thta4=67.5*3.1415926*2/360

;print,sin(thta2)

r1=fltarr(180,90)
r2=fltarr(180,90)
r3=fltarr(180,90)
r4=fltarr(180,90)

x10=sin(thta1)*cos(y1)
x20=sin(y1)*cos(y2)
y10=sin(thta1)*sin(y1)
y20=sin(y1)*sin(y2)

x11=sin(thta1)*cos(y1)
x21=sin(y1)*cos(y2)
y11=sin(thta1)*sin(y1)
y21=sin(y1)*sin(y2)

x12=sin(thta1)*cos(y1)
x22=sin(y1)*cos(y2)
y12=sin(thta1)*sin(y1)
y22=sin(y1)*sin(y2)

x13=sin(thta1)*cos(y1)
x23=sin(y1)*cos(y2)
y13=sin(thta1)*sin(y1)
y23=sin(y1)*sin(y2)

for i=0,179 do begin
for j=0,89 do begin
r1(i,j)=x10(i)*x20(j)+y10(i)*y20(j)+cos(y20(j))*cos(thta1)
r2(i,j)=x11(i)*x21(j)+y11(i)*y21(j)+cos(y21(j))*cos(thta2)
r3(i,j)=x12(i)*x22(j)+y12(i)*y22(j)+cos(y22(j))*cos(thta3)
r4(i,j)=x13(i)*x23(j)+y13(i)*y23(j)+cos(y23(j))*cos(thta4)

endfor
endfor
help,r1

s1=acos(r1)*360/2/3.1415926
s2=acos(r2)*360/2/3.1415926
s3=acos(r3)*360/2/3.1415926
s4=acos(r4)*360/2/3.1415926

;print,acos(r1)*360/2/3.1415926
s11=180-s1
s22=180-s2
s33=180-s3
s44=180-s4

;help,x1
;print,s1
!P.background=255
!P.multi=[0,2,2]

CONTOUR,s11,rangle,vangle,xrange=[0,180],yrange=[0,90],ytitl e='Viewing
zenith angle',xtitle='Relative azimuth angle',/C_LABELS
CONTOUR,s22,rangle,vangle,xrange=[0,180],yrange=[0,90],ytitl e='Viewing
zenith angle',xtitle='Relative azimuth angle',/C_LABELS
CONTOUR,s33,rangle,vangle,xrange=[0,180],yrange=[0,90],ytitl e='Viewing
zenith angle',xtitle='Relative azimuth angle',/C_LABELS
CONTOUR,s44,rangle,vangle,xrange=[0,180],yrange=[0,90],ytitl e='Viewing
zenith angle',xtitle='Relative azimuth angle',/C_LABELS

end

Report message to a moderator

Previous Topic:	Re: heap_free fatal error
Next Topic:	animation_image_doc.pro

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed Oct 08 19:32:54 PDT 2025

Total time taken to generate the page: 0.00460 seconds