<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Lena,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I'm CC-ing ncl-talk again. Thanks for providing all your files offline.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I was able to run your script, and got the same results as you. The issue that you have 8700 points, but 8400 of them are duplicates. When you plot this, it's plotting all 8700 values, but since 8400 of them have the same x,y pair, they get plotted in the same location.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You may want to recheck your calculations for these variables.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I've attached a version of your script that draws the XY plot, with all the duplicate points colored.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 9, 2014 at 4:57 PM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Lena,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It is hard to debug this code without knowing what "data1" looks like.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Can you write data1 to a NetCDF file and send that file to me?</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">To write data1 to a NetCDF file, you can add this code to the end of your script:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">fout = addfile("<a href="http://formary.nc" target="_blank">formary.nc</a>","c")</div><div class="gmail_default" style="font-size:small">fout->data1 = data1</div><div class="gmail_extra"><br></div><div class="gmail_extra"><div class="gmail_default" style="font-size:small">and then send me "<a href="http://formary.nc" target="_blank">formary.nc</a>". Hopefully it's not too large and you can just email it to me. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You can do this offline (i.e., don't include <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a> on the response), and then when I response, I will include ncl-talk again.</div><div class="gmail_extra"><br></div><div class="gmail_default" style="font-size:small">Thanks,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"></div><br><div class="gmail_quote"><div><div class="h5">On Wed, Oct 8, 2014 at 2:48 AM, Lena Frey <span dir="ltr"><<a href="mailto:Frey_L@web.de" target="_blank">Frey_L@web.de</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div><div style="font-family:Verdana;font-size:12.0px"><div>Hi all!</div>
<div> </div>
<div>I have made a scatter plot. I get no error message, but in my plots are too less data points.</div>
<div>The data dimensions are time,lat,lon with time=348. I would like to have the data points for all time steps of a specific region in the plot.</div>
<div>I de-seasonalized the data and then used the ndtooned command to reshape the data for the scatter plot.</div>
<div>I get 8700 values, but I can't see all the values in the plot.</div>
<div> </div>
<div>Have I used the wrong command to reshape the data?</div>
<div> </div>
<div>Here you can see the important part of my script. You can find attached the plot.</div>
<div> </div>
<div>
<div><br>
;*********************************************************************************<br>
; remove annual cycle and create monthly anomalies<br>
;*********************************************************************************</div>
<div>cf_ann = rmMonAnnCycTLL(cf)<br>
albedo_ann = rmMonAnnCycTLL(albedo)<br>
</div>
<div>printVarSummary(cf_ann)<br>
printMinMax(cf_ann,0)<br>
printVarSummary(albedo_ann)<br>
printMinMax(albedo_ann,0)</div>
<div> </div>
<div>;********************************************************************************<br>
; de-seasonalized data<br>
;********************************************************************************</div>
<div>cf_season = cf-cf_ann<br>
albedo_season = albedo-albedo_ann</div>
<div><br>
printVarSummary(cf_season)<br>
printMinMax(cf_season,0)<br>
printVarSummary(albedo_season)<br>
printMinMax(albedo_season,0)</div>
<div> </div>
<div>
<div>;********************************************************************************<br>
; reshape<br>
;********************************************************************************</div>
<div>cf_reshape = ndtooned(cf_season)<br>
albedo_reshape = ndtooned(albedo_season)</div>
<div><br>
printVarSummary(cf_reshape)<br>
printMinMax(cf_reshape,0)<br>
printVarSummary(albedo_reshape)<br>
printMinMax(albedo_reshape)</div>
<div> </div>
<div>
<div>;********************************************************************************<br>
; regression line<br>
;********************************************************************************</div>
<div>rc = regline(cf_reshape,albedo_reshape)</div>
<div>print(rc)</div>
<div> </div>
<div>
<div>;********************************************************************************<br>
; create array<br>
;********************************************************************************</div>
<div>x = cf_reshape <br>
y = albedo_reshape</div>
<div>printVarSummary(x)<br>
printVarSummary(y)</div>
<div> </div>
<div><br>
data1 = new((/2,dimsizes(y)/), typeof(y))<br>
data1(0,:) = y<br>
data1(1,:) = rc*(x-rc@xave) + rc@yave</div>
<div> </div>
<div>
<div>;*********************************************************************************<br>
; plotting parameters<br>
;*********************************************************************************</div>
<div>wks = gsn_open_wks("pdf","scatter_Australian2")</div>
<div>res = True<br>
res@gsnMaximize = True</div>
<div>
<div>res@xyMarkLineModes = (/"Markers","Lines"/)<br>
res@xyMarkers = 5<br>
res@xyMarkerColor = "black"<br>
res@xyMarkerSizeF = 0.0025<br>
res@xyDashPatterns = 1<br>
res@xyLineThicknesses = (/1,2/)</div>
<div> </div>
<div>res@tiMainString = "Australian"<br>
res@tiXAxisString = "cloud fraction"<br>
res@tiYAxisString = "albedo"</div>
<div> </div>
<div>
<div>plot = gsn_csm_xy (wks,x,data1,res)</div>
<div><br>
end</div>
</div>
<div> </div>
<div> </div>
</div>
<div> </div>
</div>
<div>Thanks.</div>
<div> </div>
<div>Best regards,</div>
<div>Lena</div>
</div>
<div> </div>
<div> </div>
</div>
<div> </div>
<div> </div>
<div> </div>
<div> </div>
</div>
<div> </div>
</div></div></div><br></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></div>
</blockquote></div><br></div>