본문 바로가기

IDL

라디오존데 데이터 읽어서 plot 그리기...

사용자 삽입 이미지
pro aaa

data=fltarr(4,9)
openr, lun, 'c:\raob.txt',/get_lun
readf, lun, data
free_lun, lun

mean1=fltarr(9)
min1=mean1
max1=mean1
s=mean1

for i=0, 8 do begin
mean1(i)=mean(data(1:3,i))
min1(i)=min(data(1:3,i))
max1(i)=max(data(1:3,i))
s(i)=STDDEV(data(1:3,i))
endfor

x=[1000,925,850,700,500,400,300,200,100]

device, decomposed=0
loadct, 0

window, xsize=400, ysize=700
plot, mean1,x, color=255, /ylog, yrange=[1000,100]

for i=0, 8 do begin
oplot, [min1(i),max1(i)],[x(i),x(i)], color=200
endfor
USERSYM, [0,0],[-0.7,0.8]
for i=0, 8 do begin
oplot, [min1(i)],[x(i)], color=200, psym=8
endfor

for i=0, 8 do begin
oplot, [max1(i)],[x(i)], color=200, psym=8
endfor

oplot, min1, x, color=140
oplot, max1, x, color=140

end



최저. 최고값을 |--------------| 로 처리해봤다. psym을 이렇게 써먹는 거구만~

IDL은 배우면 배울수록 신기한 거 같다... ㄷㄷㄷ