后使用快捷导航没有帐号？
只需一步, 快速开始
只需一步，快速开始
查看: 6798|回复: 26
ncep资料求湿位涡GS程序运行错误，求高手指点·
在线时间 小时
登录后查看更多精彩内容~
才可以下载或查看，没有帐号？
'sdfopen d:\ncep\air.2011.nc'
'sdfopen d:\ncep\rhum.2011.nc'
'sdfopen d:\ncep\uwnd.2011.nc'
'sdfopen d:\ncep\vwnd.2011.nc'
'set mpdset cnworld'
'set lat 15 45'
'set lon 95 145'
*'set mproj lambert'
'set ylevs 995 900 800 700 600 500 400 300'
*选lev和z都会出现cannot contour grid-all undefined values的错误
*'set lev 995 300'
'set z 2 20'
'set zlog on'
while(i&=979)
'set t 'i''
& &&&'q time'
& &&&res = subwrd(result,3)
& &&&year= substr(res,1,12)
& &&&'q time'
& &&&res = subwrd(result,3)
& &&&hour = substr(res,1,2)
& &&&say hour
& &&&'q time'
& &&&res = subwrd(result,3)
& &&&day=substr(res,4,2)
& &&&say day
& &&&'q time'
& &&&res = subwrd(result,3)
& &&&month=substr(res,6,3)
& &&&say month
*****计算假相当位温*****
'set grads off'
'define prs=lev*100'
'define es=100*(6.112*exp(17.67*（air-273.15)/(air-29.65)))'
'define qs=0.62197*es/(prs-es)'
'define qv=rhum.2*qs/100'
'define e=(prs/100.)*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+.5*log(air)-log(e)-4.775)'
'undefine e'
'define theta=air*pow((100000./prs),(0.-0.28*qv)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
*****计算湿位涡3*****
'define vo=hcurl(uwnd.3,vwnd.4)'
'define f=2*7.292*sin(lat*3.)*0.00001'
'define g=9.8'
'define dp=100*(lev(z-1)-lev(z+1))'
'define deqt=eqt(z-1)-eqt(z+1)'
'define du=uwnd.3(z-1)-uwnd.3(z+1)'
'define dv=vwnd.4(z-1)-vwnd.4(z+1)'
'define dx=2.0**cos(lat*3.)*3.'
'define dy=2.0**3.'
'define dtx=cdiff(eqt,x)'
'define dty=cdiff(eqt,y)'
'define pv1=-g*(vo+f)*deqt/dp'
'define pv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define pv=pv1+pv2'
'set lon 121'
'set lat 21 37'
'set gxout contour'
'set csmooth on'
'set cthick 7'
'set clopts -1 -1 0.10'
'd pv1*10e5'
*'d pv2*10e6'
*'d pv*10e5'
*'disable print'
res=subwrd(result,3)
'draw title 'res'& & 'pv
'printim d:\ncep\pv\pv1.'month'.'day'.'hour'.500.gif white x'
密码修改失败请联系qq：
在线时间 小时
错误就是前面注释的那里么，有没有不保密的数据测试下呢，楼主最好把ctl和出错的详细信息说一下啦
密码修改失败请联系qq：
在线时间 小时
即使你 'set z 2 20'，也不能这样做垂直差分
'define dp=100*(lev(z-1)-lev(z+1))'
'define deqt=eqt(z-1)-eqt(z+1)'
'define du=uwnd.3(z-1)-uwnd.3(z+1)'
'define dv=vwnd.4(z-1)-vwnd.4(z+1)'
必须对垂直层循环，一层一层的计算，不要想一次搞定，如下：
'set gxout fwrite'
'set fwrite mpv_ver_'year'.grd'
while(zz&=3)
level=1000-zz*25
lev1=1000-(zz-1)*25
lev2=1000-(zz+1)*25
'set lev 'level''
'define vo=hcurl(UGRD.3,VGRD.4)'
'define f=2*7.292*sin(lat*3.)*0.00001'
'define g=9.8'
'define dp=5000.0'
'define dthse=thse(lev='lev1')-thse(lev='lev2')'
'define du=UGRD.3(lev='lev1')-UGRD.3(lev='lev2')'
'define dv=VGRD.4(lev='lev1')-VGRD.4(lev='lev2')'
'define dx=2.0**cos(lat*3.)*3.'
'define dy=2.0**3.'
'define dtx=cdiff(thse,x)'
'define dty=cdiff(thse,y)'
'define mpv1=-g*(vo+f)*dthse/dp'&&
'define mpv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define mpv=mpv1+mpv2'
while(zz&=17)
level=1000-zz*50
lev1=1000-(zz-1)*50
lev2=1000-(zz+1)*50
'set lev 'level''
'define vo=hcurl(UGRD.3,VGRD.4)'
'define f=2*7.292*sin(lat*3.)*0.00001'
'define g=9.8'
'define dp=10000.0'
'define dthse=thse(lev='lev1')-thse(lev='lev2')'
'define du=UGRD.3(lev='lev1')-UGRD.3(lev='lev2')'
'define dv=VGRD.4(lev='lev1')-VGRD.4(lev='lev2')'
'define dx=2.0**cos(lat*3.)*3.'
'define dy=2.0**3.'
'define dtx=cdiff(thse,x)'
'define dty=cdiff(thse,y)'
'define mpv1=-g*(vo+f)*dthse/dp'&&
'define mpv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define mpv=mpv1+mpv2'
密码修改失败请联系qq：
在线时间 小时
数据是ncep2.5*2.5资料！
密码修改失败请联系qq：
在线时间 小时
mofangbao 发表于
错误就是前面注释的那里么，有没有不保密的数据测试下呢，楼主最好把ctl和出错的详细信息说一下啦
ncep 2.5*2.5资料！
密码修改失败请联系qq：
在线时间 小时
本帖最后由 abd 于
09:32 编辑
请问楼主这个问题解决了吗？麻烦贴出来共享一下好吗？我也在纠结这个湿位涡的计算了....谢谢！
密码修改失败请联系qq：
在线时间 小时
abd 发表于
请问楼主这个问题解决了吗？麻烦贴出来共享一下好吗？我也在纠结这个湿位涡的计算了....谢谢！
没有呢!不知哪位高人能指点一下!!!
密码修改失败请联系qq：
在线时间 小时
abd 发表于
请问楼主这个问题解决了吗？麻烦贴出来共享一下好吗？我也在纠结这个湿位涡的计算了....谢谢！
没有呢!不知哪位高人能指点一下!!!
密码修改失败请联系qq：
在线时间 小时
abd 发表于
请问楼主这个问题解决了吗？麻烦贴出来共享一下好吗？我也在纠结这个湿位涡的计算了....谢谢！
没有呢!不知哪位高人能指点一下!!!
密码修改失败请联系qq：
在线时间 小时
abd 发表于
请问楼主这个问题解决了吗？麻烦贴出来共享一下好吗？我也在纠结这个湿位涡的计算了....谢谢！
没有呢!不知哪位高人能指点一下!!!
密码修改失败请联系qq：手机签到经验翻倍！快来扫一扫！
都说作图用A卡，是这样么？用来P图和CAD的电脑，求高手给看下配置够用不？
1627浏览 / 25回复
配置品牌型号数量当时的单价现在的单价商家数量CPU1?1489? 1489主板1?899? 899内存1?999? 999硬盘1?775? 775显卡1?1899? 1899机箱1?599? 599电源1?399? 399散热器1?109? 109显示器1?2960? 2960合计金额：10128 元
主板可以用微星B85M-E45，显卡换丽台K2000，电源换海韵520
内存用1600频率的就可以，骇客8G 1600*2
主板B85就好了，电源选游戏悍将不怕BOOM？
其实你用不着显卡，不信你先别买显卡就用其他的装起来试试。另外你这 2133 的内存也只能跑在 1600，除非你把 CPU 换成带K的。
靠核显就能作图？ 有点怀疑
好贴必须顶！
那个是苹果用a卡！微软a卡，我去了个～它！貌似xxx认证显卡民用级560ti、760、770…
想要做图快，请用N制专业卡
想要做图快，请用N制专业卡
顶顶顶顶顶
顶顶顶顶顶
顶顶顶顶顶
路过顶楼主了！
顶顶顶顶顶
如果只作图不玩游戏，那么可以考虑作图专业级显卡，A卡有蓝宝石Firepro系列，N卡有丽台Quadro系列。但如果仅仅是简单的CAD和PS，那么其实CPU自带核显就足够了~~~~除非你有更复杂的应用，或者每天工作量非常大，那么需要专业显卡来帮你渲染加速提升工作效率。我亲戚开公司，他们公司的美工，每天都在PS，用的就是CPU核显，完全没问题。至于AutoCAD，据我所知也有不少人在用CPU核显跑，关键看你具体应用了。
顶顶顶顶顶
非专业显卡对作图没多少帮助的
支持键盘翻页 ( 左右 )&
用户名/注册邮箱/注册手机号
其他第三方号登录Photoshop经典教程10000例设计和作图高手_百度文库
您的浏览器Javascript被禁用，需开启后体验完整功能，
赠送免券下载特权
10W篇文档免费专享
部分付费文档8折起
每天抽奖多种福利
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
Photoshop经典教程10000例设计和作图高手
阅读已结束，下载本文需要
想免费下载本文？
定制HR最喜欢的简历
下载文档到电脑，同时保存到云知识，更方便管理
加入VIP
还剩57页未读，
定制HR最喜欢的简历
你可能喜欢关于子窗口GDI+画图和相应OnPaint()消息的问题，求高手详细解答！在线等...
[问题点数：100分，结帖人whereisrxy]
关于子窗口GDI+画图和相应OnPaint()消息的问题，求高手详细解答！在线等...
[问题点数：100分，结帖人whereisrxy]
不显示删除回复
显示所有回复
显示星级回复
显示得分回复
只显示楼主
2012年12月 总版技术专家分月排行榜第一2012年10月 总版技术专家分月排行榜第一2012年9月 总版技术专家分月排行榜第一2012年8月 总版技术专家分月排行榜第一
2012年11月 总版技术专家分月排行榜第三
2012年12月 总版技术专家分月排行榜第一2012年10月 总版技术专家分月排行榜第一2012年9月 总版技术专家分月排行榜第一2012年8月 总版技术专家分月排行榜第一
2012年11月 总版技术专家分月排行榜第三
匿名用户不能发表回复！|初中物理题，入射角反射角法线作图啥的有木有高手解答一下_百度知道
初中物理题，入射角反射角法线作图啥的有木有高手解答一下
我有更好的答案
啊你在不，还有几个题
😂😂
不要跑啊噗
采纳率：78%
来自团队：
为您推荐：
其他类似问题
换一换
回答问题，赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。