<div>Hi <span style="line-height: 1.5;">NCL users,</span></div><div><span style="line-height: 1.5;"><br></span></div><div>While I get no response, I have finished it by myself after one-day continuous try/test.</div><div><br></div><div>What I think useful now is the overlay procedure instead of gsn_coordinates. The basic idea is to draw another plot showing prob only (0-0.05) and then overlay it.</div><div>I paste a little bit codes on how I coped with it here, and if someone think it is helpful he/she can search this question via this mail.</div><div><br></div><div>Thanks anyways ! It is really good to learn more about NCL ploting.</div><div>Jie</div><div><br></div><div><div>  ;=========================For p<0.05, fill using MonoPattern and Monocolor==============================</div><div>  res3 = True                                  ; res3 for probability plot</div><div>  res3@gsnDraw             = False       ; Do not draw plot</div><div>  res3@gsnFrame            = False       ; Do not advance frome</div><div>  res3@cnInfoLabelOn       = False       ; turn off info label</div><div>  res3@cnLinesOn           = False       ; do not draw contour lines</div><div>  res3@cnLineLabelsOn      = False       ; do not draw contour labels</div><div>  res3@lbLabelBarOn        = False       ; do not need lable bar for prob </div><div>  res3@cnLevelSelectionMode = "ExplicitLevels" </div><div>  res3@cnLevels = (/0.,0.05/) </div><div>  res3@cnFillOn = True </div><div>  res3@cnMonoFillPattern = True  ; only one pattern</div><div>  res3@cnFillPattern = 17        ; stipple dot fill pattern</div><div>  res3@cnMonoFillColor = True</div><div>  res3@cnFillColor   = "black"</div><div>  res3@cnFillScaleF  = 0.6     ; default is 1, larger is corse, lower is dense</div><div>  </div><div>  plot3   = gsn_csm_contour(wks,prob,res3)</div><div><br></div><div>  overlay (plot, plot3)</div></div><div><div>.........</div><div><br></div><div style="font-size: 12px;font-family: Arial Narrow;padding:2px 0 2px 0;">------------------ Original ------------------</div><div style="font-size: 12px;background:#efefef;padding:8px;"><div><b>From: </b> "393069045"<393069045@qq.com>;</div><div><b>Date: </b> Mon, Sep 10, 2018 00:40 AM</div><div><b>To: </b> "ncl-talk"<ncl-talk@ucar.edu>;<wbr></div><div></div><div><b>Subject: </b> [ncl-talk] plot the dot on the grids with significant trends</div></div><div><br></div>Hi NCL users,<div><br></div><div>I am trying to plot the spatial distribution of the trends in annual evaporation. My plan is to show trends for all grids but cover those grids with dots if its trends are significant (i.e. p<0.05).</div><div><br></div><div>I calculate rc using the NCL guide and use rc_sig to show only the sigificant trend. I have tested the way to derive rc is correct. I then tried <span style="line-height: 1.5;">the procedure </span><span style="line-height: 1.5;">gsn_coordinates but I got a very big PDF which shows all black. Does any one can help me on this?</span><span style="font-family: verdana, sans-serif; font-size: 13.3333px; line-height: 16px; background-color: rgb(255, 231, 198);"></span></div><div><span style="line-height: 1.5;"><br></span></div><div><span style="line-height: 1.5;">If what I did is not correct, could you please give me some clues for drawing such plots?</span></div><div><span style="line-height: 1.5;"><br></span></div><div><span style="line-height: 1.5;">Thanks so much</span></div><div><span style="line-height: 1.5;"><br></span></div><div><span style="line-height: 1.5;">Jie.</span></div><div><span style="line-height: 1.5;"><br></span></div><div><span style="line-height: 1.5;"><br></span></div><div><span style="line-height: 1.5;"><br></span></div><div>;----ATTACH CODES HERE-----</div><div>;-----------read data and calculate rc------</div><div><div>  Ea                 = ncfile->Ea  ; mm/month</div><div>  Ea_annual      = month_to_annual (Ea, 0)  ; monthly to annaul</div><div>  </div><div>  Ea_annual!0    = "time"</div><div>  Ea_annual&time =(/ispan(1980,2015,1)/)  </div><div><br></div><div>  Ea_annual_0  = Ea_annual (lat|:,lon|:,time|:)  ;re-order to make time as rightmost</div><div>  </div><div>  time         = Ea_annual_0&time </div><div>  rc            = regCoef(time,Ea_annual_0)    ;tends with /plat,plon/</div><div><br></div><div>  rc@long_name = "Trend of annual Ea"     ; trend</div><div>  rc@units         = "(mm yr~S~-1~N~)"  </div><div>  </div><div>  printVarSummary(rc)</div><div>  printMinMax (rc, True)</div><div>  </div><div> ;-----make all data become (plat,plon)------</div><div>  tval = onedtond(rc@tval , dimsizes(rc))</div><div>  df   = onedtond(rc@nptxy, dimsizes(rc)) - 2   ;The FillValue grid will be =0 for nptxy. So df may<0</div><div>  df@_FillValue=-9999</div><div>  df   = where(df.le.0,df@_FillValue,df)</div><div>  </div><div>  b    = tval      ; b must be same size as tval (and df)</div><div>  b    = 0.5</div><div>  </div><div>;----betainc (x,a,b): 0<x<1, a>0,b>0)-----</div><div>;----ony if this value<0.05, means p<0.05----</div><div>  prob = betainc(df/(df+tval^2), df/2.0, b)    </div><div>  copy_VarCoords(Ea_annual_0(:,:,0), prob)</div><div><br></div><div> ;-----only significant rc with p<0.05-------</div><div>  rc_sig = new ((/plat,plon/),float)</div><div>  rc_sig@_FillValue=-9999</div><div>  copy_VarCoords(Ea_annual_0(:,:,0), rc_sig)</div><div>  rc_sig= where(prob.lt.0.05,rc,rc_sig@_FillValue)</div><div><br></div><div>HERE SKIP SOME LINES</div><div></div></div><div><br></div><div>;---plot------------------</div><div><div>   plot=new(1,graphic)</div><div>   plot = gsn_csm_contour_map_overlay(wks,rc,Basin,res,res2)</div><div>  </div><div>   pres                   = True</div><div>   pres@gsMarkerSizeF     = 0.007  ; default is 0.007</div><div>   gsn_coordinates(wks,plot,rc_sig,pres)</div></div><div><br></div><div>end </div></div>