关于wrf数据画温度平流的问题探讨

supersupergirl

求各位大佬指教，我在用wrf数据画温度平流的时候，参考了家园几位大神的计算公式，改动后脚本如下：
  tc_p = tk+273.166
  copy_VarCoords(tk, tc_p)

  dlon = (y(1,2)-y(1,1))*0.0174533
  dTdX = new((/141,153/),typeof(tc_p),tc_p@_FillValue)
  do nl = 0,141-1
  dX =6378388.*cos(0.0174533*x(nl,1))*dlon
  dTdX(nl:nl,:) = center_finite_diff_n(tc_p(nl:nl,:),dX,True,0,1)
  end do
  dlat = (x(2,1)-x(1,1))*0.0174533
  dTdY = new((/141,153/),typeof(tc_p),tc_p@_FillValue)
  do ml = 0,153-1  
  dY = 6378388.*dlat;
  dTdY(:,ml:ml) = center_finite_diff_n(tc_p(:,ml:ml),dY,True,0,0)
  end do
  advT=-(u_p*dTdX+v_p*dTdY)*100000
  advT=smth9(advT,0.5,0.25,False)

现在的问题是advT的值是正负无穷，查看后得知是由于dY=0导致，请问是什么原因，真的是因为wrf数据每个网格的dx和dy不同吗，但我已经使用了wrf_user_intrp3d插值过了呀。
而且此脚本也可出图，但就像有位大佬画的一样，暖平流里夹杂着冷平流，而且也无法取值。求各位大佬解答！！
另附上各参数值：
Variable: dY
Type: float
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes:   [1]
Coordinates:
(0)      0

advT最大最小值：(0)     min=-inf   max=inf
dX =6378388.*cos(0.0174533*x(nl,1))*dlon    ;[1]
  tcc = tc_p(nl:nl,:)                         ;[south_north | 1] x [west_east | 153]
  dTdX(nl:nl,:) = center_finite_diff_n(tc_p(nl:nl,:),dX,True,0,1)
  dd = dTdX(nl:nl,:)                          ;[1] x [153]
  end do
  dlat = (x(2,1)-x(1,1))*0.0174533
  dTdY = new((/141,153/),typeof(tc_p),tc_p@_FillValue)
  do ml = 0,153-1  
dY = 6378388.*dlat                  ;[1]
tcc1 = tc_p(:,ml:ml)                ; [south_north | 141] x [west_east | 1]
dTdY(:,ml:ml) = center_finite_diff_n(tc_p(:,ml:ml),dY,True,0,0)
dd1 = dTdY(:,ml:ml)                 ;[141] x [1]

以及参考的家园伙伴们的公式，如下：
ncl 画WRF的温度平流，质疑
http://bbs.06climate.com/forum.p ... &fromuid=110106
(出处: 气象家园)

温度平流随时间的剖面图。
http://bbs.06climate.com/forum.p ... &fromuid=110106
(出处: 气象家园)

计算温度平流
http://bbs.06climate.com/forum.p ... &fromuid=110106
(出处: 气象家园)

supersupergirl

如果这种方法不行 那求问各位大佬该用哪种方法计算温度平流呢

supersupergirl

下图为所绘制的单个时刻700hpa温度平流图

supersupergirl

此时艾特下我的大神好友@小其其格

supersupergirl

气象家园真的没人了吗 这些大神们都不搞科研了吗

supersupergirl

@river 跪求大神帮忙

river

给你个公式吧
'define tt=TMPprs-273.16'
******温度平流Temperature advection单位℃/s*******
'define dtx = cdiff(tt,x)'
'define dty = cdiff(tt,y)'
'define dx = cdiff(lon,x)*3.1416/180'
'define dy = cdiff(lat,y)*3.1416/180'
'define Tadv= -1*( (u*dtx)/(cos(lat*3.1416/180)*dx) + v*dty/dy )/6.37e6'

supersupergirl

river 发表于 2019-8-25 14:47
给你个公式吧
'define tt=TMPprs-273.16'
******温度平流Temperature advection单位℃/s*******

感谢大神 我就是按照这个公式来的 还有你这个是grads的公式吧 ncl好像不需要/dy这种吧？现在的问题是好像物理量设置有问题，因为是常数 我看其他帖子上也有这样的疑问 但是最后怎么解决也没说 求大神指教

river

supersupergirl 发表于 2019-8-25 14:51
感谢大神 我就是按照这个公式来的 还有你这个是grads的公式吧 ncl好像不需要/dy这种吧？现在的问题是好像 ...

是GrADS的公式，主要是我不会NCL，你应该请教NCL大神@尽头的尽头 @兰溪之水

supersupergirl

river 发表于 2019-8-25 17:44
是GrADS的公式，主要是我不会NCL，你应该请教NCL大神@尽头的尽头

好的 还是依旧感谢river大神 貌似已经解决了 虽然还在研究中