Re: [ferret_users] number of error bars

[saurabh rathore <rohitsrb2020@xxxxxxxxx>]

Hi Ansley,

sorry to say but this time i couldn't understand the logic to plot error bars. As i want to plot 1 standard deviation error bars but with some separation between two consecutive error bars, just to avoid the crowding of error bars. i didn't understand the use of sst_erb, errup, errdn, as how they will incorporate the 1 standard deviation.

kindly help me in detail. these are the steps i followed till now for plotting the time series and the respective error bars of this time series.

use sla_new11.nc                                                            !! input file
 
let gm=sla_new[x=@ave,y=@ave,d=1]                            !! mean of variable in basin
let gst=(sla_new[x=@ave,y=@ave,l=@var,d=1])^0.5       !!  one standard deviation of the time series   

set viewport a11
plot/thick/axes/set/nolab/nokey/vlimits=0:.1:.01 sla_new[x=@ave,y=@ave,l=@spz:11,d=1]        !! time series

go error_bars polygon/col=black/nolab/nokey/ov gm gst                           !! error bars of 1 standard deviation

regards, saurabh

On Thu, Dec 22, 2016 at 7:36 AM, Ansley C. Manke <ansley.b.manke@xxxxxxxx> wrote:

Hi,
You can define the arguments to the script so that they are on a subset of the axis locations.  For instance, looking at the example at the start of the script,

yes? go/help error_bars
 Full path name: /home/users/tmap/ferret/rhel6-x86_64-linux/go/error_bars.jnl
\cancel mode verify
! error_bars.jnl
! 6/99 *sh*
! 7/01 *acm* update /xlimits to /hlimits

! Description: Overlay error bars on a plot
...

! examples
...

! 2) Use a data set, but create some dummy error numbers
!
! Note that if using transformations to create a 1-dimensional variable
! from a multi-dimensional variable, the limits must be made explicitly.
! sst[x=@ave,y=@ave] would cause this script to fail where it uses the
! `var,return=shape` syntax.

!       yes? use coads_climatology
!       yes? let sst_ave = sst[x=0:360@ave,y=-60:60@ave]
!       yes? let errup = 0.1 + 0*sst_ave
!       yes? let errdn = 0.2 + 0*sst_ave
!
!       yes? plot sst_ave
!       yes? GO error_bars poly/line=2 sst_ave, errup, errdn

to plot half as many error bars, define the arguments to error_bars with a subset of the L index:

yes? use coads_climatology
yes? let sst_ave = sst[x=0:360@ave,y=-60:60@ave]
yes? plot sst_ave

yes? let sst_erb = sst[x=0:360@ave,y=-60:60@ave,L=1:12:2]
yes? let errup = 0.1 + 0*sst_erb
yes? let errdn = 0.2 + 0*sst_erb

yes? GO error_bars poly/line=2 sst_erb, errup, errdn

On 12/20/2016 5:20 PM, saurabh rathore wrote:

Dear Ferret users,

i am struggling with the number of error bars in the plot, to many error bars making the plot to congested. i want to plot few error bars so that i looks nice. i am attaching my figure. please suggest me how to manipulate the number of error bars and their positions.

let gm=sla_new[x=@ave,y=@ave,d=1]
let gst=(sla_new[x=@ave,y=@ave,l=@var,d=1])^0.5
go error_bars polygon/col=black/nolab/nokey gm gst

​regards, saurabh​

--

REGARDS

Saurabh Rathore

Research Scholar (PhD.)

Centre For Oceans, Rivers, Atmosphere & Land Science Technology

Indian Institute Of Technology, Kharagpur

contact :- 91- 8345984434

--

REGARDS

Saurabh Rathore

Research Scholar (PhD.)

Centre For Oceans, Rivers, Atmosphere & Land Science Technology

Indian Institute Of Technology, Kharagpur

contact :- 91- 8345984434