WRF初始场中尺度滤波试验出错求助！

schliezer

yeah... 发表于 2017-2-17 09:41
哦哦哦 谢谢您 那请问您有计算响应曲线吗？就是那个r和拉姆打的关系

计算过，也画了图，但图放不上来。

schliezer

yeah... 发表于 2017-2-17 09:42
还有我想请教下 这个Barnes滤波能实现高通滤波吗？

找到上图的方法了，Barnes高通低通应该都可以吧，主要根据你所要研究的对象的尺度来调整参数。但我觉得高通滤波的效果不是很好，不知道是不是我选的参数大小有问题，滤出来的波动在200-700km，还是太长了。

yeah...

schliezer 发表于 2017-2-17 16:25
找到上图的方法了，Barnes高通低通应该都可以吧，主要根据你所要研究的对象的尺度来调整参数。但我觉得高 ...

谢谢回复啦啦啦 还想请教下这个您是求了R0  还是根据r 和拉姆大 的关系式画出来的呀？方便问下是用什么画的吗？   另外我在读程序的时候 r=distance(-90.0+(j-1)*0.75,-90.0+(k-1)*0.75,(i-1)*0.75,(l-1)*0.75)计算距离的时候  那个为什么是-90.0  不好意思哦 是菜鸟 所以很不懂  谢谢您的指导哈

yeah...

schliezer 发表于 2017-2-17 16:25
找到上图的方法了，Barnes高通低通应该都可以吧，主要根据你所要研究的对象的尺度来调整参数。但我觉得高 ...

请问高通滤波是原始的序列减去低频的吗？

schliezer

yeah... 发表于 2017-2-17 19:13
请问高通滤波是原始的序列减去低频的吗？

高通滤波就是留下高频率，短波长的波啊。

schliezer

yeah... 发表于 2017-2-17 16:48
谢谢回复啦啦啦 还想请教下这个您是求了R0  还是根据r 和拉姆大 的关系式画出来的呀？方便问下是用什么画 ...

function r=SphericalDistance(lon1,lat1,lon2,lat2)
R=6371;
r=R*acos(cos(lat1*pi/180)*cos(lat2*pi/180)*cos((lon1-lon2)*pi/180)+sin(lat1*pi/180)*sin(lat2*pi/180));
end

carlye2002

可以研究看看

carlye2002

感谢楼主，好东西，值得一顶！

yeah...

schliezer 发表于 2017-2-19 08:55
function r=SphericalDistance(lon1,lat1,lon2,lat2)
R=6371;
r=R*acos(cos(lat1*pi/180)*cos(lat2*pi/ ...

谢谢谢谢 我还想请教您一下 我看到的一个程序里面 原始数据经纬度格点数 是41*31 但是最后滤波以后变成了31*21 我想请教下这样滤波的话是会失去一些数据吗

schliezer

yeah... 发表于 2017-3-15 10:13
谢谢谢谢 我还想请教您一下 我看到的一个程序里面 原始数据经纬度格点数 是41*31 但是最后滤波以后变成了 ...

太久了，遗忘比较严重，你查查文献吧