任意两点间斜剖面做法总结【利用GRADS自带函数】

麦田_smile

river 发表于 2014-5-2 19:43
理论上是可以的吧，不过想起来还是比较复杂啊，感觉fortran实现起来可能更方便一些

是的 比较复杂 fortran的话也不容易吧 毕竟是倾斜的平行四边形

river

麦田_smile 发表于 2014-5-3 09:59
是的 比较复杂 fortran的话也不容易吧 毕竟是倾斜的平行四边形

都不容易，但感觉fortran 要比grads 容易实现，而且好懂，因为grads涉及到插值

dream

river 发表于 2014-5-2 19:55
本来就是，我可没说是什么新的地形之类的······

恩，我那时候连什么是orog都不知道

river

dream 发表于 2014-5-7 21:59
恩，我那时候连什么是orog都不知道

谁光看了这几个字母就知道它是地形文件，那我服了她了

一个不小心

学习一下。最近一直在搜

气象vivi

'reinit'
'open  d:\grads\ctl\2012.ctl'
'set lon 240 270'
'set z 1 21'  
'set zlog on'
'set x 1'
'set y 1'
lon1=-90
lon2=-60
lat1=55
lat2=20
'set t 1'
lon=lon1
'collect 1 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(VVELprs,'lon','lat')'
lon=lon+1
endwhile
'd coll2gr(1,-u)'
'printim d:\grads\pic1\xiepoumian\1.bmp white'
'c'
运行的时候总是报错，没明白是什么原因，求楼主解答。。。

                               
登录/注册后可看大图

气象vivi

'reinit'
'open  d:\grads\ctl\2012.ctl'
'set lon 240 270'
'set z 1 21'  
'set zlog on'
'set x 1'
'set y 1'
lon1=-90
lon2=-60
lat1=55
lat2=20
'set t 1'
lon=lon1
'collect 1 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(VVELprs,'lon','lat')'
lon=lon+1
endwhile
'd coll2gr(1,-u)'
'printim d:\grads\pic1\xiepoumian\1.bmp white'
'c'
运行的时候总是报错，没明白是什么原因，求楼主解答。。。

                               
登录/注册后可看大图

气象vivi

我知道我的问题了，后面的维数没有设。。。

成长的快乐

求下载！点赞！

river

成长的快乐 发表于 2014-5-19 22:57
求下载！点赞！

等你升级到小雨就能下载了，或者可以在论坛的资料站下载，那个没有限制