<div dir="ltr"><div>Hi,</div><div><br></div><div>I'm trying to reproduce with <i>center_finite_difference</i> function the computation of wind divergence, to compare with <i>uv2dv_<span class="" style="" id=":244.53" tabindex="-1">cfd</span></i> function. However, I'm getting some strange values for divergence.<br></div><div><br></div><div>For <i><span class="" style="" id=":244.55" tabindex="-1">Dv</span>/Dy</i>, I'm using:</div><div><br></div><div><span class="" style="" id=":244.56" tabindex="-1">dv</span>_<span class="" style="" id=":244.57" tabindex="-1">dy</span> = center_finite_diff_n( v850, lat*g2r, False, 0, 0 )/<span class="" style="" id=":244.58" tabindex="-1">rTerra</span><br></div><div><br></div><div>For Du/<span class="" style="" id=":244.59" tabindex="-1">Dx</span>, I'm using:</div><div><br></div><div><i><span class="" style="" id=":244.60" tabindex="-1">dlon</span> = (<span class="" style="" id=":244.61" tabindex="-1">lon</span>(2)-<span class="" style="" id=":244.62" tabindex="-1">lon</span>(1))*g2r<br><span class="" style="" id=":244.63" tabindex="-1">du</span>_<span class="" style="" id=":244.64" tabindex="-1">dx</span> = new( (/<span class="" style="" id=":244.65" tabindex="-1">dimsizes</span>(u850)/), <span class="" style="" id=":244.66" tabindex="-1">typeof</span>( u850 ), u850@_<span class="" style="" id=":244.67" tabindex="-1">FillValue</span> )<br>term3 = new( (/<span class="" style="" id=":244.68" tabindex="-1">dimsizes</span>(v850)/), <span class="" style="" id=":244.69" tabindex="-1">typeof</span>(v850), v850@_<span class="" style="" id=":244.70" tabindex="-1">FillValue</span> )<br>do i = 0, <span class="" style="" id=":244.71" tabindex="-1">dimsizes</span>(lat)-1<br>      <span class="" style="" id=":244.72" tabindex="-1">dx</span> = <span class="" style="" id=":244.73" tabindex="-1">rTerra</span> * cos( g2r*lat(i) ) * <span class="" style="" id=":244.74" tabindex="-1">dlon</span><br>      <span class="" style="" id=":244.75" tabindex="-1">du</span>_<span class="" style="" id=":244.76" tabindex="-1">dx</span>(i,:) = center_finite_diff( u850(i,:), <span class="" style="" id=":244.77" tabindex="-1">dx</span>, True, 0 )</i></div><div><i>      term3(i,:) = ( v850(i,:)/<span class="" style="" id=":244.78" tabindex="-1">rTerra</span> ) * tan( lat(i)*g2r )<br>  end do</i></div><div><i><br></i></div><div>with <i>rTerra</i> being the mean radius of the Earth, and finally, for divergence (using the same equation calculated by <i>uv2dv_cfd</i>):</div><div><br></div><div><i>div = du_dx + dv_dy - term3</i></div><div><br></div><div>With <i>uv2dv_cfd</i>, I'm using:</div><div><br></div><div><i>divV = uv2dv_cfd( u850, v850, lat, lon, 3 )</i></div><div><i><br></i></div><div>Below, a comparison between the results for latitude 90, for both calculations:</div><div><br></div><div>(0)       div manual      div function<br>(0)    -0.435886     8.54325e-07<br>(1)        0.282385     8.54325e-07<br>(2) 0.641552     8.54325e-07<br>(3) 0.640924     8.54325e-07<br>(4) -3.11505     8.54325e-07</div><div>...</div><div>(138)        0.0877819     8.54325e-07<br>(139)      0.806053     8.54325e-07<br>(140)       -2.59076     8.54325e-07<br>(141)       -1.87243     8.54325e-07<br>(142)       2.96093     8.54325e-07<br>(143)        -0.795052     8.54325e-07</div><div><br></div><div>And below, the comparison for latitude 87.5, for both calculations:</div><div><br></div><div>(0)       div manual      div function<br>(0)    2.54982e-06     2.54933e-06<br>(1)      3.44865e-06     3.44799e-06<br>(2)      4.34874e-06     4.34791e-06<br>(3)      1.30406e-06     1.30381e-06<br>(4)      2.92225e-06     2.9217e-06</div><div>...</div><div>(138)      2.74009e-06     2.73957e-06<br>(139)    1.3016e-06     1.30136e-06<br>(140)     3.98753e-06     3.98678e-06<br>(141)    7.39285e-06     7.39144e-06<br>(142)    2.72861e-06     2.72809e-06<br>(143)    2.54953e-06     2.54905e-06</div><div><br></div><div>I've noted this problem only in the poles and seems to be related with <i>cos</i> function in the computation of gradient of zonal wind over longitude, since it doesn't give zero (or a very small value) for <i>cos(90)</i> or <i>cos(-90).</i></div><div><i></i><br></div><div>I'm using NCEP Reanalysis 1 data and NCL 6.4.0.</div><div><br></div><div>Am I doing something wrong in these calculations?</div><div><br></div><div>Thanks,</div><div><br></div><div>Mateus</div><div><br></div><div><br></div><div><i></i></div><div><i></i></div></div>