grads计算散度注意事项及其它
2010-01-25 23:14
387 查看
'reinit'
'sdfopen e:/data/ncep/data/uwnd.nc'
'sdfopen e:/data/ncep/data/vwnd.nc'
'set gxout fwrite'
'set fwrite e:/data/ncep/data/divmav.grd'
'set lev 1000'
k=1967
while(k<=1998)
'set lon 10 240'
'set lat 20 50'
'define a=ave(uwnd.1,time=01may'k',time=01sep'k')'
'define b=ave(vwnd.2,time=01may'k',time=01sep'k')'
'define hh=hdivg(a,b)'
'set lon 97.5 102.5'
'set lat 37.5 42.5'
'd hh'
k=k+1
endwhile
'disable fwrite'
;
注意:这是一个求散度的grads程序,蓝色为最终要求的区域(的散度),注意其位置,红色为自行设定的区域,考虑到中央差的特性,这个自行设定的区域必须大于蓝色范围。
此外grads 中time表示具体的年月日时间,t只表示时次。
如
k=1950
while(k<=1999)
'define a=ave(uwnd,time=01may'k',time=01sep'k')
'd a'
k=k+1
endwhile
此处的time不能写成t.
'sdfopen e:/data/ncep/data/uwnd.nc'
'sdfopen e:/data/ncep/data/vwnd.nc'
'set gxout fwrite'
'set fwrite e:/data/ncep/data/divmav.grd'
'set lev 1000'
k=1967
while(k<=1998)
'set lon 10 240'
'set lat 20 50'
'define a=ave(uwnd.1,time=01may'k',time=01sep'k')'
'define b=ave(vwnd.2,time=01may'k',time=01sep'k')'
'define hh=hdivg(a,b)'
'set lon 97.5 102.5'
'set lat 37.5 42.5'
'd hh'
k=k+1
endwhile
'disable fwrite'
;
注意:这是一个求散度的grads程序,蓝色为最终要求的区域(的散度),注意其位置,红色为自行设定的区域,考虑到中央差的特性,这个自行设定的区域必须大于蓝色范围。
此外grads 中time表示具体的年月日时间,t只表示时次。
如
k=1950
while(k<=1999)
'define a=ave(uwnd,time=01may'k',time=01sep'k')
'd a'
k=k+1
endwhile
此处的time不能写成t.
相关文章推荐
- 关于浮点数计算注意事项
- ArcGis 字段计算表达式写法注意事项
- Bitmap那些事之内存占用计算和加载注意事项
- Bitmap那些事之内存占用计算和加载注意事项
- Bitmap那些事之内存占用计算和加载注意事项
- 计算Cell高度的注意事项(我总会出错)
- linux主机load average的概念&&计算过程&&注意事项
- 索引列上计算引起的索引失效及优化措施以及注意事项
- 使用NSAttributedString来计算行高注意事项
- C# 根据当前时间计算周、月、季度、年的换算及注意事项
- Bitmap那些事之内存占用计算和加载注意事项
- iOS不得姐项目--精华模块上拉下拉的注意事项,日期显示,重构子控制器,计算cell的高度(只计算一次),图片帖子的显示
- matlab中画折线/柱状图(横坐标为字符串、及其它注意事项)
- 计算网站有效PV注意事项
- UIDocumentInterRactionController 使用时的注意事项(其它应用打开)
- 计算UCS2编码的长度以及注意事项(C编程实现)
- Bitmap那些事之内存占用计算和加载注意事项
- Bitmap那些事之内存占用计算和加载注意事项
- Bitmap那些事之内存占用计算和加载注意事项