<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><div class=""></div><div class="">The red line shows the difference <font face="Menlo" class="">w_cloud-w_envir</font> which is, as expected, positive meaning that updrafts are moister and warmer than the surrounding. </div><div class=""><br class=""></div><div class="">Coming to Python I decided to use a similar approach with numpy masked arrays</div><div class=""><br class=""></div><div class=""><div class=""><font face="Menlo" class="">mse_cloud_3D=np.ma.masked_where( ((w<1) | (qc<1E-04)), mse)</font></div><div class=""><font face="Menlo" class="">mse_cloud=np.ma.mean(mse_cloud_3D, axis=(0,2))</font></div><div class=""><font face="Menlo" class=""><br class=""></font></div><div class=""><font face="Menlo" class="">mse_envir_3D=np.ma.masked_where( ((w>=1) & (qc>=1E-04)), mse)</font></div><div class=""><font face="Menlo" class="">mse_envir=np.ma.mean(mse_envir_3D, axis=(0,2))</font></div></div><div class=""><br class=""></div><div class="">For the average MSE the profile looks almost exactly the same (good!) but the difference instead is completely different, also considering the order of magnitude.</div><div class=""></div></body></html>