<div dir="ltr">Dear All,<div><br></div><div>I tried to shade <span style="color:rgb(0,0,0);font-family:Tahoma">where only significant grid points and wrote this script:</span></div><div><span style="color:rgb(0,0,0);font-family:Tahoma"><br></span></div><div><span style="color:rgb(0,0,0);font-family:Tahoma"><br></span></div><div><div style=""><br></div><div style=""><font color="#000000" face="Tahoma">sig2D!0=&quot;lat&quot;</font></div><div style=""><font color="#000000" face="Tahoma">sig2D!1=&quot;lon&quot;</font></div><div style=""><font color="#000000" face="Tahoma">;Coordinate</font></div><div style=""><font color="#000000" face="Tahoma">sig2D&amp;lat=lat</font></div><div style=""><font color="#000000" face="Tahoma">sig2D&amp;lon=lon</font></div><div style=""><br></div><div style=""><font color="#000000" face="Tahoma">alb_clr_der_summer2D!0=&quot;lat&quot;</font></div><div style=""><font color="#000000" face="Tahoma">alb_clr_der_summer2D!1=&quot;lon&quot;</font></div><div style=""><font color="#000000" face="Tahoma">;Coordinate</font></div><div style=""><font color="#000000" face="Tahoma">alb_clr_der_summer2D&amp;lat=lat</font></div><div style=""><font color="#000000" face="Tahoma">alb_clr_der_summer2D&amp;lon=lon</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">a = addfile(&quot;$NCARG_ROOT/lib/ncarg/data/cdf/<a href="http://landsea.nc">landsea.nc</a>&quot;,&quot;r&quot;)</font></div><div style=""><font color="#000000" face="Tahoma">lsdata = a-&gt;LSMASK</font></div><div style=""><font color="#000000" face="Tahoma">lsm  = landsea_mask(lsdata,alb_clr_der_summer2D&amp;lat,alb_clr_der_summer2D&amp;lon)</font></div><div style=""><font color="#000000" face="Tahoma">alb_clr_der_summer2D = mask(alb_clr_der_summer2D,lsm.eq.1,False)</font></div><div style=""><font color="#000000" face="Tahoma">;alb_clr_der_summer2D = mask(alb_clr_der_summer2D,lsm.eq.3,False)</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">lsm  = landsea_mask(lsdata,sig2D&amp;lat,sig2D&amp;lon)</font></div><div style=""><font color="#000000" face="Tahoma">sig2D = mask(sig2D,lsm.eq.1,False)</font></div><div style=""><font color="#000000" face="Tahoma">;sig2D = mask(sig2D,lsm.eq.3,False)</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">print(&quot;Plot&quot;)</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">wks = gsn_open_wks(&quot;eps&quot;,&quot;1st_der_test&quot;)</font></div><div style=""><font color="#000000" face="Tahoma">gsn_define_colormap(wks,&quot;BlWhRe&quot;)   ; choose color</font></div><div style=""><font color="#000000" face="Tahoma">;gsn_define_colormap(wks,&quot;BlAqGrYeOrRe&quot;) </font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</font></div><div style=""><font color="#000000" face="Tahoma">  res                     = True               ; plot mods desired</font></div><div style=""><font color="#000000" face="Tahoma">  res@cnFillOn            = True              ; turn on color   </font></div><div style=""><font color="#000000" face="Tahoma">  res@gsnSpreadColors     = True               ; use full color map</font></div><div style=""><font color="#000000" face="Tahoma">  res@cnLinesOn           = False              ; no contour lines</font></div><div style=""><font color="#000000" face="Tahoma">  res@cnLineLabelsOn      = False              ; no line labels</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">  res@mpMaxLatF = 29.5                         ; specify the plot domain</font></div><div style=""><font color="#000000" face="Tahoma">  res@mpMinLatF = 0.5                         ;                         </font></div><div style=""><font color="#000000" face="Tahoma">  res@mpMinLonF = -44.5                        ;</font></div><div style=""><font color="#000000" face="Tahoma">  res@mpMaxLonF = 15.5     </font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">  res@gsnAddCyclic        = False </font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma"> ; res@cnLevelSelectionMode =  &quot;AutomaticLevels&quot;   ; AutomaticLevels</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">  res@cnLevelSelectionMode =  &quot;ManualLevels&quot;   ; AutomaticLevels</font></div><div style=""><font color="#000000" face="Tahoma">  res@cnMinLevelValF       = -1.              ; min level</font></div><div style=""><font color="#000000" face="Tahoma">  res@cnMaxLevelValF       =  1.               ; max level</font></div><div style=""><font color="#000000" face="Tahoma">  res@cnLevelSpacingF      =  0.1               ; interval</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</font></div><div style=""><font color="#000000" face="Tahoma">res@tiMainString        = &quot;1st der test&quot;</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">res@gsnDraw              = False           ; Do not draw plot</font></div><div style=""><font color="#000000" face="Tahoma">res@gsnFrame             = False           ; Do not advance frame</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">plot= gsn_csm_contour_map_ce(wks,alb_clr_der_summer2D(:,:),res)</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">;delete(wks)</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">; ========================= PLOT 2 ==============================</font></div><div style=""><font color="#000000" face="Tahoma"> res2 = True                            ; res2 probability plots</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">  res2@gsnDraw              = False           ; Do not draw plot</font></div><div style=""><font color="#000000" face="Tahoma">  res2@gsnFrame             = False           ; Do not advance frame</font></div><div style=""><font color="#000000" face="Tahoma">  res2@cnLevelSelectionMode = &quot;ManualLevels&quot; ; set manual contour levels</font></div><div style=""><font color="#000000" face="Tahoma">  res2@cnMinLevelValF      = 0.00        ; set min contour level</font></div><div style=""><font color="#000000" face="Tahoma">  res2@cnMaxLevelValF      = 1.05        ; set max contour level</font></div><div style=""><font color="#000000" face="Tahoma">  res2@cnLevelSpacingF     = 0.05        ; set contour spacing</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">  res2@cnInfoLabelOn       = False       ; turn off info label</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">  res2@cnLinesOn           = False       ; do not draw contour lines</font></div><div style=""><font color="#000000" face="Tahoma">  res2@cnLineLabelsOn      = False       ; do not draw contour labels</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">  res2@cnFillScaleF        = 0.6         ; add extra density</font></div><div style=""><font color="#000000" face="Tahoma">;  delete(alpha@long_name)</font></div><div style=""><font color="#000000" face="Tahoma">;  delete(alpha@units)</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">  plot2   = gsn_csm_contour(wks,gsn_add_cyclic_point(sig2D(:,:)), res2) </font></div><div style=""><font color="#000000" face="Tahoma">  plot2   = ShadeGtContour(plot2, 0.95, 17)  ; shade all areas less than the</font></div><div style=""><font color="#000000" face="Tahoma">                                             ; 0.05 contour level</font></div><div style=""><font color="#000000" face="Tahoma">  overlay (plot, plot2)</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style=""><font color="#000000" face="Tahoma">  draw (plot)</font></div><div style=""><font color="#000000" face="Tahoma">  frame(wks)</font></div><div style=""><font color="#000000" face="Tahoma"><br></font></div><div style="color:rgb(0,0,0);font-family:Tahoma"><br></div></div><div style="color:rgb(0,0,0);font-family:Tahoma"><br></div><div style="color:rgb(0,0,0);font-family:Tahoma"><br></div><div style="color:rgb(0,0,0);font-family:Tahoma">but I got this warning:</div><div style="color:rgb(0,0,0);font-family:Tahoma"><br></div><div style="color:rgb(0,0,0);font-family:Tahoma"><br></div><div style="color:rgb(0,0,0);font-family:Tahoma">







