load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin

 dir1 = "/cygdrive/d/FNL2000-2011/200006/"
....................
 print(filepath1)
 n1 = dimsizes(filepath1)
 print(n1)
................................
 high = 100
 u1 = f1[:]->U_GRD_3_ISBL(:,{high},::-1,:)
 v1 = f1[:]->V_GRD_3_ISBL(:,{high},::-1,:)
 temp1 = f1[:]->TMP_3_ISBL(:,{high},::-1,:)
 rh1 = f1[:]->R_H_3_ISBL(:,{high},::-1,:)
 p4 = f1[:]->lv_ISBL6(:,{high})
 q1 = temp1
 qu1 = temp1
 qv1 = temp1
 q1 = mixhum_ptrh(conform(rh1,p4,0),temp1,rh1,2)
 qu1 = q1 * u1
 qv1 = q1 * v1

 qduv = uv2dvF_Wrap(qu1,qv1)
 scale = 1.e10
 qduv = qduv*scale
 print(qduv(40:43,90,90))
 qduvnew = reshape(qduv,(/61,4,181,360/))
 qqduv = dim_avg_n(qduvnew,1)
 printVarSummary(qqduv)
 
 system("rm -rf 2000mo0607div100hpa-day.bin")
 path = "./2000mo0607div100hpa-day.bin"
 fbindirwrite(path,qqduv) 
end

http://blog.sina.com.cn/s/blog_d8f6ec6b0102vcsk.html   -yangcan