水汽通量直接用u*shum+v*shum行吗?

568679112

'd uwnd*shum.5+vwnd.2*shum.5'

做个霸气的木头

默哥之前有详细的帖子，你在论坛搜搜吧~

陌小默的

我之前想试过，你可以尝试下，两种算法都画个图对比下，也不难，就是改个表达式

568679112

做个霸气的木头 发表于 2012-6-3 14:29

                               
登录/注册后可看大图

默哥之前有详细的帖子，你在论坛搜搜吧~

好  我去瞅瞅他的  3q

568679112

陌小默的 发表于 2012-6-3 15:53

                               
登录/注册后可看大图

我之前想试过，你可以尝试下，两种算法都画个图对比下，也不难，就是改个表达式

汗。。我这个小菜鸟一开始就是想那样画来着，结果没图出来。  然后这样画有图了，发现有图出来了

陌小默的

568679112 发表于 2012-6-3 03:52

                               
登录/注册后可看大图

汗。。我这个小菜鸟一开始就是想那样画来着，结果没图出来。  然后这样画有图了，发现有图出来了

额，好吧

568679112

陌小默的 发表于 2012-6-3 18:32

                               
登录/注册后可看大图

额，好吧

还要请教你下 我第一步这么写对不对？就想求一个时次的水汽通量
'reinit'
'sdfopen f:/teacher_xiao/uwnd.2011.nc'
'sdfopen f:/teacher_xiao/vwnd.2011.nc'
'sdfopen f:/teacher_xiao/shum.2011.nc'
'set fwrite f:/teacher_xiao/quv.grd'
'set gxout fwrite'
'set lon 50 150'
'set lat 0 70'
'set t 770'
nz=1
while(nz<=6)
'set z 'nz''
'define qu=uwnd*shum.5/9.8'
'd qu'
nz=nz+1
endwhile

nz=1
while(nz<=6)
'set z 'nz''
'define qv=vwnd.2*shum.5/9.8'
'd qv'
nz=nz+1
endwhile

nz=1
while(nz<=6)
'set z 'nz''
'define quv=hdivg(uwnd*shum.5,vwnd.2*shum.5)/9.8'
'd quv'
nz=nz+1
endwhile

'disable fwrite'
'reinit'
;
得到的quv.grd的ctl文件能不能这么写？open它时说打不开二进制文件。。。   我现在真的是一头雾水了。。。
dset quv.grd
undef -999.0
title quv
xdef 40 linear 50 2.5
ydef 28 linear 0 2.5
zdef 6  linear 1000 50
tdef 1  linear jan2011 1mo
vars 3
qu 6 99 a
qv 6 99 b
quv 6 99 c
endvars

陌小默的

首先，你ctl中dset quv.grd，quv。grd的路径呢？！！你怎么没写文件路径，这个还是要写吧，不是fortran啊

568679112

陌小默的 发表于 2012-6-3 23:29

                               
登录/注册后可看大图

首先，你ctl中dset quv.grd，quv。grd的路径呢？！！你怎么没写文件路径，这个还是要写吧，不是fortran啊

真是不胜感激涕零啊  把路径写全了能打开了~~~~可下面显示整层水汽通量时又不行了  不知道怎么做了 grads打开下面的gs文件 全国山河一片红啊 全是错~~~~
'reinit'
'open f:/teacher_xiao/quv.ctl'
'sdfopen f:/teacher_xiao/pres.sfc.2011.nc'
'set fwrite f:/teacher_xiao/quv_all.grd'
'set gxout fwrite'
'set gxout shaded'
'set lon 50 140'
'set lat 0 67'
'set grads off'
'set mpdset cnworld'

'set z 1'
'define quall=vint(pres.2,qu,500)*9.8'
'define qvall=vint(pres.2,qv,500)*9.8'
'd quall'
'd qvall'
'd hdivg(quall,qvall)'
'disable fwrite'
'reinit'
;

陌小默的

=你第六行为什么要set gxout shaded.不红才怪