求助 水汽通量计算

chincy09

近来需要计算水汽通量，用matlab2010a读取nc文件，读出来八层比湿，气压，17层风速，不晓得之后如何操作，请大侠指点。
代码如下：
%------------ 读取所需nc文件
shum = netcdf.open('G:\NCEP\water vapor need\shum.mon.mean.nc','nc_nowrite');
[varname, xtype, dimids, atts] = netcdf.inqVar(shum,0);
varid = netcdf.inqVarID(shum,varname);
data1 = netcdf.getVar(shum,varid)
whos data1

pres = netcdf.open('G:\NCEP\water vapor need\pres.sfc.mon.mean.nc','nc_nowrite');
[varname, xtype, dimids, atts] = netcdf.inqVar(pres,0);
varid = netcdf.inqVarID(pres,varname);
data2 = netcdf.getVar(pres,varid)
whos data2

uwnd = netcdf.open('G:\NCEP\water vapor need\uwnd.mon.mean.nc','nc_nowrite');
[varname, xtype, dimids, atts] = netcdf.inqVar(uwnd,0);
varid = netcdf.inqVarID(uwnd,varname);
data3 = netcdf.getVar(uwnd,varid)
whos data3

vwnd = netcdf.open('G:\NCEP\water vapor need\vwnd.mon.mean.nc','nc_nowrite');
[varname, xtype, dimids, atts] = netcdf.inqVar(vwnd,0);
varid = netcdf.inqVarID(vwnd,varname);
data4 = netcdf.getVar(vwnd,varid)
whos data4

尽头的尽头

不会matlab。帮忙顶下

飞天猫

#---- Calculate
prs = 700
g = 9.8
es = 6.112*DataMath.Exp(17.67*(air-273.16)/(air-29.65))
qs = 0.62197*es/(prs-0.378*es)
q = qs*rhum/100
qhdivg = DataMath.Hdivg(q*uwnd/g,q*vwnd/g)
qv = rhum*es/100
uv = DataMath.Magnitude(uwnd, vwnd)
uvq = uv*qv/(9.8*1000)
这个是 Meteoroinfo 里的计算 程序     可以参考一下

太阳雨

不会matlab 顶一下

wlzhongouc

找水汽通量的计算公式  http://wenku.baidu.com/view/b4acb50f7cd184254b3535b8.html