请高手指教

zjhzjh

求助 fwrite出错。以下是我写的计算水汽通量散度的gs，直接出图是对的，但fwrite写出后写ctl打开画图，数据显示是-9e+08 to 0.请教高手指点
'reinit'
'set gxout fwrite'
'set fwrite d:\zjh\30pj\UP_ave_30y_0712qvf1.grd'
'sdfopen F:\ERA5_07m_30y\1981-2010_07_12.nc'
'set dfile 1'
'set lon 100 120'
'set lat 30 40'
'set lev 850'
'define vq1=ave(q*1000*mag(u,v)/9.8,t=1,t=240,8)'
'define vqj1=sqrt(ave(pow((q*1000*mag(u,v)/9.8)-vq1,2),t=1,t=240,8))'
'define div1=ave(hdivg(u*q*1000/9.8,v*q*1000/9.8)*1000000,t=1,t=240,8)'
'define divj1=sqrt(ave(pow(hdivg(u*q*1000/9.8,v*q*1000/9.8)*1000000-div1,2),t=1,t=240,8))'

'define vq2=ave(q*1000*mag(u,v)/9.8,t=2,t=240,8)'
'define vqj2=sqrt(ave(pow((q*1000*mag(u,v)/9.8)-vq1,2),t=2,t=240,8))'
'define div2=ave(hdivg(u*q*1000/9.8,v*q*1000/9.8)*1000000,t=2,t=240,8)'
'define divj2=sqrt(ave(pow(hdivg(u*q*1000/9.8,v*q*1000/9.8)*1000000-div1,2),t=2,t=240,8))'
'd vq1'
'd vqj1'
'd div1'
'd divj1'
'd vq2'
'd vqj2'
'd div2'
'd divj2'
'disable fwrite'

'set gxout contour'
'set lev 850'
'd div1'
'printim  d:/zjh/30pj/qf1.png white'
c
'd div2'
'printim  d:/zjh/30pj/qf2.png white'
c
'd divj1'
'printim  d:/zjh/30pj/qfj1.png white'
c
'd divj2'
'printim  d:/zjh/30pj/qfj2.png white'
；