<div dir="ltr"><div><div>What to say? <br><br>NCL's epflux function matches that at :<br> <a href="http://www.esrl.noaa.gov/psd/data/epflux/">http://www.esrl.noaa.gov/psd/data/epflux/</a><br><br></div>That attached plot for 7 Dec 2016 also includes contours of EP-divergence (dashed are negative).<br><br>---<br><br></div>The epflux_1.ncl code contains ....<br><br><pre>show_accel = 1 ; <=== plot acceleration (not EP divergence)<br></pre><pre>[SNIP]<br> if (show_accel .eq. 1) then
dudt@_FillValue = -999.0
dudt(0,:) = dudt@_FillValue ; Hide the 1000 mb level
res_con@cnLevelSpacingF = 5. ; Contour level Spacing
plotvec2 = gsn_contour(wks,dudt,res_con) ; Creates plot for du/dt = div(F)/(a*cos(phi))
overlay(plotvec,plotvec2)
else
EPdiv@_FillValue = -999.0
EPdiv(0,:) = EPdiv@_FillValue ; Hide the 1000 mb level
res_con@cnLevelSpacingF = 200. ; Contour level Spacing
plotvec3 = gsn_contour(wks,EPdiv,res_con) ; Creates plot for div(F)
overlay(plotvec,plotvec3)
end if<br><br></pre><pre>is plotting 'acceleration' (dudt) not EP-divergence.<br><br>====<br></pre><pre>Attached is the script used. Also, for your perusal, <br>is the epflux function currently scheduled to be part <br>of the 6.4.0. Please offer suggestions for code improvements<br>or corrections. <br></pre><pre>The data used may readily be downloaded from<br><br><a href="http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.pressure.html">http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.pressure.html</a><br></pre><pre><a href="http://air.2016.nc">air.2016.nc</a>, <a href="http://uwnd.2016.nc">uwnd.2016.nc</a>, <a href="http://vwnd.2016.nc">vwnd.2016.nc</a><br><br>Regards<br></pre>D<br><div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Dec 9, 2016 at 9:37 AM, <span dir="ltr"><<a href="mailto:hshp@mail.iap.ac.cn" target="_blank">hshp@mail.iap.ac.cn</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear NCL group,<br><br> Based on many previous studies, the basic divergence of EP-flux (especially in the troposphere at mid-latitudes) should be negative. However, the results calculated by your program [<a href="https://www.ncl.ucar.edu/Applications/EPflux.shtml" target="_blank">https://www.ncl.ucar.edu/<wbr>Applications/EPflux.shtml</a>] are inconsistent the previous ones.<br><br> Also, the divergence displayed at the website [<a href="http://www.esrl.noaa.gov/psd/data/epflux/" target="_blank">http://www.esrl.noaa.gov/psd/<wbr>data/epflux/</a>] is mainly positive in the troposphere.<br><br> Have you ever checked or thought about this issue?<br>Looking forward to hearing from you.<br>Best regards.<br><span><br><br>Dr. Shengping He, Research Assistant<br>Climate Change Research Center, Chinese Academy of Sciences<br>Tel: <a href="tel:%2B861082995191" value="+861082995191" target="_blank">+861082995191</a><br>E-mail：<a href="mailto:hshp@mail.iap.ac.cn" target="_blank">hshp@mail.iap.ac.cn</a><br><br></span><br><br><br><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>