<div dir="ltr"><div><div><div>Hi, <br><br></div>I have plotted a panel of two bias maps. The plot, input files and the script is attached. Problem is although my data is only on land points, on the first bias map, the ocean is colored as orange, which should not be the case. And the maps show few values of the contours on the land points irregularly, but I do not want any values to be shown. Can you please help me to do these two corections? <br><br>Thank you for your time and kind efforts.<br><br></div>Best regards,<br></div>Ruksana <br><div><div><div><br><br><br>begin<br><br>; read in data<br><br>f3=addfile(&quot;.........../WAH_APHRO_1998_2007_ANN_sum_precip_bias.nc&quot;, &quot;r&quot;);<br>f4=addfile(&quot;.........../WAH_GPCC_1998_2007_ANN_sum_precip_bias.nc&quot;, &quot;r&quot;);<br><br>ppt3 = f3-&gt;precip<br>time = f3-&gt;time<br>lat = f3-&gt;lat<br>lon = f3-&gt;lon<br><br>ppt4 = f4-&gt;precip<br>time = f4-&gt;time<br>lat = f4-&gt;lat<br>lon = f4-&gt;lon<br><br>wks = gsn_open_wks(&quot;eps&quot;,&quot;ppt1&quot;)  <br>plot = new(2,graphic)   <br><br>printVarSummary(ppt3)<br>printVarSummary(ppt4)<br><br> gsn_define_colormap(wks,&quot;percent_11lev&quot;)      ; choose colormap for bias map<br>  ;gsn_reverse_colormap(wks)<br>  res                          = True<br>  res@gsnDraw                  = False<br>  res@gsnFrame                   = False<br>  res@cnInfoLabelOn            = False <br>  res@cnFillOn                 = True                 ; turn on color<br>  res@cnLinesOn                = False             ; turn off contour lines<br>  res@cnLevelSpacingF          = 0.10              ; contour interval<br>  res@cnFillDrawOrder          = &quot;PreDraw&quot;         ; draw contours first<br>  res@lbLabelStride            = 2                 ; stride on label bar<br>  res@gsnSpreadColors          = True              ; use full colormap<br>  res@gsnSpreadColorEnd        = -3                ; -3 don&#39;t use land color<br>  res@gsnAddCyclic             = False             ; regional data <br>  res@gsnStringFontHeightF     = 0.01              <br>  res@tiMainString             = &quot;&quot;<br>  res@gsnRightString           = &quot;&quot;<br>  ;res@mpProjection            = &quot;LambertConformal&quot;<br>  ;res@gsnMaskLambertConformal = True<br>  res@gsnLeftString            = &quot;&quot;<br>  res@mpMaxLatF                = 30;max(lat)  ; zoom in on region<br>  res@mpMinLatF                = 15;min(lat)<br>  res@mpMinLonF                = 83;min(lon)<br>  res@mpMaxLonF                = 97; max(lon)<br>  res@mpCenterLonF             = 10                ; def is zero<br>  res@mpLandFillColor          = &quot;Transparent&quot;<br>  res@cnLevelSelectionMode= &quot;ManualLevels&quot;<br>  <br>  res@cnMinLevelValF      =  -2150.00       ;for bias maps <br>  res@cnMaxLevelValF      =  4450;100<br>  res@cnLevelSpacingF     = 150 ; 150<br> <br>  res@lbLabelBarOn         = False<br>  res@lbLabelAutoStride   = True<br>  ppt3@_FillValue    = 0.0<br>  ppt3@_FillValue    = -999<br><br>plot(0) = gsn_csm_contour_map(wks,ppt3(0,:,:),res)<br>plot(1) = gsn_csm_contour_map(wks,ppt4(0,:,:),res)<br><br>;***********************************************************<br>; create panel<br>;***********************************************************<br><br>resP                          = True <br>resP@gsnFrame                 = False <br>resP@gsnMaximize              = True<br>resP@gsnPanelLabelBar         = True<br>resP@lbLabelFontHeightF       = 0.01<br>gsn_panel(wks,plot,(/2,2/),resP)          ; now draw as one plot <br><br>end<br><br></div></div></div></div>