Hi, The error-bars script is just doing the plotting step, given the time steps, and the length of the bars. So you're sending the script these arguments: ! Usage: $1 $2 $3 $4 $5 ! yes? GO error_bars plot_command base errors [err_down] ["transpose"] ! $1 - the plot command with qualifiers (eg PLOT/VS/OVER/LINE=2) ! $2 - underlying points to be marked with error bars ! $3 - error values associated with points (or error limit above) ! $4 - [optional] error limit below (if different that $4) Argument 2 is a function defining points along the line where there will be a bar drawn. Argument 3 is a function defining the plus-and-minus values for the bars at each of the points (or argument 3 is the error above and argument 4 is the error limit below) So argument 3 might be something that varies along the line you're plotting, or in your case it's just a single value, so the bars represent the value plus and minus the standard deviation at each place where you draw it. Plot the entire time series at its original time resolution let gm = sla_new[x=@ave,y=@ave,d=1] !! mean of variable in basin plot gm !! time series Then for instance, at every 10th point, the error bars, plus or minus the value gst. let gst=(sla_new[x=@ave,y=@ave,l=@var,d=1])^0.5 !! one standard deviation of the time series let nt = `gm,return=Lsize` let gm_10 = gm[L=1:`nt`:10] !! The function plotted above, at every 10th point let gm_10_with_std = 0*gm_10 + `gst` !! The standard-deviation value repeated at every 10th time point go error_bars polygon/col=black/nolab/nokey/ov gm_10 `gst`+0*gm_10 !! error bars of 1 standard deviation -Ansley On 12/21/2016 2:26 PM, saurabh rathore
wrote:
|