- 积分
- 4593
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-10-2
- 最后登录
- 1970-1-1

|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
pro SOILMOISTURE_mpdi_36
year=2005
day1=julday(1,1,year)-julday(1,1,year)+1 ;对应实际资料第二天
; print,day1
; stop
day2=julday(12,31,year)-julday(1,1,year)+1
row_xj0=532 ;xinjiang area
row_xj1=603
line_xj0=328
line_xj1=435
row_xj=row_xj1-row_xj0+1
line_xj=line_xj1-line_xj0+1
mpdi1 =fltarr(row_xj,line_xj,day2-day1+1)
TAO_VWC=fltarr(row_xj,line_xj,day2-day1+1)
mpdimin=fltarr(row_xj,line_xj)
openr,lun,'J:\soilmoisture\PRODUCT\2005mpdi_minnew.txt',/get_lun
readf,lun,mpdimin
close,lun
free_lun,lun
for day=day1, day2 do begin
dday=STRTRIM(string(day),1)
if day lt 100 then dday='0'+STRTRIM(string(day),1)
if day lt 10 then dday='00'+STRTRIM(string(day),1)
orbit="d"
row=721
line=721
l=72
n=108
channel=uintarr(row,line,6)
mpdimin=fltarr(l,n)
;****************************************
i=0
filename='NL'+STRTRIM(string(year),1)+dday+orbit+'.36v'
openr,lun,filepath(filename,root_dir='J:\amsr\'+STRTRIM(string(year),1)),/get_lun
TB=uintarr(row,line)
readu,lun,TB
channel(*,*,i)=TB
close,lun
free_lun,lun
i++
这个程序是读取AMSR-E数据中的一部分,但是每次都在红色字体那报错,“OPENR:Error opening file. Unit 100, File:E:\data\2006.8.29\2006\2006001d.36v”
请问各位IDL的高手,程序错在哪里,要怎么改正?
|
|