输出:
system("rm -f /public/home/sunxiaoyun/datadir/datasets/Q1Q2/ncep_daily/"+iyear+"/airtemp.bin")
do i = 0,file_num-1
do j = 0,size_p-1
do k = 0,size_lat-1
do m = 0,size_lon-1
fbindirwrite("/public/home/sunxiaoyun/datadir/datasets/Q1Q2/ncep_daily/"+iyear+"/airtemp.bin",t(i,j,k,m))
end do
end do
end do
end do
输入:
fili="/home/ncep_daily/"+iyear+"/airtemp.bin"
setfileoption("bin","ReadByteOrder","LittleEndian")
t_daily=fbindirread(fili,0,(/368,26,181,360/),"float")
lon_ori=fspan(0,359,360)
lat_ori=fspan(-90,90,181);
level=(/1000,975,950,925,900,850,800,750,700,650,600,550,500,450,400,350,300,250,200,150,100,70,50,30,20,10/)
t_daily!0="time"
t_daily!1="level"
t_daily!2="lat"
t_daily!3="lon"
t_daily@_FillValue =32766
t_daily&lat = lat_ori
t_daily&lon = lon_ori
t_daily&level = level