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>