輸出：
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
