问

925450204

我的程序为什么我的for循环进行不下去？？？？

;data=make_array(144,73,17,718,/float)
cdfid=ncdf_open('atmprfzhuan.nc')
varid2=ncdf_varid(cdfid,'refrac')
varid1=ncdf_varid(cdfid,'alt_refrac')
ncdf_varget,cdfid,varid2,atm1
ncdf_varget,cdfid,varid1,atm2
;ncdf_attget,cdfid,varid,'scale_factor',a   ;系数
;ncdf_attget,cdfid,varid,'add_offset',b      ;偏值
ncdf_close,cdfid
help,atm1,atm2
;openw,lun,'d:\1.dat',/get_lun
;writeu,lun,data1,data2
;free_lun,lun
;print,data[*,*,1,1]
;print,atm1[0],atm2[0]
;print,atm2,ecm2
;a=replicate(0.000,50)
;for i=1L,50L,1L do begin
;a(i)=i*1000
;endfor
;print,a
cdfid=ncdf_open('ucar1.nc')
varid2=ncdf_varid(cdfid,'refrac')
varid1=ncdf_varid(cdfid,'alt_refrac')
ncdf_varget,cdfid,varid2,ropp1
ncdf_varget,cdfid,varid1,ropp2
ncdf_close,cdfid
help,atm1,atm2
a=[0,2,4,6,8,10,12,14,16,18,20,22,24,26,28,30,32,34,36,38,40,42,44,46,48,50]
atmrec=interpol(atm1,atm2/1000,a)
ropprec=interpol(ropp1,ropp2/1000,a)
;print,atmrec,ecmrec
bias=atmrec-ropprec
print,bias
;print,b
plot,bias,a,xtitle="bias",ytitle="alt(km)",xrange=[-10,10],yrange=[0,50],$
title='refrac COSMIC-ROPP',ycharsize=0.85,xcharsize=0.85
plotS,[0,0],[0,50],LineStyle=2,color='0000FF'XL
Window, size=400,xsize=400
LoadCT, 3
image = LoadData(7)
TVImage, image, Position=[0.1, 0.1, 0.75, 0.9]
Colorbar, Position=[0.8, 0.1, 0.86, 0.9], /Right, $
/Vertical, Division=5, Format='(F5.1)'
image = TVRD()
image = BytArr(3, 400, 400)
TVLCT, r, g, b, /Get
image3D[0, *,*]=r[image]
image3D[1, *,*] =g [image]
image3D[2, *, *] = b[image]
image3D=TVRD(true=1)
    Device, Decomposed=1
    !P.Background='FFFFFF'XL
    !P.Color='000000'XL
Write_JPEG,' COSMIC-ROPP.jpg', image3D, true=1, quality=75
end

尽头的尽头

以后发帖把题目写清楚一点吧。。。

chenhongfengye

你把for都注释了，怎么进行循环啊？

925450204

chenhongfengye 发表于 2012-12-28 09:47

                               
登录/注册后可看大图

你把for都注释了，怎么进行循环啊？

我搞定了，现在就是如何把画的图保存为图片格式？？有程序吗？

chenhongfengye

我们一般都保存为.ps格式

chenhongfengye

然后用Gsview看图

王比比

有程序吗

925450204

chenhongfengye 发表于 2012-12-29 12:42

                               
登录/注册后可看大图

然后用Gsview看图

求idl保存为ps格式图片的程序

chenhongfengye

outfile=path_out_dir+'cmp_precipitable-3_0000.ps'

set_plot,'ps'
;device,/color,/Helvetica

device,filename=outfile,/portrait,Ysize=12,Yoffset=-0.7,Xsize=24,Xoffset=-2.3

!x.thick=1.5
!x.charsize=1
!p.charthick=2
!y.style=1
!y.thick=1.5
!p.charsize=1.75
!p.thick=1.5
!y.charsize=1
!p.font=1
!x.style=1

xxr

请问我输入mage = LoadData(7）显示没有loaddata这个文件是怎么回事？