<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>