cdtcurl文档
cdtcurl计算旋度的z分量对椭圆地球网格向量。
参见:cdtgradient,cdtdivergence和埃克曼。
内容
语法
Cz = cdtcurl(纬度,经度,U, V)
描述
Cz = cdtcurl(纬度,经度,U, V)使用cdtdim估计中的每个网格单元的尺寸纬度、经度网格的网格,然后计算旋度向量U, V。
例子
加载这个示例10米的风u和v。注意,当我们风数据加载时,我们只是以平均风速为这个例子。
文件名=“ERA_Interim_2017.nc”;u =意味着(ncread(文件名,“u10”),3);v =意味着(ncread(文件名,“v10”),3);lat =双(ncread(文件名,“纬度”));朗=双(ncread(文件名,“经”));(纬度、经度)= meshgrid(纬度、经度);
原始数据是在一个网格,从0到360个经度。我宁愿把中间的本初子午线的地图,让我们重定位的网格。这一步不是必须的,但它是一个偏好:
(纬度,经度,u, v) =回到中心位置(纬度,经度,u, v);
这里有风在全球地图:
图earthimage举行在q = quiversc(经度、纬度、u, v,“密度”,75);传奇(q,“风速”)
说明使用cdtcurl函数,我们可以重建图S11.03aTalley et al。第六版的书,描述物理海洋学。
首先使用windstress函数将10米风速矢量风应力在N / m ^ 2:
[Taux, Tauy] = windstress (u, v);
现在使用cdtcurl全球的风应力旋度:
C = cdtcurl(纬度,经度,Taux Tauy);
使用岛面具出网格细胞相应的土地,和后Talley et al .,南半球卷发乘以- 1:
%面具土地:土地=岛(纬度、经度);C(土地)=南;Taux(土地)=南;Tauy(土地)=南;%南半球旋度乘以- 1:C = C *标志(Lat);
现在我们准备创建Talley et al。”年代地图相关的风应力风应力旋度:
图pcolor(经度、纬度、C)阴影平持有在%覆盖各大洲:边界(“国家”,“facecolor”rgb (“灰色”),“edgecolor”,“没有”);%覆盖风应力向量:quiversc(经度、纬度、Taux Tauy,“k”,“密度”,100)%飞机从来都不是一个好主意,但是我们会公布的模拟图:colormap (jet (20))%设置颜色限制轴和colorbar:2 caxis ([2] * 1 e) cb = colorbar (“位置”,“southoutside”);包含(cb、的风应力旋度x标志(lat) (\ times10 ^ {7} N / m ^ 3)”)%设置地图范围:轴([10 150 -60 30])%的格式标记标签度:xtickformat (“度”)ytickformat (“度”)
如果你好奇,这是相同的图更合适,发散colormapcmocean(Thyng et al ., 2016)
cmocean落下帷幕
作者信息
这个函数的一部分气候数据为Matlab工具箱。的功能和支持文档是乍得a·格林写的金宝app德州大学奥斯丁分校。