请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 175|回复: 5

[程序设计] 计算风应力旋度

[复制链接]

新浪微博达人勋

发表于 2024-9-26 18:36:43 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
在家园里看到有说用[curlz,~]=curl(lon,lat,taux,tauy)求出旋度,经纬度要转换成距离。但是我用距离求不出来,直接用lon,lat才能求出值,和另一种方法求出的一样,麻烦大家帮忙看看我的代码哪里有错,谢谢!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2024-9-26 18:38:49 | 显示全部楼层
clc,clear,close all;
lon=ncread('fengchang40_year.nc','lon');
lat=ncread('fengchang40_year.nc','lat');
time=ncread('fengchang40_year.nc','time');
taux40_year=ncread('fengchang40_year.nc','taux40_year');
tauy40_year=ncread('fengchang40_year.nc','tauy40_year');
%方法1
curlmean_year=NaN(40,1);
for i=1:40
    taux1=taux40_year(:,:,i);
    tauy1=tauy40_year(:,:,i);
    [mm,nn]=curl(taux1,tauy1);
    curlmean_year(i)=nanmean(mm(:));
end
% 方法1(距离)
%网格距离
[lon,lat]=meshgrid(lon,lat);
d=2*pi*6371000/360;%一个纬度网格的距离,单位m,(地球半径为6371 km),大约为111km
d_lon=lon*d.*cosd(lat);
d_lat=lat*d;
curlz_year=NaN(40,1);
for i=1:40
    taux1=taux40_year(:,:,i);
    tauy1=tauy40_year(:,:,i);
    [curlz,~]=curl(d_lon',d_lat',taux1,tauy1);
    curlz_year(i)=nanmean(curlz(:));
    disp(i)
end
% 方法2
[dtauy_dx, ~] = gradient(tauy40_year);%求偏导
[~, dtaux_dy] = gradient(taux40_year);
wind_stress_curl =dtauy_dx-dtaux_dy;

curl_year=NaN(40,1);
for i=1:40
    wind_stress_curl11=wind_stress_curl(:,:,i);

   curl_year(i)=nanmean(wind_stress_curl11(:));

end
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2024-9-26 18:42:53 | 显示全部楼层
wywen 发表于 2024-9-26 18:38
clc,clear,close all;
lon=ncread('fengchang40_year.nc','lon');
lat=ncread('fengchang40_year.nc','la ...

求指导{:eb511:}
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2024-9-26 18:57:01 | 显示全部楼层
wywen 发表于 2024-9-26 18:38
clc,clear,close all;
lon=ncread('fengchang40_year.nc','lon');
lat=ncread('fengchang40_year.nc','la ...

经纬度范围比较小
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 7 天前 | 显示全部楼层
计算风应力旋度的程序家园里面是有的,看看这个:http://bbs.06climate.com/forum.p ... hlight=%D0%FD%B6%C8
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 7 天前 | 显示全部楼层
MengXC 发表于 2024-9-27 09:29
计算风应力旋度的程序家园里面是有的,看看这个:http://bbs.06climate.com/forum.php?mod=viewthread&tid= ...

好的,谢谢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表