<div dir="ltr"><div><div><div>Hi ,<br></div>Could someone point out the error in this script? It plots no values. Many times i checked, the values are proper. I think graphical part has issues. <br><br>;****************************************************<br>;****************************************************<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>;***************************************************<br><br>;***************************************************<br>begin<br>;***************************************************<br>;prec*.25x.25cpc1948-present-usa grid<br>;others*2.25x2.25cpc1948-present-globalgrid<br>f1    = addfile(&quot;tas_Amon_CCSM4_historical_r3i1p1_185001-200512.nc&quot;, &quot;r&quot;)   ; note the &quot;s&quot; of addfile<br>u1 = f1-&gt;tas<br>printVarSummary(u1)                            ; [time| 720]x[lat| 91]x[lon| 180]<br>;;;;;;1961-2005;;;;;;;;;;;;;;;;;<br>ua1=u1(lat|:,lon|:,time|1332:)<br>ua11=u1(time|1332:,lat|:,lon|:)<br>printVarSummary(ua1)                            ; [time| 720]x[lat| 91]x[lon| 180]<br>   ua11      = ua11-273.15<br>   ua11@units= &quot;degC&quot;<br>;return<br>   yyyymm = cd_calendar(ua11&amp;time, -1)<br>   yyyy   = yyyymm/100<br>   dimx = dimsizes(ua11)<br>   ntim = dimx(0)                ; all years and months<br><br>   year  = ispan(yyyy(0), yyyy(ntim-2), 1)<br>print(year)<br>   nyrs  = dimsizes(year)<br>print(nyrs)<br>;return<br>;delete(u)<br>xann1  = month_to_annual(ua11 , 1)                 ; [year| 60]x[lat| 91]x[lon| 180]<br>printVarSummary(xann1)<br>xann1&amp;year  = year<br>rc1           = regCoef(year,xann1(lat|:,lon|:,year|:))<br>copy_VarCoords(ua11(0,:,:), rc1)<br>rc1=rc1*45<br>tval1 = onedtond(rc1@tval , dimsizes(rc1))<br>printVarSummary(tval1)<br>df1   = onedtond(rc1@nptxy, dimsizes(rc1)) - 2<br>printVarSummary(df1)<br>b = 0.5<br>t4=onedtond(0.5,dimsizes(rc1))<br>prob1  = (1-betainc(df1/(df1+tval1^2),df1/2.0,t4))<br>print(prob1)<br>printVarSummary(prob1)<br>copy_VarCoords(rc1,prob1)<br>printVarSummary(prob1)<br>;return<br>print(prob1)<br>print(min(prob1))<br>print(max(prob1))<br>printVarSummary(prob1)<br>;**************************************************<br>; plot parameters<br>;**************************************************<br><br>  wks = gsn_open_wks(&quot;x11&quot;,&quot;NorthernPlains_Significant_AirTemp_1961-2005&quot;)                  ; Open ps file<br>  gsn_define_colormap(wks,&quot;gui_default&quot;)           ; assign colormap<br>plot = new(1,graphic)                       ; create a plot array<br><br>  res                       = True<br>res@gsnDraw             = False           ; don&#39;t draw<br>res@gsnFrame            = False           ; don&#39;t advance frame<br>res@cnInfoLabelOn       = False           ; turn off cn info label<br> res@cnLinesOn       = False                   ; turn off contour lines<br>   res@cnFillOn        = True               ; color plot desired<br>   res@cnLineLabelsOn  = False<br> res@lbLabelBarOn        = False           ; turn off individual cb&#39;s<br>res@txFontHeightF = 0.030                   ; adjust font height of gsnLeftString/gsnRightString titles<br>res@cnLevelSelectionMode = &quot;ExplicitLevels&quot;                     ; set explicit contour levels<br>;  res@cnLevels    = (/ 10,20,30,40,50,60,70,80,90,100/)<br>  res@cnLevels    = (/ 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1/)<br>res@vpWidthF              = 0.90                ; change aspect ratio of plot<br>  res@vpHeightF             = 0.90<br>  res@vpXF                  = 0.05<br>  res@vpYF                  = 0.95<br><br>; map resources for projection<br>res@mpGridAndLimbOn       = False<br>res@gsnTickMarksOn = False<br>;***************************************************<br>res@mpLimitMode           = &quot;LatLon&quot;<br>  res@mpMinLatF             =   40.<br>  res@mpMaxLatF             =   49.<br>  res@mpMinLonF             = 242.<br>  res@mpMaxLonF             =  264.<br> res@mpCenterLonF           = 252<br><br>   res@mpPerimOn              = True                    ; draw box around map<br><br>res@mpGeophysicalLineThicknessF = 3.0<br>res@mpGeophysicalLineColor = &quot;Black&quot;; (/22/)<br>;res@mpNationalLineThicknessF  = 3.0<br>res@mpOutlineBoundarySets = &quot;GeophysicalAndUSStates&quot; ; add state boundaries<br>res@mpNationalLineColor  = res@mpGeophysicalLineColor<br>res@mpUSStateLineThicknessF = 3.0<br>res@mpUSStateLineColor  = res@mpGeophysicalLineColor<br><br>res@gsnAddCyclic            = False<br>; to take advantage of NCL&#39;s automatic labeling, we need to assign som attributes:<br><br><br>plot(0) = gsn_csm_contour_map_ce(wks,prob1,res)<br>;return<br>;************************************************<br>; create panel<br>;************************************************<br>  resP                  = True                   ; modify the panel plot<br>  resP@gsnFrame         = False                  ; don&#39;t advance the frame, so we can use gsn_text_ndc<br>  resP@gsnPanelLabelBar = True                   ; add common colorbar<br>  resP@lbLabelFontHeightF = 0.015                ; set font height of Label Bar labels<br>  resP@gsnPanelBottom   = 0.2                    ; shrink panel plot by setting bottom edge of plot<br> resP@gsnPanelTop      = 0.9                    ; shrink panel plot by setting top edge of plot<br> resP@pmLabelBarWidthF  = 0.7                 ; make thinner<br>  resP@pmLabelBarHeightF = 0.09<br>  resP@gsnPanelYWhiteSpacePercent = 5.           ; increase spacing along Y-axis between panel plots<br>  resP@txString   = &quot;NorthernPlains_AirTemperature_Significant_1961-2005&quot;<br>  gsn_panel(wks,plot,(/1,1/),resP)               ; now draw as one plot<br><br><br><br>end<br><br></div>Thank you,<br></div>Adv<br><div><div><br></div></div></div>