<div dir="ltr"><div><div><div><div><div><div>One issue is that that when the WMO criteria  for detecting a 'tropopause break' was not found, the function correctly assigns a values of -999.0. This should have been returned as an attribute (_FillValue) of the returned variable. This was not done. A bug report has been submitted: NCL-2730.<br><br>---<br></div>The fix is 'trivial'. Manually set the value after return.<br><br>    ptrop = trop_wmo(ta&level, ta(time|:,lat|:,lon|:,level|:),0,False)<br>    ptrop@long_name = "tropopause"<br>    ptrop@units     = ta&level@units<br>    ptrop@_FillValue= -999.0   <====<br><br>---<br><br></div>More important, I think that the order you did the computations should be altered.<br><br></div>[a] You computed a temperature climatology. Then you computed the tropopause levels  using the temperature climatology.<br><br></div>[b] I suggest that the tropopause height  for each year-month be computed first. Then compute the climatology of the tropopause heights.<br><br></div>I made some changes to your code [ :-) ].<br><br></div><div>[a] calculate_tropopause.ncl_v0: your 'updated' code<br></div>[b] calculate_tropopause.ncl: I think this is the better approach<br><br><div>Note: the differences are not 'visually' significant except for Dec.<br><br></div><div>Good Luck<br></div><div>D<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Mar 1, 2018 at 8:14 AM, 黄涵 <span dir="ltr"><<a href="mailto:hedgehogh@pku.edu.cn" target="_blank">hedgehogh@pku.edu.cn</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi all,<br><br>  I have been calculating tropopause height recently, but I found the results from the function "trop_wmo" is quite strange.  <br>  Attached are the code of using "trop_wmo" to calculate the tropopause height and corresponding images.<br><br>  The dataset is from NCEP2, you can download from: <a href="https://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBSearch.pl?&Dataset=NCEP/DOE+AMIP-II+Reanalysis+(Reanalysis-2)+Monthly+Averages&Variable=Air+Temperature" target="_blank">https://www.esrl.noaa.gov/psd/<wbr>cgi-bin/db_search/DBSearch.pl?<wbr>&Dataset=NCEP/DOE+AMIP-II+<wbr>Reanalysis+(Reanalysis-2)+<wbr>Monthly+Averages&Variable=Air+<wbr>Temperature</a> <br><br>  The version of NCL is 6.4.0<br><br>  Thank you very much.<br><br>Han<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>