求水汽通量的变化

duan2008

请问各位高手：我的目的是求两个不同时段的夏季 水汽通量之差，分别是2011-1996与1980-1995两个时段：程序如下：
'reinit'
'sdfopen E:\downloads\uwnd.mon.mean.nc'
'sdfopen E:\downloads\vwnd.mon.mean.nc'
'sdfopen E:\downloads\rhum.mon.mean.nc'
'sdfopen E:\downloads\air.mon.mean.nc'
'sdfopen E:\downloads\pres.sfc.mon.mean.nc'
'set dfile 1'
'set lon 110 170'
'set lat 20 60'
'set z 1 8'
'set t 385 768'
'define prs=lev'
  'define g=9.8'
  'define es=(6.112*exp((17.67*(air.4-273.16))/(air.4-29.65)))'
  'define qs=(0.62197*es/(prs-0.378*es))'
  'define qv=rhum.3*qs/100'
  'define uq=qv*uwnd.1*1000'
  'define vq=qv*vwnd.2*1000'
'define uqs=vint(pres.5(lev=0),uq,300)'
  'define vqs=vint(pres.5(lev=0),vq,300)'
   
'set z 1'
'define sq=mag(uqs,vqs)'
'set t 6 8'
'define sq1=ave(sq,t+384,t=576,12)'
'modify sq1 seasonal'
'define sq2=sum(sq1,t=390,t=392)'
'set t 6 8'
'define sq3=ave(sq,t+576,t=768,12)'
'modify sq3 seasonal'
'define sq4=sum(sq3,t=582,t=585)'
'define sq5=(sq4-sq3)'

'd sq5'
'run cbarn'

但是运行结果等值线出现缺失，其中'sdfopen E:\downloads\uwnd.mon.mean.nc'
'sdfopen E:\downloads\vwnd.mon.mean.nc'
'sdfopen E:\downloads\rhum.mon.mean.nc'
'sdfopen E:\downloads\air.mon.mean.nc'
'sdfopen E:\downloads\pres.sfc.mon.mean.nc' 等这些数据我已经确认是下的每一个月的在分析数据，所以想请高手们帮我看一下程序，十分谢谢！！

陌小默的

不是资料问题，这套资料我也用过，应该是你程序的问题你的set dfile 1去掉试试，感觉可能会影响点，上次兰溪给我指出过。而且你运行的时候一点错误没提示？

陌小默的

不是资料问题，这套资料我也用过，应该是你程序的问题，你的set dfile 1去掉试试，感觉可能会影响点，上次兰溪给我指出过。而且你运行的时候一点错误没提示？

duan2008

陌小默的 发表于 2012-5-13 21:55

                               
登录/注册后可看大图

不是资料问题，这套资料我也用过，应该是你程序的问题，你的set dfile 1去掉试试，感觉可能会影响点，上次兰 ...

好像没有错误提示，不知道是什么原因，还请版主帮忙看一下，不知道您当时算的是什么程序，可以参考一下吗，谢谢。

陌小默的

duan2008 发表于 2012-5-14 09:14

                               
登录/注册后可看大图

好像没有错误提示，不知道是什么原因，还请版主帮忙看一下，不知道您当时算的是什么程序，可以参考一下吗 ...

我用的程序已经传到论坛了，跟你的大致差不多哈

chengongyan

你没有设置输出等值线？你最后那个run cbarn是运行分色图色标，你都没有设置输出分色图干嘛要运行这个？

seanmichael

陌小默的 发表于 2012-5-13 21:55
不是资料问题，这套资料我也用过，应该是你程序的问题，你的set dfile 1去掉试试，感觉可能会影响点，上次 ...

您好 可以给我这份资料的连接吗小白本科生不胜感激