[ncl-talk] Complicated nd to 1d to nd

[Dennis Shea]

Is this what you want?

nlat   = dimsizes(lat)
mlon = dimsizes(lon)
nmdl = dimsizes(model)
npct  = dimsizes(percentile)
nst4  = dimsizes(statfour)
ngrd = nlat*mlon
percentileStatsT = new((/ngrd,nmdl,npct,nst4/),"double")

pST_1d = ndtooned(percentileStatsT)
pST_5d = onedtond(pST_1d,(/nlat,mlon,nmdl,npct,nst4/))
delete([/percentileStatsT, pST_1d/] )         ; <== only if memory is a
problem

On Tue, May 11, 2021 at 1:36 AM Buzan, Jonathan via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi NCL-Talk,
>
> My script is taking forever to execute because of loops.
> I am attempting to remedy this situation by using ndtooned, then executing
> the commands, followed by rebuilding the array back to nd.
>
> In my example below percentileStatsT is an array 4d array that is
> percentileStatsT(lat*lon,model,percentile,statistics), but I need to
> convert the lat*lon to lat, lon aka a 5d array:
> nd_T_percentileStats(lat,lon,model,percentile,statistics)
>
> However, I am not sure how to convert lat*lon into lat,lon. Is there a
> command in ncl that will do this?
>
> Cheers,
> -Jonathan
>
>
> Code snippet:
> ex.
> T(time,lat,lon)  ;my data of temperatures
>
> ;;;create new array to store statistics.;;;;;
> model = ispan(0,0,1) ; this will be the unlimited domain
> model!0 = "model"
> model&model = model
> model at long_name = "model"
>
> percentile = ispan(1,99,1)
> percentile!0 = "percentile"
> percentile&percentile = percentile
> percentile at long_name = "percentile"
>
> statfour = ispan(0,3,1)
> statfour!0 = "statfour"
> statfour&statfour = statfour
> statfour at long_name = "Mean, SD, Skewness, Kurtosis"
>
> percentileStatsT  =
> new((/dimsizes(lat)*dimsizes(lon),dimsizes(model),dimsizes(percentile),dimsizes(statfour)/),
> "double")
>
>
> ;;;;; manipulate data to calculate statistics;;;;;;
> T_1d = ndtooned(T(lat|:,lon|:,time|:) ; reorder dimensions and convert to
> 1d
>
> ;;;;; operate on data;;;;;;;
> ;;;;; breaks;;;;;;;
> do k = 0, (dimsizes(lat) * dimsizes(lon) -1)
>         index_time_st = dimsizes(time) * k
>         index_time_end = dimsizes(time) * k + (dimsizes(time) -1)
>         onedT = T_1d(index_time_st:index_time_end)
>         workT_Index = dim_pqsort_n(onedT,2,0)
>         manip_T = onedT(workT_Index)    ; reorder data from lowest to
> highest
>         percentile = new((/99/),"integer”)
>         do l = 0, dimsizes(percentile) -1 ; time is 29200 data points.
> dividing the array into percentiles
>                 starti = 292*l + 146
>                 endi = 292*(l+1) + 146 -1
>                 percentileStatsT(k,:,l,:) = dim_stat4(onedT(starti:endi))
>         end do
> end do
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210511/b8e3db0b/attachment.html>