中央差分算500hPa涡度

lt5201314

    最近涉及到计算涡度，用center_finite_diff_n和uv2vr_cfd分别算了500hPa的涡度，都是使用的中央差分，结果一致，来分享下center_finite_diff_n写的程序，但是对比其他前辈画的，值偏小，不知是否是因为海陆差异，还需要进一步检验。

;**********************************计算涡度***********************************************
dlon = (lon(2)-lon(1)) * 0.0174533
dv500dX = new((/731, 71, 41/), typeof(v500), v500@_FillValue)

do n1=0,70
    dX = 6378388. * cos(0.0174533*lat(n1)) * dlon
    dv500dX(:,n1:n1,:) = center_finite_diff_n(v500(:,n1:n1,:), dX, False, 0, 2)  ; False 边界点采用单边差分
end do

;**********************************************************************************
dlat = (lat(2)-lat(1)) * 0.0174533
du500dY = new((/731, 71, 41/), typeof(u500), u500@_FillValue)

do n2=0,40
    dY = 6378388.  * dlat
    du500dY(:,:,n2:n2) = center_finite_diff_n(u500(:,:,n2:n2), dY, False, 0, 1)
end do

vor500 = dv500dX - du500dY

瑾珂

你好，我想问一下dlon = (lon(2)-lon(1)) * 0.0174533里面的 0.0174533是怎么得到的

lt5201314

瑾珂 发表于 2021-09-23 15:34
你好，我想问一下dlon = (lon(2)-lon(1)) * 0.0174533里面的 0.0174533是怎么得到的

1°=π/180≈0.174533

来自苹果APP客户端

瑾珂

lt5201314 发表于 2021-9-23 16:31
1°=π/180≈0.174533

谢谢，我照着你这个写了一个程序，结果总是出问题

lt5201314

瑾珂 发表于 2021-9-24 16:08
谢谢，我照着你这个写了一个程序，结果总是出问题

你根据报错信息调整看看，图上说的是64行左右维度不相等，可以printVarSummary输出看一下

瑾珂

lt5201314 发表于 2021-9-24 16:40
你根据报错信息调整看看，图上说的是64行左右维度不相等，可以printVarSummary输出看一下

试了一下结果又出现新的问题

瑾珂

总是出问题，不知道哪里需要改进

Lu_uup

瑾珂 发表于 2021-9-24 16:08
谢谢，我照着你这个写了一个程序，结果总是出问题

我也出现了这个问题，求问解决了吗？