[ncl-talk] Plot Wind Speed significant more than 90%

[???????? ???????]

Excuse me.

My name is Apiwat from University of Phayao

I want to plot map significant same to fig 1.  This fig 1 is plot Light and Dark shading are significant at 90% and 95%. But my script (show fig 2.) Contour is not smooth same to fig 1. Can you help me check my script that where script is missing?

Thank you very much.

Apiwat



My Script.

probt = 100.*(1. -ttest(wspost,wspoststd^2,ntimu,wspre,wsprestd^2,ntimv,True,False)) ; all areas whose value > 95 are significant

    probt at _FillValue = -999

    probt  = where(probt.lt.90.,probt at _FillValue,probt)

    copy_VarCoords(wspre,probt)

    pvalue  = dim_avg_Wrap(probt(lat|:,lon|:,month|i_season(1,:)))

    copy_VarCoords(wspre_avg,pvalue)

    pvalue  = where(wsdiff.lt.0,pvalue*(-1),pvalue)

    pvalue at _FillValue = 0.

;=====================================;

; Plot Composite

;=====================================;

    wks = gsn_open_wks("png","TestWindSEA")               ; send graphics to PNG file

    ;gsn_define_colormap(wks,"GMT_no_green")

    res = True

    res at gsnDraw   = False

    res at gsnFrame  = False

    res at cnFillOn = True ;à»Ô´¡ÒÃàµÔÁÊÕ

    res at cnLinesOn    = False   ; »Ô´¡ÒÃáÊ´§ contour

    res at cnFillPalette  = (/"(/0.125490,0.376471,1.000000/)","(/0.125490,0.623529,1.000000/)", \

                                  "(/0.125490,0.749020,1.000000/)","(/0.000000,0.811765,1.000000/)","(/0.164706,1.000000,1.000000/)", \

                                  "(/0.333333,1.000000,1.000000/)","(/0.498039,1.000000,1.000000/)","white",  \

                                  "white","(/1.000000,0.941176,0.000000/)","(/1.000000,0.749020,0.000000/)", \

                                  "(/1.000000,0.658824,0.000000/)","(/1.000000,0.541176,0.000000/)","(/1.000000,0.439216,0.000000/)", \

                                  "(/1.000000,0.301961,0.000000/)","(/1.000000,0.000000,0.000000/)"/)

    res at cnLineLabelsOn  = False

    res at mpLandFillColor = "transparent"

    res at mpGeophysicalLineThicknessF = 3.0

    res at mpCenterLonF    = 180.0

    ;res at tiMainString = "500hPa Wind  1986-2014 minus 1957-1985"

    res at gsnRightString = "m/s"

    res at gsnLeftString = "(a) 500hPa Wind"

    res at cnLevelSelectionMode  = "ExplicitLevels"

    ;res at cnLevels  = (/-95,-94,-93,-92,-91,-90,0,90,91,92,93,94,95/)

    ;res at cnMinLevelValF    = 95

    ;res at cnMaxLevelValF    = 99

    ;res at cnLevelSpacingF = 0.5

    res at mpMinLatF            = 15         ; zoom in on map

    res at mpMaxLatF            = 65

    res at mpMinLonF            = 70

    res at mpMaxLonF            = 150

    res at lbLabelBarOn = False     ; »Ô´á¶ºáÊ´§¤èÒ

    plot  = gsn_csm_contour_map_ce(wks,pvalue,res)

    plot = ZeroNegDashLineContour(plot)



[Inline image OWAPstImg360541][Inline image OWAPstImg675377]

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171224/99205982/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Outlook-Inline ima.png
Type: image/png
Size: 143395 bytes
Desc: Outlook-Inline ima.png
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171224/99205982/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Outlook-Inline ima.png
Type: image/png
Size: 101068 bytes
Desc: Outlook-Inline ima.png
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171224/99205982/attachment-0003.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MyScript.ncl
Type: application/octet-stream
Size: 6956 bytes
Desc: MyScript.ncl
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171224/99205982/attachment-0001.obj>