我已经计算了美元∂θ_ \ varepsilon /∂z美元“高度”、“纬度”“经度”维度。
从wrf进口numpy np进口xarray如xr进口(getvar)乙= np。阵列(getvar (ncfile,“乙”,timeidx =)) z = np。阵列(getvar (ncfile、“z”timeidx =时间,单位=“公里”))deth = np。梯度(eth轴= 0)dz = np。梯度(z轴= 0)deth_dz = deth / dz dthdz = xr。数据Array(data=deth_dz, dims=['Height', 'Latitude', 'Longitude'])我想画出水平(公里)的深度和厚度dthdz。这有点类似这的问题。但是当我做dthdz.sum(暗=“高度”),我得到非常低的深度值。我只是想知道我做错了。谢谢你！
