[ncl-talk] removing ensemble mean from original data

[Rashed Mahmood]

Please add ncl-talk in subsequent messages...

In that case, my suggestion in previous message should work.

Rashed.

On Fri, Jun 12, 2020 at 1:55 PM Sri nandini <bax8609 at uni-hamburg.de> wrote:

> Thank you.
>
> I wish to remove the ensemble mean from each corresponding time of each
> ensemble.
>
> Sri
> On 12.06.20 13:48, Rashed Mahmood wrote:
>
> If you just want to remove ensemble mean then you could do the following:
>
> x_new      = x
> ens_mean = dim_avg_n_Wrap(x,1)
> do n=0,nens-1
>    x_new(:,n,:,:) = x(:,n,:,:) - ens_mean
> end do
>
> However, please make sure that this is what you want!, it is up to you to
> decide and know what you are trying to accomplish here. Your variable names
> suggest that there could be some confusion...
>
> Cheers,
> Rashed
>
>
>
>
>
> On Fri, Jun 12, 2020 at 12:50 PM Sri nandini via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Hello dear ncl-users,
>>
>> I want to know the best way to remove the ensemble mean from the
>> corresponding time of each ensemble member. Can someone tell me if the
>> following code is correct? i do not get any error but i wanted to be
>> sure of the method.
>>
>> Thanx in advance
>>
>> Sri
>>
>> ;==================================================================
>> ;  Remove the ensemble mean from each ensemble member
>> ;==================================================================
>>     dimx = dimsizes(x)
>>     ntim = dimx(0)          ; 240
>>     nens = dimx(1)          ; 100
>>     nlat = dimx(2)          ; 45
>>     mlon = dimx(3)          ; 90
>>
>>     nmos = 12
>>     nyrs = ntim/nmos         ; 21
>>
>>     printVarSummary(x)       ;[time | 240] x [ens | 100] x [lat | 45] x
>> [lon | 90]
>>     ens_mean=dim_avg_n_Wrap(x,1)
>>     printVarSummary(ens_mean);[time | 240]  x [lat | 45] x [lon | 90]
>>
>>      detrended_x = x
>>      do nt=0,ntim-1
>>        nmo = nt%12
>>       do ne=0,nens-1
>>          detrended_x(nt,ne,:,:) = (/ x(nt,ne,:,:) - ens_mean(nmo,:,:) /)
>> ; detrended array
>>       end do
>>     end do
>>     detrended_x at long_name = “Ensemble mean removed”
>>
>>     printVarSummary(detrended_x)
>>    [time | 240] x [ens | 100] x [lat | 45] x [lon | 90]
>>
>> _______________________________________________
>> 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/20200612/020d5a7c/attachment.html>