【分享贴】5%pdf间隔的子函数

立方光年的团子

没贡献来赚一点的。。
气象家园坚挺

正文：
有时候大家会觉得contour画图间隔怎么取都不得劲，不能很好凸显分布特征
有一种方法是用作出pdf分布后，按5%的间隔作为contour的间隔
于是就编了一套以5%为间隔的ncl子程序
可自行修改nbin和bin更改n%的pdf

补充：pdf就是[size=13.3333px]probability，概率分布

使用方法：

1. load “get_pdf_5percent.ncl”
   
2. begin
       ;...此处有一堆乱码请忽视，我删不掉，也不知怎么出现的
       colors = (/0,1,2, 13, 18, 34, 50, 66,\
   
3.                 82, 98,114, 128, 0,    \
   
4.                 130,135,146,162,178,\
   
5.                 194,210,226,242,250,255/)
   
6. 
   
7.     ;凑够21个color（此处可以随意更改）
   
8. 
   
9.     c_colrbin = toint(get_pdf_5percent(varname))
   
10.     ncolr = dimsizes(c_colrbin)+3    wks = gsn_open_wks("png",picname)
    
11.     gsn_define_colormap(wks, colors(:ncolr-1))
    
12.     ;...
    
13.     res@cnLevelSelectionMode    = "ExplicitLevels"              
    
14.     res@cnLevels                = c_colrbin
    
15.     ;...
    
16.     plot = gsn_csm_contour_map_ce(wks,varname , res)
    
17. end
    

load “get_pdf_5percent.ncl”<br /> <div align=\"left\">begin</div><div align=\"left\">    ;...此处有一堆乱码请忽视，我删不掉，也不知怎么出现的</div>    colors = (/0,1,2, 13, 18, 34, 50, 66,\<br />                 82, 98,114, 128, 0,    \<br />                 130,135,146,162,178,\<br />                 194,210,226,242,250,255/)<br /> <br />     ;凑够21个color（此处可以随意更改）<br /> <br />     c_colrbin = toint(get_pdf_5percent(varname))<br />     ncolr = dimsizes(c_colrbin)+3    wks = gsn_open_wks("png",picname)<br />     gsn_define_colormap(wks, colors(:ncolr-1))<br />     ;...<br />     res@cnLevelSelectionMode    = "ExplicitLevels"              <br />     res@cnLevels                = c_colrbin<br />     ;...<br />     plot = gsn_csm_contour_map_ce(wks,varname , res)<br /> end<br />

chongzika

加一个效果图更好