valueo3.bak

用IDL编程绘制等高线。希望对大家有所帮助。

pro valueo3
str=''
y=IntArr(1);年
m=intarr(1);月
d=intarr(1);日
hour=intArr(1);小时
minute=intarr(1);分钟
degrees=fltarr(1);太阳高度角
lat_data=fltArr(1);纬度
lon_data=fltArr(1);经度
profile_data=fltArr(65);臭氧剖面值

;!P.Multi=[0,1,4,0,0];指定一页画1列2行图，即一页画2个图
TvLCT,255,0,0,5
thisDevice=!D.Name
Set_Plot,'ps'
Device,filename='2002.ps',XSize=14,YSize=15,Xoffset=0.03,Yoffset=0.01,/Encapsulated,/Preview ;,decomposed=0


clrtble1=intarr(12)
clrtble2=intarr(12)
clrtble1=[1,2,3,4,5,6,7,8,9,10,11,12]
clrtble2=[1,2,3,4,5,6,7,8,9,10,11,12]

red=  [0,205,0,  100,45, 0,  0,  152,255,255,255,255,167]
green=[0,205,0,  100,255,180,255,255,255,150,106,0,  34]
blue= [0,205,255,255,234,4,  0,  13, 0,  1,  106,0,  27]

backgound=[0,0,0]

TVLCT,red,green,blue

openr,lun,'pO3 1-15.txt',/Get_lun

k=0;记录数
while(eof(lun)ne 1)do begin
readF,lun,str,y,m,d,hour,minute,degrees,lat_data,lon_data,profile_data,format='(a,/,I4,2X,I2,2X,I2,2x,I2,2X,I2,2X,f6.2,2X,f7.2,2X,f7.2,/,11F9.2,/,11F9.2,/,11F9.2,/,11F9.2,/,11F9.2,/,10f9.2)'
k=k+1
endwhile

;lon=fltarr(k)
;lat=fltarr(k)
value=fltarr(k,3)

free_lun,lun
openr,lun,'pO3 1-15.txt',/Get_lun
i=0
while(eof(lun)ne 1)do begin
readF,lun,str,y,m,d,hour,minute,degrees,lat_data,lon_data,profile_data,format='(a,/,I4,2X,I2,2X,I2,2x,I2,2X,I2,2X,f6.2,2X,f7.2,2X,f7.2,/,11F9.2,/,11F9.2,/,11F9.2,/,11F9.2,/,11F9.2,/,10f9.2)'

;lon[i]=lon_data
;lat[i]=lat_data
value(i,0)=lon_data
value(i,1)=lat_data
value(i,2)=profile_data(9)
i=i+1
endwhile

contour,value(*,2),value(*,0),value(*,1),Xstyle=1,Ystyle=1,Title='200275-105 28.5km',XTitle='lon',Ytitle='lat',$
    Nlevels=7,Xticks=10
device,Preview=0,/close_file
end