<p class=""><span class="">(0)     gsn_add_cyclic: Warning: The range of your longitude data is not 360.</span></p>
<p class=""><span class="">(0)     You may want to set gsnAddCyclic to False to avoid a warning</span></p>
<p class=""><span class="">(0)     message from the spline function.</span></p>
<p class=""><span class="">warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for X axis failed: consider adjusting trXTensionF value</span></p>
<p class=""><span class="">warning:IrTransInitialize: error creating spline approximation for trXCoordPoints; defaulting to linear</span></p><p class=""><span class=""><br></span></p><p class="">lat ranges from 0.5 to 29.5 and lon ranges from -45.5 to 14.5 and the resolution is 1*1 degree.</p><p class=""><br></p><p class="">Thanks,</p><p class="">Ana</p></div><div style="color:rgb(0,0,0);font-family:Tahoma"><br></div><div style="color:rgb(0,0,0);font-family:Tahoma"><br></div><div style="color:rgb(0,0,0);font-family:Tahoma"><br></div><div style="color:rgb(0,0,0);font-family:Tahoma"><br></div><div style="color:rgb(0,0,0);font-family:Tahoma"><br></div><div style="color:rgb(0,0,0);font-family:Tahoma"><br></div><br style="color:rgb(0,0,0);font-family:Tahoma"></div>