<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">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">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">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">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><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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>_______________________________________________<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>