[ncl-talk] calculating intra standard deviation

[Sri nandini]

Hello

I want to calculate intra standard deviation between 100 ensemble 
members using this code:

hist_an1=[year|156,ens|100,lat,lon]

ntim=156

nens=100

do ne=0,nens-1                         ; loop over each ensemble member
      do n=0,ntim-1                       ; loop over each year for 
current member
         work = hist_an1(n,ne,:,:)     ; convenience and efficiency
         hist_std(n,ne,:,:) = dim_stddev_n(work,0) ; ensemble std
      end do
    end do
   printVarSummary(hist_std)

What i want is intra annual ensemble Standard deviation in the format: 
[year|156,ens|100,lat,lon] but i run into error below:

fatal:Dimension size mismatch on subscript #2, left-hand and right-hand 
side dimensions do not match

pointing to the line:         hist_std(n,ne,:,:) = dim_stddev_n(work,0) 
; ensemble std

Can someone please help me with this?

Sri