[ncl-talk] Fwd: stat_desp code

Ehsan Taghizadeh ehsantaghizadeh at yahoo.com
Sun May 19 01:21:13 MDT 2024


Dear Ali, The community has successfully received your email. I believe that if someone has the answer to your query, they will respond accordingly.

-----------------------------------------------------------------------------Sincerely,Ehsan Taghizade

    On Tuesday, May 7, 2024 at 04:37:02 PM PDT, ali mughal via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:  
 
 Hi Everyone
I sent the following email a while back not sure if it reached the  NCL community. 
Kindly let me know if someone has a solution to it. 
Also will the support for NCL be decommissioned prior to complete migration to NCL-python ?
Regards 

---------- Forwarded message ---------
From: ali mughal <mughalali655 at gmail.com>
Date: Mon, Apr 29, 2024 at 2:49 PM
Subject: stat_desp code
To: <ncl-talk at mailman.ucar.edu>, Ncl-talk <ncl-talk at ucar.edu>


Hi My following little code is not providing value of the percentiles probably because of missing values 

dir_agcd = "/data2/WRF_ERA5_out/SWWA/AGCD_regrid_mon/means/"
dir_wrf = "/data2/WRF_ERA5_out/SWWA/wrf_seasmean/"
ens = "R3"
yr_start = "1980"
yr_end = "2019"
nlon = 523
nlat = 346
opt = True
system("mkdir -p " + diro)
f_agcd_tmax = addfile(dir_agcd + "agcd_v1-0-1_" + "tmax_" + yr_start + "_" + yr_end + "_seasmean.nc", "r")
agcd_tmax = f_agcd_tmax->tmax
delete(f_agcd_tmax)
f_wrf_tasmax = addfile(dir_wrf + "tasmax_" + ens + "_" + yr_start + "_" + yr_end + "_seasmean.nc", "r")
wrf_tasmax = (f_wrf_tasmax->tasmax)
wrf_tasmax = wrf_tasmax-273
delete(f_wrf_tasmax)

do i = 0, nlat - 1, 1
    do j = 0, nlon - 1, 1
        agcd_tmax_out = agcd_tmax(:, i, j)
        if all(ismissing(agcd_tmax_out)) .eq. False then
            wrf_tasmax_out = wrf_tasmax(:, i, j)
            stat_agcd_tmax = stat_dispersion(agcd_tmax_out, opt)
            stat_wrf_tasmax = stat_dispersion(wrf_tasmax_out, opt)
        else
            delete(agcd_tmax_out)
        end if
    end do
end do
print(stat_agcd_tmax)
print(stat_wrf_tasmax)

Output Variable: stat_agcd_tmax
Type: float
Total Size: 120 bytes
            30 values
Number of Dimensions: 1
Dimensions and sizes:   [30]
Coordinates:
Number Of Attributes: 2
  long_name :   Robust Dispersion Statistics
  _FillValue :  1e+20
(0)     31.07983
(1)     0.768555
(2)     29.00586
(3)     29.96567
(4)     30.07379
(5)     30.28866
(6)     30.6044
(7)     30.78744
(8)     31.18414
(9)     31.38833
(10)    31.58258
(11)    31.85179
(12)    31.96011
(13)    32.05429
(14)    32.88403
(15)    3.878166
(16)    5.046049
(17)    0.7660876
(18)    156
(19)    156
(20)     0
(21)     0
(22)    1e+20
(23)    1e+20
(24)    1e+20
(25)    1e+20
(26)    1e+20
(27)    1e+20
(28)    -0.2533386
(29)    -0.2925659


Variable: stat_wrf_tasmax
Type: float
Total Size: 120 bytes
            30 values
Number of Dimensions: 1
Dimensions and sizes:   [30]
Coordinates:
Number Of Attributes: 3
  units :       K
  long_name :   Robust Dispersion Statistics: Daily Maximum Near-Surface Air Temperature
  _FillValue :  1e+20
(0)     25.65089
(1)     0.8997155
(2)     23.26544
(3)     24.28622
(4)     24.36398
(5)     24.49762
(6)     25.16293
(7)     25.38058
(8)     25.80746
(9)     26.04553
(10)    26.27335
(11)    26.43805
(12)    26.66354
(13)    26.73364
(14)    28.01169
(15)    4.746246
(16)    5.275275
(17)    0.8968271
(18)    156
(19)    156
(20)     0
(21)     0
(22)    1e+20
(23)    1e+20
(24)    1e+20
(25)    1e+20
(26)    1e+20
(27)    1e+20
(28)    -0.4017957
(29)    -0.3656675


_______________________________________________
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/20240519/5b7f0518/attachment.htm>


More information about the ncl-talk mailing list