<div dir="ltr"><div><div>Dear NCL user,<br><br></div>I have plotted the lag correlation between 2 variables. I also want to highlight those parts which are significance at 95 % with black lines, over the attached plot (or as an overlay plot). So it can easily be vizulise that which part of the given region is more significant with black lines or dot points.<br></div><div><br>please provide some suggestion on that.<br></div><div><br></div>I have data of one month with 31 days for each variables. The scipt I used is given below <br><div><br>                <img src="cid:ii_j27vukc10_15bca65d8f87f421" style="margin-right: 0px;" width="227" height="227"><br>​<br><br>;***********************************************<br>  in1  = addfile(&quot;/home/kunal/<a href="http://anomaly_du.nc">anomaly_du.nc</a>&quot;,&quot;r&quot;)<br>  in2  = addfile(&quot;/home/kunal/<a href="http://anomaly_ch.nc">anomaly_ch.nc</a>&quot;,&quot;r&quot;)<br> <br>  <br>  tmp2 = in1-&gt;DUDP001<br>  tmp1 = in2-&gt;ch<br>;************************************************<br>; reorder to get time as right most dimension<br>;***********************************************<br>  ts1 = tmp1(lat|:,lon|:,time|:)<br>  ts2 = tmp2(lat|:,lon|:,time|:)<br>;************************************************<br>; calculate cross correlations<br>;************************************************<br>  maxlag = 4                                 ; set lag<br>; note, the max lag should not be more than N/4<br>  <br>  ccr = esccr(ts2,ts1,maxlag)                   ; calc cross correlations<br><br>; copy meta data and coordinate variables using contributed functions<br>  copy_VarAtts(ts1, ccr)                    <br>  copy_VarCoords_1(ts2,ccr)<br>;************************************************<br>; plot the correlations<br>;************************************************<br>;************************************************<br> ;---Open workstation and change color map<br>  wks_type = &quot;png&quot;                       <br>  wks        = gsn_open_wks(wks_type, &quot;du-ch&quot;)<br>;************************************************<br><br>  res@mpMinLonF = 59.5<br>  res@mpMaxLonF = 69.8<br>  res@mpMinLatF = 16.5<br>  res@mpMaxLatF = 25.0<br>  <br> plot = gsn_csm_contour_map(wks,ccr(:,:,lag),res) <br><br></div><div>Thank You. <br></div><div><div><br clear="all"><div><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Regards<br></div><div>Kunal<br></div><div><p style="margin:0px;border-collapse:collapse;font-family:tahoma,verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div></div></div></div>
</div></div></div></div>