python wrfout求垂直梯度

赵新宇

请问一下，wrfout变量量例如风场ua，求垂直梯度。我之前的方法是：Dz = np.arange(getvar(fn,'z'),axis=0);du_dz = np.arange(getvar(fn,'ua'),axis=0)/Dz
后来问我师姐，说是要插值到等高面上，然后再求梯度：level = np.arange(1000,20000+500,500) ;dz = 500
但是计算发现，由于地形问题，低层有许多nan值，在计算其他诊断量的时候，导致中高层也出现许多nan值。
所以想请问一下，应该如何计算垂直梯度？？？