<html><head></head><body><div class="ydpa3f1b47cyahoo-style-wrap" style="font-family:bookman old style, new york, times, serif;font-size:16px;"><div><div dir="ltr" data-setdir="false">Dear Ali,</div><div dir="ltr" data-setdir="false"> <span>The community has successfully received your email. I believe that if someone has the answer to your query, they will respond accordingly.</span><br></div><div><br></div><div class="ydpa3f1b47csignature"><div dir="ltr" style="font-family:old times, serif;font-size:16px;"><div><div dir="ltr" style="color:rgb(0, 0, 0);font-family:old times, serif;font-size:16px;"><i><font size="2">-----------------------------------------------------------------------------</font></i></div><div dir="ltr" style="color:rgb(0, 0, 0);font-family:old times, serif;font-size:16px;"><i><font size="2">Sincerely,</font></i></div><div dir="ltr" style="color:rgb(0, 0, 0);font-family:old times, serif;font-size:16px;"><i><font size="2">Ehsan Taghizade<br></font></i></div></div></div></div></div><div><br></div>
        
        </div><div id="yahoo_quoted_6583820445" class="yahoo_quoted">
            <div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
                
                <div>
                        On Tuesday, May 7, 2024 at 04:37:02 PM PDT, ali mughal via ncl-talk <ncl-talk@mailman.ucar.edu> wrote:
                    </div>
                    <div><br></div>
                    <div><br></div>
                
                
                <div><div id="yiv0089781123"><div><div dir="ltr">Hi Everyone<div><br clear="none"></div><div>I sent the following email a while back not sure if it reached the  NCL community. </div><div><br clear="none"></div><div>Kindly let me know if someone has a solution to it. </div><div><br clear="none"></div><div>Also will the support for NCL be decommissioned prior to complete migration to NCL-python ?</div><div><br clear="none"></div><div>Regards <br clear="none"><br clear="none"><div id="yiv0089781123yqt39706" class="yiv0089781123yqt3995451247"><div class="yiv0089781123gmail_quote"><div dir="ltr" class="yiv0089781123gmail_attr">---------- Forwarded message ---------<br clear="none">From: <strong class="yiv0089781123gmail_sendername">ali mughal</strong> <span><<a rel="nofollow noopener noreferrer" shape="rect" ymailto="mailto:mughalali655@gmail.com" target="_blank" href="mailto:mughalali655@gmail.com">mughalali655@gmail.com</a>></span><br clear="none">Date: Mon, Apr 29, 2024 at 2:49 PM<br clear="none">Subject: stat_desp code<br clear="none">To:  <<a rel="nofollow noopener noreferrer" shape="rect" ymailto="mailto:ncl-talk@mailman.ucar.edu" target="_blank" href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>>, Ncl-talk <<a rel="nofollow noopener noreferrer" shape="rect" ymailto="mailto:ncl-talk@ucar.edu" target="_blank" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br clear="none"></div><br clear="none"><br clear="none"><div dir="ltr">Hi <div>My following little code is not providing value of the percentiles probably because of missing values <br clear="none"><br clear="none">dir_agcd = "/data2/WRF_ERA5_out/SWWA/AGCD_regrid_mon/means/"<br clear="none">dir_wrf = "/data2/WRF_ERA5_out/SWWA/wrf_seasmean/"<br clear="none">ens = "R3"<br clear="none">yr_start = "1980"<br clear="none">yr_end = "2019"<br clear="none">nlon = 523<br clear="none">nlat = 346<br clear="none">opt = True<br clear="none">system("mkdir -p " + diro)<br clear="none">f_agcd_tmax = addfile(dir_agcd + "agcd_v1-0-1_" + "tmax_" + yr_start + "_" + yr_end + "_<a rel="nofollow noopener noreferrer" shape="rect" target="_blank" href="http://seasmean.nc">seasmean.nc</a>", "r")<br clear="none">agcd_tmax = f_agcd_tmax->tmax<br clear="none">delete(f_agcd_tmax)<br clear="none">f_wrf_tasmax = addfile(dir_wrf + "tasmax_" + ens + "_" + yr_start + "_" + yr_end + "_<a rel="nofollow noopener noreferrer" shape="rect" target="_blank" href="http://seasmean.nc">seasmean.nc</a>", "r")<br clear="none">wrf_tasmax = (f_wrf_tasmax->tasmax)<br clear="none">wrf_tasmax = wrf_tasmax-273<br clear="none">delete(f_wrf_tasmax)<br clear="none"><br clear="none">do i = 0, nlat - 1, 1<br clear="none">    do j = 0, nlon - 1, 1<br clear="none">        agcd_tmax_out = agcd_tmax(:, i, j)<br clear="none">        if all(ismissing(agcd_tmax_out)) .eq. False then<br clear="none">            wrf_tasmax_out = wrf_tasmax(:, i, j)<br clear="none">            stat_agcd_tmax = stat_dispersion(agcd_tmax_out, opt)<br clear="none">            stat_wrf_tasmax = stat_dispersion(wrf_tasmax_out, opt)<br clear="none">        else<br clear="none">            delete(agcd_tmax_out)<br clear="none">        end if<br clear="none">    end do<br clear="none">end do<br clear="none">print(stat_agcd_tmax)<br clear="none">print(stat_wrf_tasmax)<br clear="none"></div><div><br clear="none"></div><div>Output </div><div>Variable: stat_agcd_tmax<br clear="none">Type: float<br clear="none">Total Size: 120 bytes<br clear="none">            30 values<br clear="none">Number of Dimensions: 1<br clear="none">Dimensions and sizes:   [30]<br clear="none">Coordinates:<br clear="none">Number Of Attributes: 2<br clear="none">  long_name :   Robust Dispersion Statistics<br clear="none">  _FillValue :  1e+20<br clear="none">(0)     31.07983<br clear="none">(1)     0.768555<br clear="none">(2)     29.00586<br clear="none">(3)     29.96567<br clear="none">(4)     30.07379<br clear="none">(5)     30.28866<br clear="none">(6)     30.6044<br clear="none">(7)     30.78744<br clear="none">(8)     31.18414<br clear="none">(9)     31.38833<br clear="none">(10)    31.58258<br clear="none">(11)    31.85179<br clear="none">(12)    31.96011<br clear="none">(13)    32.05429<br clear="none">(14)    32.88403<br clear="none">(15)    3.878166<br clear="none">(16)    5.046049<br clear="none">(17)    0.7660876<br clear="none">(18)    156<br clear="none">(19)    156<br clear="none">(20)     0<br clear="none">(21)     0<br clear="none">(22)    1e+20<br clear="none">(23)    1e+20<br clear="none">(24)    1e+20<br clear="none">(25)    1e+20<br clear="none">(26)    1e+20<br clear="none">(27)    1e+20<br clear="none">(28)    -0.2533386<br clear="none">(29)    -0.2925659<br clear="none"><br clear="none"><br clear="none">Variable: stat_wrf_tasmax<br clear="none">Type: float<br clear="none">Total Size: 120 bytes<br clear="none">            30 values<br clear="none">Number of Dimensions: 1<br clear="none">Dimensions and sizes:   [30]<br clear="none">Coordinates:<br clear="none">Number Of Attributes: 3<br clear="none">  units :       K<br clear="none">  long_name :   Robust Dispersion Statistics: Daily Maximum Near-Surface Air Temperature<br clear="none">  _FillValue :  1e+20<br clear="none">(0)     25.65089<br clear="none">(1)     0.8997155<br clear="none">(2)     23.26544<br clear="none">(3)     24.28622<br clear="none">(4)     24.36398<br clear="none">(5)     24.49762<br clear="none">(6)     25.16293<br clear="none">(7)     25.38058<br clear="none">(8)     25.80746<br clear="none">(9)     26.04553<br clear="none">(10)    26.27335<br clear="none">(11)    26.43805<br clear="none">(12)    26.66354<br clear="none">(13)    26.73364<br clear="none">(14)    28.01169<br clear="none">(15)    4.746246<br clear="none">(16)    5.275275<br clear="none">(17)    0.8968271<br clear="none">(18)    156<br clear="none">(19)    156<br clear="none">(20)     0<br clear="none">(21)     0<br clear="none">(22)    1e+20<br clear="none">(23)    1e+20<br clear="none">(24)    1e+20<br clear="none">(25)    1e+20<br clear="none">(26)    1e+20<br clear="none">(27)    1e+20<br clear="none">(28)    -0.4017957<br clear="none">(29)    -0.3656675<br clear="none"></div><div><br clear="none"><br clear="none"></div></div>
</div></div></div></div>
</div></div><div class="yqt3995451247" id="yqt35086">_______________________________________________<br clear="none">ncl-talk mailing list<br clear="none"><a shape="rect" ymailto="mailto:ncl-talk@mailman.ucar.edu" href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a><br clear="none">List instructions, subscriber options, unsubscribe:<br clear="none"><a shape="rect" href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br clear="none"></div></div>
            </div>
        </div></body></html>