[ncl-talk] dim_stat_n calculation over time,lev dimension

[Sri nandini]

Hello ncl-users,

I am computing for skewness and kurtosis over the time and lev 
dimensions in my data. I want to use this output to plot a scatter plot 
of skewness versus kurtosis and i am using the ncl dim_stat_n function 
to compute this.

I am stuck in my code at this function, see the entire code below:

f     = addfile ("Hist_monthly_corrected_pi_trend.nc", "r")
   hist_anom    = f->hist_trend({1192152:1366728},:,:,:)

   hist_anom=hist_anom*100
   printVarSummary(hist_anom)
   dimx = dimsizes(hist_anom)
   ntim = dimx(0)          ; 240
   nens = dimx(1)          ;100
   nlat = dimx(2)           ; 45
   mlon = dimx(3)         ; 90
   nmos = 12
   nyrs = ntim/nmos        ; 20

   hist_ensavg=dim_avg_n(hist_anom,1)
   printVarSummary(hist_ensavg)
   hist_anom1 = hist_anom
   do n=0,nens-1
   hist_anom1(:,n,:,:) = hist_anom(:,n,:,:) - hist_ensavg
   end do
   printVarSummary (hist_anom1)

  do nmo=0,nmos-1
       work := reshape(hist_anom1(nmo::nmos,:,:,:) ,(/nyrs*nens,nlat,mlon/))
           printVarSummary(work)                            ; (2000,45,90)
    end do
      copy_VarCoords(hist_anom1(0,0,:,:), work(0,:,:))
         printVarSummary(work)

   TTstat11 = dim_stat4_n(work, (/1,2/))
    printVarSummary(TTstat11)

This results in :

Variable: TTstat11
Type: float
Total Size: 0 bytes
             0 values
Number of Dimensions: 3
Dimensions and sizes:    [4] x [2000] x [0]
Coordinates:
Number Of Attributes: 1
   _FillValue :    -9.96921e+36

How can i compute this to result in skewness (2000) to plot as a scatter 
plot?

Best

Sri