求助，用ncl输出含多变量不同维数数据的grd文件

寒枫沁雪

    fbinrecwrite(grdout,-1,cu)   
     fbinrecwrite(grdout,-1,cv)  
     fbinrecwrite(grdout,-1,cw)      
     fbinrecwrite(grdout,-1,cqv)   
     fbinrecwrite(grdout,-1,cqc)  
     fbinrecwrite(grdout,-1,cqr)  
     fbinrecwrite(grdout,-1,cqi)
     fbinrecwrite(grdout,-1,cqs)
     fbinrecwrite(grdout,-1,cqg)
     fbinrecwrite(grdout,-1,nc)  
     fbinrecwrite(grdout,-1,cnr)  
     fbinrecwrite(grdout,-1,cni)
     fbinrecwrite(grdout,-1,ns)
     fbinrecwrite(grdout,-1,ng)      
     fbinrecwrite(grdout,-1,ccldfra)  
     fbinrecwrite(grdout,-1,cgeopt)  
     fbinrecwrite(grdout,-1,ch)
     fbinrecwrite(grdout,-1,ctc)
     fbinrecwrite(grdout,-1,crh)
     fbinrecwrite(grdout,-1,cdbz)
     fbinrecwrite(grdout,-1,rssr)
     fbinrecwrite(grdout,-1,cpsfc)  
     fbinrecwrite(grdout,-1,chgt)  
     fbinrecwrite(grdout,-1,ctsk)
     fbinrecwrite(grdout,-1,crainc)
     fbinrecwrite(grdout,-1,crainnc)
     fbinrecwrite(grdout,-1,cxland)
     fbinrecwrite(grdout,-1,cslp)
     fbinrecwrite(grdout,-1,cmax_dbz)
     fbinrecwrite(grdout,-1,ccband)  
     fbinrecwrite(grdout,-1,cvil)  
     fbinrecwrite(grdout,-1,cvisl)   
     fbinrecwrite(grdout,-1,cttop)   
     fbinrecwrite(grdout,-1,cztop)  
     fbinrecwrite(grdout,-1,ctbot)  
     fbinrecwrite(grdout,-1,czbot)  

我用fbinrecwrite函数把这二三十个变量写进同一个grd文件中，这些变量的维数有的不同。结果用如下ctl打开，发现前几个变量可以打开，后面的变量提示错误。low level i/o Error   请问是不是写入的方式不对
options sequential template
TITLE WRF  DATA
undef 1.e30
xdef  752 linear  111.51984   0.01351351
ydef  627 linear   25.87175   0.01351351
zdef   19 levels  
1000.00000
950.00000
900.00000
850.00000
800.00000
750.00000
700.00000
650.00000
600.00000
550.00000
500.00000
450.00000
400.00000
350.00000
300.00000
250.00000
200.00000
150.00000
100.00000
tdef    49 linear 00Z13MAR2016   60MN      
VARS   36
U             19  0  x-wind component (m s-1)
V             19  0  y-wind component (m s-1)
W             19  0  z-wind component (m s-1)
qvapor        19  0  Water vapor mixing ratio (kg kg-1)
qcloud        19  0  Cloud water mixing ratio (kg kg-1)
qrain         19  0  Rain water mixing ratio (kg kg-1)
qice          19  0  Ice mixing ratio (kg kg-1)
qsnow         19  0  Snow mixing ratio (kg kg-1)
qgraup        19  0  Graupel mixing ratio (kg kg-1)
qncloud       19  0  cloud water Number concentration (  kg(-1))
qnrain        19  0  Rain Number concentration (  kg(-1))
qnice         19  0  Ice Number concentration (  kg-1)
qnsnow        19  0  Snow Number concentration (  kg(-1))
qngraupel     19  0  Graupel Number concentration (  kg(-1))
cldfra        19  0  CLOUD FRACTION (-)
geopt         19  0  Geopotential (m2/s2)
height        19  0  Model height (km)
tc            19  0  Temperature (C)
rh            19  0  Relative Humidity (%)
dbz           19  0  Reflectivity (-)
rssr          19  0  relative super saturation ratio
psfc           1  0  SFC PRESSURE (Pa)
hgt            1  0  Terrain Height (m)
tsk            1  0  SURFACE SKIN TEMPERATURE (K)
rainc          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
rainnc         1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
xland          1  0  LAND MASK (1 FOR LAND, 2 FOR WATER) (-)
slp            1  0  Sea Levelp Pressure (hPa)
max_dbz        1  0  Max Reflectivity (-)
cband         0 0  cloud band (mm)
vil           0 0  vertical integration of liquid water (mm)
visl          0 0  vertical integration of supercooled water (mm)
cloudtopt     0 0 cloudtop temprature (degree)
cloudtoph     0 0 cloudtop height (km)
cloudbott     0 0 cloudbottom temprature (degree)
cloudboth     0 0 cloudbottom height (km)
endvars