From Anieklal at cas.iitd.ac.in  Thu May  2 03:25:45 2024
From: Anieklal at cas.iitd.ac.in (Anie K Lal)
Date: Thu, 02 May 2024 14:55:45 +0530
Subject: [ncl-talk] Calculation of Cloud Liquid water path from WRF out file
Message-ID: <fa3409af4f5a5d74986d9689ac915ca6@cas.iitd.ac.in>

Hi all

I am tring to get CLWP from WRF out file.

I tried prcwater_dp as follows:

q1 = f1->QCLOUD[:,:,:,:]; kg/kg

dp_a = wrf_user_getvar(f1,"p",-1)
dp1 = dp_a[:,:,:,:]; Pa

q1!0 = "time"
q1!1 = "lev"
q1!2 = "lat"
q1!3 = "lon"

dp1!0 = "time"
dp1!1 = "lev"
dp1!2 = "lat"
dp1!3 = "lon"

lwp_a = prcwater_dp (q1(time|:,lat|:,lon|:,lev|: ), 
dp1(time|:,lat|:,lon|:,lev|: ) ; kg/m2

Here I am getting values like
27.09562 45.50217 61.3766 68.90879 67.90661 66.62698 64.72494 67.85168 
64.0595 62.04931 55.23169 42.21736 23.00407... in kg/m2

I also tried wrf_user_getvar(f1,"pw",-1). Here I am getting values like

5.871577 6.066606 6.3876 6.375202 6.273459 5.973719 6.050189 6.115168 
6.225941 6.275452 6.347611 6.305155 6.381238 ... in kg/m2

But the realistic values are very small compared to what I am getting, 
based on previous literature:
To what extents do urbanization and air pollution affect fog? states 
values ~0-160 g/m2
ACP - Radiative forcing estimates of sulfate aerosol in coupled 
climate-chemistry models with emphasis on the role of the temporal 
variability states values ~0-240 g/m2 and many examples are there. They 
all show clwp in g/m2 which is about 0-300 g/m2 range.

So why I am getting very high values. Are the units in kg/m2 correct for 
the functions 'pw' and prcwater_dp?Or, am I doing something wrong!

Looking forward to any kind of help.

Thank you
Anie

From Anieklal at cas.iitd.ac.in  Thu May  2 08:16:42 2024
From: Anieklal at cas.iitd.ac.in (Anie K Lal)
Date: Thu, 02 May 2024 19:46:42 +0530
Subject: [ncl-talk] Calculation of Cloud Liquid water path from WRF out
 file
In-Reply-To: <fa3409af4f5a5d74986d9689ac915ca6@cas.iitd.ac.in>
References: <fa3409af4f5a5d74986d9689ac915ca6@cas.iitd.ac.in>
Message-ID: <2204081badbfd6292e9c747a0d303d8b@cas.iitd.ac.in>

Do I need to use dpres_hybrid_ccm for finding LWP? If so, how do I find 
hyai, and  hybi from wrf output.

Looking forward to any kind of help.




-------- Original Message --------
Subject: Calculation of Cloud Liquid water path from WRF out file
Date: 02-05-2024 14:55
 From: Anie K Lal <Anieklal at cas.iitd.ac.in>
To: ncl-talk at mailman.ucar.edu

Hi all

I am tring to get CLWP from WRF out file.

I tried prcwater_dp as follows:

q1 = f1->QCLOUD[:,:,:,:]; kg/kg

dp_a = wrf_user_getvar(f1,"p",-1)
dp1 = dp_a[:,:,:,:]; Pa

q1!0 = "time"
q1!1 = "lev"
q1!2 = "lat"
q1!3 = "lon"

dp1!0 = "time"
dp1!1 = "lev"
dp1!2 = "lat"
dp1!3 = "lon"

lwp_a = prcwater_dp (q1(time|:,lat|:,lon|:,lev|: ),
dp1(time|:,lat|:,lon|:,lev|: ) ; kg/m2

Here I am getting values like
27.09562 45.50217 61.3766 68.90879 67.90661 66.62698 64.72494 67.85168
64.0595 62.04931 55.23169 42.21736 23.00407... in kg/m2

I also tried wrf_user_getvar(f1,"pw",-1). Here I am getting values like

5.871577 6.066606 6.3876 6.375202 6.273459 5.973719 6.050189 6.115168
6.225941 6.275452 6.347611 6.305155 6.381238 ... in kg/m2

But the realistic values are very small compared to what I am getting,
based on previous literature:
To what extents do urbanization and air pollution affect fog? states
values ~0-160 g/m2
ACP - Radiative forcing estimates of sulfate aerosol in coupled
climate-chemistry models with emphasis on the role of the temporal
variability states values ~0-240 g/m2 and many examples are there.
They all show clwp in g/m2 which is about 0-300 g/m2 range.

So why I am getting very high values. Are the units in kg/m2 correct
for the functions 'pw' and prcwater_dp?Or, am I doing something wrong!

Looking forward to any kind of help.

Thank you
Anie

From setareh.rahimi at gmail.com  Fri May  3 09:53:56 2024
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Fri, 3 May 2024 19:23:56 +0330
Subject: [ncl-talk] Calculation of Cloud Liquid water path from WRF out
 file
In-Reply-To: <fa3409af4f5a5d74986d9689ac915ca6@cas.iitd.ac.in>
References: <fa3409af4f5a5d74986d9689ac915ca6@cas.iitd.ac.in>
Message-ID: <CAN+DFgyjUxpxW=jZh=XoDPkCDgp-uZzaELF6v_fV_7hDr+6R2g@mail.gmail.com>

Hi,

I think it would be better to ask this question from the WRF-Chem forum.
And please have a look at this link.

https://bb.cgd.ucar.edu/cesm/threads/which-variable-is-the-cloud-liquid-water-path.2465/


On Thu, May 2, 2024 at 12:55?PM Anie K Lal via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi all
>
> I am tring to get CLWP from WRF out file.
>
> I tried prcwater_dp as follows:
>
> q1 = f1->QCLOUD[:,:,:,:]; kg/kg
>
> dp_a = wrf_user_getvar(f1,"p",-1)
> dp1 = dp_a[:,:,:,:]; Pa
>
> q1!0 = "time"
> q1!1 = "lev"
> q1!2 = "lat"
> q1!3 = "lon"
>
> dp1!0 = "time"
> dp1!1 = "lev"
> dp1!2 = "lat"
> dp1!3 = "lon"
>
> lwp_a = prcwater_dp (q1(time|:,lat|:,lon|:,lev|: ),
> dp1(time|:,lat|:,lon|:,lev|: ) ; kg/m2
>
> Here I am getting values like
> 27.09562 45.50217 61.3766 68.90879 67.90661 66.62698 64.72494 67.85168
> 64.0595 62.04931 55.23169 42.21736 23.00407... in kg/m2
>
> I also tried wrf_user_getvar(f1,"pw",-1). Here I am getting values like
>
> 5.871577 6.066606 6.3876 6.375202 6.273459 5.973719 6.050189 6.115168
> 6.225941 6.275452 6.347611 6.305155 6.381238 ... in kg/m2
>
> But the realistic values are very small compared to what I am getting,
> based on previous literature:
> To what extents do urbanization and air pollution affect fog? states
> values ~0-160 g/m2
> ACP - Radiative forcing estimates of sulfate aerosol in coupled
> climate-chemistry models with emphasis on the role of the temporal
> variability states values ~0-240 g/m2 and many examples are there. They
> all show clwp in g/m2 which is about 0-300 g/m2 range.
>
> So why I am getting very high values. Are the units in kg/m2 correct for
> the functions 'pw' and prcwater_dp?Or, am I doing something wrong!
>
> Looking forward to any kind of help.
>
> Thank you
> Anie
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20240503/059cef44/attachment.htm>

From Anieklal at cas.iitd.ac.in  Fri May  3 21:59:18 2024
From: Anieklal at cas.iitd.ac.in (Anie K Lal)
Date: Sat, 04 May 2024 09:29:18 +0530
Subject: [ncl-talk] Calculation of Cloud Liquid water path from WRF out
 file
In-Reply-To: <CAN+DFgyjUxpxW=jZh=XoDPkCDgp-uZzaELF6v_fV_7hDr+6R2g@mail.gmail.com>
References: <fa3409af4f5a5d74986d9689ac915ca6@cas.iitd.ac.in>
 <CAN+DFgyjUxpxW=jZh=XoDPkCDgp-uZzaELF6v_fV_7hDr+6R2g@mail.gmail.com>
Message-ID: <c077c41e162efdcdf47620b7aa191ee8@cas.iitd.ac.in>

Hi

Thank you for your response. I have already asked this in the WRF-Forum. 
Link to my query:
https://forum.mmm.ucar.edu/threads/calculation-of-cloud-liquid-water-path-from-wrf-out-file-using-ncl.14950/

Looking forward to any kind of help or advice.

Thank you




-------- Original Message --------
Subject: Re: [ncl-talk] Calculation of Cloud Liquid water path from WRF 
out file
Date: 03-05-2024 21:23
 From: Setareh Rahimi <setareh.rahimi at gmail.com>
To: Anie K Lal <Anieklal at cas.iitd.ac.in>, ncl-talk at ucar.edu

Hi,

I think it would be better to ask this question from the WRF-Chem
forum. And please have a look at this link.

https://bb.cgd.ucar.edu/cesm/threads/which-variable-is-the-cloud-liquid-water-path.2465/

On Thu, May 2, 2024 at 12:55?PM Anie K Lal via ncl-talk
<ncl-talk at mailman.ucar.edu> wrote:

> Hi all
> 
> I am tring to get CLWP from WRF out file.
> 
> I tried prcwater_dp as follows:
> 
> q1 = f1->QCLOUD[:,:,:,:]; kg/kg
> 
> dp_a = wrf_user_getvar(f1,"p",-1)
> dp1 = dp_a[:,:,:,:]; Pa
> 
> q1!0 = "time"
> q1!1 = "lev"
> q1!2 = "lat"
> q1!3 = "lon"
> 
> dp1!0 = "time"
> dp1!1 = "lev"
> dp1!2 = "lat"
> dp1!3 = "lon"
> 
> lwp_a = prcwater_dp (q1(time|:,lat|:,lon|:,lev|: ),
> dp1(time|:,lat|:,lon|:,lev|: ) ; kg/m2
> 
> Here I am getting values like
> 27.09562 45.50217 61.3766 68.90879 67.90661 66.62698 64.72494
> 67.85168
> 64.0595 62.04931 55.23169 42.21736 23.00407... in kg/m2
> 
> I also tried wrf_user_getvar(f1,"pw",-1). Here I am getting values
> like
> 
> 5.871577 6.066606 6.3876 6.375202 6.273459 5.973719 6.050189
> 6.115168
> 6.225941 6.275452 6.347611 6.305155 6.381238 ... in kg/m2
> 
> But the realistic values are very small compared to what I am
> getting,
> based on previous literature:
> To what extents do urbanization and air pollution affect fog? states
> 
> values ~0-160 g/m2
> ACP - Radiative forcing estimates of sulfate aerosol in coupled
> climate-chemistry models with emphasis on the role of the temporal
> variability states values ~0-240 g/m2 and many examples are there.
> They
> all show clwp in g/m2 which is about 0-300 g/m2 range.
> 
> So why I am getting very high values. Are the units in kg/m2 correct
> for
> the functions 'pw' and prcwater_dp?Or, am I doing something wrong!
> 
> Looking forward to any kind of help.
> 
> Thank you
> Anie
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
S.Rahimi

CAUTION: This email originated from outside of IIT Delhi. Do not click
links or open attachments unless you recognize the sender and know the
content is safe.

From mughalali655 at gmail.com  Tue May  7 17:36:36 2024
From: mughalali655 at gmail.com (ali mughal)
Date: Wed, 8 May 2024 07:36:36 +0800
Subject: [ncl-talk] Fwd: stat_desp code
In-Reply-To: <CAB5RZE1KukShPiE9RbWcneDC_bJqMhbnWDyiYTcr5QCtkcVYhA@mail.gmail.com>
References: <CAB5RZE1KukShPiE9RbWcneDC_bJqMhbnWDyiYTcr5QCtkcVYhA@mail.gmail.com>
Message-ID: <CAB5RZE31dxK8xH+VRXgkWRePqgUQv12njnqKi_Q+OTyyXL+ttQ@mail.gmail.com>

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20240508/c9527bb9/attachment.htm>

From mughalali655 at gmail.com  Tue May  7 17:36:36 2024
From: mughalali655 at gmail.com (ali mughal)
Date: Wed, 8 May 2024 07:36:36 +0800
Subject: [ncl-talk] Fwd: stat_desp code
In-Reply-To: <CAB5RZE1KukShPiE9RbWcneDC_bJqMhbnWDyiYTcr5QCtkcVYhA@mail.gmail.com>
References: <CAB5RZE1KukShPiE9RbWcneDC_bJqMhbnWDyiYTcr5QCtkcVYhA@mail.gmail.com>
Message-ID: <CAB5RZE31dxK8xH+VRXgkWRePqgUQv12njnqKi_Q+OTyyXL+ttQ@mail.gmail.com>

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20240508/c9527bb9/attachment-0001.htm>

From setareh.rahimi at gmail.com  Fri May 10 07:59:28 2024
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Fri, 10 May 2024 17:29:28 +0330
Subject: [ncl-talk] dpres_hybrid_ccm
Message-ID: <CAN+DFgy4fhrybNeoEmbbjaqtrhte3mqBxrFTXmdmtqfO1QujbQ@mail.gmail.com>

Dear all,

I need to use "*dpres_hybrid_ccm"* function available at:
https://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_hybrid_ccm.shtml
.
But I do not know what "*hyai*" and "*hybi*" are.

Could you please explain what they are *and how I can find them from the
WRF-Chem model output*?

Many thanks in advance.
Best wishes,
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20240510/dd03437c/attachment.htm>

From setareh.rahimi at gmail.com  Wed May 15 04:23:48 2024
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Wed, 15 May 2024 13:53:48 +0330
Subject: [ncl-talk] How to plot two data set completely similar to
 each-other ?
Message-ID: <CAN+DFgxrj0TRWNGTpYwkyk0xDeYn+s0bn2F6yxuYT-QNWe9LTg@mail.gmail.com>

Dear NCL users,
I would like to make a comparison between simulated precipitation from
WRF-Chem and GPM data.

I plotted precipitation from the WRF-Chem and the GPM (attached). However,
the plotted precipitation from the GPM is not completely similar to the one
for the WRF-Chem (label bars, and zooming). So I need your help to make the
plotted precipitation from the GPM completely similar to the plotted
precipitation from the WRF-Chem.

Please kindly advise me to sort this issue out.
Thanks.
Kind regards,


-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20240515/e6cba5a1/attachment-0001.htm>
-------------- next part --------------
begin
;
; Open file 
; Read U10 and V10, Cumulus (rinc) and Non-cumulus (rainnc) prc
;

  ;f       = addfile ("wrfout_d01_2021-03-10_13:00:00","r")
  f       = addfile ("wrfout_d01_2021-03-06_19:00:00","r")
 
  rainc  = f->RAINC                 ; (Time, south_north, west_east)
  rainnc  = f->RAINNC  
   
  u10     = f->U10(148:173,:,:)                   ; (Time, south_north, west_east)
  v10     = f->V10(148:173,:,:)         



 u = dim_avg_n_Wrap (u10, 0)
 v = dim_avg_n_Wrap (v10,0) 
 ;-------------------------------------------------------
;---------------------------------------------------------
  times   = wrf_user_getvar(f,"times",-1)
  ntim    = dimsizes(times)          ; # time steps
  print(times)

  ;---Set the two lat/lon corners that we want to zoom in on.
  
  minlat =   34.6
  maxlat =   36.6
  minlon =   49.7
  maxlon =   53.6


  lats = (/ minlat, maxlat /)
  lons = (/ minlon, maxlon /)
  loc  = wrf_user_ll_to_xy(f, lons, lats, True)

; Use NCL operator > to make sure all values >=0.0
; Sum components and assign attributes
;
  rainc  = rainc  > 0.0         
  rainnc  = rainnc > 0.0

  
  rainTot = rainc + rainnc
  
 
 ;printVarSummary(rainTot)
  
      
  rainTot at description = "Total Precipitation"
  rainTot at units       =  rainc at units

  wks = gsn_open_wks("png","prc-13")

           
  
  res                       = True             ; plot mods desired
  res at gsnMaximize           = True             ; maximize size
  res at gsnScalarContour      = True               ; contours desired
  res at gsnLeftString         = "Totall precipitation (mm)"
  res at gsnRightString        = "2021/03/13"

  res at cnFillPalette         = "BlAqGrYeOrReVi200" ; select color map
  res at cnFillOn              = True             ; color plot desired
  ;res at cnFillPalette         = colors           ; define colors for contour plot
  res at cnLinesOn             = False            ; turn off contour lines
  res at cnLineLabelsOn        = False            ; turn off contour labels
  res at cnFillMode            = "AreaFill"     
  ;res at cnLevelSelectionMode  = "ExplicitLevels" ; explicit [unequal] cn levels
  ;res at cnLevels              = (/0,28,30,32,34,36,38,40,42,44,46,48,50,52,54,56,58,60,62,64,66,68,70,72,74,76,78,80,82/)
  ;res at cnLevels              = (/0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20/)  

  res at cnLevelSelectionMode = "ManualLevels"     ; set manual contour levels
  res at cnMinLevelValF       =  0             ; set min contour level
  res at cnMaxLevelValF       = 20  ; set max contour level
  res at cnLevelSpacingF      = 1         ; set contour spacing




  res at vcGlyphStyle          = "WindBarb"
  res at vcRefLengthF          = 0.025            ; ref vec length
  res at vcMinDistanceF        = 0.025            ; larger means sparser
  res at vcWindBarbTickLengthF = 0.4              ; default 0.3
  res at vcRefAnnoOn           = False
  
  res = wrf_map_resources(f,res)
  res at gsnAddCyclic          = False            ; regional data: not cyclic 
  res at tfDoNDCOverlay        = True

  res at mpFillOn                    = False
  res at mpGeophysicalLineColor      = "black"    ; wrf_map_resources uses "gray"
  res at mpUSStateLineColor          = "black"
  res at mpGeophysicalLineThicknessF = 2.0        ; wrf_map_resources uses 0.5
  res at mpUSStateLineThicknessF     = 2.0

;;;----------------------------------------------------------

  res1        = res         ; Make copy of common resource list

  res1 at ZoomIn = True        ; These five resources are required
  res1 at Xstart = loc(0,0)    ; when zooming in on WRF data and
  res1 at Xend   = loc(0,1)    ; keeping the same map projection.
  res1 at Ystart = loc(1,0)
  res1 at Yend   = loc(1,1)


  res1 = wrf_map_resources(f, res1)
  res1 at tfDoNDCOverlay              = True            ; Tells NCL this is a native projection
; res1 at tfDoNDCOverlay              = "NDCViewport"   ; NCL V6.5.0 or later

;---Overwrite some of the resources set by wrf_map_resources.
  res1 at mpUSStateLineColor          = "black"
  res1 at mpNationalLineColor         = "black"
  res1 at mpGeophysicalLineColor      = "black"
  res1 at mpUSStateLineThicknessF     = 2.0
  res1 at mpNationalLineThicknessF    = 2.0
  res1 at mpGeophysicalLineThicknessF = 2.0
;;;---------------------------------------------------------- 


  h = (rainTot(173,:,:) - rainTot(148,:,:)) 
 
  printVarSummary(h)

 plot = gsn_csm_vector_scalar_map(wks,u,v,h(loc(1,0):loc(1,1),loc(0,0):loc(0,1)),res1) 

  shapefile_name1 = "Ostan.shp"
  shapefile_name2 = "City-Poligon6.shp"
  lnres = True
  lnres at gsLineThicknessF = 6.0
  lnres at gsLineColor = "Black"
  shape1 = gsn_add_shapefile_polylines(wks,plot,shapefile_name1,lnres)
  shape2 = gsn_add_shapefile_polylines(wks,plot,shapefile_name2,lnres)



  res at gsnDraw = False
  res at gsnFrame= False


draw(plot)
frame(wks)



end

-------------- next part --------------
A non-text attachment was scrubbed...
Name: GPM-12.png
Type: image/png
Size: 90132 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20240515/e6cba5a1/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: WRF-12.png
Type: image/png
Size: 212887 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20240515/e6cba5a1/attachment-0003.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gpm
Type: application/octet-stream
Size: 1860 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20240515/e6cba5a1/attachment-0001.obj>

From setareh.rahimi at gmail.com  Sat May 18 15:48:26 2024
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Sun, 19 May 2024 01:18:26 +0330
Subject: [ncl-talk] ****Removing some unnecessary descriptions on plot ****
Message-ID: <CAN+DFgzeO59VVyE3F5gWrYCjBtnQ5qmbHhbmZMORsbD6Eay5gQ@mail.gmail.com>

Dear all,

I was using this script:
https://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/LEVELS_INTERP/wrf_PressureLevel.php
But after running the script I noticed some descriptions above the plots,
which made it somehow unpleasant. I want to put the plots next to each
other, but they are not exactly the same size. I want to plot all the
images at the same size so I need to remove those descriptions.
Would you please kindly advise me on how to remove those descriptions?
Thanks
-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20240519/86296302/attachment.htm>

From ehsantaghizadeh at yahoo.com  Sun May 19 01:21:13 2024
From: ehsantaghizadeh at yahoo.com (Ehsan Taghizadeh)
Date: Sun, 19 May 2024 07:21:13 +0000 (UTC)
Subject: [ncl-talk] Fwd: stat_desp code
In-Reply-To: <CAB5RZE31dxK8xH+VRXgkWRePqgUQv12njnqKi_Q+OTyyXL+ttQ@mail.gmail.com>
References: <CAB5RZE1KukShPiE9RbWcneDC_bJqMhbnWDyiYTcr5QCtkcVYhA@mail.gmail.com>
 <CAB5RZE31dxK8xH+VRXgkWRePqgUQv12njnqKi_Q+OTyyXL+ttQ@mail.gmail.com>
Message-ID: <40333386.896260.1716103273143@mail.yahoo.com>

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>

From ehsantaghizadeh at yahoo.com  Sun May 19 02:10:30 2024
From: ehsantaghizadeh at yahoo.com (Ehsan Taghizadeh)
Date: Sun, 19 May 2024 08:10:30 +0000 (UTC)
Subject: [ncl-talk] ncl build
In-Reply-To: <CAJ0V9jrTKt6rjLmgwEt84ecg0iOD=bB59S2xpW8RnLr3dPa-QA@mail.gmail.com>
References: <CAJ0V9jrTKt6rjLmgwEt84ecg0iOD=bB59S2xpW8RnLr3dPa-QA@mail.gmail.com>
Message-ID: <1274241765.903815.1716106230048@mail.yahoo.com>

Dear Debasish, 
I suggest installing the precompiled binary designed for CentOS 7.6, which you can obtain from the provided link:https://www.earthsystemgrid.org/dataset/ncl.662.dap/file.html

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

 

    On Tuesday, March 19, 2024 at 03:11:08 AM PDT, Debasish Hazra via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:  
 
 Hi,Is there any precompiled latest version of NCL in Rocky 8 ? 

ThanksDebasish
_______________________________________________
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/d4f2b15c/attachment.htm>