[ncl-talk] Question about the built-in function trop_wmo

[Dennis Shea]

Source variable:  air(time, level, lat, lon)

Then you do

tf   = addfile(dir + "air.mon.mean.nc", "r")
temp = tf->air
temp = temp + 273.15                               ; units conversion (from
degC to K)
ta   = dim_avg_n_Wrap(temp, 0)               ; compute climatology

===
Thie climatology is calculated over all months and year yielding

  ta(level,lat,lon)

===

I have *not* run the script or looked at the source  fortran code for the
'trop_wmo' function.

The temporal averaging over all years and months has smeared out the
vertical profiles at each grid point.

I *speculate* that the 'vertical change criterion' used by the function to
detect the tropopause has been 'smoothed away'.

You might try a single year month for testing. Keep in mind that the
monthly averages are also averages.


Good luck


On Wed, Sep 20, 2017 at 8:15 AM, 王娜 <wangna at mail.iap.ac.cn> wrote:

> Hi all,
>
>
> I am trying to calculate the tropopause pressure with the built-in
> function trop_wmo, but got some abnormal values (equal to 85 hPa or -999
> hPa) when using the monthly temperature from NCEP/NCAR reanalysis 1 (17
> levels between 1000 hPa and 10 hPa).  It seems like that this problem is
> related to the vertical resolution of the input temperature. If I use the
> ERA-interim reanalysis data (32 levels between 1000 hPa and 10 hPa), the
> abnormal values are non-existent.  Does the function trop_wmo not apply to
> the renanlysis data with coarse resolution? Any help on how to fix this
> problem?
>
> Attached are my ncl script and input data.
>
>
> Thanks
>
> Wang Na
>
>
> ------------------------------
> *从mail.iap.ac.cn <http://mail.iap.ac.cn>发来的超大附件*
> air.mon.mean.nc
> <https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjqIu3sKBT90UTs0mTbKm3VUL3srL3Zt1Sn2LjkKCTyCu3sKBT98dantETs0mTbKm3VUL3srL3Ztdan7EryDuonECzcqpS4CWSeIkasKAwZEXaWCm3VUL3srL3ZtMfeYpzsFWSnKvSbKyTgk8Ss2maUrUUUUrjnCgTuDn0ka6r47ur4kAr47Xw17KFXF6SZ0LaZt1Tsvk1q3Ur7jAUn5U-7jjUefUI7jOU88Ur7jlUekUrJjYUePU2DjOUefU77jjUerUID0mUK1BOg3=&code=4enzgir3> (305MB,
> 2017年10月05日 22:15 到期)
> 查看下载信息
> <https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjqIu3sKBT90UTs0mTbKm3VUL3srL3Zt1Sn2LjkKCTyCu3sKBT98dantETs0mTbKm3VUL3srL3Ztdan7EryDuonECzcqpS4CWSeIkasKAwZEXaWCm3VUL3srL3ZtMfeYpzsFWSnKvSbKyTgk8Ss2maUrUUUUrjnCgTuDn0ka6r47ur4kAr47Xw17KFXF6SZ0LaZt1Tsvk1q3Ur7jAUn5U-7jjUefUI7jOU88Ur7jlUekUrJjYUePU2DjOUefU77jjUerUID0mUK1BOg3=&code=4enzgir3>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170920/c34449ed/attachment.html>