[ncl-talk] problem probability plotting
Rick Brownrigg
brownrig at ucar.edu
Mon Apr 1 15:06:11 MDT 2019
Hi,
I'm not sure what you find wrong with your plot -- I'm not an atmospheric
scientist, so I don't know by just looking.
However, you have these lines for plotting the markers:
llat = probt&lat
llon = probt&lon
plot2=gsn_add_polymarker(wks,plot1,llon,llat,res3)
I don't see where "probt" is defined anywhere, and I would expect that this
should give you an error. Did you mean simply "prob&lat", "prob%lon"?
Rick
On Tue, Mar 26, 2019 at 4:43 PM sunmin park <mireiyue at gmail.com> wrote:
> Dear NCL users
>
> I try to plot probability>0.95% on a map
> Used regcoef to get trend and betainc to get probability. Attached the
> figure from the following script.
>
> I do appreciate any help!
>
> Best,
> Sun-
> ************************************
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
>
> begin
>
> diri =
> "/Users/spark/03_Drought/Season_data/04_TOTALprecipitaton_season_mmperday_trend_21cen/"
> fili = "PRE_MAMmean4Dmmperday_21CEN.nc"
> filenc = addfile (diri+fili+".nc", "r")
> varname = "pre"
> pme = filenc->$varname$
> time = filenc->time
>
> pmee = dim_avg_n_Wrap(pme,0)
>
> rc = regCoef(time,pmee(lat|:,lon|:,time|:))
> rc at long_name = "regression coefficient"
> rc at units = "mm/day/century"
> copy_VarCoords(pmee(0,:,:),rc)
>
> rcc = rc*100
> copy_VarCoords(rc,rcc)
>
> tval = onedtond(rc at tval , dimsizes(rc))
> df = onedtond(rc at nptxy-2, dimsizes(rc))
> df at _FillValue=-9999
> df = where(df.le.0,df at _FillValue,df)
>
> b=tval
> b=0.5
> prob = betainc(df/(df+tval^2),df/2.0,b)
> prob at long_name = "probability"
> copy_VarCoords(rcc,prob)
>
> do i=0, dimsizes(lat)-1
> do j =0, dimsizes(lon)-1
> if(prob(i,j) .ltt. 0.05) then
> print("lat "+lat(i)+" lon "+lon(j))
> end if
> end do
> end do
>
> wks = gsn_open_wks("eps","PRE_US_trend_MAM")
> gsn_define_colormap(wks,"BlRe") ;WhiteYellowOrangeRed") ;WhiteBlue")
> ; choose colormap
>
> res = True ; plot mods desired
> res at gsnDraw = False ; don't draw yet
> res at gsnFrame = False ; don't advance frame yet
> res at gsnAddCyclic = False
>
> res at cnFillOn = True ; turn on color
> res at gsnSpreadColors = True ; use full color map
> res at cnLinesOn = False ; no contour lines
> res at cnLineLabelsOn = False ; no line labels
>
> res at mpMaxLatF = 60 ; choose subregion
> res at mpMinLatF = 20
> res at mpMaxLonF = -75
> res at mpMinLonF = -150
>
> res at cnLevelSelectionMode = "ManualLevels"
> res at cnMinLevelValF = -0.004 ; min level
> res at cnMaxLevelValF = 0.004 ; max level
> res at cnLevelSpacingF = 0.0001 ; interval
> res at tmXBLabelFontHeightF = 0.014 ; adjust some font heights
> res at tmYLLabelFontHeightF = 0.014
> res at tiMainFontHeightF = 0.022
> res at txFontHeightF = 0.017
>
> res at lbLabelBarOn = False
> res at mpOutlineBoundarySets = "GeophysicalAndUSStates"
> res at mpGeophysicalLineThicknessF = 1.0
> plot1 = gsn_csm_contour_map_ce(wks,rcc,res)
>
> res3 = True
> res3 at cnInfoLabelOn = False
> res3 at cnLinesOn = False
> res3 at cnLineLabelsOn = False
> res3 at cnFillDotSizeF = 0.004
> res3 at cnLevelSelectionMode = "ManualLevels"
> res3 at cnMinLevelValF = 0.00
> res3 at cnMaxLevelValF = 1.05
> res3 at cnLevelSpacingF = 0.05
> llat = probt&lat
> llon = probt&lon
> plot2=gsn_add_polymarker(wks,plot1,llon,llat,res3)
>
> ppres = True
> ppres at gsnPanelLabelBar = True ; common label bar
> gsn_panel(wks,plot1,(/1,1/),ppres)
>
> delete(wks)
> end
> _______________________________________________
> 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/20190401/1638af20/attachment.html>
More information about the ncl-talk
mailing list