QuickSCAT海表风速和风应力详细下载方法 以及 计算风应力和风应力旋度matlab程序

血嗜子木

zane 发表于 2017-6-1 09:34
我算得没有出错  你什么问题？

谢谢解决了

血嗜子木

zane 发表于 2017-6-1 09:34
我算得没有出错  你什么问题？

谢谢解决了

vuchiyuan

mark~~~~~~~~~~~~~~

moriaty

第一次需要使用风应力资料，谢谢楼主的帮助

shou_ju

很好啊，楼主很棒，受教了

龙井

大神你好，有点问题想问一下
（1）计算风应力旋度的时候，程序里为什么要使用前差后差，中差，没看懂，可以解释一下吗
（2）运行时结果出现了点问题
首先，源程序有点小问题，第75行，原来为

1. long(ii, jj)=lon(jj)*111176*cos(lat(ii)*rad);

复制代码

改为了

1. long(ii, jj)=lon(ii)*111176*cos(lat(jj)*rad);

复制代码

，还有程序的最后一行

1. curlZ(lt, 1)=curlZ(lt, lt-1)

复制代码

，是错误的，我改为这样是否合理

1. curlZ(lt, 1)=curlZ(lt, ln-1);

复制代码

。修改之后，运行显示，出现了inf和-inf，这是为什么呀

zane

龙井 发表于 2017-6-18 17:17
大神你好，有点问题想问一下
（1）计算风应力旋度的时候，程序里为什么要使用前差后差，中差，没看懂，可 ...

第一，我不知道为什么层主为什么会问这个问题，如果不用到计算流体力学里面的前插后插，如何用matlab计算出数学公式里面的差分；程序里面先用的是中心差，后面的四个循环是计算中心差不能计算到的周边插值点，用到的是前后差。
第二，程序第75行是在计算在同一纬度上，不同格点在x方向的距离，所以是cos(lat(ii)*rad)没错。 不过我现在看，怎发现么程序写错了，应该是cos(lat(ii))*rad？  anyway,你可以试试看。
还有最后一行的代码，它跟99行有一样的作用，就是直接把前后差计算不了的格点值，把它等于倒与之相近插值点的值。我反倒觉得应该是curlZ(lt, 1)=curlZ(lt, 2); 不确定。
另外，函数下载原出处：
http://cn.mathworks.com/matlabce ... curl-from-wind-data

龙井

zane 发表于 2017-6-19 23:45
第一，我不知道为什么层主为什么会问这个问题，如果不用到计算流体力学里面的前插后插，如何用matlab计算 ...

有点傻了，那是边界的问题，搞懂了，谢谢

毛毛雨915983428

好东西  好材料

maying0391

有计算风应力的吗？