计算温度平流

ABCD

通过一天的折腾，发现楼主你的程序应该是不对的，dX，dY理应进行中央差分，但你这里却是一个常数

342901wzx

function  TempAdvection(T,u0,v0,tau,lev,lat,lon) ;计算温度平流
begin
       ntau   = dimsizes(tau)   
       nlon   = dimsizes(lon)
       nlat   = dimsizes(lat)
       nlev   = dimsizes(lev)
       Tadv   = T
       dlon = (lon(2)-lon(1))*0.0174533 ; convert to radians
       dTdX = new((/ntau,nlev,nlat,nlon/), typeof(T), T@_FillValue)

       do nl=0,nlat-1                      ; loop over each latitude
         dX = 6378388.*cos(0.0174533*lat(nl))*dlon  ; constant at this latitude
         dTdX(:,:,nl:nl,:) = center_finite_diff_n (T(:,:,nl:nl,:), dX ,False,0,3)
       end do

       dlat = (lat(2)-lat(1))*0.0174533 ; convert to radians
       dTdY = new((/ntau,nlev,nlat,nlon/), typeof(T), T@_FillValue)
       do ml=0,nlon-1
          dY = 6378388.*dlat ; constant at this latitude
          dTdY(:,:,:,ml:ml) = center_finite_diff_n(T(:,:,:,ml:ml),dY,False,0,2)
       end do
       Tadv  =  -(/ u0*dTdX + v0*dTdY /)  
       Tadv  =   (/ Tadv  * 1e5       /)
       return(Tadv)
end

a123456

ABCD 发表于 2020-4-21 20:51
通过一天的折腾，发现楼主你的程序应该是不对的，dX，dY理应进行中央差分，但你这里却是一个常数

你好，请问你这里弄懂了吗？我也没想懂，想请教一下，拜托！

y-yang

码温度平流

沐..

342901wzx 发表于 2020-5-6 09:32
function  TempAdvection(T,u0,v0,tau,lev,lat,lon) ;计算温度平流
begin
       ntau   = dimsizes(tau) ...

哇 这个是温度平流子程序么   

rainbow12138

请问楼主： 用此代码算温度平流时， 温度需要从南到北排序吗？

dh5002

342901wzx 发表于 2020-5-6 09:32
function  TempAdvection(T,u0,v0,tau,lev,lat,lon) ;计算温度平流
begin
       ntau   = dimsizes(tau) ...

请问一下，tau是什么变量呢

wsx

验证过，和ncl函数advect_variable结果一致，算法可以作为求偏导的参考

林登万

感谢楼主！

duoci01

wsx 发表于 2022-3-22 19:12
验证过，和ncl函数advect_variable结果一致，算法可以作为求偏导的参考

我用两个方法算了一遍，结果对不上啊。以哪个为主。求解答