From debasish.hazra5 at gmail.com  Wed Jun  1 17:38:20 2022
From: debasish.hazra5 at gmail.com (Debasish Hazra)
Date: Wed, 1 Jun 2022 19:38:20 -0400
Subject: [ncl-talk] seasonality remove and trend
In-Reply-To: <CAOF1d_7drT6A97WmN=TD2grDGEqw16Fch6pyYhVtU4fAxC5x+w@mail.gmail.com>
References: <CAJ0V9joyqTud8W1DzSxmekdeP3LOBv52Sn9h3zCsMUER05d6Kg@mail.gmail.com>
 <CAOF1d_7drT6A97WmN=TD2grDGEqw16Fch6pyYhVtU4fAxC5x+w@mail.gmail.com>
Message-ID: <CAJ0V9jrbxOGTcSG6eyX+797y7cGBCOhpU=SztXKzmdDnA-DtDw@mail.gmail.com>

Thanks Dennis.

Deba

On Mon, May 30, 2022 at 9:39 AM Dennis Shea <shea at ucar.edu> wrote:

> One approach:
>
> [1] compute the climatological monthly means and then calculate the
> monthly anomalies: x(time,lat,lon)
>
>   xMonClm = *clmMonTLL* <https://www.ncl.ucar.edu/Document/Functions/Contributed/clmMonTLL.shtml> (x)                  ; monthly climatology   (12,lat,lon)
>   xMonAnom = *calcMonAnomTLL* <https://www.ncl.ucar.edu/Document/Functions/Contributed/calcMonAnomTLL.shtml> (x,xMonClm)    ; (time,lat,lon);  (0,1,2)
>
>
> [2]* xDtrend = dtrend_n
> <https://www.ncl.ucar.edu/Document/Functions/Built-in/dtrend_n.shtml>*(xMonAnom,True,0)
> ; trend returned as attribute
>      printVarSummary(xDtrend)
>
> ======
> NCL is not a statistics package. It does have some simple regression
> functions.
>
> *regline*
> <https://www.ncl.ucar.edu/Document/Functions/Built-in/regline.shtml>,
> *regline_stats*
> <https://www.ncl.ucar.edu/Document/Functions/Contributed/regline_stats.shtml>,
> *regCoef*
> <https://www.ncl.ucar.edu/Document/Functions/Built-in/regCoef.shtml>,
> *regCoef_n*
> <https://www.ncl.ucar.edu/Document/Functions/Built-in/regCoef_n.shtml>,
> *reg_multlin_stats*
> <https://www.ncl.ucar.edu/Document/Functions/Contributed/reg_multlin_stats.shtml>
>
> *bootstrap_regcoef*
> <https://www.ncl.ucar.edu/Document/Functions/Bootstrap/bootstrap_regcoef.shtml>
>
>
> <https://www.ncl.ucar.edu/Document/Functions/Built-in/equiv_sample_size.shtml>
> <https://www.ncl.ucar.edu/Document/Functions/Built-in/rtest.shtml>
>
> On Sun, May 29, 2022 at 2:30 PM Debasish Hazra via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Hi,
>> I am analyzing 25 years monthly data (300 months) of a variable and would
>> like to measure the trend of the variable from beginning to the end of the
>> time. Questions are :
>> 1. How do I remove the seasonal variability by using a 12-month running
>> average prior to computing the linear trend. What function in NCL to be
>> used.
>> 2. Is there any weighted least square linear regression available in NCL
>> to be used for regression and trend analysis ?
>>
>> Appreciate your suggestions.
>> Thanks,
>> Deba
>> _______________________________________________
>> 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/20220601/b6cc74fb/attachment.htm>

From najibgal at yahoo.com  Thu Jun  2 04:35:28 2022
From: najibgal at yahoo.com (Najib Yusuf)
Date: Thu, 2 Jun 2022 10:35:28 +0000 (UTC)
Subject: [ncl-talk] Assigned plot is replaced by unassigned plot
References: <23082715.6212652.1654166128079.ref@mail.yahoo.com>
Message-ID: <23082715.6212652.1654166128079@mail.yahoo.com>

Dear NCL experts and colleagues,
I have a challenge that I need a guide about please, I assigned ncl to plot 4 plots (plot0,plot1,plot2,plot3) though I ended up having the 4 plots but with plot3 (the last) not displayed but replaced by another plot that I did not assign as the first in the /2,2/ panel. What may likely cause this please. I hope my question is clear.
Thank you?

------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/8011c72c/attachment.htm>

From ciarlo.james at gmail.com  Thu Jun  2 05:12:48 2022
From: ciarlo.james at gmail.com (James Ciarlo')
Date: Thu, 2 Jun 2022 13:12:48 +0200
Subject: [ncl-talk] Fwd:  plotting rotated grid - good overlay or bad box
In-Reply-To: <CAKhgdeTcOrc=FpeuAKEbGa72sk1TmJbMFHa4CU=jwTXmhV3vZA@mail.gmail.com>
References: <CAKhgdeTcOrc=FpeuAKEbGa72sk1TmJbMFHa4CU=jwTXmhV3vZA@mail.gmail.com>
Message-ID: <CAKhgdeSQf5=RPCMgvCz=xAzYCGgn-AmRN_+vLP6yzSFJfH5YCg@mail.gmail.com>

Dear NCL users,

I am trying to plot some data in rotated coordinates. I followed the
instructions on the website (
https://www.ncl.ucar.edu/Applications/rotatedltln.shtml) and I got to a
situation where I either have to switch off the terrain map, or cut out
part of the data.

Specifically, the example shows that res at tfDoNDCOverlay should be True, but
I get the result in the image "plot1rot_true.png", when I switch it to
False, the map aligns well, but you can see that part of the Mediterranean
on the East got cut out.

I am not sure what else to do. I have attached my script and a data file.

The file is found here:
https://www.dropbox.com/sh/8l0a3b5e3pfpq13/AACiICNlYCsAdvjP8X2xpAxea?dl=0

Thanks
James
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/31dbe2dd/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot1rot_true.png
Type: image/png
Size: 179875 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/31dbe2dd/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot1rot_false.png
Type: image/png
Size: 169701 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/31dbe2dd/attachment-0003.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot1rot.ncl
Type: application/octet-stream
Size: 2363 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/31dbe2dd/attachment-0001.obj>

From asphilli at ucar.edu  Thu Jun  2 10:55:45 2022
From: asphilli at ucar.edu (Adam Phillips)
Date: Thu, 2 Jun 2022 10:55:45 -0600
Subject: [ncl-talk] Assigned plot is replaced by unassigned plot
In-Reply-To: <23082715.6212652.1654166128079@mail.yahoo.com>
References: <23082715.6212652.1654166128079.ref@mail.yahoo.com>
 <23082715.6212652.1654166128079@mail.yahoo.com>
Message-ID: <CACvWuqBOcxp_yhu8xKU=bdKTRNWPN31j4TMbJmc1szhrs+gw8Q@mail.gmail.com>

Hi Najib,
Without seeing your script, I do not believe that anyone on ncl-talk can
assist you to solve the issue. Please send your script to ncl-talk.
Adam

On Thu, Jun 2, 2022 at 4:35 AM Najib Yusuf via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear NCL experts and colleagues,
>
> I have a challenge that I need a guide about please, I assigned ncl to
> plot 4 plots (plot0,plot1,plot2,plot3) though I ended up having the 4 plots
> but with plot3 (the last) not displayed but replaced by another plot that I
> did not assign as the first in the /2,2/ panel. What may likely cause this
> please. I hope my question is clear.
>
> Thank you
>
> *------------------------------------*
>
> *Najib Yusuf Galadanci Ph.D.*
>
> *Assistant Director (R & D)*
>
> *Lower Atmospheric Dynamics(Aerosol radiative forcing, air quality
> observations and modeling)Centre for Atmospheric ResearchNational Space
> Research and Development Agency. Kogi StateUniversity Campus, Anyigba. Kogi
> State. Nigeria+234 80 3871 4158------------------------------------Know
> that Victory comes with Patience, Relief with Affliction, and Ease with
> Hardship.Prophet Muhammad Rasulullah *
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/997a301e/attachment.htm>

From asphilli at ucar.edu  Thu Jun  2 12:00:47 2022
From: asphilli at ucar.edu (Adam Phillips)
Date: Thu, 2 Jun 2022 12:00:47 -0600
Subject: [ncl-talk] Fwd: plotting rotated grid - good overlay or bad box
In-Reply-To: <CAKhgdeSQf5=RPCMgvCz=xAzYCGgn-AmRN_+vLP6yzSFJfH5YCg@mail.gmail.com>
References: <CAKhgdeTcOrc=FpeuAKEbGa72sk1TmJbMFHa4CU=jwTXmhV3vZA@mail.gmail.com>
 <CAKhgdeSQf5=RPCMgvCz=xAzYCGgn-AmRN_+vLP6yzSFJfH5YCg@mail.gmail.com>
Message-ID: <CACvWuqAByWEsQboNfAn_dAV3n=N=hnT0a56jTXSnFNTNsDNOzQ@mail.gmail.com>

Hi James,
It looks like you took your code from example #1
<https://www.ncl.ucar.edu/Applications/Scripts/rotatedltln_1.ncl> in the
rotated lat/lon examples. That example demonstrated 2 methods to plot
rotated data:
1 - Set tfDoNDCOverlay=True and do not transform the data onto the map
projection.
2 - Attach the latitude/longitude coordinates (lat2d,lon2d) to the array to
be plotted, do not use tfDoNDCOverlay thus allowing NCL to transform the
data to the map projection.

It looks like from your script that you're trying to use method #2, but
were setting tfDoNDCOverlay = True. I have modified your script, and
tweaked the map settings so that the entire data map is shown.
Hope that helps!
Adam

On Thu, Jun 2, 2022 at 5:13 AM James Ciarlo' via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

>
> Dear NCL users,
>
> I am trying to plot some data in rotated coordinates. I followed the
> instructions on the website (
> https://www.ncl.ucar.edu/Applications/rotatedltln.shtml) and I got to a
> situation where I either have to switch off the terrain map, or cut out
> part of the data.
>
> Specifically, the example shows that res at tfDoNDCOverlay should be True,
> but I get the result in the image "plot1rot_true.png", when I switch it to
> False, the map aligns well, but you can see that part of the Mediterranean
> on the East got cut out.
>
> I am not sure what else to do. I have attached my script and a data file.
>
> The file is found here:
> https://www.dropbox.com/sh/8l0a3b5e3pfpq13/AACiICNlYCsAdvjP8X2xpAxea?dl=0
>
> Thanks
> James
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/428d37d9/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot1rot.png
Type: image/png
Size: 176441 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/428d37d9/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot1rot.ncl
Type: application/octet-stream
Size: 2460 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/428d37d9/attachment.obj>

From najibgal at yahoo.com  Thu Jun  2 13:20:30 2022
From: najibgal at yahoo.com (Najib Yusuf)
Date: Thu, 2 Jun 2022 19:20:30 +0000 (UTC)
Subject: [ncl-talk] Assigned plot is replaced by unassigned plot
In-Reply-To: <CACvWuqBOcxp_yhu8xKU=bdKTRNWPN31j4TMbJmc1szhrs+gw8Q@mail.gmail.com>
References: <23082715.6212652.1654166128079.ref@mail.yahoo.com>
 <23082715.6212652.1654166128079@mail.yahoo.com>
 <CACvWuqBOcxp_yhu8xKU=bdKTRNWPN31j4TMbJmc1szhrs+gw8Q@mail.gmail.com>
Message-ID: <380316284.6458657.1654197630537@mail.yahoo.com>

Hi Adam,
Thank you. Kindly find attached the script.
Najib

------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

 

    On Thursday, June 2, 2022, 05:56:12 PM GMT+1, Adam Phillips <asphilli at ucar.edu> wrote:  
 
 Hi Najib,Without seeing your script, I do not believe that anyone on ncl-talk can assist you to solve the?issue. Please send your script to ncl-talk.Adam
On Thu, Jun 2, 2022 at 4:35 AM Najib Yusuf via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:

Dear NCL experts and colleagues,
I have a challenge that I need a guide about please, I assigned ncl to plot 4 plots (plot0,plot1,plot2,plot3) though I ended up having the 4 plots but with plot3 (the last) not displayed but replaced by another plot that I did not assign as the first in the /2,2/ panel. What may likely cause this please. I hope my question is clear.
Thank you?

------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

_______________________________________________
ncl-talk mailing list
ncl-talk at mailman.ucar.edu
List instructions, subscriber options, unsubscribe:
https://mailman.ucar.edu/mailman/listinfo/ncl-talk



-- 
Adam Phillips 
Associate Scientist,? Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/?? 303-497-1726 
  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/8d52ee20/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: najiby_script.ncl
Type: application/octet-stream
Size: 8107 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/8d52ee20/attachment.obj>

From kunal.bali9 at gmail.com  Thu Jun  2 14:09:51 2022
From: kunal.bali9 at gmail.com (Kunal Bali)
Date: Thu, 2 Jun 2022 12:09:51 -0800
Subject: [ncl-talk] add_topography
Message-ID: <CAOiHCtTP55L+q+vxjJVuyZxb3SQko1ozhEJ6scb5MVB7Crichg@mail.gmail.com>

Hi NCL users,

Could you please give me an idea of how I can add the topography to the
attached plot?


Thanks

---
Kunal Bali
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/2b5393d1/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.pdf
Type: application/pdf
Size: 69482 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/2b5393d1/attachment-0001.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.ncl
Type: application/octet-stream
Size: 7837 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220602/2b5393d1/attachment-0001.obj>

From brownrig at ucar.edu  Thu Jun  2 19:21:42 2022
From: brownrig at ucar.edu (Rick Brownrigg)
Date: Thu, 2 Jun 2022 19:21:42 -0600
Subject: [ncl-talk] Assigned plot is replaced by unassigned plot
In-Reply-To: <380316284.6458657.1654197630537@mail.yahoo.com>
References: <23082715.6212652.1654166128079.ref@mail.yahoo.com>
 <23082715.6212652.1654166128079@mail.yahoo.com>
 <CACvWuqBOcxp_yhu8xKU=bdKTRNWPN31j4TMbJmc1szhrs+gw8Q@mail.gmail.com>
 <380316284.6458657.1654197630537@mail.yahoo.com>
Message-ID: <CAGKRhbGJx0Tx=iwVziHqmSx=wo3RrSBftAOJ83snQ5pwuX+7MQ@mail.gmail.com>

Hi Najib,

In your calls to gsn_csm_contour_map_ce(), you are plotting the same
variable, O3diff1, for each of plot1, plot2, plot3. Did you mean for plots
2 and 3 to depict variables O3diff2 & O3diff3?

Also, I'm unsure of your intent with respect to the resources -- it *looks*
like you make the calls to gsn_csm_contour_map_ce, and then set the
resources (in particular, note the "dangling" resources after the creation
of plot3).  Of course any resource settings should be done before the
creation of the plot they are intended to affect.

Hope that helps...
Rick


On Thu, Jun 2, 2022 at 1:20 PM Najib Yusuf via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi Adam,
>
> Thank you. Kindly find attached the script.
>
> Najib
>
> *------------------------------------*
>
> *Najib Yusuf Galadanci Ph.D.*
>
> *Assistant Director (R & D)*
>
> *Lower Atmospheric Dynamics(Aerosol radiative forcing, air quality
> observations and modeling)Centre for Atmospheric ResearchNational Space
> Research and Development Agency. Kogi StateUniversity Campus, Anyigba. Kogi
> State. Nigeria+234 80 3871 4158------------------------------------Know
> that Victory comes with Patience, Relief with Affliction, and Ease with
> Hardship.Prophet Muhammad Rasulullah *
>
>
>
> On Thursday, June 2, 2022, 05:56:12 PM GMT+1, Adam Phillips <
> asphilli at ucar.edu> wrote:
>
>
> Hi Najib,
> Without seeing your script, I do not believe that anyone on ncl-talk can
> assist you to solve the issue. Please send your script to ncl-talk.
> Adam
>
> On Thu, Jun 2, 2022 at 4:35 AM Najib Yusuf via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
> Dear NCL experts and colleagues,
>
> I have a challenge that I need a guide about please, I assigned ncl to
> plot 4 plots (plot0,plot1,plot2,plot3) though I ended up having the 4 plots
> but with plot3 (the last) not displayed but replaced by another plot that I
> did not assign as the first in the /2,2/ panel. What may likely cause this
> please. I hope my question is clear.
>
> Thank you
>
> *------------------------------------*
>
> *Najib Yusuf Galadanci Ph.D.*
>
> *Assistant Director (R & D)*
>
> Lower Atmospheric Dynamics
>
> (Aerosol radiative forcing, air quality observations and modeling)
>
> Centre for Atmospheric Research
>
> National Space Research and Development Agency. Kogi State
>
> University Campus, Anyigba. Kogi State. Nigeria
>
> +234 80 3871 4158
> ------------------------------------
> Know that Victory comes with Patience,
> Relief with Affliction, and Ease with Hardship.
>
> Prophet Muhammad Rasulullah
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
> _______________________________________________
> 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/20220602/4638f211/attachment.htm>

From ciarlo.james at gmail.com  Fri Jun  3 01:43:53 2022
From: ciarlo.james at gmail.com (James Ciarlo')
Date: Fri, 3 Jun 2022 09:43:53 +0200
Subject: [ncl-talk] Fwd: plotting rotated grid - good overlay or bad box
In-Reply-To: <CACvWuqAByWEsQboNfAn_dAV3n=N=hnT0a56jTXSnFNTNsDNOzQ@mail.gmail.com>
References: <CAKhgdeTcOrc=FpeuAKEbGa72sk1TmJbMFHa4CU=jwTXmhV3vZA@mail.gmail.com>
 <CAKhgdeSQf5=RPCMgvCz=xAzYCGgn-AmRN_+vLP6yzSFJfH5YCg@mail.gmail.com>
 <CACvWuqAByWEsQboNfAn_dAV3n=N=hnT0a56jTXSnFNTNsDNOzQ@mail.gmail.com>
Message-ID: <CAKhgdeTFG_jaA5oBYn0QDArf3=wxs_hyoa54k623-+cNk4+-cA@mail.gmail.com>

Hi Adam,

Thanks for your reply. To be honest, I had tried that system, but I was
trying to avoid that blank space and the inclined box.

I was hoping to get something like the method #1, that focuses only on the
box with the data.

On Thu, 2 Jun 2022 at 20:01, Adam Phillips <asphilli at ucar.edu> wrote:

> Hi James,
> It looks like you took your code from example #1
> <https://www.ncl.ucar.edu/Applications/Scripts/rotatedltln_1.ncl> in the
> rotated lat/lon examples. That example demonstrated 2 methods to plot
> rotated data:
> 1 - Set tfDoNDCOverlay=True and do not transform the data onto the map
> projection.
> 2 - Attach the latitude/longitude coordinates (lat2d,lon2d) to the array
> to be plotted, do not use tfDoNDCOverlay thus allowing NCL to transform the
> data to the map projection.
>
> It looks like from your script that you're trying to use method #2, but
> were setting tfDoNDCOverlay = True. I have modified your script, and
> tweaked the map settings so that the entire data map is shown.
> Hope that helps!
> Adam
>
> On Thu, Jun 2, 2022 at 5:13 AM James Ciarlo' via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>>
>> Dear NCL users,
>>
>> I am trying to plot some data in rotated coordinates. I followed the
>> instructions on the website (
>> https://www.ncl.ucar.edu/Applications/rotatedltln.shtml) and I got to a
>> situation where I either have to switch off the terrain map, or cut out
>> part of the data.
>>
>> Specifically, the example shows that res at tfDoNDCOverlay should be True,
>> but I get the result in the image "plot1rot_true.png", when I switch it to
>> False, the map aligns well, but you can see that part of the Mediterranean
>> on the East got cut out.
>>
>> I am not sure what else to do. I have attached my script and a data file.
>>
>> The file is found here:
>> https://www.dropbox.com/sh/8l0a3b5e3pfpq13/AACiICNlYCsAdvjP8X2xpAxea?dl=0
>>
>> Thanks
>> James
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220603/310daca8/attachment.htm>

From najibgal at yahoo.com  Fri Jun  3 04:45:27 2022
From: najibgal at yahoo.com (Najib Yusuf)
Date: Fri, 3 Jun 2022 10:45:27 +0000 (UTC)
Subject: [ncl-talk] Assigned plot is replaced by unassigned plot
In-Reply-To: <CAGKRhbGJx0Tx=iwVziHqmSx=wo3RrSBftAOJ83snQ5pwuX+7MQ@mail.gmail.com>
References: <23082715.6212652.1654166128079.ref@mail.yahoo.com>
 <23082715.6212652.1654166128079@mail.yahoo.com>
 <CACvWuqBOcxp_yhu8xKU=bdKTRNWPN31j4TMbJmc1szhrs+gw8Q@mail.gmail.com>
 <380316284.6458657.1654197630537@mail.yahoo.com>
 <CAGKRhbGJx0Tx=iwVziHqmSx=wo3RrSBftAOJ83snQ5pwuX+7MQ@mail.gmail.com>
Message-ID: <555542295.6737661.1654253127546@mail.yahoo.com>

Hi Adam,
Oh yes, thank you for this important observation. I?meant to plot O3diff0 as plot0 and O3diff1 as plot1, O3diff2 and O3diff3 as plot2 and plot3 respectively.
I moved the resources;?
;--Panel all three plots in a m x n configurationpres = True
; add lat long info to plot? ? ? ? ? ? ? ? ; add lat long info to plot? ? txres? ? ? ? ? ? ? ?= True? ? txres at txFontHeightF = 0.017
? ?panres = True? ?panres at gsnFrame = False? ?panres at gsnMaximize = True? ? panres at gsnPanelTop = 0.96? ? ?panres at gsnPanelBottom = 0.05;? ? ?gsn_panel (wks, (/plot0,plot1,plot2,plot3/),(/2,2/),panres)
? ? frame(wks);-------------------------------------------------------------------------?delete(res)
above the?calls to?gsn_csm_contour_map_ce but somehow it didn't plot.
Please see the panel plots am getting after correcting the observation and before moving the above resources.



Thank you
------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

 

    On Friday, June 3, 2022, 02:21:54 AM GMT+1, Rick Brownrigg <brownrig at ucar.edu> wrote:  
 
 Hi Najib,
In your calls to gsn_csm_contour_map_ce(), you are plotting the same variable, O3diff1, for each of plot1, plot2, plot3. Did you mean for plots 2 and 3 to depict variables O3diff2 & O3diff3?
Also, I'm unsure of your intent with respect to the resources -- it *looks* like you make the calls to?gsn_csm_contour_map_ce, and then set the resources (in particular, note the "dangling" resources after the creation of plot3).? Of course any resource settings should be done before the creation of the plot they are intended to affect.
Hope that helps...Rick

On Thu, Jun 2, 2022 at 1:20 PM Najib Yusuf via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:

Hi Adam,
Thank you. Kindly find attached the script.
Najib

------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

 

    On Thursday, June 2, 2022, 05:56:12 PM GMT+1, Adam Phillips <asphilli at ucar.edu> wrote:  
 
 Hi Najib,Without seeing your script, I do not believe that anyone on ncl-talk can assist you to solve the?issue. Please send your script to ncl-talk.Adam
On Thu, Jun 2, 2022 at 4:35 AM Najib Yusuf via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:

Dear NCL experts and colleagues,
I have a challenge that I need a guide about please, I assigned ncl to plot 4 plots (plot0,plot1,plot2,plot3) though I ended up having the 4 plots but with plot3 (the last) not displayed but replaced by another plot that I did not assign as the first in the /2,2/ panel. What may likely cause this please. I hope my question is clear.
Thank you?

------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

_______________________________________________
ncl-talk mailing list
ncl-talk at mailman.ucar.edu
List instructions, subscriber options, unsubscribe:
https://mailman.ucar.edu/mailman/listinfo/ncl-talk



-- 
Adam Phillips 
Associate Scientist,? Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/?? 303-497-1726 
  _______________________________________________
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/20220603/018107f7/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: najiby_plot.png
Type: image/png
Size: 121881 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220603/018107f7/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 1654253097209blob.jpg
Type: image/png
Size: 121881 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220603/018107f7/attachment-0001.png>

From brownrig at ucar.edu  Fri Jun  3 11:10:23 2022
From: brownrig at ucar.edu (Rick Brownrigg)
Date: Fri, 3 Jun 2022 11:10:23 -0600
Subject: [ncl-talk] Assigned plot is replaced by unassigned plot
In-Reply-To: <555542295.6737661.1654253127546@mail.yahoo.com>
References: <23082715.6212652.1654166128079.ref@mail.yahoo.com>
 <23082715.6212652.1654166128079@mail.yahoo.com>
 <CACvWuqBOcxp_yhu8xKU=bdKTRNWPN31j4TMbJmc1szhrs+gw8Q@mail.gmail.com>
 <380316284.6458657.1654197630537@mail.yahoo.com>
 <CAGKRhbGJx0Tx=iwVziHqmSx=wo3RrSBftAOJ83snQ5pwuX+7MQ@mail.gmail.com>
 <555542295.6737661.1654253127546@mail.yahoo.com>
Message-ID: <CAGKRhbFB7k0w0rpbygjwgobmP482w-ZsKDnHYpsVn2OWqL-F4Q@mail.gmail.com>

I was not referring to the paneling resources. In your original script, you
had this:

 plot0 = gsn_csm_contour_map_ce(wks,O3diff0,res)
         res at gsnRightString = "O3 ppb"
         res at gsnLeftString= "DJF_dry cold"
         res at gsnCenterString = "Full minus Zero Nigeria"
         res at txFontHeightF   = 0.011
 plot1 = gsn_csm_contour_map_ce(wks,O3diff1,res)
          res at gsnRightString = "O3 ppb"
          res at gsnLeftString= "MAM_dry warm"
          res at gsnCenterString = "Full minus Zero Nigeria"
          res at txFontHeightF   = 0.011
 plot2 = gsn_csm_contour_map_ce(wks,O3diff1,res)
          res at gsnRightString = "O3 ppb"
          res at gsnLeftString= "JJA_monsoon"
          res at gsnCenterString = "Full minus Zero Nigeria"
          res at txFontHeightF   = 0.011
 plot3 = gsn_csm_contour_map_ce(wks,O3diff1,res)
          res at gsnRightString = "O3 ppb"
          res at gsnLeftString= "SON_post_monsoon"
          res at gsnCenterString = "Full minus Zero Nigeria"
          res at txFontHeightF   = 0.011

The resources get applied to the *next* call to gsn_csm_contour_map_ce().
Notice that the resources in blue are never used -- is this what you
intended?  More to the point, these 4 groups of resource settings differ by
only the gsnLeftString -- are the four plots properly labelled with the
original script?

Rick

On Fri, Jun 3, 2022 at 4:45 AM Najib Yusuf <najibgal at yahoo.com> wrote:

> Hi Adam,
>
> Oh yes, thank you for this important observation. I meant to plot O3diff0
> as plot0 and O3diff1 as plot1, O3diff2 and O3diff3 as plot2 and plot3
> respectively.
>
> I moved the resources;
>
> ;--Panel all three plots in a m x n configuration
> pres = True
>
> ; add lat long info to plot
>                 ; add lat long info to plot
>     txres               = True
>     txres at txFontHeightF = 0.017
>
>    panres = True
>    panres at gsnFrame = False
>    panres at gsnMaximize = True
>     panres at gsnPanelTop = 0.96
>      panres at gsnPanelBottom = 0.05
> ;     gsn_panel (wks, (/plot0,plot1,plot2,plot3/),(/2,2/),panres)
>
>     frame(wks)
> ;-------------------------------------------------------------------------
>  delete(res)
>
> above the calls to gsn_csm_contour_map_ce but somehow it didn't plot.
>
> Please see the panel plots am getting after correcting the observation and
> before moving the above resources.
>
> [image: Inline image]
>
>
> Thank you
>
> *------------------------------------*
>
> *Najib Yusuf Galadanci Ph.D.*
>
> *Assistant Director (R & D)*
>
> *Lower Atmospheric Dynamics(Aerosol radiative forcing, air quality
> observations and modeling)Centre for Atmospheric ResearchNational Space
> Research and Development Agency. Kogi StateUniversity Campus, Anyigba. Kogi
> State. Nigeria+234 80 3871 4158------------------------------------Know
> that Victory comes with Patience, Relief with Affliction, and Ease with
> Hardship.Prophet Muhammad Rasulullah *
>
>
>
> On Friday, June 3, 2022, 02:21:54 AM GMT+1, Rick Brownrigg <
> brownrig at ucar.edu> wrote:
>
>
> Hi Najib,
>
> In your calls to gsn_csm_contour_map_ce(), you are plotting the same
> variable, O3diff1, for each of plot1, plot2, plot3. Did you mean for plots
> 2 and 3 to depict variables O3diff2 & O3diff3?
>
> Also, I'm unsure of your intent with respect to the resources -- it
> *looks* like you make the calls to gsn_csm_contour_map_ce, and then set the
> resources (in particular, note the "dangling" resources after the creation
> of plot3).  Of course any resource settings should be done before the
> creation of the plot they are intended to affect.
>
> Hope that helps...
> Rick
>
>
> On Thu, Jun 2, 2022 at 1:20 PM Najib Yusuf via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
> Hi Adam,
>
> Thank you. Kindly find attached the script.
>
> Najib
>
> *------------------------------------*
>
> *Najib Yusuf Galadanci Ph.D.*
>
> *Assistant Director (R & D)*
>
> Lower Atmospheric Dynamics
>
> (Aerosol radiative forcing, air quality observations and modeling)
>
> Centre for Atmospheric Research
>
> National Space Research and Development Agency. Kogi State
>
> University Campus, Anyigba. Kogi State. Nigeria
>
> +234 80 3871 4158
> ------------------------------------
> Know that Victory comes with Patience,
> Relief with Affliction, and Ease with Hardship.
>
> Prophet Muhammad Rasulullah
>
>
>
>
>
> On Thursday, June 2, 2022, 05:56:12 PM GMT+1, Adam Phillips <
> asphilli at ucar.edu> wrote:
>
>
> Hi Najib,
> Without seeing your script, I do not believe that anyone on ncl-talk can
> assist you to solve the issue. Please send your script to ncl-talk.
> Adam
>
> On Thu, Jun 2, 2022 at 4:35 AM Najib Yusuf via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
> Dear NCL experts and colleagues,
>
> I have a challenge that I need a guide about please, I assigned ncl to
> plot 4 plots (plot0,plot1,plot2,plot3) though I ended up having the 4 plots
> but with plot3 (the last) not displayed but replaced by another plot that I
> did not assign as the first in the /2,2/ panel. What may likely cause this
> please. I hope my question is clear.
>
> Thank you
>
> *------------------------------------*
>
> *Najib Yusuf Galadanci Ph.D.*
>
> *Assistant Director (R & D)*
>
> Lower Atmospheric Dynamics
>
> (Aerosol radiative forcing, air quality observations and modeling)
>
> Centre for Atmospheric Research
>
> National Space Research and Development Agency. Kogi State
>
> University Campus, Anyigba. Kogi State. Nigeria
>
> +234 80 3871 4158
> ------------------------------------
> Know that Victory comes with Patience,
> Relief with Affliction, and Ease with Hardship.
>
> Prophet Muhammad Rasulullah
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
> _______________________________________________
> 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/20220603/443cbcb9/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 1654253097209blob.jpg
Type: image/jpeg
Size: 121881 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220603/443cbcb9/attachment-0001.jpg>

From najibgal at yahoo.com  Fri Jun  3 11:24:17 2022
From: najibgal at yahoo.com (Najib Yusuf)
Date: Fri, 3 Jun 2022 17:24:17 +0000 (UTC)
Subject: [ncl-talk] Assigned plot is replaced by unassigned plot
In-Reply-To: <CAGKRhbFB7k0w0rpbygjwgobmP482w-ZsKDnHYpsVn2OWqL-F4Q@mail.gmail.com>
References: <23082715.6212652.1654166128079.ref@mail.yahoo.com>
 <23082715.6212652.1654166128079@mail.yahoo.com>
 <CACvWuqBOcxp_yhu8xKU=bdKTRNWPN31j4TMbJmc1szhrs+gw8Q@mail.gmail.com>
 <380316284.6458657.1654197630537@mail.yahoo.com>
 <CAGKRhbGJx0Tx=iwVziHqmSx=wo3RrSBftAOJ83snQ5pwuX+7MQ@mail.gmail.com>
 <555542295.6737661.1654253127546@mail.yahoo.com>
 <CAGKRhbFB7k0w0rpbygjwgobmP482w-ZsKDnHYpsVn2OWqL-F4Q@mail.gmail.com>
Message-ID: <1383840092.432726.1654277057392@mail.yahoo.com>

Oh I see now, I will rearrange and revert.? Yes, three out of the four plots appears well properly as labelled but the last one is missing/replace with different labelled.?
Thank you
Najib

Sent from Yahoo Mail on Android 
 
  On Fri, Jun 3, 2022 at 6:11 PM, Rick Brownrigg<brownrig at ucar.edu> wrote:   I was not referring to the paneling resources. In your original script, you had this:
?plot0 = gsn_csm_contour_map_ce(wks,O3diff0,res)
? ? ? ? ?res at gsnRightString = "O3 ppb"
? ? ? ? ?res at gsnLeftString= "DJF_dry cold"
? ? ? ? ?res at gsnCenterString = "Full minus Zero Nigeria"
? ? ? ? ?res at txFontHeightF ? = 0.011
?plot1 = gsn_csm_contour_map_ce(wks,O3diff1,res)
? ? ? ? ? res at gsnRightString = "O3 ppb"
? ? ? ? ? res at gsnLeftString= "MAM_dry warm"
? ? ? ? ? res at gsnCenterString = "Full minus Zero Nigeria"
? ? ? ? ? res at txFontHeightF ? = 0.011
?plot2 = gsn_csm_contour_map_ce(wks,O3diff1,res)
? ? ? ? ? res at gsnRightString = "O3 ppb"
? ? ? ? ? res at gsnLeftString= "JJA_monsoon"
? ? ? ? ? res at gsnCenterString = "Full minus Zero Nigeria"
? ? ? ? ? res at txFontHeightF ? = 0.011
?plot3 = gsn_csm_contour_map_ce(wks,O3diff1,res)
? ? ? ? ? res at gsnRightString = "O3 ppb"
? ? ? ? ? res at gsnLeftString= "SON_post_monsoon"
? ? ? ? ? res at gsnCenterString = "Full minus Zero Nigeria"
? ? ? ? ? res at txFontHeightF ? = 0.011
The resources get applied to the *next* call to gsn_csm_contour_map_ce(). Notice that the resources in blue are never used -- is this what you intended?? More to the point, these 4 groups of resource settings differ by only the gsnLeftString -- are the four plots properly labelled with the original script?

Rick

On Fri, Jun 3, 2022 at 4:45 AM Najib Yusuf <najibgal at yahoo.com> wrote:

Hi Adam,
Oh yes, thank you for this important observation. I?meant to plot O3diff0 as plot0 and O3diff1 as plot1, O3diff2 and O3diff3 as plot2 and plot3 respectively.
I moved the resources;?
;--Panel all three plots in a m x n configurationpres = True
; add lat long info to plot? ? ? ? ? ? ? ? ; add lat long info to plot? ? txres? ? ? ? ? ? ? ?= True? ? txres at txFontHeightF = 0.017
? ?panres = True? ?panres at gsnFrame = False? ?panres at gsnMaximize = True? ? panres at gsnPanelTop = 0.96? ? ?panres at gsnPanelBottom = 0.05;? ? ?gsn_panel (wks, (/plot0,plot1,plot2,plot3/),(/2,2/),panres)
? ? frame(wks);-------------------------------------------------------------------------?delete(res)
above the?calls to?gsn_csm_contour_map_ce but somehow it didn't plot.
Please see the panel plots am getting after correcting the observation and before moving the above resources.



Thank you
------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

 

    On Friday, June 3, 2022, 02:21:54 AM GMT+1, Rick Brownrigg <brownrig at ucar.edu> wrote:  
 
 Hi Najib,
In your calls to gsn_csm_contour_map_ce(), you are plotting the same variable, O3diff1, for each of plot1, plot2, plot3. Did you mean for plots 2 and 3 to depict variables O3diff2 & O3diff3?
Also, I'm unsure of your intent with respect to the resources -- it *looks* like you make the calls to?gsn_csm_contour_map_ce, and then set the resources (in particular, note the "dangling" resources after the creation of plot3).? Of course any resource settings should be done before the creation of the plot they are intended to affect.
Hope that helps...Rick

On Thu, Jun 2, 2022 at 1:20 PM Najib Yusuf via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:

Hi Adam,
Thank you. Kindly find attached the script.
Najib

------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

 

    On Thursday, June 2, 2022, 05:56:12 PM GMT+1, Adam Phillips <asphilli at ucar.edu> wrote:  
 
 Hi Najib,Without seeing your script, I do not believe that anyone on ncl-talk can assist you to solve the?issue. Please send your script to ncl-talk.Adam
On Thu, Jun 2, 2022 at 4:35 AM Najib Yusuf via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:

Dear NCL experts and colleagues,
I have a challenge that I need a guide about please, I assigned ncl to plot 4 plots (plot0,plot1,plot2,plot3) though I ended up having the 4 plots but with plot3 (the last) not displayed but replaced by another plot that I did not assign as the first in the /2,2/ panel. What may likely cause this please. I hope my question is clear.
Thank you?

------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

_______________________________________________
ncl-talk mailing list
ncl-talk at mailman.ucar.edu
List instructions, subscriber options, unsubscribe:
https://mailman.ucar.edu/mailman/listinfo/ncl-talk



-- 
Adam Phillips 
Associate Scientist,? Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/?? 303-497-1726 
  _______________________________________________
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/20220603/ee63e595/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 1654253097209blob.jpg
Type: image/jpeg
Size: 121881 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220603/ee63e595/attachment.jpg>

From najibgal at yahoo.com  Fri Jun  3 12:21:42 2022
From: najibgal at yahoo.com (Najib Yusuf)
Date: Fri, 3 Jun 2022 18:21:42 +0000 (UTC)
Subject: [ncl-talk] Assigned plot is replaced by unassigned plot
In-Reply-To: <CAGKRhbFB7k0w0rpbygjwgobmP482w-ZsKDnHYpsVn2OWqL-F4Q@mail.gmail.com>
References: <23082715.6212652.1654166128079.ref@mail.yahoo.com>
 <23082715.6212652.1654166128079@mail.yahoo.com>
 <CACvWuqBOcxp_yhu8xKU=bdKTRNWPN31j4TMbJmc1szhrs+gw8Q@mail.gmail.com>
 <380316284.6458657.1654197630537@mail.yahoo.com>
 <CAGKRhbGJx0Tx=iwVziHqmSx=wo3RrSBftAOJ83snQ5pwuX+7MQ@mail.gmail.com>
 <555542295.6737661.1654253127546@mail.yahoo.com>
 <CAGKRhbFB7k0w0rpbygjwgobmP482w-ZsKDnHYpsVn2OWqL-F4Q@mail.gmail.com>
Message-ID: <1009409870.6929433.1654280502930@mail.yahoo.com>

Dear Rick and Adam,
Thank you so much for your help. It works.
I really appreciate.

------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

 

    On Friday, June 3, 2022, 06:11:25 PM GMT+1, Rick Brownrigg <brownrig at ucar.edu> wrote:  
 
 I was not referring to the paneling resources. In your original script, you had this:
?plot0 = gsn_csm_contour_map_ce(wks,O3diff0,res)
? ? ? ? ?res at gsnRightString = "O3 ppb"
? ? ? ? ?res at gsnLeftString= "DJF_dry cold"
? ? ? ? ?res at gsnCenterString = "Full minus Zero Nigeria"
? ? ? ? ?res at txFontHeightF ? = 0.011
?plot1 = gsn_csm_contour_map_ce(wks,O3diff1,res)
? ? ? ? ? res at gsnRightString = "O3 ppb"
? ? ? ? ? res at gsnLeftString= "MAM_dry warm"
? ? ? ? ? res at gsnCenterString = "Full minus Zero Nigeria"
? ? ? ? ? res at txFontHeightF ? = 0.011
?plot2 = gsn_csm_contour_map_ce(wks,O3diff1,res)
? ? ? ? ? res at gsnRightString = "O3 ppb"
? ? ? ? ? res at gsnLeftString= "JJA_monsoon"
? ? ? ? ? res at gsnCenterString = "Full minus Zero Nigeria"
? ? ? ? ? res at txFontHeightF ? = 0.011
?plot3 = gsn_csm_contour_map_ce(wks,O3diff1,res)
? ? ? ? ? res at gsnRightString = "O3 ppb"
? ? ? ? ? res at gsnLeftString= "SON_post_monsoon"
? ? ? ? ? res at gsnCenterString = "Full minus Zero Nigeria"
? ? ? ? ? res at txFontHeightF ? = 0.011
The resources get applied to the *next* call to gsn_csm_contour_map_ce(). Notice that the resources in blue are never used -- is this what you intended?? More to the point, these 4 groups of resource settings differ by only the gsnLeftString -- are the four plots properly labelled with the original script?

Rick

On Fri, Jun 3, 2022 at 4:45 AM Najib Yusuf <najibgal at yahoo.com> wrote:

Hi Adam,
Oh yes, thank you for this important observation. I?meant to plot O3diff0 as plot0 and O3diff1 as plot1, O3diff2 and O3diff3 as plot2 and plot3 respectively.
I moved the resources;?
;--Panel all three plots in a m x n configurationpres = True
; add lat long info to plot? ? ? ? ? ? ? ? ; add lat long info to plot? ? txres? ? ? ? ? ? ? ?= True? ? txres at txFontHeightF = 0.017
? ?panres = True? ?panres at gsnFrame = False? ?panres at gsnMaximize = True? ? panres at gsnPanelTop = 0.96? ? ?panres at gsnPanelBottom = 0.05;? ? ?gsn_panel (wks, (/plot0,plot1,plot2,plot3/),(/2,2/),panres)
? ? frame(wks);-------------------------------------------------------------------------?delete(res)
above the?calls to?gsn_csm_contour_map_ce but somehow it didn't plot.
Please see the panel plots am getting after correcting the observation and before moving the above resources.



Thank you
------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

 

    On Friday, June 3, 2022, 02:21:54 AM GMT+1, Rick Brownrigg <brownrig at ucar.edu> wrote:  
 
 Hi Najib,
In your calls to gsn_csm_contour_map_ce(), you are plotting the same variable, O3diff1, for each of plot1, plot2, plot3. Did you mean for plots 2 and 3 to depict variables O3diff2 & O3diff3?
Also, I'm unsure of your intent with respect to the resources -- it *looks* like you make the calls to?gsn_csm_contour_map_ce, and then set the resources (in particular, note the "dangling" resources after the creation of plot3).? Of course any resource settings should be done before the creation of the plot they are intended to affect.
Hope that helps...Rick

On Thu, Jun 2, 2022 at 1:20 PM Najib Yusuf via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:

Hi Adam,
Thank you. Kindly find attached the script.
Najib

------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

 

    On Thursday, June 2, 2022, 05:56:12 PM GMT+1, Adam Phillips <asphilli at ucar.edu> wrote:  
 
 Hi Najib,Without seeing your script, I do not believe that anyone on ncl-talk can assist you to solve the?issue. Please send your script to ncl-talk.Adam
On Thu, Jun 2, 2022 at 4:35 AM Najib Yusuf via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:

Dear NCL experts and colleagues,
I have a challenge that I need a guide about please, I assigned ncl to plot 4 plots (plot0,plot1,plot2,plot3) though I ended up having the 4 plots but with plot3 (the last) not displayed but replaced by another plot that I did not assign as the first in the /2,2/ panel. What may likely cause this please. I hope my question is clear.
Thank you?

------------------------------------

Najib Yusuf Galadanci Ph.D.

Assistant Director (R & D)

Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.
Prophet Muhammad Rasulullah

?

_______________________________________________
ncl-talk mailing list
ncl-talk at mailman.ucar.edu
List instructions, subscriber options, unsubscribe:
https://mailman.ucar.edu/mailman/listinfo/ncl-talk



-- 
Adam Phillips 
Associate Scientist,? Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/?? 303-497-1726 
  _______________________________________________
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/20220603/cd9ea1e7/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 1654253097209blob.jpg
Type: image/jpeg
Size: 121881 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220603/cd9ea1e7/attachment-0001.jpg>

From dmjalal90 at gmail.com  Mon Jun  6 05:00:10 2022
From: dmjalal90 at gmail.com (Md. Jalal Uddin)
Date: Mon, 6 Jun 2022 19:00:10 +0800
Subject: [ncl-talk] How to access example data for computing the
 Penman-Monteith reference evapotranspiration
Message-ID: <CAN5hBbsiZ5tK-iTs7d=8Zdnqx9A=EtS4BhDU=S5oCjAbXhU9gw@mail.gmail.com>

Dear all,

I would like to compute the Penman-Monteith reference evapotranspiration
via FAO-56 methodology using
https://www.ncl.ucar.edu/Applications/Scripts/fao56_1.ncl example. But I
could not access the data used in this example.

Could anyone share the source of data to compute evapotranspiration?

Regards,
Jalal
-- 
*Md. Jalal Uddin*
Founder and Director of Research Society
PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
M.Sc. in Applied Meteorology (NUIST)
B.Sc. in Disaster Management (PSTU)
Email: founder-and-director at researchsociety20.org, 20205103002 at nuist.edu.cn
<dmjalal90 at nuist.edu.cn>
Website: *https://researchsociety20.org/founder-and-director/
<https://researchsociety20.org/founder-and-director/>*
ResearchGate profile: https://www.researchgate.net/profile/Md_Uddin125
YouTube channel:https://www.youtube.com/jalalsdream
Cell: +8613260859092, +8801792052662
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220606/95f18d22/attachment.htm>

From shea at ucar.edu  Mon Jun  6 07:47:42 2022
From: shea at ucar.edu (Dennis Shea)
Date: Mon, 6 Jun 2022 07:47:42 -0600
Subject: [ncl-talk] How to access example data for computing the
 Penman-Monteith reference evapotranspiration
In-Reply-To: <CAN5hBbsiZ5tK-iTs7d=8Zdnqx9A=EtS4BhDU=S5oCjAbXhU9gw@mail.gmail.com>
References: <CAN5hBbsiZ5tK-iTs7d=8Zdnqx9A=EtS4BhDU=S5oCjAbXhU9gw@mail.gmail.com>
Message-ID: <CAOF1d_79VKBi17+QpZnn7SzRiSwoB-fu2ngu_KZgd6ZYsLUtdw@mail.gmail.com>

Unfortunately, the BCL* fieles are no longer available on systems
accessible to me.

Sorry.

On Mon, Jun 6, 2022 at 5:00 AM Md. Jalal Uddin via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear all,
>
> I would like to compute the Penman-Monteith reference evapotranspiration
> via FAO-56 methodology using
> https://www.ncl.ucar.edu/Applications/Scripts/fao56_1.ncl example. But I
> could not access the data used in this example.
>
> Could anyone share the source of data to compute evapotranspiration?
>
> Regards,
> Jalal
> --
> *Md. Jalal Uddin*
> Founder and Director of Research Society
> PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
> M.Sc. in Applied Meteorology (NUIST)
> B.Sc. in Disaster Management (PSTU)
> Email: founder-and-director at researchsociety20.org,
> 20205103002 at nuist.edu.cn <dmjalal90 at nuist.edu.cn>
> Website: *https://researchsociety20.org/founder-and-director/
> <https://researchsociety20.org/founder-and-director/>*
> ResearchGate profile: https://www.researchgate.net/profile/Md_Uddin125
> YouTube channel:https://www.youtube.com/jalalsdream
> Cell: +8613260859092, +8801792052662
> _______________________________________________
> 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/20220606/800fdb2b/attachment.htm>

From asphilli at ucar.edu  Mon Jun  6 14:23:42 2022
From: asphilli at ucar.edu (Adam Phillips)
Date: Mon, 6 Jun 2022 14:23:42 -0600
Subject: [ncl-talk] How to access example data for computing the
 Penman-Monteith reference evapotranspiration
In-Reply-To: <CAOF1d_79VKBi17+QpZnn7SzRiSwoB-fu2ngu_KZgd6ZYsLUtdw@mail.gmail.com>
References: <CAN5hBbsiZ5tK-iTs7d=8Zdnqx9A=EtS4BhDU=S5oCjAbXhU9gw@mail.gmail.com>
 <CAOF1d_79VKBi17+QpZnn7SzRiSwoB-fu2ngu_KZgd6ZYsLUtdw@mail.gmail.com>
Message-ID: <CACvWuqCF+S_7vfd77MSkgTKagm7N5cQfybKuRw=Tmves=Q=hRg@mail.gmail.com>

Hi Jalal,
I was able to retrieve the data from one of Dennis's directories. I will
email you offline with instructions on how to access it.
Adam

On Mon, Jun 6, 2022 at 7:48 AM Dennis Shea via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Unfortunately, the BCL* fieles are no longer available on systems
> accessible to me.
>
> Sorry.
>
> On Mon, Jun 6, 2022 at 5:00 AM Md. Jalal Uddin via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Dear all,
>>
>> I would like to compute the Penman-Monteith reference evapotranspiration
>> via FAO-56 methodology using
>> https://www.ncl.ucar.edu/Applications/Scripts/fao56_1.ncl example. But I
>> could not access the data used in this example.
>>
>> Could anyone share the source of data to compute evapotranspiration?
>>
>> Regards,
>> Jalal
>> --
>> *Md. Jalal Uddin*
>> Founder and Director of Research Society
>> PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
>> M.Sc. in Applied Meteorology (NUIST)
>> B.Sc. in Disaster Management (PSTU)
>> Email: founder-and-director at researchsociety20.org,
>> 20205103002 at nuist.edu.cn <dmjalal90 at nuist.edu.cn>
>> Website: *https://researchsociety20.org/founder-and-director/
>> <https://researchsociety20.org/founder-and-director/>*
>> ResearchGate profile: https://www.researchgate.net/profile/Md_Uddin125
>> YouTube channel:https://www.youtube.com/jalalsdream
>> Cell: +8613260859092, +8801792052662
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220606/4efe1ac5/attachment.htm>

From asphilli at ucar.edu  Mon Jun  6 17:02:25 2022
From: asphilli at ucar.edu (Adam Phillips)
Date: Mon, 6 Jun 2022 17:02:25 -0600
Subject: [ncl-talk] Fwd: plotting rotated grid - good overlay or bad box
In-Reply-To: <CAKhgdeTFG_jaA5oBYn0QDArf3=wxs_hyoa54k623-+cNk4+-cA@mail.gmail.com>
References: <CAKhgdeTcOrc=FpeuAKEbGa72sk1TmJbMFHa4CU=jwTXmhV3vZA@mail.gmail.com>
 <CAKhgdeSQf5=RPCMgvCz=xAzYCGgn-AmRN_+vLP6yzSFJfH5YCg@mail.gmail.com>
 <CACvWuqAByWEsQboNfAn_dAV3n=N=hnT0a56jTXSnFNTNsDNOzQ@mail.gmail.com>
 <CAKhgdeTFG_jaA5oBYn0QDArf3=wxs_hyoa54k623-+cNk4+-cA@mail.gmail.com>
Message-ID: <CACvWuqDhKSHwk6TSt4QLtQmH2+y=LpmCA=tFnkhSNmPAiA_P2Q@mail.gmail.com>

Hi James,
I have altered your script to go with method 1 from the 1st rotated lat/lon
example. However, I had to alter the mpCenterLatF/mpCenterLonF settings a
bit from the projection origins. I do not know why the alteration was
required, as usually one can set the projection origins and all is well.
Thus, you might want to alter the settings a bit on your own (lines 41-42)
to see if you can get a better fit of the data to the map.
Good luck,
Adam

On Fri, Jun 3, 2022 at 1:44 AM James Ciarlo' <ciarlo.james at gmail.com> wrote:

> Hi Adam,
>
> Thanks for your reply. To be honest, I had tried that system, but I was
> trying to avoid that blank space and the inclined box.
>
> I was hoping to get something like the method #1, that focuses only on the
> box with the data.
>
> On Thu, 2 Jun 2022 at 20:01, Adam Phillips <asphilli at ucar.edu> wrote:
>
>> Hi James,
>> It looks like you took your code from example #1
>> <https://www.ncl.ucar.edu/Applications/Scripts/rotatedltln_1.ncl> in the
>> rotated lat/lon examples. That example demonstrated 2 methods to plot
>> rotated data:
>> 1 - Set tfDoNDCOverlay=True and do not transform the data onto the map
>> projection.
>> 2 - Attach the latitude/longitude coordinates (lat2d,lon2d) to the array
>> to be plotted, do not use tfDoNDCOverlay thus allowing NCL to transform the
>> data to the map projection.
>>
>> It looks like from your script that you're trying to use method #2, but
>> were setting tfDoNDCOverlay = True. I have modified your script, and
>> tweaked the map settings so that the entire data map is shown.
>> Hope that helps!
>> Adam
>>
>> On Thu, Jun 2, 2022 at 5:13 AM James Ciarlo' via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>>
>>> Dear NCL users,
>>>
>>> I am trying to plot some data in rotated coordinates. I followed the
>>> instructions on the website (
>>> https://www.ncl.ucar.edu/Applications/rotatedltln.shtml) and I got to a
>>> situation where I either have to switch off the terrain map, or cut out
>>> part of the data.
>>>
>>> Specifically, the example shows that res at tfDoNDCOverlay should be True,
>>> but I get the result in the image "plot1rot_true.png", when I switch it to
>>> False, the map aligns well, but you can see that part of the Mediterranean
>>> on the East got cut out.
>>>
>>> I am not sure what else to do. I have attached my script and a data file.
>>>
>>> The file is found here:
>>> https://www.dropbox.com/sh/8l0a3b5e3pfpq13/AACiICNlYCsAdvjP8X2xpAxea?dl=0
>>>
>>> Thanks
>>> James
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at mailman.ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>
>>
>> --
>> Adam Phillips
>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
>> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>>
>> <http://www.cgd.ucar.edu/staff/asphilli>
>>
>

-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220606/587a8ce6/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot1rot.ncl
Type: application/octet-stream
Size: 2591 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220606/587a8ce6/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot1rot.png
Type: image/png
Size: 197613 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220606/587a8ce6/attachment.png>

From ciarlo.james at gmail.com  Tue Jun  7 06:57:39 2022
From: ciarlo.james at gmail.com (James Ciarlo')
Date: Tue, 7 Jun 2022 14:57:39 +0200
Subject: [ncl-talk] Fwd: plotting rotated grid - good overlay or bad box
In-Reply-To: <CACvWuqDhKSHwk6TSt4QLtQmH2+y=LpmCA=tFnkhSNmPAiA_P2Q@mail.gmail.com>
References: <CAKhgdeTcOrc=FpeuAKEbGa72sk1TmJbMFHa4CU=jwTXmhV3vZA@mail.gmail.com>
 <CAKhgdeSQf5=RPCMgvCz=xAzYCGgn-AmRN_+vLP6yzSFJfH5YCg@mail.gmail.com>
 <CACvWuqAByWEsQboNfAn_dAV3n=N=hnT0a56jTXSnFNTNsDNOzQ@mail.gmail.com>
 <CAKhgdeTFG_jaA5oBYn0QDArf3=wxs_hyoa54k623-+cNk4+-cA@mail.gmail.com>
 <CACvWuqDhKSHwk6TSt4QLtQmH2+y=LpmCA=tFnkhSNmPAiA_P2Q@mail.gmail.com>
Message-ID: <CAKhgdeQJJR6oDi6SapovTG_6uc84wJai_1cp-rj0zoC6MaFhKQ@mail.gmail.com>

Hi Adam,

Yes that fixed the problem thank you very much. It didn't cross my mind to
move the centre point.

It's very unfortunate that the defined centre points in the data need to be
corrected manually for the best fit but I'm really happy to have neat plots
again.

Thanks once again.

Best regards
James

On Tue, 7 Jun 2022 at 01:02, Adam Phillips <asphilli at ucar.edu> wrote:

> Hi James,
> I have altered your script to go with method 1 from the 1st rotated
> lat/lon example. However, I had to alter the mpCenterLatF/mpCenterLonF
> settings a bit from the projection origins. I do not know why the
> alteration was required, as usually one can set the projection origins and
> all is well. Thus, you might want to alter the settings a bit on your own
> (lines 41-42) to see if you can get a better fit of the data to the map.
> Good luck,
> Adam
>
> On Fri, Jun 3, 2022 at 1:44 AM James Ciarlo' <ciarlo.james at gmail.com>
> wrote:
>
>> Hi Adam,
>>
>> Thanks for your reply. To be honest, I had tried that system, but I was
>> trying to avoid that blank space and the inclined box.
>>
>> I was hoping to get something like the method #1, that focuses only on
>> the box with the data.
>>
>> On Thu, 2 Jun 2022 at 20:01, Adam Phillips <asphilli at ucar.edu> wrote:
>>
>>> Hi James,
>>> It looks like you took your code from example #1
>>> <https://www.ncl.ucar.edu/Applications/Scripts/rotatedltln_1.ncl> in
>>> the rotated lat/lon examples. That example demonstrated 2 methods to plot
>>> rotated data:
>>> 1 - Set tfDoNDCOverlay=True and do not transform the data onto the map
>>> projection.
>>> 2 - Attach the latitude/longitude coordinates (lat2d,lon2d) to the array
>>> to be plotted, do not use tfDoNDCOverlay thus allowing NCL to transform the
>>> data to the map projection.
>>>
>>> It looks like from your script that you're trying to use method #2, but
>>> were setting tfDoNDCOverlay = True. I have modified your script, and
>>> tweaked the map settings so that the entire data map is shown.
>>> Hope that helps!
>>> Adam
>>>
>>> On Thu, Jun 2, 2022 at 5:13 AM James Ciarlo' via ncl-talk <
>>> ncl-talk at mailman.ucar.edu> wrote:
>>>
>>>>
>>>> Dear NCL users,
>>>>
>>>> I am trying to plot some data in rotated coordinates. I followed the
>>>> instructions on the website (
>>>> https://www.ncl.ucar.edu/Applications/rotatedltln.shtml) and I got to
>>>> a situation where I either have to switch off the terrain map, or cut out
>>>> part of the data.
>>>>
>>>> Specifically, the example shows that res at tfDoNDCOverlay should be
>>>> True, but I get the result in the image "plot1rot_true.png", when I switch
>>>> it to False, the map aligns well, but you can see that part of the
>>>> Mediterranean on the East got cut out.
>>>>
>>>> I am not sure what else to do. I have attached my script and a data
>>>> file.
>>>>
>>>> The file is found here:
>>>> https://www.dropbox.com/sh/8l0a3b5e3pfpq13/AACiICNlYCsAdvjP8X2xpAxea?dl=0
>>>>
>>>> Thanks
>>>> James
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at mailman.ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>
>>>
>>> --
>>> Adam Phillips
>>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
>>> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>>>
>>> <http://www.cgd.ucar.edu/staff/asphilli>
>>>
>>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220607/af482961/attachment.htm>

From dmjalal90 at gmail.com  Wed Jun  8 03:42:13 2022
From: dmjalal90 at gmail.com (Md. Jalal Uddin)
Date: Wed, 8 Jun 2022 17:42:13 +0800
Subject: [ncl-talk] Need a clarification for "getMidMonDay" function in NCL
Message-ID: <CAN5hBbthSq-dubxj7Bd3-CyodvuOg=EikJXY5-T_uArE82eVqg@mail.gmail.com>

Hi all,

The fao56_1.ncl (https://www.ncl.ucar.edu/Applications/Scripts/fao56_1.ncl)
script used "getMidMonDay" function to get the middle day of the month, but
it seems that this function is not exist.

How to get the middle day of the month to compute evapotranspiration? Is it
similar to the "day_of_year" function (
https://www.ncl.ucar.edu/Document/Functions/Built-in/day_of_year.shtml)?

Regards,
Jalal
-- 
*Md. Jalal Uddin*
Founder and Director of Research Society
PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
M.Sc. in Applied Meteorology (NUIST)
B.Sc. in Disaster Management (PSTU)
Email: founder-and-director at researchsociety20.org, 20205103002 at nuist.edu.cn
<dmjalal90 at nuist.edu.cn>
Website: *https://researchsociety20.org/founder-and-director/
<https://researchsociety20.org/founder-and-director/>*
ResearchGate profile: https://www.researchgate.net/profile/Md_Uddin125
YouTube channel:https://www.youtube.com/jalalsdream
Cell: +8613260859092, +8801792052662
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220608/6a775fb4/attachment.htm>

From shea at ucar.edu  Wed Jun  8 09:23:07 2022
From: shea at ucar.edu (Dennis Shea)
Date: Wed, 8 Jun 2022 09:23:07 -0600
Subject: [ncl-talk] Need a clarification for "getMidMonDay" function in
 NCL
In-Reply-To: <CAN5hBbthSq-dubxj7Bd3-CyodvuOg=EikJXY5-T_uArE82eVqg@mail.gmail.com>
References: <CAN5hBbthSq-dubxj7Bd3-CyodvuOg=EikJXY5-T_uArE82eVqg@mail.gmail.com>
Message-ID: <CAOF1d_4yEqkqLe=NKO7W8FREQL0HuZuQkgc3KSTzb4sj0cdxJg@mail.gmail.com>

Oops. My bad!  I did not document the function.
getMidMonDay is a utility function contained within the 'crop' library
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/crop.ncl"
====
The code is 'simple'

undef("getMidMonDay")
function getMidMonDay(time:numeric)
;
; utility function
; The input 'time' (seconds/hours/days since ...") has the middle of the
month
; This returns the (*approximate*) 'middle day' of the current month as an
integer
; "Approximate" because day 15 is used for all months ...
;      rather than 14 (Feb), 14.5 (Leap year Feb) or 15.5 (31 day months).
; This is more than adequate.

local date, yyyy, mm, dd, jday
begin
    date = cd_calendar(time, 0)
    yyyy = toint(date(:,0))
    mm   = toint(date(:,1))

    dd   = conform(mm, 15, -1)    ; approx mid month

    jday = day_of_year(yyyy, mm, dd)

    copy_VarCoords(time, jday)
    jday at long_name = "mid-month day of current year"
    return(jday)
end

On Wed, Jun 8, 2022 at 3:43 AM Md. Jalal Uddin via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi all,
>
> The fao56_1.ncl (https://www.ncl.ucar.edu/Applications/Scripts/fao56_1.ncl)
> script used "getMidMonDay" function to get the middle day of the month, but
> it seems that this function is not exist.
>
> How to get the middle day of the month to compute evapotranspiration? Is
> it similar to the "day_of_year" function (
> https://www.ncl.ucar.edu/Document/Functions/Built-in/day_of_year.shtml)?
>
> Regards,
> Jalal
> --
> *Md. Jalal Uddin*
> Founder and Director of Research Society
> PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
> M.Sc. in Applied Meteorology (NUIST)
> B.Sc. in Disaster Management (PSTU)
> Email: founder-and-director at researchsociety20.org,
> 20205103002 at nuist.edu.cn <dmjalal90 at nuist.edu.cn>
> Website: *https://researchsociety20.org/founder-and-director/
> <https://researchsociety20.org/founder-and-director/>*
> ResearchGate profile: https://www.researchgate.net/profile/Md_Uddin125
> YouTube channel:https://www.youtube.com/jalalsdream
> Cell: +8613260859092, +8801792052662
> _______________________________________________
> 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/20220608/876999d7/attachment.htm>

From dmjalal90 at gmail.com  Wed Jun  8 17:59:13 2022
From: dmjalal90 at gmail.com (Md. Jalal Uddin)
Date: Thu, 9 Jun 2022 07:59:13 +0800
Subject: [ncl-talk] Need a clarification for "getMidMonDay" function in
 NCL
In-Reply-To: <CAOF1d_4yEqkqLe=NKO7W8FREQL0HuZuQkgc3KSTzb4sj0cdxJg@mail.gmail.com>
References: <CAN5hBbthSq-dubxj7Bd3-CyodvuOg=EikJXY5-T_uArE82eVqg@mail.gmail.com>
 <CAOF1d_4yEqkqLe=NKO7W8FREQL0HuZuQkgc3KSTzb4sj0cdxJg@mail.gmail.com>
Message-ID: <CAN5hBbtKws0ezGnm9ikTh7TUFpXOM8jS-e7cts37j3p9Je8Fdg@mail.gmail.com>

Thank you so much, sir. I am using NCL 6.6.2 on Windows 10. I could not
load the ?crop.ncl? file. Could I have this file?

I am sorry for the inconvenience caused.

??:?? PM ? ???, ? ???, ???? ? Dennis Shea <shea at ucar.edu> ???????:

> Oops. My bad!  I did not document the function.
> getMidMonDay is a utility function contained within the 'crop' library
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/crop.ncl"
> ====
> The code is 'simple'
>
> undef("getMidMonDay")
> function getMidMonDay(time:numeric)
> ;
> ; utility function
> ; The input 'time' (seconds/hours/days since ...") has the middle of the
> month
> ; This returns the (*approximate*) 'middle day' of the current month as an
> integer
> ; "Approximate" because day 15 is used for all months ...
> ;      rather than 14 (Feb), 14.5 (Leap year Feb) or 15.5 (31 day months).
> ; This is more than adequate.
>
> local date, yyyy, mm, dd, jday
> begin
>     date = cd_calendar(time, 0)
>     yyyy = toint(date(:,0))
>     mm   = toint(date(:,1))
>
>     dd   = conform(mm, 15, -1)    ; approx mid month
>
>     jday = day_of_year(yyyy, mm, dd)
>
>     copy_VarCoords(time, jday)
>     jday at long_name = "mid-month day of current year"
>     return(jday)
> end
>
> On Wed, Jun 8, 2022 at 3:43 AM Md. Jalal Uddin via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Hi all,
>>
>> The fao56_1.ncl (
>> https://www.ncl.ucar.edu/Applications/Scripts/fao56_1.ncl) script
>> used "getMidMonDay" function to get the middle day of the month, but it
>> seems that this function is not exist.
>>
>> How to get the middle day of the month to compute evapotranspiration? Is
>> it similar to the "day_of_year" function (
>> https://www.ncl.ucar.edu/Document/Functions/Built-in/day_of_year.shtml)?
>>
>> Regards,
>> Jalal
>> --
>> *Md. Jalal Uddin*
>> Founder and Director of Research Society
>> PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
>> M.Sc. in Applied Meteorology (NUIST)
>> B.Sc. in Disaster Management (PSTU)
>> Email: founder-and-director at researchsociety20.org,
>> 20205103002 at nuist.edu.cn <dmjalal90 at nuist.edu.cn>
>> Website: *https://researchsociety20.org/founder-and-director/
>> <https://researchsociety20.org/founder-and-director/>*
>> ResearchGate profile: https://www.researchgate.net/profile/Md_Uddin125
>> YouTube channel:https://www.youtube.com/jalalsdream
>> Cell: +8613260859092, +8801792052662
>>
> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
> --
*Md. Jalal Uddin*
Founder and Director of Research Society
PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
M.Sc. in Applied Meteorology (NUIST)
B.Sc. in Disaster Management (PSTU)
Email: founder-and-director at researchsociety20.org, 20205103002 at nuist.edu.cn
<dmjalal90 at nuist.edu.cn>
Website: *https://researchsociety20.org/founder-and-director/
<https://researchsociety20.org/founder-and-director/>*
ResearchGate profile: https://www.researchgate.net/profile/Md_Uddin125
YouTube channel:https://www.youtube.com/jalalsdream
Cell: +8613260859092, +8801792052662
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220609/85e4c03a/attachment.htm>

From dmjalal90 at gmail.com  Wed Jun  8 19:33:09 2022
From: dmjalal90 at gmail.com (Md. Jalal Uddin)
Date: Thu, 9 Jun 2022 09:33:09 +0800
Subject: [ncl-talk] Need a clarification for "getMidMonDay" function in
 NCL
In-Reply-To: <CAN5hBbtKws0ezGnm9ikTh7TUFpXOM8jS-e7cts37j3p9Je8Fdg@mail.gmail.com>
References: <CAN5hBbthSq-dubxj7Bd3-CyodvuOg=EikJXY5-T_uArE82eVqg@mail.gmail.com>
 <CAOF1d_4yEqkqLe=NKO7W8FREQL0HuZuQkgc3KSTzb4sj0cdxJg@mail.gmail.com>
 <CAN5hBbtKws0ezGnm9ikTh7TUFpXOM8jS-e7cts37j3p9Je8Fdg@mail.gmail.com>
Message-ID: <CAN5hBbumRQw76gFpfC2v6b_fHjv6LRbG7F8aBqo3mpEewGkriQ@mail.gmail.com>

I got the ?crop.ncl? file. Thank you so much for your continuous support.

On Thu, 9 Jun 2022 at 07:59, Md. Jalal Uddin <dmjalal90 at gmail.com> wrote:

> Thank you so much, sir. I am using NCL 6.6.2 on Windows 10. I could not
> load the ?crop.ncl? file. Could I have this file?
>
> I am sorry for the inconvenience caused.
>
> ??:?? PM ? ???, ? ???, ???? ? Dennis Shea <shea at ucar.edu> ???????:
>
>> Oops. My bad!  I did not document the function.
>> getMidMonDay is a utility function contained within the 'crop' library
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/crop.ncl"
>> ====
>> The code is 'simple'
>>
>> undef("getMidMonDay")
>> function getMidMonDay(time:numeric)
>> ;
>> ; utility function
>> ; The input 'time' (seconds/hours/days since ...") has the middle of the
>> month
>> ; This returns the (*approximate*) 'middle day' of the current month as
>> an integer
>> ; "Approximate" because day 15 is used for all months ...
>> ;      rather than 14 (Feb), 14.5 (Leap year Feb) or 15.5 (31 day months).
>> ; This is more than adequate.
>>
>> local date, yyyy, mm, dd, jday
>> begin
>>     date = cd_calendar(time, 0)
>>     yyyy = toint(date(:,0))
>>     mm   = toint(date(:,1))
>>
>>     dd   = conform(mm, 15, -1)    ; approx mid month
>>
>>     jday = day_of_year(yyyy, mm, dd)
>>
>>     copy_VarCoords(time, jday)
>>     jday at long_name = "mid-month day of current year"
>>     return(jday)
>> end
>>
>> On Wed, Jun 8, 2022 at 3:43 AM Md. Jalal Uddin via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>> Hi all,
>>>
>>> The fao56_1.ncl (
>>> https://www.ncl.ucar.edu/Applications/Scripts/fao56_1.ncl) script
>>> used "getMidMonDay" function to get the middle day of the month, but it
>>> seems that this function is not exist.
>>>
>>> How to get the middle day of the month to compute evapotranspiration? Is
>>> it similar to the "day_of_year" function (
>>> https://www.ncl.ucar.edu/Document/Functions/Built-in/day_of_year.shtml)?
>>>
>>> Regards,
>>> Jalal
>>> --
>>> *Md. Jalal Uddin*
>>> Founder and Director of Research Society
>>> PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
>>> M.Sc. in Applied Meteorology (NUIST)
>>> B.Sc. in Disaster Management (PSTU)
>>> Email: founder-and-director at researchsociety20.org,
>>> 20205103002 at nuist.edu.cn <dmjalal90 at nuist.edu.cn>
>>> Website: *https://researchsociety20.org/founder-and-director/
>>> <https://researchsociety20.org/founder-and-director/>*
>>> ResearchGate profile: https://www.researchgate.net/profile/Md_Uddin125
>>> YouTube channel:https://www.youtube.com/jalalsdream
>>> Cell: +8613260859092, +8801792052662
>>>
>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at mailman.ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>> --
> *Md. Jalal Uddin*
> Founder and Director of Research Society
> PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
> M.Sc. in Applied Meteorology (NUIST)
> B.Sc. in Disaster Management (PSTU)
> Email: founder-and-director at researchsociety20.org,
> 20205103002 at nuist.edu.cn <dmjalal90 at nuist.edu.cn>
> Website: *https://researchsociety20.org/founder-and-director/
> <https://researchsociety20.org/founder-and-director/>*
> ResearchGate profile: https://www.researchgate.net/profile/Md_Uddin125
> YouTube channel:https://www.youtube.com/jalalsdream
> Cell: +8613260859092, +8801792052662
>


-- 
*Md. Jalal Uddin*
Founder and Director of Research Society
PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
M.Sc. in Applied Meteorology (NUIST)
B.Sc. in Disaster Management (PSTU)
Email: founder-and-director at researchsociety20.org, 20205103002 at nuist.edu.cn
<dmjalal90 at nuist.edu.cn>
Website: *https://researchsociety20.org/founder-and-director/
<https://researchsociety20.org/founder-and-director/>*
ResearchGate profile: https://www.researchgate.net/profile/Md_Uddin125
YouTube channel:https://www.youtube.com/jalalsdream
Cell: +8613260859092, +8801792052662
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220609/16685fcc/attachment.htm>

From najibgal at yahoo.com  Thu Jun  9 07:22:48 2022
From: najibgal at yahoo.com (Najib Yusuf)
Date: Thu, 9 Jun 2022 13:22:48 +0000 (UTC)
Subject: [ncl-talk] Multiple_Timeser_annavg_O3CONOX
References: <1098861008.1464498.1654780969019.ref@mail.yahoo.com>
Message-ID: <1098861008.1464498.1654780969019@mail.yahoo.com>

Dear NCL experts and colleagues,
I am trying to plot multiple time series annual averages(monthly files) of three variables from same files using xy2 example, but I am getting error of dimension subscripting, when I tried using only One File it plotted. Kindly assist. See attached the script, and below is the VarSummaries. The error message is?fatal:VarVarWrite: Number of dimensions on left hand side does not match right hand side.?
data(0,:) = var_ann_flip(55,:,:)data(1,:) = var_ann_flip(55,:,:)data(2,:) = var_ann_flip(55,:,:)

Thank you for your support always

Variable: f_all
Type: stringTotal Size: 96 bytes? ? ? ? ? ? 12 valuesNumber of Dimensions: 1Dimensions and sizes:? ?[12]Coordinates:(0)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-01.nc(1)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-02.nc(2)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-03.nc(3)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-04.nc(4)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-05.nc(5)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-06.nc(6)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-07.nc(7)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-08.nc(8)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-09.nc(9)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-10.nc(10)? ? /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-11.nc(11)? ? /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-12.nc(0)? ? ?1.016496e-05(0)? ? ?3.577813e-09
Variable: var_annType: floatTotal Size: 12386304 bytes? ? ? ? ? ? 3096576 valuesNumber of Dimensions: 3Dimensions and sizes:? ?[lev | 56] x [lat | 192] x [lon | 288]Coordinates:? ? ? ? ? ? lev: [1.867879997007549..992.5000106104562]? ? ? ? ? ? lat: [ -90..? 90]? ? ? ? ? ? lon: [? ?0..358.75]Number Of Attributes: 5? cell_methods :? ? ? ? time: mean? long_name :? ?O3 concentration? units :? ? ? ?mol/mol? mdims :? ? ? ?1? average_op_ncl :? ? ? dim_avg_n over dimension(s): time
Variable: var_ann_flipType: floatTotal Size: 12386304 bytes? ? ? ? ? ? 3096576 valuesNumber of Dimensions: 3Dimensions and sizes:? ?[lev | 56] x [lat | 192] x [lon | 288]Coordinates:? ? ? ? ? ? lev: [1.867879997007549..992.5000106104562]? ? ? ? ? ? lat: [ -90..? 90]? ? ? ? ? ? lon: [-180..178.75]Number Of Attributes: 6? cell_methods :? ? ? ? time: mean? long_name :? ?O3 concentration? units :? ? ? ?mol/mol? mdims :? ? ? ?1? average_op_ncl :? ? ? dim_avg_n over dimension(s): time? lonFlip :? ? ?longitude coordinate variable has been reordered via lonFlip
Variable: dataType: floatTotal Size: 2304 bytes? ? ? ? ? ? 576 valuesNumber of Dimensions: 2Dimensions and sizes:? ?[3] x [192]Coordinates:Number Of Attributes: 1? _FillValue :? 9.96921e+36fatal:VarVarWrite: Number of dimensions on left hand side does not match right hand sidefatal:["Execute.c":8637]:Execute: Error occurred at or near line 50 in file Timeser_annavg_O3CONOX.ncl



------------------------------------
Najib Yusuf Galadanci PhDAssistant Director (R & D)
Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.Prophet Muhammad Rasulullah (PBUH)
?


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220609/ea328a31/attachment-0001.htm>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: Najiby_script3.txt
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220609/ea328a31/attachment-0001.txt>

From brownrig at ucar.edu  Thu Jun  9 07:41:32 2022
From: brownrig at ucar.edu (Rick Brownrigg)
Date: Thu, 9 Jun 2022 07:41:32 -0600
Subject: [ncl-talk] Multiple_Timeser_annavg_O3CONOX
In-Reply-To: <1098861008.1464498.1654780969019@mail.yahoo.com>
References: <1098861008.1464498.1654780969019.ref@mail.yahoo.com>
 <1098861008.1464498.1654780969019@mail.yahoo.com>
Message-ID: <CAGKRhbG7G_F5yZuY4Zxw9+KySzW_26Cmn3RwrMNsEoyHjFN2pQ@mail.gmail.com>

Hi,

The problem is just what the error message says; in these lines:

data(0,:) = var_ann_flip(55,:,:)
data(1,:) = var_ann_flip(55,:,:)
data(2,:) = var_ann_flip(55,:,:)

data is [3] x [192]
whereas var_ann_flip is [lev | 56] x [lat | 192] x [lon | 288]

I would speculate that the line that creates data should be something like:

data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)
/),float)

Hope that helps...
Rick


On Thu, Jun 9, 2022 at 7:23 AM Najib Yusuf via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear NCL experts and colleagues,
>
> I am trying to plot multiple time series annual averages(monthly files) of
> three variables from same files using xy2 example, but I am getting error
> of dimension subscripting, when I tried using only One File it plotted.
> Kindly assist. See attached the script, and below is the VarSummaries. The
> error message is fatal:VarVarWrite: Number of dimensions on left hand
> side does not match right hand side.
>
> data(0,:) = var_ann_flip(55,:,:)
> data(1,:) = var_ann_flip(55,:,:)
> data(2,:) = var_ann_flip(55,:,:)
>
>
> Thank you for your support always
>
>
> Variable: f_all
> Type: string
> Total Size: 96 bytes
>             12 values
> Number of Dimensions: 1
> Dimensions and sizes:   [12]
> Coordinates:
> (0)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-01.nc
> (1)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-02.nc
> (2)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-03.nc
> (3)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-04.nc
> (4)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-05.nc
> (5)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-06.nc
> (6)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-07.nc
> (7)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-08.nc
> (8)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-09.nc
> (9)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-10.nc
> (10)    /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-11.nc
> (11)    /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-12.nc
> (0)     1.016496e-05
> (0)     3.577813e-09
>
> Variable: var_ann
> Type: float
> Total Size: 12386304 bytes
>             3096576 values
> Number of Dimensions: 3
> Dimensions and sizes:   [lev | 56] x [lat | 192] x [lon | 288]
> Coordinates:
>             lev: [1.867879997007549..992.5000106104562]
>             lat: [ -90..  90]
>             lon: [   0..358.75]
> Number Of Attributes: 5
>   cell_methods :        time: mean
>   long_name :   O3 concentration
>   units :       mol/mol
>   mdims :       1
>   average_op_ncl :      dim_avg_n over dimension(s): time
>
> Variable: var_ann_flip
> Type: float
> Total Size: 12386304 bytes
>             3096576 values
> Number of Dimensions: 3
> Dimensions and sizes:   [lev | 56] x [lat | 192] x [lon | 288]
> Coordinates:
>             lev: [1.867879997007549..992.5000106104562]
>             lat: [ -90..  90]
>             lon: [-180..178.75]
> Number Of Attributes: 6
>   cell_methods :        time: mean
>   long_name :   O3 concentration
>   units :       mol/mol
>   mdims :       1
>   average_op_ncl :      dim_avg_n over dimension(s): time
>   lonFlip :     longitude coordinate variable has been reordered via
> lonFlip
>
> Variable: data
> Type: float
> Total Size: 2304 bytes
>             576 values
> Number of Dimensions: 2
> Dimensions and sizes:   [3] x [192]
> Coordinates:
> Number Of Attributes: 1
>   _FillValue :  9.96921e+36
> fatal:VarVarWrite: Number of dimensions on left hand side does not match
> right hand side
> fatal:["Execute.c":8637]:Execute: Error occurred at or near line 50 in
> file Timeser_annavg_O3CONOX.ncl
>
>
>
> *------------------------------------*
> *Najib Yusuf Galadanci PhD*
> *Assistant Director (R & D)*
>
> *Lower Atmospheric Dynamics(Aerosol radiative forcing, air quality
> observations and modeling)Centre for Atmospheric ResearchNational Space
> Research and Development Agency. Kogi StateUniversity Campus, Anyigba. Kogi
> State. Nigeria+234 80 3871 4158------------------------------------Know
> that Victory comes with Patience, Relief with Affliction, and Ease with
> Hardship.Prophet Muhammad Rasulullah (PBUH) *
>
>
> _______________________________________________
> 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/20220609/c4e2ec3b/attachment.htm>

From najibgal at yahoo.com  Thu Jun  9 07:55:45 2022
From: najibgal at yahoo.com (Najib Yusuf)
Date: Thu, 9 Jun 2022 13:55:45 +0000 (UTC)
Subject: [ncl-talk] Multiple_Timeser_annavg_O3CONOX
In-Reply-To: <CAGKRhbG7G_F5yZuY4Zxw9+KySzW_26Cmn3RwrMNsEoyHjFN2pQ@mail.gmail.com>
References: <1098861008.1464498.1654780969019.ref@mail.yahoo.com>
 <1098861008.1464498.1654780969019@mail.yahoo.com>
 <CAGKRhbG7G_F5yZuY4Zxw9+KySzW_26Cmn3RwrMNsEoyHjFN2pQ@mail.gmail.com>
Message-ID: <1553082357.1460104.1654782945859@mail.yahoo.com>

Hi Rick,
Yes,?data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)?/),float)
data(0,:) = var_ann_flip(55,:,:)
data(1,:) = var_ann_flip(55,:,:)
data(2,:) = var_ann_flip(55,:,:)left hand side is 2 dimension and right side is 3 dimension but the error is:Number of subscripts on left-hand-side do not match number of dimensions of variable: (2),? Subscripts used: (3)
It is not clear to me how I can make left hand side subscripting to 3 because the dimension is 2. Kindly clear my confusion.
Thank you sir


------------------------------------
Najib Yusuf Galadanci PhDAssistant Director (R & D)
Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.Prophet Muhammad Rasulullah (PBUH)
?

 

    On Thursday, June 9, 2022, 02:41:44 PM GMT+1, Rick Brownrigg <brownrig at ucar.edu> wrote:  
 
 Hi,
The problem is just what the error message says; in these lines:
data(0,:) = var_ann_flip(55,:,:)
data(1,:) = var_ann_flip(55,:,:)
data(2,:) = var_ann_flip(55,:,:)
data is [3] x [192] 
whereas var_ann_flip is [lev | 56] x [lat | 192] x [lon | 288]
I would speculate that the line that creates data should be something like:
data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon) /),float)
Hope that helps...Rick

On Thu, Jun 9, 2022 at 7:23 AM Najib Yusuf via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:

Dear NCL experts and colleagues,
I am trying to plot multiple time series annual averages(monthly files) of three variables from same files using xy2 example, but I am getting error of dimension subscripting, when I tried using only One File it plotted. Kindly assist. See attached the script, and below is the VarSummaries. The error message is?fatal:VarVarWrite: Number of dimensions on left hand side does not match right hand side.?
data(0,:) = var_ann_flip(55,:,:)data(1,:) = var_ann_flip(55,:,:)data(2,:) = var_ann_flip(55,:,:)

Thank you for your support always

Variable: f_all
Type: stringTotal Size: 96 bytes? ? ? ? ? ? 12 valuesNumber of Dimensions: 1Dimensions and sizes:? ?[12]Coordinates:(0)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-01.nc(1)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-02.nc(2)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-03.nc(3)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-04.nc(4)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-05.nc(5)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-06.nc(6)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-07.nc(7)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-08.nc(8)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-09.nc(9)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-10.nc(10)? ? /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-11.nc(11)? ? /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-12.nc(0)? ? ?1.016496e-05(0)? ? ?3.577813e-09
Variable: var_annType: floatTotal Size: 12386304 bytes? ? ? ? ? ? 3096576 valuesNumber of Dimensions: 3Dimensions and sizes:? ?[lev | 56] x [lat | 192] x [lon | 288]Coordinates:? ? ? ? ? ? lev: [1.867879997007549..992.5000106104562]? ? ? ? ? ? lat: [ -90..? 90]? ? ? ? ? ? lon: [? ?0..358.75]Number Of Attributes: 5? cell_methods :? ? ? ? time: mean? long_name :? ?O3 concentration? units :? ? ? ?mol/mol? mdims :? ? ? ?1? average_op_ncl :? ? ? dim_avg_n over dimension(s): time
Variable: var_ann_flipType: floatTotal Size: 12386304 bytes? ? ? ? ? ? 3096576 valuesNumber of Dimensions: 3Dimensions and sizes:? ?[lev | 56] x [lat | 192] x [lon | 288]Coordinates:? ? ? ? ? ? lev: [1.867879997007549..992.5000106104562]? ? ? ? ? ? lat: [ -90..? 90]? ? ? ? ? ? lon: [-180..178.75]Number Of Attributes: 6? cell_methods :? ? ? ? time: mean? long_name :? ?O3 concentration? units :? ? ? ?mol/mol? mdims :? ? ? ?1? average_op_ncl :? ? ? dim_avg_n over dimension(s): time? lonFlip :? ? ?longitude coordinate variable has been reordered via lonFlip
Variable: dataType: floatTotal Size: 2304 bytes? ? ? ? ? ? 576 valuesNumber of Dimensions: 2Dimensions and sizes:? ?[3] x [192]Coordinates:Number Of Attributes: 1? _FillValue :? 9.96921e+36fatal:VarVarWrite: Number of dimensions on left hand side does not match right hand sidefatal:["Execute.c":8637]:Execute: Error occurred at or near line 50 in file Timeser_annavg_O3CONOX.ncl



------------------------------------
Najib Yusuf Galadanci PhDAssistant Director (R & D)
Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.Prophet Muhammad Rasulullah (PBUH)
?


_______________________________________________
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/20220609/16043a7c/attachment.htm>

From brownrig at ucar.edu  Thu Jun  9 08:33:50 2022
From: brownrig at ucar.edu (Rick Brownrigg)
Date: Thu, 9 Jun 2022 08:33:50 -0600
Subject: [ncl-talk] Multiple_Timeser_annavg_O3CONOX
In-Reply-To: <1553082357.1460104.1654782945859@mail.yahoo.com>
References: <1098861008.1464498.1654780969019.ref@mail.yahoo.com>
 <1098861008.1464498.1654780969019@mail.yahoo.com>
 <CAGKRhbG7G_F5yZuY4Zxw9+KySzW_26Cmn3RwrMNsEoyHjFN2pQ@mail.gmail.com>
 <1553082357.1460104.1654782945859@mail.yahoo.com>
Message-ID: <CAGKRhbEYMqpuTmXjBksb5aLHQ+18E4dBDgXymGguJTyLbW=R6Q@mail.gmail.com>

I don't really understand what you are trying to do conceptually, but a few
comments:

- In assigning one array to another, the size and shapes must be the same.
It looks like you are trying to take 3 slices of var_ann_flip at level 55
and assign it to the variable data.  Those slices would each be a 2D
subarray of dimension lat x lon.  In that case, data needs to be
dimensioned as I suggested:

    data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)
 /),float)

and then the assignments to data become:

    data(0,:,:) = var_ann_flip(55,:,:)
    data(1,:,:) = var_ann_flip(55,:,:)
    data(2,:,:) = var_ann_flip(55,:,:)

- That takes care of the error message, but do you really mean for data
0..2 to contain the same slice (55) of var_ann_flip?  Should they be
different slices?

- All of this is contained inside a do-loop of three iterations. Is this
correct? Were you trying to build up the data and plots iteratively?

- Further down in the script, there is this:

   plot0  = gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)
   plot1  = gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)
   plot2  = gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)

Note that the three plots are of the same data(0,:,:)  -- should that be
data(0,:,:)  data(1,:,:) data(2,:,:)?

Rick


On Thu, Jun 9, 2022 at 7:55 AM Najib Yusuf <najibgal at yahoo.com> wrote:

> Hi Rick,
>
> Yes, data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)
>  /),float)
>
> data(0,:) = var_ann_flip(55,:,:)
> data(1,:) = var_ann_flip(55,:,:)
> data(2,:) = var_ann_flip(55,:,:)
> left hand side is 2 dimension and right side is 3 dimension but the error
> is
> :Number of subscripts on left-hand-side do not match number of dimensions
> of variable: (2),  Subscripts used: (3)
>
> It is not clear to me how I can make left hand side subscripting to 3
> because the dimension is 2. Kindly clear my confusion.
>
> Thank you sir
>
>
> *------------------------------------*
> *Najib Yusuf Galadanci PhD*
> *Assistant Director (R & D)*
>
> *Lower Atmospheric Dynamics(Aerosol radiative forcing, air quality
> observations and modeling)Centre for Atmospheric ResearchNational Space
> Research and Development Agency. Kogi StateUniversity Campus, Anyigba. Kogi
> State. Nigeria+234 80 3871 4158------------------------------------Know
> that Victory comes with Patience, Relief with Affliction, and Ease with
> Hardship.Prophet Muhammad Rasulullah (PBUH) *
>
>
>
> On Thursday, June 9, 2022, 02:41:44 PM GMT+1, Rick Brownrigg <
> brownrig at ucar.edu> wrote:
>
>
> Hi,
>
> The problem is just what the error message says; in these lines:
>
> data(0,:) = var_ann_flip(55,:,:)
> data(1,:) = var_ann_flip(55,:,:)
> data(2,:) = var_ann_flip(55,:,:)
>
> data is [3] x [192]
> whereas var_ann_flip is [lev | 56] x [lat | 192] x [lon | 288]
>
> I would speculate that the line that creates data should be something like:
>
> data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)
> /),float)
>
> Hope that helps...
> Rick
>
>
> On Thu, Jun 9, 2022 at 7:23 AM Najib Yusuf via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
> Dear NCL experts and colleagues,
>
> I am trying to plot multiple time series annual averages(monthly files) of
> three variables from same files using xy2 example, but I am getting error
> of dimension subscripting, when I tried using only One File it plotted.
> Kindly assist. See attached the script, and below is the VarSummaries. The
> error message is fatal:VarVarWrite: Number of dimensions on left hand
> side does not match right hand side.
>
> data(0,:) = var_ann_flip(55,:,:)
> data(1,:) = var_ann_flip(55,:,:)
> data(2,:) = var_ann_flip(55,:,:)
>
>
> Thank you for your support always
>
>
> Variable: f_all
> Type: string
> Total Size: 96 bytes
>             12 values
> Number of Dimensions: 1
> Dimensions and sizes:   [12]
> Coordinates:
> (0)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-01.nc
> (1)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-02.nc
> (2)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-03.nc
> (3)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-04.nc
> (4)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-05.nc
> (5)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-06.nc
> (6)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-07.nc
> (7)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-08.nc
> (8)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-09.nc
> (9)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-10.nc
> (10)    /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-11.nc
> (11)    /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-12.nc
> (0)     1.016496e-05
> (0)     3.577813e-09
>
> Variable: var_ann
> Type: float
> Total Size: 12386304 bytes
>             3096576 values
> Number of Dimensions: 3
> Dimensions and sizes:   [lev | 56] x [lat | 192] x [lon | 288]
> Coordinates:
>             lev: [1.867879997007549..992.5000106104562]
>             lat: [ -90..  90]
>             lon: [   0..358.75]
> Number Of Attributes: 5
>   cell_methods :        time: mean
>   long_name :   O3 concentration
>   units :       mol/mol
>   mdims :       1
>   average_op_ncl :      dim_avg_n over dimension(s): time
>
> Variable: var_ann_flip
> Type: float
> Total Size: 12386304 bytes
>             3096576 values
> Number of Dimensions: 3
> Dimensions and sizes:   [lev | 56] x [lat | 192] x [lon | 288]
> Coordinates:
>             lev: [1.867879997007549..992.5000106104562]
>             lat: [ -90..  90]
>             lon: [-180..178.75]
> Number Of Attributes: 6
>   cell_methods :        time: mean
>   long_name :   O3 concentration
>   units :       mol/mol
>   mdims :       1
>   average_op_ncl :      dim_avg_n over dimension(s): time
>   lonFlip :     longitude coordinate variable has been reordered via
> lonFlip
>
> Variable: data
> Type: float
> Total Size: 2304 bytes
>             576 values
> Number of Dimensions: 2
> Dimensions and sizes:   [3] x [192]
> Coordinates:
> Number Of Attributes: 1
>   _FillValue :  9.96921e+36
> fatal:VarVarWrite: Number of dimensions on left hand side does not match
> right hand side
> fatal:["Execute.c":8637]:Execute: Error occurred at or near line 50 in
> file Timeser_annavg_O3CONOX.ncl
>
>
>
> *------------------------------------*
> *Najib Yusuf Galadanci PhD*
> *Assistant Director (R & D)*
>
> Lower Atmospheric Dynamics
>
> (Aerosol radiative forcing, air quality observations and modeling)
>
> Centre for Atmospheric Research
>
> National Space Research and Development Agency. Kogi State
>
> University Campus, Anyigba. Kogi State. Nigeria
>
> +234 80 3871 4158
> ------------------------------------
> *Know that Victory comes with Patience, *
> *Relief with Affliction, and Ease with Hardship.*
> *Prophet Muhammad Rasulullah (PBUH)*
>
>
>
>
> _______________________________________________
> 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/20220609/ab1b992c/attachment.htm>

From najibgal at yahoo.com  Fri Jun 10 05:07:07 2022
From: najibgal at yahoo.com (Najib Yusuf)
Date: Fri, 10 Jun 2022 11:07:07 +0000 (UTC)
Subject: [ncl-talk] Multiple_Timeser_annavg_O3CONOX
In-Reply-To: <CAGKRhbFLAzuTn9JyW_WXvfnyK_EtqpN5RKouSzha38WMq5-4zg@mail.gmail.com>
References: <1098861008.1464498.1654780969019.ref@mail.yahoo.com>
 <1098861008.1464498.1654780969019@mail.yahoo.com>
 <CAGKRhbG7G_F5yZuY4Zxw9+KySzW_26Cmn3RwrMNsEoyHjFN2pQ@mail.gmail.com>
 <1553082357.1460104.1654782945859@mail.yahoo.com>
 <CAGKRhbEYMqpuTmXjBksb5aLHQ+18E4dBDgXymGguJTyLbW=R6Q@mail.gmail.com>
 <1350611259.1494282.1654786531443@mail.yahoo.com>
 <CAGKRhbFLAzuTn9JyW_WXvfnyK_EtqpN5RKouSzha38WMq5-4zg@mail.gmail.com>
Message-ID: <666573137.1916224.1654859227328@mail.yahoo.com>

Hi Rick,
Thank you for your help. It is so muddled for me. Now I got another error even though var_ann_flip and data both have 3 dimensions but leftmost dimension size is not same. How can I take care of this please. Thank you?(0)? ? ?gsn_csm_xy: Fatal: X and Y must have the same dimensions sizes, or one must be one-dimensional and both have the same rightmost dimension.

Variable: var_ann_flipType: floatTotal Size: 12386304 bytes? ? ? ? ? ? 3096576 valuesNumber of Dimensions: 3Dimensions and sizes:? ?[lev | 56] x [lat | 192] x [lon | 288]Coordinates:? ? ? ? ? ? lev: [1.867879997007549..992.5000106104562]? ? ? ? ? ? lat: [ -90..? 90]? ? ? ? ? ? lon: [-180..178.75]Number Of Attributes: 6? cell_methods :? ? ? ? time: mean? long_name :? ?O3 concentration? units :? ? ? ?mol/mol? mdims :? ? ? ?1? average_op_ncl :? ? ? dim_avg_n over dimension(s): time? lonFlip :? ? ?longitude coordinate variable has been reordered via lonFlip
Variable: dataType: floatTotal Size: 663552 bytes? ? ? ? ? ? 165888 valuesNumber of Dimensions: 3Dimensions and sizes:? ?[3] x [192] x [288]Coordinates:Number Of Attributes: 1? _FillValue :? 9.96921e+36
Variable: data (subsection)Type: floatTotal Size: 221184 bytes? ? ? ? ? ? 55296 valuesNumber of Dimensions: 2Dimensions and sizes:? ?[lat | 192] x [lon | 288]Coordinates:? ? ? ? ? ? lat: [ -90..? 90]? ? ? ? ? ? lon: [-180..178.75]Number Of Attributes: 8? lev : 992.5000106104562? _FillValue :? 9.96921e+36? cell_methods :? ? ? ? time: mean? long_name :? ?O3 concentration? units :? ? ? ?mol/mol? mdims :? ? ? ?1? average_op_ncl :? ? ? dim_avg_n over dimension(s): time? lonFlip :? ? ?longitude coordinate variable has been reordered via lonFlip(0)? ? ?gsn_csm_xy: Fatal: X and Y must have the same dimensions sizes, or one must be one-dimensional and both have the same rightmost dimension.



------------------------------------
Najib Yusuf Galadanci PhDAssistant Director (R & D)
Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.Prophet Muhammad Rasulullah (PBUH)
?

 

    On Thursday, June 9, 2022, 05:11:51 PM GMT+1, Rick Brownrigg <brownrig at ucar.edu> wrote:  
 
 Hi Najib,
OK, I see now that each time through the loop, a different variable is pulled from the files into the variable "var", from which the variable "var_ann_flip" gets constructed.? If so, then these statements don't make sense:
data(0,:) = var_ann_flip(55,:,:)
data(1,:) = var_ann_flip(55,:,:)
data(2,:) = var_ann_flip(55,:,:)
each time through the loop, all three data 0..2 get overwritten by the same "var". You probably want something like:
?? data(i,:,:) = var_ann_flip(55,:,:)
But that mean the plotting is wrong too; rather than trying to plot all three plots *each time through the loop*, you probably want to generate just one:
??? plot? = gsn_csm_xy (wks,var_ann&lat,data(i,:,:),res)
This all has implications for labelling and colors -- they too should likely be based off of indexing into an arrays of labels/colors based upon "i".
The logic is kind of muddled. Try to think about what needs to change each time through the loop and what values/variables you have in hand at each iteration. 

Hope that helps...Rick


On Thu, Jun 9, 2022 at 8:54 AM Najib Yusuf <najibgal at yahoo.com> wrote:


Hi Rick again,
I am trying to plot annual average of O3, CO and NOX this is what I wanted to represented with this variable 0,1,2.
? ? data(0,:,:)?
??? data(1,:,:)?
??? data(2,:,:)
all the three variables are surface level 55.?
So my annual average is?var_ann_flip, I thought to represent each variable as?
? ? data(0,:,:) = var_ann_flip(55,:,:)
??? data(1,:,:) = var_ann_flip(55,:,:)
??? data(2,:,:) = var_ann_flip(55,:,:)
I don't really understand what you are trying to do conceptually, but a few comments:
- In assigning one array to another, the size and shapes must be the same. It looks like you are trying to take 3 slices of var_ann_flip at level 55 and assign it to the variable data.? Those slices would each be a 2D subarray of dimension lat x lon.? In that case, data needs to be dimensioned as I suggested:
Yes
????data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)?/),float)
and then the assignments to data become:
??? data(0,:,:) = var_ann_flip(55,:,:)
??? data(1,:,:) = var_ann_flip(55,:,:)
??? data(2,:,:) = var_ann_flip(55,:,:)
- That takes care of the error message, but do you really mean for data 0..2 to contain the same slice (55) of var_ann_flip?? Should they be different slices??Yes0..2 are the variables O3, CO? and NOX all at the surface level 55.
- All of this is contained inside a do-loop of three iterations. Is this correct? Were you trying to build up the data and plots iteratively?
? ?Yes I tried to loop.- Further down in the script, there is this:
?? plot0 ?= gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)
?? plot1 ?= gsn_csm_xy (wks,var_ann&lat,data(1,:,:),res)
?? plot2 ?= gsn_csm_xy (wks,var_ann&lat,data(2,:,:),res)
Note that the three plots are of the same data(0,:,:)? -- should that be data(0,:,:)? data(1,:,:) data(2,:,:)?

Thank you. They should be as you pointed out; data(0,:,:) for O3, data(1,:,:) for CO and data(2,;,;) for NOX.
Please help. Thank you so much
Najib



------------------------------------
Najib Yusuf Galadanci PhDAssistant Director (R & D)
Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.Prophet Muhammad Rasulullah (PBUH)
?

 

    On Thursday, June 9, 2022, 03:34:02 PM GMT+1, Rick Brownrigg <brownrig at ucar.edu> wrote:  
 
 I don't really understand what you are trying to do conceptually, but a few comments:
- In assigning one array to another, the size and shapes must be the same. It looks like you are trying to take 3 slices of var_ann_flip at level 55 and assign it to the variable data.? Those slices would each be a 2D subarray of dimension lat x lon.? In that case, data needs to be dimensioned as I suggested:
??? data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)?/),float)
and then the assignments to data become:
??? data(0,:,:) = var_ann_flip(55,:,:)
??? data(1,:,:) = var_ann_flip(55,:,:)
??? data(2,:,:) = var_ann_flip(55,:,:)
- That takes care of the error message, but do you really mean for data 0..2 to contain the same slice (55) of var_ann_flip?? Should they be different slices?
- All of this is contained inside a do-loop of three iterations. Is this correct? Were you trying to build up the data and plots iteratively?
- Further down in the script, there is this:
?? plot0 ?= gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)
?? plot1 ?= gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)
?? plot2 ?= gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)
Note that the three plots are of the same data(0,:,:)? -- should that be data(0,:,:)? data(1,:,:) data(2,:,:)?

Rick


On Thu, Jun 9, 2022 at 7:55 AM Najib Yusuf <najibgal at yahoo.com> wrote:

Hi Rick,
Yes,?data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)?/),float)
data(0,:) = var_ann_flip(55,:,:)
data(1,:) = var_ann_flip(55,:,:)
data(2,:) = var_ann_flip(55,:,:)left hand side is 2 dimension and right side is 3 dimension but the error is:Number of subscripts on left-hand-side do not match number of dimensions of variable: (2),? Subscripts used: (3)
It is not clear to me how I can make left hand side subscripting to 3 because the dimension is 2. Kindly clear my confusion.
Thank you sir


------------------------------------
Najib Yusuf Galadanci PhDAssistant Director (R & D)
Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.Prophet Muhammad Rasulullah (PBUH)
?

 

    On Thursday, June 9, 2022, 02:41:44 PM GMT+1, Rick Brownrigg <brownrig at ucar.edu> wrote:  
 
 Hi,
The problem is just what the error message says; in these lines:
data(0,:) = var_ann_flip(55,:,:)
data(1,:) = var_ann_flip(55,:,:)
data(2,:) = var_ann_flip(55,:,:)
data is [3] x [192] 
whereas var_ann_flip is [lev | 56] x [lat | 192] x [lon | 288]
I would speculate that the line that creates data should be something like:
data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon) /),float)
Hope that helps...Rick

On Thu, Jun 9, 2022 at 7:23 AM Najib Yusuf via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:

Dear NCL experts and colleagues,
I am trying to plot multiple time series annual averages(monthly files) of three variables from same files using xy2 example, but I am getting error of dimension subscripting, when I tried using only One File it plotted. Kindly assist. See attached the script, and below is the VarSummaries. The error message is?fatal:VarVarWrite: Number of dimensions on left hand side does not match right hand side.?
data(0,:) = var_ann_flip(55,:,:)data(1,:) = var_ann_flip(55,:,:)data(2,:) = var_ann_flip(55,:,:)

Thank you for your support always

Variable: f_all
Type: stringTotal Size: 96 bytes? ? ? ? ? ? 12 valuesNumber of Dimensions: 1Dimensions and sizes:? ?[12]Coordinates:(0)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-01.nc(1)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-02.nc(2)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-03.nc(3)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-04.nc(4)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-05.nc(5)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-06.nc(6)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-07.nc(7)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-08.nc(8)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-09.nc(9)? ? ?/glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-10.nc(10)? ? /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-11.nc(11)? ? /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/f.e20.FCSD.cam.h0.2013-12.nc(0)? ? ?1.016496e-05(0)? ? ?3.577813e-09
Variable: var_annType: floatTotal Size: 12386304 bytes? ? ? ? ? ? 3096576 valuesNumber of Dimensions: 3Dimensions and sizes:? ?[lev | 56] x [lat | 192] x [lon | 288]Coordinates:? ? ? ? ? ? lev: [1.867879997007549..992.5000106104562]? ? ? ? ? ? lat: [ -90..? 90]? ? ? ? ? ? lon: [? ?0..358.75]Number Of Attributes: 5? cell_methods :? ? ? ? time: mean? long_name :? ?O3 concentration? units :? ? ? ?mol/mol? mdims :? ? ? ?1? average_op_ncl :? ? ? dim_avg_n over dimension(s): time
Variable: var_ann_flipType: floatTotal Size: 12386304 bytes? ? ? ? ? ? 3096576 valuesNumber of Dimensions: 3Dimensions and sizes:? ?[lev | 56] x [lat | 192] x [lon | 288]Coordinates:? ? ? ? ? ? lev: [1.867879997007549..992.5000106104562]? ? ? ? ? ? lat: [ -90..? 90]? ? ? ? ? ? lon: [-180..178.75]Number Of Attributes: 6? cell_methods :? ? ? ? time: mean? long_name :? ?O3 concentration? units :? ? ? ?mol/mol? mdims :? ? ? ?1? average_op_ncl :? ? ? dim_avg_n over dimension(s): time? lonFlip :? ? ?longitude coordinate variable has been reordered via lonFlip
Variable: dataType: floatTotal Size: 2304 bytes? ? ? ? ? ? 576 valuesNumber of Dimensions: 2Dimensions and sizes:? ?[3] x [192]Coordinates:Number Of Attributes: 1? _FillValue :? 9.96921e+36fatal:VarVarWrite: Number of dimensions on left hand side does not match right hand sidefatal:["Execute.c":8637]:Execute: Error occurred at or near line 50 in file Timeser_annavg_O3CONOX.ncl



------------------------------------
Najib Yusuf Galadanci PhDAssistant Director (R & D)
Lower Atmospheric Dynamics

(Aerosol radiative forcing, air quality observations and modeling)

Centre for Atmospheric Research

National Space Research and Development Agency. Kogi State

University Campus, Anyigba.?Kogi State. Nigeria

+234 80 3871 4158
------------------------------------Know that Victory comes with Patience,?Relief with Affliction, and Ease with Hardship.Prophet Muhammad Rasulullah (PBUH)
?


_______________________________________________
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/20220610/361f3775/attachment.htm>

From brownrig at ucar.edu  Fri Jun 10 14:35:52 2022
From: brownrig at ucar.edu (Rick Brownrigg)
Date: Fri, 10 Jun 2022 14:35:52 -0600
Subject: [ncl-talk] Multiple_Timeser_annavg_O3CONOX
In-Reply-To: <666573137.1916224.1654859227328@mail.yahoo.com>
References: <1098861008.1464498.1654780969019.ref@mail.yahoo.com>
 <1098861008.1464498.1654780969019@mail.yahoo.com>
 <CAGKRhbG7G_F5yZuY4Zxw9+KySzW_26Cmn3RwrMNsEoyHjFN2pQ@mail.gmail.com>
 <1553082357.1460104.1654782945859@mail.yahoo.com>
 <CAGKRhbEYMqpuTmXjBksb5aLHQ+18E4dBDgXymGguJTyLbW=R6Q@mail.gmail.com>
 <1350611259.1494282.1654786531443@mail.yahoo.com>
 <CAGKRhbFLAzuTn9JyW_WXvfnyK_EtqpN5RKouSzha38WMq5-4zg@mail.gmail.com>
 <666573137.1916224.1654859227328@mail.yahoo.com>
Message-ID: <CAGKRhbGwbwdVBKW7ygBxxecsW7T20sYVR-fps=sg=j2SofZwXw@mail.gmail.com>

Hi,

In your original script, you had something like this:

   plot0  = gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)

data(0,:,:) is a slice into a 3D variable, and thus a 2D variable. Plotting
latitude vs. a 2D field doesn't make sense for an XY plot. Do you want a
contour plot instead?  Or somehow  you'll need to reduce data(0,:,:) to a
1D variable.

Rick

On Fri, Jun 10, 2022 at 5:06 AM Najib Yusuf <najibgal at yahoo.com> wrote:

> Hi Rick,
>
> Thank you for your help. It is so muddled for me. Now I got another error
> even though var_ann_flip and data both have 3 dimensions but leftmost
> dimension size is not same. How can I take care of this please. Thank you
> (0)     gsn_csm_xy: Fatal: X and Y must have the same dimensions sizes, or
> one must be one-dimensional and both have the same rightmost dimension.
>
>
> Variable: var_ann_flip
> Type: float
> Total Size: 12386304 bytes
>             3096576 values
> Number of Dimensions: 3
> Dimensions and sizes:   [lev | 56] x [lat | 192] x [lon | 288]
> Coordinates:
>             lev: [1.867879997007549..992.5000106104562]
>             lat: [ -90..  90]
>             lon: [-180..178.75]
> Number Of Attributes: 6
>   cell_methods :        time: mean
>   long_name :   O3 concentration
>   units :       mol/mol
>   mdims :       1
>   average_op_ncl :      dim_avg_n over dimension(s): time
>   lonFlip :     longitude coordinate variable has been reordered via
> lonFlip
>
> Variable: data
> Type: float
> Total Size: 663552 bytes
>             165888 values
> Number of Dimensions: 3
> Dimensions and sizes:   [3] x [192] x [288]
> Coordinates:
> Number Of Attributes: 1
>   _FillValue :  9.96921e+36
>
> Variable: data (subsection)
> Type: float
> Total Size: 221184 bytes
>             55296 values
> Number of Dimensions: 2
> Dimensions and sizes:   [lat | 192] x [lon | 288]
> Coordinates:
>             lat: [ -90..  90]
>             lon: [-180..178.75]
> Number Of Attributes: 8
>   lev : 992.5000106104562
>   _FillValue :  9.96921e+36
>   cell_methods :        time: mean
>   long_name :   O3 concentration
>   units :       mol/mol
>   mdims :       1
>   average_op_ncl :      dim_avg_n over dimension(s): time
>   lonFlip :     longitude coordinate variable has been reordered via
> lonFlip
> (0)     gsn_csm_xy: Fatal: X and Y must have the same dimensions sizes, or
> one must be one-dimensional and both have the same rightmost dimension.
>
>
>
> *------------------------------------*
> *Najib Yusuf Galadanci PhD*
> *Assistant Director (R & D)*
>
> *Lower Atmospheric Dynamics(Aerosol radiative forcing, air quality
> observations and modeling)Centre for Atmospheric ResearchNational Space
> Research and Development Agency. Kogi StateUniversity Campus, Anyigba. Kogi
> State. Nigeria+234 80 3871 4158------------------------------------Know
> that Victory comes with Patience, Relief with Affliction, and Ease with
> Hardship.Prophet Muhammad Rasulullah (PBUH) *
>
>
>
> On Thursday, June 9, 2022, 05:11:51 PM GMT+1, Rick Brownrigg <
> brownrig at ucar.edu> wrote:
>
>
> Hi Najib,
>
> OK, I see now that each time through the loop, a different variable is
> pulled from the files into the variable "var", from which the variable
> "var_ann_flip" gets constructed.  If so, then these statements don't make
> sense:
>
> data(0,:) = var_ann_flip(55,:,:)
> data(1,:) = var_ann_flip(55,:,:)
> data(2,:) = var_ann_flip(55,:,:)
>
> each time through the loop, all three data 0..2 get overwritten by the
> same "var". You probably want something like:
>
>    data(i,:,:) = var_ann_flip(55,:,:)
>
> But that mean the plotting is wrong too; rather than trying to plot all
> three plots *each time through the loop*, you probably want to generate
> just one:
>
>     plot  = gsn_csm_xy (wks,var_ann&lat,data(i,:,:),res)
>
> This all has implications for labelling and colors -- they too should
> likely be based off of indexing into an arrays of labels/colors based upon
> "i".
>
> The logic is kind of muddled. Try to think about what needs to change each
> time through the loop and what values/variables you have in hand at each
> iteration.
>
> Hope that helps...
> Rick
>
>
>
> On Thu, Jun 9, 2022 at 8:54 AM Najib Yusuf <najibgal at yahoo.com> wrote:
>
>
> *Hi Rick again,*
>
> *I am trying to plot annual average of O3, CO and NOX this is what I
> wanted to represented with this variable 0,1,2.*
>
>
>
> *    data(0,:,:)     data(1,:,:)     data(2,:,:)*
>
> *all the three variables are surface level 55. *
>
> *So my annual average is var_ann_flip, I thought to represent each
> variable as *
>
>
>
> *    data(0,:,:) = var_ann_flip(55,:,:)    data(1,:,:) =
> var_ann_flip(55,:,:)    data(2,:,:) = var_ann_flip(55,:,:)*
>
> I don't really understand what you are trying to do conceptually, but a
> few comments:
>
> - In assigning one array to another, the size and shapes must be the same.
> It looks like you are trying to take 3 slices of var_ann_flip at level 55
> and assign it to the variable data.  Those slices would each be a 2D
> subarray of dimension lat x lon.  In that case, data needs to be
> dimensioned as I suggested:
>
> *Yes*
>
>     data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)
>  /),float)
>
> and then the assignments to data become:
>
>     data(0,:,:) = var_ann_flip(55,:,:)
>     data(1,:,:) = var_ann_flip(55,:,:)
>     data(2,:,:) = var_ann_flip(55,:,:)
>
> - That takes care of the error message, but do you really mean for data
> 0..2 to contain the same slice (55) of var_ann_flip?  Should they be
> different slices?
> *Yes*
> *0..2 are the variables O3, CO  and NOX all at the surface level 55.*
>
> - All of this is contained inside a do-loop of three iterations. Is this
> correct? Were you trying to build up the data and plots iteratively?
>
>   * Yes I tried to loop.*
> - Further down in the script, there is this:
>
>    plot0  = gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)
>    plot1  = gsn_csm_xy (wks,var_ann&lat,data(1,:,:),res)
>    plot2  = gsn_csm_xy (wks,var_ann&lat,data(2,:,:),res)
>
> Note that the three plots are of the same data(0,:,:)  -- should that be
> data(0,:,:)  data(1,:,:) data(2,:,:)?
>
> *Thank you. They should be as you pointed out; data(0,:,:) for O3,
> data(1,:,:) for CO and data(2,;,;) for NOX.*
>
> *Please help. Thank you so much*
>
> *Najib*
>
>
>
> *------------------------------------*
> *Najib Yusuf Galadanci PhD*
> *Assistant Director (R & D)*
>
> Lower Atmospheric Dynamics
>
> (Aerosol radiative forcing, air quality observations and modeling)
>
> Centre for Atmospheric Research
>
> National Space Research and Development Agency. Kogi State
>
> University Campus, Anyigba. Kogi State. Nigeria
>
> +234 80 3871 4158
> ------------------------------------
> *Know that Victory comes with Patience, *
> *Relief with Affliction, and Ease with Hardship.*
> *Prophet Muhammad Rasulullah (PBUH)*
>
>
>
>
>
> On Thursday, June 9, 2022, 03:34:02 PM GMT+1, Rick Brownrigg <
> brownrig at ucar.edu> wrote:
>
>
> I don't really understand what you are trying to do conceptually, but a
> few comments:
>
> - In assigning one array to another, the size and shapes must be the same.
> It looks like you are trying to take 3 slices of var_ann_flip at level 55
> and assign it to the variable data.  Those slices would each be a 2D
> subarray of dimension lat x lon.  In that case, data needs to be
> dimensioned as I suggested:
>
>     data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)
>  /),float)
>
> and then the assignments to data become:
>
>     data(0,:,:) = var_ann_flip(55,:,:)
>     data(1,:,:) = var_ann_flip(55,:,:)
>     data(2,:,:) = var_ann_flip(55,:,:)
>
> - That takes care of the error message, but do you really mean for data
> 0..2 to contain the same slice (55) of var_ann_flip?  Should they be
> different slices?
>
> - All of this is contained inside a do-loop of three iterations. Is this
> correct? Were you trying to build up the data and plots iteratively?
>
> - Further down in the script, there is this:
>
>    plot0  = gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)
>    plot1  = gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)
>    plot2  = gsn_csm_xy (wks,var_ann&lat,data(0,:,:),res)
>
> Note that the three plots are of the same data(0,:,:)  -- should that be
> data(0,:,:)  data(1,:,:) data(2,:,:)?
>
> Rick
>
>
> On Thu, Jun 9, 2022 at 7:55 AM Najib Yusuf <najibgal at yahoo.com> wrote:
>
> Hi Rick,
>
> Yes, data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)
>  /),float)
>
> data(0,:) = var_ann_flip(55,:,:)
> data(1,:) = var_ann_flip(55,:,:)
> data(2,:) = var_ann_flip(55,:,:)
> left hand side is 2 dimension and right side is 3 dimension but the error
> is
> :Number of subscripts on left-hand-side do not match number of dimensions
> of variable: (2),  Subscripts used: (3)
>
> It is not clear to me how I can make left hand side subscripting to 3
> because the dimension is 2. Kindly clear my confusion.
>
> Thank you sir
>
>
> *------------------------------------*
> *Najib Yusuf Galadanci PhD*
> *Assistant Director (R & D)*
>
> Lower Atmospheric Dynamics
>
> (Aerosol radiative forcing, air quality observations and modeling)
>
> Centre for Atmospheric Research
>
> National Space Research and Development Agency. Kogi State
>
> University Campus, Anyigba. Kogi State. Nigeria
>
> +234 80 3871 4158
> ------------------------------------
> *Know that Victory comes with Patience, *
> *Relief with Affliction, and Ease with Hardship.*
> *Prophet Muhammad Rasulullah (PBUH)*
>
>
>
>
>
> On Thursday, June 9, 2022, 02:41:44 PM GMT+1, Rick Brownrigg <
> brownrig at ucar.edu> wrote:
>
>
> Hi,
>
> The problem is just what the error message says; in these lines:
>
> data(0,:) = var_ann_flip(55,:,:)
> data(1,:) = var_ann_flip(55,:,:)
> data(2,:) = var_ann_flip(55,:,:)
>
> data is [3] x [192]
> whereas var_ann_flip is [lev | 56] x [lat | 192] x [lon | 288]
>
> I would speculate that the line that creates data should be something like:
>
> data = new((/3,dimsizes(var_ann_flip&lat),dimsizes(var_ann_flip&lon)
> /),float)
>
> Hope that helps...
> Rick
>
>
> On Thu, Jun 9, 2022 at 7:23 AM Najib Yusuf via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
> Dear NCL experts and colleagues,
>
> I am trying to plot multiple time series annual averages(monthly files) of
> three variables from same files using xy2 example, but I am getting error
> of dimension subscripting, when I tried using only One File it plotted.
> Kindly assist. See attached the script, and below is the VarSummaries. The
> error message is fatal:VarVarWrite: Number of dimensions on left hand
> side does not match right hand side.
>
> data(0,:) = var_ann_flip(55,:,:)
> data(1,:) = var_ann_flip(55,:,:)
> data(2,:) = var_ann_flip(55,:,:)
>
>
> Thank you for your support always
>
>
> Variable: f_all
> Type: string
> Total Size: 96 bytes
>             12 values
> Number of Dimensions: 1
> Dimensions and sizes:   [12]
> Coordinates:
> (0)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-01.nc
> (1)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-02.nc
> (2)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-03.nc
> (3)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-04.nc
> (4)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-05.nc
> (5)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-06.nc
> (6)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-07.nc
> (7)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-08.nc
> (8)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-09.nc
> (9)     /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-10.nc
> (10)    /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-11.nc
> (11)    /glade/work/najiby/f.e20.FCSD/f.e20.FCSD/atm/hist/
> f.e20.FCSD.cam.h0.2013-12.nc
> (0)     1.016496e-05
> (0)     3.577813e-09
>
> Variable: var_ann
> Type: float
> Total Size: 12386304 bytes
>             3096576 values
> Number of Dimensions: 3
> Dimensions and sizes:   [lev | 56] x [lat | 192] x [lon | 288]
> Coordinates:
>             lev: [1.867879997007549..992.5000106104562]
>             lat: [ -90..  90]
>             lon: [   0..358.75]
> Number Of Attributes: 5
>   cell_methods :        time: mean
>   long_name :   O3 concentration
>   units :       mol/mol
>   mdims :       1
>   average_op_ncl :      dim_avg_n over dimension(s): time
>
> Variable: var_ann_flip
> Type: float
> Total Size: 12386304 bytes
>             3096576 values
> Number of Dimensions: 3
> Dimensions and sizes:   [lev | 56] x [lat | 192] x [lon | 288]
> Coordinates:
>             lev: [1.867879997007549..992.5000106104562]
>             lat: [ -90..  90]
>             lon: [-180..178.75]
> Number Of Attributes: 6
>   cell_methods :        time: mean
>   long_name :   O3 concentration
>   units :       mol/mol
>   mdims :       1
>   average_op_ncl :      dim_avg_n over dimension(s): time
>   lonFlip :     longitude coordinate variable has been reordered via
> lonFlip
>
> Variable: data
> Type: float
> Total Size: 2304 bytes
>             576 values
> Number of Dimensions: 2
> Dimensions and sizes:   [3] x [192]
> Coordinates:
> Number Of Attributes: 1
>   _FillValue :  9.96921e+36
> fatal:VarVarWrite: Number of dimensions on left hand side does not match
> right hand side
> fatal:["Execute.c":8637]:Execute: Error occurred at or near line 50 in
> file Timeser_annavg_O3CONOX.ncl
>
>
>
> *------------------------------------*
> *Najib Yusuf Galadanci PhD*
> *Assistant Director (R & D)*
>
> Lower Atmospheric Dynamics
>
> (Aerosol radiative forcing, air quality observations and modeling)
>
> Centre for Atmospheric Research
>
> National Space Research and Development Agency. Kogi State
>
> University Campus, Anyigba. Kogi State. Nigeria
>
> +234 80 3871 4158
> ------------------------------------
> *Know that Victory comes with Patience, *
> *Relief with Affliction, and Ease with Hardship.*
> *Prophet Muhammad Rasulullah (PBUH)*
>
>
>
>
> _______________________________________________
> 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/20220610/a5792c27/attachment.htm>

From ehsantaghizadeh at yahoo.com  Sun Jun 12 04:16:12 2022
From: ehsantaghizadeh at yahoo.com (Ehsan Taghizadeh)
Date: Sun, 12 Jun 2022 10:16:12 +0000 (UTC)
Subject: [ncl-talk] add_topography
In-Reply-To: <CAOiHCtTP55L+q+vxjJVuyZxb3SQko1ozhEJ6scb5MVB7Crichg@mail.gmail.com>
References: <CAOiHCtTP55L+q+vxjJVuyZxb3SQko1ozhEJ6scb5MVB7Crichg@mail.gmail.com>
Message-ID: <136304501.1472185.1655028972749@mail.yahoo.com>

Dear KunalDid you see "NCL Graphics: Topographic maps"? I suggest download topography data from "Topography and Digital Terrain Data | NCEI", in NetCDF format and then add plot them on your map.

| 
| 
| 
|  |  |

 |

 |
| 
|  | 
Topography and Digital Terrain Data | NCEI

National Geophysical Data Center

The National Geophysical Data Center is involved in several projects to collect public domain digital elevation ...
 |

 |

 |




-----------------------------------------------------------------------------Sincerely,Ehsan Taghizadeh

 

    On Friday, June 3, 2022, 12:40:10 AM GMT+4:30, Kunal Bali via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:  
 
 Hi NCL users,
Could you please give me an idea of how I can?add the topography to the attached plot??

Thanks
---
Kunal Bali






_______________________________________________
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/20220612/f808e231/attachment.htm>

From s1073755 at pm.univpm.it  Mon Jun 13 13:20:43 2022
From: s1073755 at pm.univpm.it (Mauro Morichetti)
Date: Mon, 13 Jun 2022 21:20:43 +0200
Subject: [ncl-talk] (no subject)
Message-ID: <83E924BC-E075-41A3-9484-68B6E759D9A3@pm.univpm.it>

Hi all,

I have a problem with the script attached. I can not figure out, with the text doesn?t show on the scatter plot. I think it is  stupid question, but I going crazy, also because I used with other data without problems.

Thanks


---------

Ph.D. Mauro Morichetti
Department of Industrial Engineering and Mathematical Science
Marche Polytechnic University, Ancona, Italy
m.morichetti at pm.univpm.it
-------------- next part --------------
A non-text attachment was scrubbed...
Name: DAILY.csv
Type: text/csv
Size: 277296 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220613/da22029c/attachment-0001.csv>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: REGRESSION_WRF.ncl
Type: application/octet-stream
Size: 4619 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220613/da22029c/attachment-0001.obj>

From barry.h.lynn at gmail.com  Mon Jun 13 22:23:06 2022
From: barry.h.lynn at gmail.com (Barry Lynn)
Date: Tue, 14 Jun 2022 07:23:06 +0300
Subject: [ncl-talk] (no subject)
In-Reply-To: <83E924BC-E075-41A3-9484-68B6E759D9A3@pm.univpm.it>
References: <83E924BC-E075-41A3-9484-68B6E759D9A3@pm.univpm.it>
Message-ID: <CAN7Bv24ri+Enr3rqZ6xbCM_F9SjgPxQKntRbaXxd_+XkDo_W+g@mail.gmail.com>

Hello Mauro:

Could you please describe your problem in more detail.  What text does not
show up?  I have a program to make regression plots, so I could just send
it to you.

Barry

On Mon, Jun 13, 2022 at 10:20 PM Mauro Morichetti via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi all,
>
> I have a problem with the script attached. I can not figure out, with the
> text doesn?t show on the scatter plot. I think it is  stupid question, but
> I going crazy, also because I used with other data without problems.
>
> Thanks
>
>
> ---------
>
> Ph.D. Mauro Morichetti
> Department of Industrial Engineering and Mathematical Science
> Marche Polytechnic University, Ancona, Italy
> m.morichetti at pm.univpm.it
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 

Barry H. Lynn, Ph.D
Senior Scientist, Lecturer,
The Institute of Earth Sciences,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

Weather It Is, LTD
Weather and Climate Focus
https://weather-it-is.com <http://weather-it-is.com>
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220614/060af5f1/attachment.htm>

From barry.h.lynn at gmail.com  Tue Jun 14 08:11:33 2022
From: barry.h.lynn at gmail.com (Barry Lynn)
Date: Tue, 14 Jun 2022 17:11:33 +0300
Subject: [ncl-talk] (no subject)
In-Reply-To: <CAN7Bv24ri+Enr3rqZ6xbCM_F9SjgPxQKntRbaXxd_+XkDo_W+g@mail.gmail.com>
References: <83E924BC-E075-41A3-9484-68B6E759D9A3@pm.univpm.it>
 <CAN7Bv24ri+Enr3rqZ6xbCM_F9SjgPxQKntRbaXxd_+XkDo_W+g@mail.gmail.com>
Message-ID: <CAN7Bv25o2U4h_besmQtrUkLDjSOEWZ7i2v2jh3fkQU76deA0cw@mail.gmail.com>

Hi Mauro:

Just letting the group that the problem was solved as indicated.

On Tue, Jun 14, 2022 at 7:23 AM Barry Lynn <barry.h.lynn at gmail.com> wrote:

> Hello Mauro:
>
> Could you please describe your problem in more detail.  What text does not
> show up?  I have a program to make regression plots, so I could just send
> it to you.
>
> Barry
>
> On Mon, Jun 13, 2022 at 10:20 PM Mauro Morichetti via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Hi all,
>>
>> I have a problem with the script attached. I can not figure out, with the
>> text doesn?t show on the scatter plot. I think it is  stupid question, but
>> I going crazy, also because I used with other data without problems.
>>
>> Thanks
>>
>>
>> ---------
>>
>> Ph.D. Mauro Morichetti
>> Department of Industrial Engineering and Mathematical Science
>> Marche Polytechnic University, Ancona, Italy
>> m.morichetti at pm.univpm.it
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>
> --
>
> Barry H. Lynn, Ph.D
> Senior Scientist, Lecturer,
> The Institute of Earth Sciences,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> Weather It Is, LTD
> Weather and Climate Focus
> https://weather-it-is.com <http://weather-it-is.com>
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
>
>

-- 

Barry H. Lynn, Ph.D
Senior Scientist, Lecturer,
The Institute of Earth Sciences,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

Weather It Is, LTD
Weather and Climate Focus
https://weather-it-is.com <http://weather-it-is.com>
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220614/e04d8bcb/attachment.htm>

From s1073755 at pm.univpm.it  Tue Jun 14 08:39:52 2022
From: s1073755 at pm.univpm.it (Mauro Morichetti)
Date: Tue, 14 Jun 2022 16:39:52 +0200
Subject: [ncl-talk] (no subject)
In-Reply-To: <CAN7Bv25S3+sBffsd3ebS5hEMeAEUS0ohKK6kXpjYmU9XM=crCg@mail.gmail.com>
References: <83E924BC-E075-41A3-9484-68B6E759D9A3@pm.univpm.it>
 <CAN7Bv24ri+Enr3rqZ6xbCM_F9SjgPxQKntRbaXxd_+XkDo_W+g@mail.gmail.com>
 <0CF2097B-3B56-4847-85E4-6A7DD255B308@pm.univpm.it>
 <CAN7Bv25S3+sBffsd3ebS5hEMeAEUS0ohKK6kXpjYmU9XM=crCg@mail.gmail.com>
Message-ID: <C9358737-16AA-499D-8B9D-F5BB7F2260C7@pm.univpm.it>

Hi,

yes the error it is exactly the range. 

Thanks, and I am sorry for dumb question.

Regards


---------

Ph.D. Mauro Morichetti
Department of Industrial Engineering and Mathematical Science
Marche Polytechnic University, Ancona, Italy
m.morichetti at pm.univpm.it

> On 14Jun, 2022, at 13:17, Barry Lynn <barry.h.lynn at gmail.com> wrote:
> 
> Mauro:
> 
> My use of this command is:
> 
> text = gsn_add_text(wks,map,label,lon_cen,lat_loc,txres)
> 
> yours:
> gsn_add_text(wks,plot,"N pt. = "+nxy, 100,30,txres)
> 
> The most obvious thing I can think of is that your points (100,30) need to have a range of 0 to 1.0.  
> 
> Can you try this?
> 
> On Tue, Jun 14, 2022 at 2:04 PM Mauro Morichetti <s1073755 at pm.univpm.it <mailto:s1073755 at pm.univpm.it>> wrote:
> Hi all,
> 
> the text on the folloewing lines:
> 
> txid(0) = gsn_add_text(wks,plot,"N pt. = "+nxy, 100,30,txres)
> txid(1) = gsn_add_text(wks,plot,"Y = "+sprintf("%5.2f", rc)+"X + "+sprintf("%5.2f", rc at yintercept), 100,15,txres)
> 
> Thanks
> 
> ---------
> 
> Ph.D. Mauro Morichetti
> Department of Industrial Engineering and Mathematical Science
> Marche Polytechnic University, Ancona, Italy
> m.morichetti at pm.univpm.it <mailto:m.morichetti at pm.univpm.it>
> 
>> On 14Jun, 2022, at 06:23, Barry Lynn <barry.h.lynn at gmail.com <mailto:barry.h.lynn at gmail.com>> wrote:
>> 
>> Hello Mauro:
>> 
>> Could you please describe your problem in more detail.  What text does not show up?  I have a program to make regression plots, so I could just send it to you.
>> 
>> Barry
>> 
>> On Mon, Jun 13, 2022 at 10:20 PM Mauro Morichetti via ncl-talk <ncl-talk at mailman.ucar.edu <mailto:ncl-talk at mailman.ucar.edu>> wrote:
>> Hi all,
>> 
>> I have a problem with the script attached. I can not figure out, with the text doesn?t show on the scatter plot. I think it is  stupid question, but I going crazy, also because I used with other data without problems.
>> 
>> Thanks
>> 
>> 
>> ---------
>> 
>> Ph.D. Mauro Morichetti
>> Department of Industrial Engineering and Mathematical Science
>> Marche Polytechnic University, Ancona, Italy
>> m.morichetti at pm.univpm.it <mailto:m.morichetti at pm.univpm.it>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu <mailto:ncl-talk at mailman.ucar.edu>
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk <https://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>> 
>> 
>> -- 
>> 
>> Barry H. Lynn, Ph.D
>> Senior Scientist, Lecturer,
>> The Institute of Earth Sciences, 
>> The Hebrew University of Jerusalem, 
>> Givat Ram, Jerusalem 91904, Israel 
>> Tel: 972 547 231 170
>> Fax: (972)-25662581
>> 
>> Weather It Is, LTD
>> Weather and Climate Focus
>> https://weather-it-is.com <http://weather-it-is.com/>
>> Jerusalem, Israel
>> Local: 02 930 9525
>> Cell: 054 7 231 170
>> Int-IS: x972 2 930 9525
> 
> 
> 
> -- 
> 
> Barry H. Lynn, Ph.D
> Senior Scientist, Lecturer,
> The Institute of Earth Sciences, 
> The Hebrew University of Jerusalem, 
> Givat Ram, Jerusalem 91904, Israel 
> Tel: 972 547 231 170
> Fax: (972)-25662581
> 
> Weather It Is, LTD
> Weather and Climate Focus
> https://weather-it-is.com <http://weather-it-is.com/>
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220614/f3cadf1d/attachment.htm>

From setareh.rahimi at gmail.com  Fri Jun 17 04:02:33 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Fri, 17 Jun 2022 14:32:33 +0430
Subject: [ncl-talk] Plotting NDVI from (AVHRR) sensor
Message-ID: <CAN+DFgwx3kZXP7A-MdUqEXfv7bJxtoOM_WN3P8ZHZULhkVu7SQ@mail.gmail.com>

Dear all,
I would like to use this example "
https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl",  available at:
https://www.ncl.ucar.edu/Applications/ndvi.shtml
to plot NDVI from Advanced Very High-Resolution Radiometer (AVHRR) sensor.
But, I got the following error:

fatal:error in statement
fatal:syntax error: line 55 in file ndvi before or near else
  else
-----^

fatal:error in statement
fatal:syntax error: line 57 in file ndvi before or near if
  end if
-------^
, "IF" statement is not complete. I do not know what condition has been
considered for the "IF" statement. Would you please kindly advise me?
Thanks in advance,
Best wishes,



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

From shea at ucar.edu  Fri Jun 17 12:19:27 2022
From: shea at ucar.edu (Dennis Shea)
Date: Fri, 17 Jun 2022 12:19:27 -0600
Subject: [ncl-talk] Plotting NDVI from (AVHRR) sensor
In-Reply-To: <CAN+DFgwx3kZXP7A-MdUqEXfv7bJxtoOM_WN3P8ZHZULhkVu7SQ@mail.gmail.com>
References: <CAN+DFgwx3kZXP7A-MdUqEXfv7bJxtoOM_WN3P8ZHZULhkVu7SQ@mail.gmail.com>
Message-ID: <CAOF1d_6rpKDAshBpUv8o9sewXfQ6Q51UH4c12Zznn=y2268RQg@mail.gmail.com>

Looks like I 'somehow'  deleted a line
====
Note the script has the following at the top:

TEST = False
=====

The published script should be:

  if (TEST) then   ; <=======  This was inadvertently deleted by
'someone'    :;-)
                                   ;
      ndvi  = short2flt(f->NDVI(:,::50,::25))         ; original unpacked values
  else
      ndvi  = short2flt(f->NDVI)                      ; original unpacked values
  end if



On Fri, Jun 17, 2022 at 4:03 AM Setareh Rahimi via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear all,
> I would like to use this example "
> https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl",  available at:
> https://www.ncl.ucar.edu/Applications/ndvi.shtml
> to plot NDVI from Advanced Very High-Resolution Radiometer (AVHRR) sensor.
> But, I got the following error:
>
> fatal:error in statement
> fatal:syntax error: line 55 in file ndvi before or near else
>   else
> -----^
>
> fatal:error in statement
> fatal:syntax error: line 57 in file ndvi before or near if
>   end if
> -------^
> , "IF" statement is not complete. I do not know what condition has been
> considered for the "IF" statement. Would you please kindly advise me?
> Thanks in advance,
> Best wishes,
>
>
>
> --
> S.Rahimi
>
> _______________________________________________
> 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/20220617/76c2aa10/attachment.htm>

From setareh.rahimi at gmail.com  Sat Jun 18 04:49:56 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Sat, 18 Jun 2022 15:19:56 +0430
Subject: [ncl-talk] Plotting NDVI from (AVHRR) sensor
In-Reply-To: <CAOF1d_6rpKDAshBpUv8o9sewXfQ6Q51UH4c12Zznn=y2268RQg@mail.gmail.com>
References: <CAN+DFgwx3kZXP7A-MdUqEXfv7bJxtoOM_WN3P8ZHZULhkVu7SQ@mail.gmail.com>
 <CAOF1d_6rpKDAshBpUv8o9sewXfQ6Q51UH4c12Zznn=y2268RQg@mail.gmail.com>
Message-ID: <CAN+DFgyLiKkcOeSz1CesyQuf4kXSUkYt+hyUizS3J4pgp2jjTA@mail.gmail.com>

Dear Dennis,
So many thanks for your help, but after running the script it said: "f is
undefined" !!
Why does this happen, please?
Best wishes,


On Fri, Jun 17, 2022 at 10:49 PM Dennis Shea <shea at ucar.edu> wrote:

> Looks like I 'somehow'  deleted a line
> ====
> Note the script has the following at the top:
>
> TEST = False
> =====
>
> The published script should be:
>
>   if (TEST) then   ; <=======  This was inadvertently deleted by 'someone'    :;-)
>                                    ;
>       ndvi  = short2flt(f->NDVI(:,::50,::25))         ; original unpacked values
>   else
>       ndvi  = short2flt(f->NDVI)                      ; original unpacked values
>   end if
>
>
>
> On Fri, Jun 17, 2022 at 4:03 AM Setareh Rahimi via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Dear all,
>> I would like to use this example "
>> https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl",  available
>> at: https://www.ncl.ucar.edu/Applications/ndvi.shtml
>> to plot NDVI from Advanced Very High-Resolution Radiometer (AVHRR)
>> sensor. But, I got the following error:
>>
>> fatal:error in statement
>> fatal:syntax error: line 55 in file ndvi before or near else
>>   else
>> -----^
>>
>> fatal:error in statement
>> fatal:syntax error: line 57 in file ndvi before or near if
>>   end if
>> -------^
>> , "IF" statement is not complete. I do not know what condition has been
>> considered for the "IF" statement. Would you please kindly advise me?
>> Thanks in advance,
>> Best wishes,
>>
>>
>>
>> --
>> S.Rahimi
>>
>> _______________________________________________
>> 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/20220618/44d94ab9/attachment.htm>

From setareh.rahimi at gmail.com  Sat Jun 18 05:01:09 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Sat, 18 Jun 2022 15:31:09 +0430
Subject: [ncl-talk] NDVI from Advanced Very High-Resolution Radiometer
 (AVHRR)
Message-ID: <CAN+DFgyqB1HBkNCNxcGQDKMMrPna8z+SgQawFVaegCyYnnhCyw@mail.gmail.com>

Dear all,

I need to calculate the monthly anomaly of the Normalized Difference
Vegetation Index (NDVI), from Advanced Very High-Resolution Radiometer
(AVHRR) sensor. I have downloaded AVHRR data (1986-2021) from the link
below:
https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
I first tried to prepare a script to combine all the data, and then
calculate the monthly average from daily data ( the link provides daily
data). After having monthly data I tried to figure out anomalies (the
script has been attached). But NCL stopped and returned nothing to me. I
just wonder what is wrong. Please kindly advise me on this issue.
Best wishes,


-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220618/c1228dd7/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NDVI anomaly
Type: application/octet-stream
Size: 778 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220618/c1228dd7/attachment.obj>

From shea at ucar.edu  Sat Jun 18 09:37:20 2022
From: shea at ucar.edu (Dennis Shea)
Date: Sat, 18 Jun 2022 09:37:20 -0600
Subject: [ncl-talk] Plotting NDVI from (AVHRR) sensor
In-Reply-To: <CAN+DFgyLiKkcOeSz1CesyQuf4kXSUkYt+hyUizS3J4pgp2jjTA@mail.gmail.com>
References: <CAN+DFgwx3kZXP7A-MdUqEXfv7bJxtoOM_WN3P8ZHZULhkVu7SQ@mail.gmail.com>
 <CAOF1d_6rpKDAshBpUv8o9sewXfQ6Q51UH4c12Zznn=y2268RQg@mail.gmail.com>
 <CAN+DFgyLiKkcOeSz1CesyQuf4kXSUkYt+hyUizS3J4pgp2jjTA@mail.gmail.com>
Message-ID: <CAOF1d_5VWGzCvXeEcsJ_9cifYWRL6G8pJ7hfPUW2aXRqVmxX5Q@mail.gmail.com>

re:   "f is undefined"

??  This indicates that the file pointer is not correct

The script loops over one or more files:

do nf=0,nfili-1
  pthi = diri+fili(nf)
  f    = addfile(pthi, "r")

f is a pointer to a file.
======================================

You should go into 'debug mode'

print(diri)               ; data path
print("=====")
print(fili)               ; print names of one or more files
print("=====")
print("nfili="+nfili)     ; number of files
print("=====")

:
:
do nf=0,nfili-1            ; loop over each file
  pthi = diri+fili(nf)

  print(pthi)
  print("=====")

  f    = addfile(pthi, "r")

  print(f)     ; equivalent to 'ncdump -h '


On Sat, Jun 18, 2022 at 4:50 AM Setareh Rahimi <setareh.rahimi at gmail.com>
wrote:

> Dear Dennis,
> So many thanks for your help, but after running the script it said: "f is
> undefined" !!
> Why does this happen, please?
> Best wishes,
>
>
> On Fri, Jun 17, 2022 at 10:49 PM Dennis Shea <shea at ucar.edu> wrote:
>
>> Looks like I 'somehow'  deleted a line
>> ====
>> Note the script has the following at the top:
>>
>> TEST = False
>> =====
>>
>> The published script should be:
>>
>>   if (TEST) then   ; <=======  This was inadvertently deleted by 'someone'    :;-)
>>                                    ;
>>       ndvi  = short2flt(f->NDVI(:,::50,::25))         ; original unpacked values
>>   else
>>       ndvi  = short2flt(f->NDVI)                      ; original unpacked values
>>   end if
>>
>>
>>
>> On Fri, Jun 17, 2022 at 4:03 AM Setareh Rahimi via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>> Dear all,
>>> I would like to use this example "
>>> https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl",  available
>>> at: https://www.ncl.ucar.edu/Applications/ndvi.shtml
>>> to plot NDVI from Advanced Very High-Resolution Radiometer (AVHRR)
>>> sensor. But, I got the following error:
>>>
>>> fatal:error in statement
>>> fatal:syntax error: line 55 in file ndvi before or near else
>>>   else
>>> -----^
>>>
>>> fatal:error in statement
>>> fatal:syntax error: line 57 in file ndvi before or near if
>>>   end if
>>> -------^
>>> , "IF" statement is not complete. I do not know what condition has been
>>> considered for the "IF" statement. Would you please kindly advise me?
>>> Thanks in advance,
>>> Best wishes,
>>>
>>>
>>>
>>> --
>>> S.Rahimi
>>>
>>> _______________________________________________
>>> 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/20220618/6adb8009/attachment.htm>

From shea at ucar.edu  Sat Jun 18 12:21:52 2022
From: shea at ucar.edu (Dennis Shea)
Date: Sat, 18 Jun 2022 12:21:52 -0600
Subject: [ncl-talk] NDVI from Advanced Very High-Resolution Radiometer
 (AVHRR)
In-Reply-To: <CAN+DFgyqB1HBkNCNxcGQDKMMrPna8z+SgQawFVaegCyYnnhCyw@mail.gmail.com>
References: <CAN+DFgyqB1HBkNCNxcGQDKMMrPna8z+SgQawFVaegCyYnnhCyw@mail.gmail.com>
Message-ID: <CAOF1d_7ZPXCrW1Jzj+BeL5QmBe1ohQWY4L66U24hKkgbXNs3fg@mail.gmail.com>

Based on the provided information, I have no idea why you would get
"returned nothing to me"

I know that NDVI files are 0.1 degree resolution (3600,7200). You have 36
years (1986-2021) of daily data.
So, 36*365=13140   [+ a few Leap Year days)

   prc(13140,3600,7200)    => daily data

A lot of memory is needed. Maybe your computer can not handle this?

After the
      xM = calculate_monthly_values (prc, "avg", 0,opt)
      printVarSummary(xM)

   xM(432,3600,7200)



On Sat, Jun 18, 2022 at 5:01 AM Setareh Rahimi via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear all,
>
> I need to calculate the monthly anomaly of the Normalized Difference
> Vegetation Index (NDVI), from Advanced Very High-Resolution Radiometer
> (AVHRR) sensor. I have downloaded AVHRR data (1986-2021) from the link
> below:
>
> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
> I first tried to prepare a script to combine all the data, and then
> calculate the monthly average from daily data ( the link provides daily
> data). After having monthly data I tried to figure out anomalies (the
> script has been attached). But NCL stopped and returned nothing to me. I
> just wonder what is wrong. Please kindly advise me on this issue.
> Best wishes,
>
>
> --
> S.Rahimi
>
> _______________________________________________
> 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/20220618/5356a81b/attachment.htm>

From setareh.rahimi at gmail.com  Tue Jun 21 05:41:46 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Tue, 21 Jun 2022 16:11:46 +0430
Subject: [ncl-talk] Plotting NDVI from (AVHRR) sensor
In-Reply-To: <CAOF1d_5VWGzCvXeEcsJ_9cifYWRL6G8pJ7hfPUW2aXRqVmxX5Q@mail.gmail.com>
References: <CAN+DFgwx3kZXP7A-MdUqEXfv7bJxtoOM_WN3P8ZHZULhkVu7SQ@mail.gmail.com>
 <CAOF1d_6rpKDAshBpUv8o9sewXfQ6Q51UH4c12Zznn=y2268RQg@mail.gmail.com>
 <CAN+DFgyLiKkcOeSz1CesyQuf4kXSUkYt+hyUizS3J4pgp2jjTA@mail.gmail.com>
 <CAOF1d_5VWGzCvXeEcsJ_9cifYWRL6G8pJ7hfPUW2aXRqVmxX5Q@mail.gmail.com>
Message-ID: <CAN+DFgwJbtBoEF4EhLyqBvbjYoEmGDsvQvUaU5p-brVYcJnUFA@mail.gmail.com>

Dear Dennis,
So many thanks for your attention. I did as you suggested and I attached
the result. Please check this out.
Best wishes,

On Sat, Jun 18, 2022 at 8:07 PM Dennis Shea <shea at ucar.edu> wrote:

> re:   "f is undefined"
>
> ??  This indicates that the file pointer is not correct
>
> The script loops over one or more files:
>
> do nf=0,nfili-1
>   pthi = diri+fili(nf)
>   f    = addfile(pthi, "r")
>
> f is a pointer to a file.
> ======================================
>
> You should go into 'debug mode'
>
> print(diri)               ; data path
> print("=====")
> print(fili)               ; print names of one or more files
> print("=====")
> print("nfili="+nfili)     ; number of files
> print("=====")
>
> :
> :
> do nf=0,nfili-1            ; loop over each file
>   pthi = diri+fili(nf)
>
>   print(pthi)
>   print("=====")
>
>   f    = addfile(pthi, "r")
>
>   print(f)     ; equivalent to 'ncdump -h '
>
>
> On Sat, Jun 18, 2022 at 4:50 AM Setareh Rahimi <setareh.rahimi at gmail.com>
> wrote:
>
>> Dear Dennis,
>> So many thanks for your help, but after running the script it said: "f
>> is undefined" !!
>> Why does this happen, please?
>> Best wishes,
>>
>>
>> On Fri, Jun 17, 2022 at 10:49 PM Dennis Shea <shea at ucar.edu> wrote:
>>
>>> Looks like I 'somehow'  deleted a line
>>> ====
>>> Note the script has the following at the top:
>>>
>>> TEST = False
>>> =====
>>>
>>> The published script should be:
>>>
>>>   if (TEST) then   ; <=======  This was inadvertently deleted by 'someone'    :;-)
>>>                                    ;
>>>       ndvi  = short2flt(f->NDVI(:,::50,::25))         ; original unpacked values
>>>   else
>>>       ndvi  = short2flt(f->NDVI)                      ; original unpacked values
>>>   end if
>>>
>>>
>>>
>>> On Fri, Jun 17, 2022 at 4:03 AM Setareh Rahimi via ncl-talk <
>>> ncl-talk at mailman.ucar.edu> wrote:
>>>
>>>> Dear all,
>>>> I would like to use this example "
>>>> https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl",  available
>>>> at: https://www.ncl.ucar.edu/Applications/ndvi.shtml
>>>> to plot NDVI from Advanced Very High-Resolution Radiometer (AVHRR)
>>>> sensor. But, I got the following error:
>>>>
>>>> fatal:error in statement
>>>> fatal:syntax error: line 55 in file ndvi before or near else
>>>>   else
>>>> -----^
>>>>
>>>> fatal:error in statement
>>>> fatal:syntax error: line 57 in file ndvi before or near if
>>>>   end if
>>>> -------^
>>>> , "IF" statement is not complete. I do not know what condition has been
>>>> considered for the "IF" statement. Would you please kindly advise me?
>>>> Thanks in advance,
>>>> Best wishes,
>>>>
>>>>
>>>>
>>>> --
>>>> S.Rahimi
>>>>
>>>> _______________________________________________
>>>> 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
>>
>>

-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220621/3d76391f/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ss.png
Type: image/png
Size: 280810 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220621/3d76391f/attachment-0001.png>

From setareh.rahimi at gmail.com  Tue Jun 21 05:48:52 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Tue, 21 Jun 2022 16:18:52 +0430
Subject: [ncl-talk] NDVI from Advanced Very High-Resolution Radiometer
 (AVHRR)
In-Reply-To: <CAOF1d_7ZPXCrW1Jzj+BeL5QmBe1ohQWY4L66U24hKkgbXNs3fg@mail.gmail.com>
References: <CAN+DFgyqB1HBkNCNxcGQDKMMrPna8z+SgQawFVaegCyYnnhCyw@mail.gmail.com>
 <CAOF1d_7ZPXCrW1Jzj+BeL5QmBe1ohQWY4L66U24hKkgbXNs3fg@mail.gmail.com>
Message-ID: <CAN+DFgwza-k80POTJSYnDMUD67Yr5oR13KPot7-M55LK729Mpw@mail.gmail.com>

Dear Dennis,
So many thanks for your advice. I run the script on a remote server with at
least *1 terabyte* of memory. I think this is not the case.
I used the below commands to check the files:
===========================================
fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
   print(fils)

   f   = addfiles (fils, "r")

   x   = short2flt(f[:]->NDVI)

   printVarSummary(x)
===========================================
, and attached the result. Please look at the image, and kindly help me
sort the issue out. I do not know what is wrong with the files I downloaded.
best regards,

On Sat, Jun 18, 2022 at 10:52 PM Dennis Shea <shea at ucar.edu> wrote:

> Based on the provided information, I have no idea why you would get
> "returned nothing to me"
>
> I know that NDVI files are 0.1 degree resolution (3600,7200). You have 36
> years (1986-2021) of daily data.
> So, 36*365=13140   [+ a few Leap Year days)
>
>    prc(13140,3600,7200)    => daily data
>
> A lot of memory is needed. Maybe your computer can not handle this?
>
> After the
>       xM = calculate_monthly_values (prc, "avg", 0,opt)
>       printVarSummary(xM)
>
>    xM(432,3600,7200)
>
>
>
> On Sat, Jun 18, 2022 at 5:01 AM Setareh Rahimi via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Dear all,
>>
>> I need to calculate the monthly anomaly of the Normalized Difference
>> Vegetation Index (NDVI), from Advanced Very High-Resolution Radiometer
>> (AVHRR) sensor. I have downloaded AVHRR data (1986-2021) from the link
>> below:
>>
>> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
>> I first tried to prepare a script to combine all the data, and then
>> calculate the monthly average from daily data ( the link provides daily
>> data). After having monthly data I tried to figure out anomalies (the
>> script has been attached). But NCL stopped and returned nothing to me. I
>> just wonder what is wrong. Please kindly advise me on this issue.
>> Best wishes,
>>
>>
>> --
>> S.Rahimi
>>
>> _______________________________________________
>> 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/20220621/79a31d5b/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 1401-03-31 at 3.47.26 PM.png
Type: image/png
Size: 387403 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220621/79a31d5b/attachment-0001.png>

From shea at ucar.edu  Tue Jun 21 11:39:31 2022
From: shea at ucar.edu (Dennis Shea)
Date: Tue, 21 Jun 2022 11:39:31 -0600
Subject: [ncl-talk] NDVI from Advanced Very High-Resolution Radiometer
 (AVHRR)
In-Reply-To: <CAN+DFgwza-k80POTJSYnDMUD67Yr5oR13KPot7-M55LK729Mpw@mail.gmail.com>
References: <CAN+DFgyqB1HBkNCNxcGQDKMMrPna8z+SgQawFVaegCyYnnhCyw@mail.gmail.com>
 <CAOF1d_7ZPXCrW1Jzj+BeL5QmBe1ohQWY4L66U24hKkgbXNs3fg@mail.gmail.com>
 <CAN+DFgwza-k80POTJSYnDMUD67Yr5oR13KPot7-M55LK729Mpw@mail.gmail.com>
Message-ID: <CAOF1d_7TKFbW1OYYX8efBsqCcCShYD+yj3Tuzk7ivO72n0S8pQ@mail.gmail.com>

I downloaded 5 files for a test.

%> ncl NDVI_anomaly.ncl

Copyright (C) 1995-2017 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 6.5.0
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.


Variable: fils
Type: string
Total Size: 40 bytes
            5 values
Number of Dimensions: 1
Dimensions and sizes:   [5]
Coordinates:
(0)     AVHRR-Land_v005_AVH13C1_NOAA-09_19860101_c20170612092103.nc
(1)     AVHRR-Land_v005_AVH13C1_NOAA-09_19860102_c20170612093403.nc
(2)     AVHRR-Land_v005_AVH13C1_NOAA-09_19861229_c20170612212717.nc
(3)     AVHRR-Land_v005_AVH13C1_NOAA-09_19861230_c20170612214313.nc
(4)     AVHRR-Land_v005_AVH13C1_NOAA-09_19861231_c20170612215936.nc

Variable: prc
Type: float
Total Size: 518400000 bytes
            129600000 values
Number of Dimensions: 3
Dimensions and sizes:   [time | 5] x [latitude | 3600] x [longitude | 7200]
Coordinates:
            time: [1826..2190]
            latitude: [89.975..-89.97499]
            longitude: [-179.975..179.975]
Number Of Attributes: 7
  long_name :   NOAA Climate Data Record of Normalized Difference
Vegetation Index
  units :       1
  grid_mapping :        crs
  standard_name :       normalized_difference_vegetation_index
  _FillValue_original : -9999
  _FillValue :  -9999
========================================================================
========================================================================
========================================================================
It seems to me that you have path or file problems.
I am not sure what else I can do.

Others?



On Tue, Jun 21, 2022 at 5:49 AM Setareh Rahimi <setareh.rahimi at gmail.com>
wrote:

> Dear Dennis,
> So many thanks for your advice. I run the script on a remote server with
> at least *1 terabyte* of memory. I think this is not the case.
> I used the below commands to check the files:
> ===========================================
> fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
>    print(fils)
>
>    f   = addfiles (fils, "r")
>
>    x   = short2flt(f[:]->NDVI)
>
>    printVarSummary(x)
> ===========================================
> , and attached the result. Please look at the image, and kindly help me
> sort the issue out. I do not know what is wrong with the files I downloaded.
> best regards,
>
> On Sat, Jun 18, 2022 at 10:52 PM Dennis Shea <shea at ucar.edu> wrote:
>
>> Based on the provided information, I have no idea why you would get
>> "returned nothing to me"
>>
>> I know that NDVI files are 0.1 degree resolution (3600,7200). You have 36
>> years (1986-2021) of daily data.
>> So, 36*365=13140   [+ a few Leap Year days)
>>
>>    prc(13140,3600,7200)    => daily data
>>
>> A lot of memory is needed. Maybe your computer can not handle this?
>>
>> After the
>>       xM = calculate_monthly_values (prc, "avg", 0,opt)
>>       printVarSummary(xM)
>>
>>    xM(432,3600,7200)
>>
>>
>>
>> On Sat, Jun 18, 2022 at 5:01 AM Setareh Rahimi via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>> Dear all,
>>>
>>> I need to calculate the monthly anomaly of the Normalized Difference
>>> Vegetation Index (NDVI), from Advanced Very High-Resolution Radiometer
>>> (AVHRR) sensor. I have downloaded AVHRR data (1986-2021) from the link
>>> below:
>>>
>>> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
>>> I first tried to prepare a script to combine all the data, and then
>>> calculate the monthly average from daily data ( the link provides daily
>>> data). After having monthly data I tried to figure out anomalies (the
>>> script has been attached). But NCL stopped and returned nothing to me. I
>>> just wonder what is wrong. Please kindly advise me on this issue.
>>> Best wishes,
>>>
>>>
>>> --
>>> S.Rahimi
>>>
>>> _______________________________________________
>>> 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/20220621/97d58106/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NDVI_anomaly.ncl
Type: application/octet-stream
Size: 811 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220621/97d58106/attachment.obj>

From coquart at cerfacs.fr  Wed Jun 22 02:58:41 2022
From: coquart at cerfacs.fr (Laure Coquart)
Date: Wed, 22 Jun 2022 10:58:41 +0200 (CEST)
Subject: [ncl-talk] NCL using ESMF for remapping
Message-ID: <748246265.61488514.1655888321217.JavaMail.zimbra@cerfacs.fr>

Hi, 

I developed a jupyter notebook to perform remapping using xesmf. 
I have a reference NCL program to see if what I did is correct, but I have differences of order(50) between the remapped file obtained with the NCL program using ESMF patch method the file obtained by my jupyter notebook using xesmf patch method. 

Did someone ever tried to use xesmf to remap data ? 
Thanks for your help, 
Best regards, 
Laure 
----------- 
L. Coquart: IE au [ http://intranet.cerfacs.fr/intranet_equipes/globc/ | CECI UMR 5318 ] ; [ http://cerfacs.fr/~coquart/ | Page perso ] 
Not working Tuesday, Thursday, Friday 


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220622/79e2ee10/attachment.htm>

From setareh.rahimi at gmail.com  Wed Jun 22 03:07:04 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Wed, 22 Jun 2022 13:37:04 +0430
Subject: [ncl-talk] NDVI from Advanced Very High-Resolution Radiometer
 (AVHRR)
In-Reply-To: <CAOF1d_7TKFbW1OYYX8efBsqCcCShYD+yj3Tuzk7ivO72n0S8pQ@mail.gmail.com>
References: <CAN+DFgyqB1HBkNCNxcGQDKMMrPna8z+SgQawFVaegCyYnnhCyw@mail.gmail.com>
 <CAOF1d_7ZPXCrW1Jzj+BeL5QmBe1ohQWY4L66U24hKkgbXNs3fg@mail.gmail.com>
 <CAN+DFgwza-k80POTJSYnDMUD67Yr5oR13KPot7-M55LK729Mpw@mail.gmail.com>
 <CAOF1d_7TKFbW1OYYX8efBsqCcCShYD+yj3Tuzk7ivO72n0S8pQ@mail.gmail.com>
Message-ID: <CAN+DFgzd9g=EKzdAE-rbUo_7uFdhi_1AK+SnuohFp7o1wzGVjw@mail.gmail.com>

Dear Dennis,
I profoundly appreciate your kind attention and the time you considered
downloading the files and checking them. Now I found out what is the
problem. so many many thanks!
Best wishes,



On Tue, Jun 21, 2022 at 10:09 PM Dennis Shea <shea at ucar.edu> wrote:

> I downloaded 5 files for a test.
>
> %> ncl NDVI_anomaly.ncl
>
> Copyright (C) 1995-2017 - All Rights Reserved
>  University Corporation for Atmospheric Research
>  NCAR Command Language Version 6.5.0
>  The use of this software is governed by a License Agreement.
>  See http://www.ncl.ucar.edu/ for more details.
>
>
> Variable: fils
> Type: string
> Total Size: 40 bytes
>             5 values
> Number of Dimensions: 1
> Dimensions and sizes:   [5]
> Coordinates:
> (0)     AVHRR-Land_v005_AVH13C1_NOAA-09_19860101_c20170612092103.nc
> (1)     AVHRR-Land_v005_AVH13C1_NOAA-09_19860102_c20170612093403.nc
> (2)     AVHRR-Land_v005_AVH13C1_NOAA-09_19861229_c20170612212717.nc
> (3)     AVHRR-Land_v005_AVH13C1_NOAA-09_19861230_c20170612214313.nc
> (4)     AVHRR-Land_v005_AVH13C1_NOAA-09_19861231_c20170612215936.nc
>
> Variable: prc
> Type: float
> Total Size: 518400000 bytes
>             129600000 values
> Number of Dimensions: 3
> Dimensions and sizes:   [time | 5] x [latitude | 3600] x [longitude | 7200]
> Coordinates:
>             time: [1826..2190]
>             latitude: [89.975..-89.97499]
>             longitude: [-179.975..179.975]
> Number Of Attributes: 7
>   long_name :   NOAA Climate Data Record of Normalized Difference
> Vegetation Index
>   units :       1
>   grid_mapping :        crs
>   standard_name :       normalized_difference_vegetation_index
>   _FillValue_original : -9999
>   _FillValue :  -9999
> ========================================================================
> ========================================================================
> ========================================================================
> It seems to me that you have path or file problems.
> I am not sure what else I can do.
>
> Others?
>
>
>
> On Tue, Jun 21, 2022 at 5:49 AM Setareh Rahimi <setareh.rahimi at gmail.com>
> wrote:
>
>> Dear Dennis,
>> So many thanks for your advice. I run the script on a remote server with
>> at least *1 terabyte* of memory. I think this is not the case.
>> I used the below commands to check the files:
>> ===========================================
>> fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
>>    print(fils)
>>
>>    f   = addfiles (fils, "r")
>>
>>    x   = short2flt(f[:]->NDVI)
>>
>>    printVarSummary(x)
>> ===========================================
>> , and attached the result. Please look at the image, and kindly help me
>> sort the issue out. I do not know what is wrong with the files I downloaded.
>> best regards,
>>
>> On Sat, Jun 18, 2022 at 10:52 PM Dennis Shea <shea at ucar.edu> wrote:
>>
>>> Based on the provided information, I have no idea why you would get
>>> "returned nothing to me"
>>>
>>> I know that NDVI files are 0.1 degree resolution (3600,7200). You have
>>> 36 years (1986-2021) of daily data.
>>> So, 36*365=13140   [+ a few Leap Year days)
>>>
>>>    prc(13140,3600,7200)    => daily data
>>>
>>> A lot of memory is needed. Maybe your computer can not handle this?
>>>
>>> After the
>>>       xM = calculate_monthly_values (prc, "avg", 0,opt)
>>>       printVarSummary(xM)
>>>
>>>    xM(432,3600,7200)
>>>
>>>
>>>
>>> On Sat, Jun 18, 2022 at 5:01 AM Setareh Rahimi via ncl-talk <
>>> ncl-talk at mailman.ucar.edu> wrote:
>>>
>>>> Dear all,
>>>>
>>>> I need to calculate the monthly anomaly of the Normalized Difference
>>>> Vegetation Index (NDVI), from Advanced Very High-Resolution Radiometer
>>>> (AVHRR) sensor. I have downloaded AVHRR data (1986-2021) from the link
>>>> below:
>>>>
>>>> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
>>>> I first tried to prepare a script to combine all the data, and then
>>>> calculate the monthly average from daily data ( the link provides daily
>>>> data). After having monthly data I tried to figure out anomalies (the
>>>> script has been attached). But NCL stopped and returned nothing to me. I
>>>> just wonder what is wrong. Please kindly advise me on this issue.
>>>> Best wishes,
>>>>
>>>>
>>>> --
>>>> S.Rahimi
>>>>
>>>> _______________________________________________
>>>> 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
>>
>>

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

From setareh.rahimi at gmail.com  Sun Jun 26 10:12:34 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Sun, 26 Jun 2022 20:42:34 +0430
Subject: [ncl-talk] NDVI from Advanced Very High-Resolution Radiometer
 (AVHRR)
In-Reply-To: <CAOF1d_7TKFbW1OYYX8efBsqCcCShYD+yj3Tuzk7ivO72n0S8pQ@mail.gmail.com>
References: <CAN+DFgyqB1HBkNCNxcGQDKMMrPna8z+SgQawFVaegCyYnnhCyw@mail.gmail.com>
 <CAOF1d_7ZPXCrW1Jzj+BeL5QmBe1ohQWY4L66U24hKkgbXNs3fg@mail.gmail.com>
 <CAN+DFgwza-k80POTJSYnDMUD67Yr5oR13KPot7-M55LK729Mpw@mail.gmail.com>
 <CAOF1d_7TKFbW1OYYX8efBsqCcCShYD+yj3Tuzk7ivO72n0S8pQ@mail.gmail.com>
Message-ID: <CAN+DFgyLXT7cBi+ARCYZMeSZhRLYtJ_3x8mZ4qSnAxB-nnsZ4w@mail.gmail.com>

Dear Dennis,
I put all NDVI files in one folder(like the attached image), and tried to
run the attached script to plot the NDVI anomaly from 1986 to 2021.
However, I faced another issue(attached). May I ask you, please to advise
me why this message appeared? How can I remove this issue, please?
Many thanks for considering my request.
Best wishes,

On Tue, Jun 21, 2022 at 10:09 PM Dennis Shea <shea at ucar.edu> wrote:

> I downloaded 5 files for a test.
>
> %> ncl NDVI_anomaly.ncl
>
> Copyright (C) 1995-2017 - All Rights Reserved
>  University Corporation for Atmospheric Research
>  NCAR Command Language Version 6.5.0
>  The use of this software is governed by a License Agreement.
>  See http://www.ncl.ucar.edu/ for more details.
>
>
> Variable: fils
> Type: string
> Total Size: 40 bytes
>             5 values
> Number of Dimensions: 1
> Dimensions and sizes:   [5]
> Coordinates:
> (0)     AVHRR-Land_v005_AVH13C1_NOAA-09_19860101_c20170612092103.nc
> (1)     AVHRR-Land_v005_AVH13C1_NOAA-09_19860102_c20170612093403.nc
> (2)     AVHRR-Land_v005_AVH13C1_NOAA-09_19861229_c20170612212717.nc
> (3)     AVHRR-Land_v005_AVH13C1_NOAA-09_19861230_c20170612214313.nc
> (4)     AVHRR-Land_v005_AVH13C1_NOAA-09_19861231_c20170612215936.nc
>
> Variable: prc
> Type: float
> Total Size: 518400000 bytes
>             129600000 values
> Number of Dimensions: 3
> Dimensions and sizes:   [time | 5] x [latitude | 3600] x [longitude | 7200]
> Coordinates:
>             time: [1826..2190]
>             latitude: [89.975..-89.97499]
>             longitude: [-179.975..179.975]
> Number Of Attributes: 7
>   long_name :   NOAA Climate Data Record of Normalized Difference
> Vegetation Index
>   units :       1
>   grid_mapping :        crs
>   standard_name :       normalized_difference_vegetation_index
>   _FillValue_original : -9999
>   _FillValue :  -9999
> ========================================================================
> ========================================================================
> ========================================================================
> It seems to me that you have path or file problems.
> I am not sure what else I can do.
>
> Others?
>
>
>
> On Tue, Jun 21, 2022 at 5:49 AM Setareh Rahimi <setareh.rahimi at gmail.com>
> wrote:
>
>> Dear Dennis,
>> So many thanks for your advice. I run the script on a remote server with
>> at least *1 terabyte* of memory. I think this is not the case.
>> I used the below commands to check the files:
>> ===========================================
>> fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
>>    print(fils)
>>
>>    f   = addfiles (fils, "r")
>>
>>    x   = short2flt(f[:]->NDVI)
>>
>>    printVarSummary(x)
>> ===========================================
>> , and attached the result. Please look at the image, and kindly help me
>> sort the issue out. I do not know what is wrong with the files I downloaded.
>> best regards,
>>
>> On Sat, Jun 18, 2022 at 10:52 PM Dennis Shea <shea at ucar.edu> wrote:
>>
>>> Based on the provided information, I have no idea why you would get
>>> "returned nothing to me"
>>>
>>> I know that NDVI files are 0.1 degree resolution (3600,7200). You have
>>> 36 years (1986-2021) of daily data.
>>> So, 36*365=13140   [+ a few Leap Year days)
>>>
>>>    prc(13140,3600,7200)    => daily data
>>>
>>> A lot of memory is needed. Maybe your computer can not handle this?
>>>
>>> After the
>>>       xM = calculate_monthly_values (prc, "avg", 0,opt)
>>>       printVarSummary(xM)
>>>
>>>    xM(432,3600,7200)
>>>
>>>
>>>
>>> On Sat, Jun 18, 2022 at 5:01 AM Setareh Rahimi via ncl-talk <
>>> ncl-talk at mailman.ucar.edu> wrote:
>>>
>>>> Dear all,
>>>>
>>>> I need to calculate the monthly anomaly of the Normalized Difference
>>>> Vegetation Index (NDVI), from Advanced Very High-Resolution Radiometer
>>>> (AVHRR) sensor. I have downloaded AVHRR data (1986-2021) from the link
>>>> below:
>>>>
>>>> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
>>>> I first tried to prepare a script to combine all the data, and then
>>>> calculate the monthly average from daily data ( the link provides daily
>>>> data). After having monthly data I tried to figure out anomalies (the
>>>> script has been attached). But NCL stopped and returned nothing to me. I
>>>> just wonder what is wrong. Please kindly advise me on this issue.
>>>> Best wishes,
>>>>
>>>>
>>>> --
>>>> S.Rahimi
>>>>
>>>> _______________________________________________
>>>> 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
>>
>>

-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220626/10b08da1/attachment-0001.htm>
-------------- next part --------------
begin
;************************************************
; Read the file
;************************************************
  ymStrt = 198601                              ; climatology start year/month
  ymLast = 202112                              ;      "      last  year/month  a      = addfile("air.mon.mean.nc","r")
  
;--------------------------------------
  fils = systemfunc("ls *.nc")
   print(fils)

   a   = addfiles (fils, "r")
;-----------------------------------------
  time   = a[:]->time
  yyyymm = cd_calendar(time, -1)
  ntStrt = ind(yyyymm.eq.ymStrt)                ; start time index
  ntLast = ind(yyyymm.eq.ymLast)                ; last  time index
  ndvi    =short2flt(a[:]->NDVI(ntStrt:ntLast,:,:))            ; (time,lat,lon)
  ntim = dimsizes(time)
 printMinMax(ndvi,0)
  months = (/"January", "February", "March", "April" \
            ,"May", "June", "July", "August"         \
            ,"September", "October", "November"      \
            ,"December" /)

;************************************************
; Compute the climatology using a function in contributed.ncl
;************************************************  
  ;prcClm   = clmMonTLL(tmp)     
  ;printVarSummary(prcClm)
  xM = calculate_monthly_values (ndvi, "avg", 0,opt)
  printVarSummary(xM)
  anom = calcMonAnomTLL(ndvi,xM)
;************************************************
; create colors
;*************************************************
  wks = gsn_open_wks("png","NDVI-anomaly_ghcn")             ; send graphics to PNG file

;************************************************
; create panel plots
;*************************************************
  plot   = new (12, graphic)                    ; create graphical array


  res                      = True               ; plot options desired
  res at gsnDraw              = False
  res at gsnFrame             = False
  res at cnFillOn             = True               ; turn on color fill
  res at cnFillPalette        = "amwg256"          ; set color map
  res at cnInfoLabelOn        = False              ; turn off contour info label
  res at cnLinesOn            = False              ; turn off contour lines
  res at cnLineLabelsOn       = False              ; turn off line labels
   

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

  ;res at mpFillOn             = False              ; turn off gray continents
  res at mpCenterLonF         = 0                  ; Centers the plot at 0
  ;res at lbLabelBarOn         = False              ; No single label bar


 res at mpMaxLatF            = 60        ; choose subregion
 res at mpMinLatF            = 0
 res at mpMaxLonF            = 70
 res at mpMinLonF            = -50
 res at mpDataBaseVersion    = "MediumRes"      ; better map outlines
 res at gsnRightString      = " mm"




  resP                     = True               ; panel options
  resP at gsnPanelMainString  = " NDVI anomaly : "+(ymStrt/100)+"-"+(ymLast/100)
  resP at gsnMaximize         = True               ; maximize image
  resP at gsnPanelLabelBar    = True               ; Add common label bar
   




  do nmo=0,2                                  ; loop over the months
     res at gsnLeftString     = months(nmo)
     plot(nmo) = gsn_csm_contour_map(wks,anom(nmo,:,:), res)  ; create plot
  end do


    do nmo=9,11                                  ; loop over the months
     res at gsnLeftString     = months(nmo)
     plot(nmo) = gsn_csm_contour_map(wks,anom(nmo,:,:), res)  ; create plot
   end do




  ;shapefile_name1 = "clip.basin3.shp"
  shapefile_name2 = "IRN_adm0.shp"
  lnres = True
  lnres at gsLineThicknessF = 5
  lnres at gsLineColor = "Black"

  ;shape1 = gsn_add_shapefile_polylines(wks,plot,shapefile_name1,lnres)
  shape2 = gsn_add_shapefile_polylines(wks,plot,shapefile_name2,lnres)
  
  ;gsn_panel(wks,plot,(/2,3/),resP)
  gsn_panel(wks,(/plot(0),plot(1),plot(2),plot(9),plot(10),plot(11)/),(/2,3/),resP)

end






-------------- next part --------------
A non-text attachment was scrubbed...
Name: 1.png
Type: image/png
Size: 220549 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220626/10b08da1/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 2.png
Type: image/png
Size: 350055 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220626/10b08da1/attachment-0003.png>

From jpclark186000 at gmail.com  Sun Jun 26 12:15:44 2022
From: jpclark186000 at gmail.com (Joseph Clark)
Date: Sun, 26 Jun 2022 14:15:44 -0400
Subject: [ncl-talk] How does dpres_plevel determine pressure weights
Message-ID: <CAD_pzv6QhMvfCu+PGP7kULKmbaiyhxauwTgG3Qqp+2SXLC-7ug@mail.gmail.com>

Hi NCL users,

Does anyone know where I can find the source code for the function
dpres_plevel? I'm confused about the exact formula used to determine
pressure weights, particularly near the surface and TOA. There clearly is a
constraint that the sum of dp's adds to the surface pressure, but I'm not
sure exactly what's going on in the function.

Joseph
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220626/58345c7f/attachment.htm>

From dave.allured at noaa.gov  Sun Jun 26 17:42:38 2022
From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate)
Date: Sun, 26 Jun 2022 17:42:38 -0600
Subject: [ncl-talk] How does dpres_plevel determine pressure weights
In-Reply-To: <CAD_pzv6QhMvfCu+PGP7kULKmbaiyhxauwTgG3Qqp+2SXLC-7ug@mail.gmail.com>
References: <CAD_pzv6QhMvfCu+PGP7kULKmbaiyhxauwTgG3Qqp+2SXLC-7ug@mail.gmail.com>
Message-ID: <CALqwTFNkLiASwhobSSOk-dPMrotLjwmGimxaMUKz=eL9-CYYgg@mail.gmail.com>

Download the 6.6.2 source code distribution from the NCL website.  Within,
that function is in ni/src/lib/nfpfort/dpres_plevel_dp.f.  It is





On Sun, Jun 26, 2022 at 12:16 PM Joseph Clark via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi NCL users,
>
> Does anyone know where I can find the source code for the function
> dpres_plevel? I'm confused about the exact formula used to determine
> pressure weights, particularly near the surface and TOA. There clearly is a
> constraint that the sum of dp's adds to the surface pressure, but I'm not
> sure exactly what's going on in the function.
>
> Joseph
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220626/edfe6170/attachment.htm>

From dave.allured at noaa.gov  Sun Jun 26 17:52:30 2022
From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate)
Date: Sun, 26 Jun 2022 17:52:30 -0600
Subject: [ncl-talk] How does dpres_plevel determine pressure weights
In-Reply-To: <CAD_pzv6QhMvfCu+PGP7kULKmbaiyhxauwTgG3Qqp+2SXLC-7ug@mail.gmail.com>
References: <CAD_pzv6QhMvfCu+PGP7kULKmbaiyhxauwTgG3Qqp+2SXLC-7ug@mail.gmail.com>
Message-ID: <CALqwTFOcQbv0xJQ9tirJZmKzf5Tj6iEMXVKQbix1vJggMP34jA@mail.gmail.com>

Download the 6.6.2 source code distribution from the NCL website.  Within,
that function is in ni/src/lib/nfpfort/dpres_plevel_dp.f.  It is called by
a C wrapper function which can be found in ni/src/lib/nfp/preshybW.c.

You can also browse the NCL source code on the NCL github site.


On Sun, Jun 26, 2022 at 12:16 PM Joseph Clark via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi NCL users,
>
> Does anyone know where I can find the source code for the function
> dpres_plevel? I'm confused about the exact formula used to determine
> pressure weights, particularly near the surface and TOA. There clearly is a
> constraint that the sum of dp's adds to the surface pressure, but I'm not
> sure exactly what's going on in the function.
>
> Joseph\
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220626/7d8a72e7/attachment.htm>

From appopson4 at gmail.com  Mon Jun 27 03:51:20 2022
From: appopson4 at gmail.com (Appo derbetini)
Date: Mon, 27 Jun 2022 10:51:20 +0100
Subject: [ncl-talk] Troubles in vertical integral using NCL(vibeta,
 dpres) and Geocat (comp.dpres_plevel)
Message-ID: <CAHVpkHgX6cauobvxZcfPYfgx_hsALm95DFxERpRs6wj5BXXVSQ@mail.gmail.com>

Dear NCL community,

I want to integrate moisture flux using ncl or python.
For this purpose i wrote ncl codes with VIBETA  and DPRES functions. I
tried the same thing in python (GEOCAT.COMP).

Surprisingly, each code is doing its own work !!!

Attached, codes and figures produced.

There is something I am doing wrong?

Any suggestion is welcome.

Best regards
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220627/4f5e6d3f/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: vert_int_vibeta.ncl
Type: application/octet-stream
Size: 2660 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220627/4f5e6d3f/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: verticalIntegral.py
Type: text/x-python
Size: 4287 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220627/4f5e6d3f/attachment.py>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: vert_int_Dpres.ncl
Type: application/octet-stream
Size: 2453 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220627/4f5e6d3f/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Vert_integral_MF_PYTHON.png
Type: image/png
Size: 85363 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220627/4f5e6d3f/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Vert_integral_Vibeta_MF_NCL.png
Type: image/png
Size: 120624 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220627/4f5e6d3f/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Vert_integral_Dpress_MF_NCL.png
Type: image/png
Size: 66685 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220627/4f5e6d3f/attachment-0002.png>

From jpclark186000 at gmail.com  Mon Jun 27 10:02:31 2022
From: jpclark186000 at gmail.com (Joseph Clark)
Date: Mon, 27 Jun 2022 12:02:31 -0400
Subject: [ncl-talk] How does dpres_plevel determine pressure weights
In-Reply-To: <CALqwTFOcQbv0xJQ9tirJZmKzf5Tj6iEMXVKQbix1vJggMP34jA@mail.gmail.com>
References: <CAD_pzv6QhMvfCu+PGP7kULKmbaiyhxauwTgG3Qqp+2SXLC-7ug@mail.gmail.com>
 <CALqwTFOcQbv0xJQ9tirJZmKzf5Tj6iEMXVKQbix1vJggMP34jA@mail.gmail.com>
Message-ID: <CAD_pzv5vi9VStFT4ieruHtFOEA6P5WEGoXy4UzPJD6H-AFzH=Q@mail.gmail.com>

Thanks Dave, very helpful.

For anyone who may want to know in the future (if such an individual
exists), here is the relevant loop that I found in the source code:

dplvl(1) = (plvl(1)+plvl(2))*0.5d0 - ptop do kl=2,klvl-1 dplvl(kl)= 0.5d0*
(plvl(kl+1) - plvl(kl-1)) end do dplvl(klvl) = psfcmx -(plvl(klvl)+plvl(klvl
-1))*0.5d0

The dp's are centered differences except at the top and bottom. At the top,
the dp's are a forward difference divided by 2, minus the pressure at the
top (which is input into the function) and then at the bottom the dp's are
equal to the surface pressure minus a forward difference divided by 2.

Joe

On Sun, Jun 26, 2022 at 7:52 PM Dave Allured - NOAA Affiliate <
dave.allured at noaa.gov> wrote:

> Download the 6.6.2 source code distribution from the NCL website.  Within,
> that function is in ni/src/lib/nfpfort/dpres_plevel_dp.f.  It is called by
> a C wrapper function which can be found in ni/src/lib/nfp/preshybW.c.
>
> You can also browse the NCL source code on the NCL github site.
>
>
> On Sun, Jun 26, 2022 at 12:16 PM Joseph Clark via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Hi NCL users,
>>
>> Does anyone know where I can find the source code for the function
>> dpres_plevel? I'm confused about the exact formula used to determine
>> pressure weights, particularly near the surface and TOA. There clearly is a
>> constraint that the sum of dp's adds to the surface pressure, but I'm not
>> sure exactly what's going on in the function.
>>
>> Joseph\
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220627/96ad1fa1/attachment-0001.htm>

From jpclark186000 at gmail.com  Mon Jun 27 10:09:42 2022
From: jpclark186000 at gmail.com (Joseph Clark)
Date: Mon, 27 Jun 2022 12:09:42 -0400
Subject: [ncl-talk] How does dpres_plevel determine pressure weights
In-Reply-To: <CAD_pzv5vi9VStFT4ieruHtFOEA6P5WEGoXy4UzPJD6H-AFzH=Q@mail.gmail.com>
References: <CAD_pzv6QhMvfCu+PGP7kULKmbaiyhxauwTgG3Qqp+2SXLC-7ug@mail.gmail.com>
 <CALqwTFOcQbv0xJQ9tirJZmKzf5Tj6iEMXVKQbix1vJggMP34jA@mail.gmail.com>
 <CAD_pzv5vi9VStFT4ieruHtFOEA6P5WEGoXy4UzPJD6H-AFzH=Q@mail.gmail.com>
Message-ID: <CAD_pzv67qUvepUXT-=hFPCY6qh_Zfzw_nh0G1SpmrC2urLrb+g@mail.gmail.com>

I think I explained the top/bottom a little off, but the equation is there.
Not really a forward difference. Sorry for any confusion.

On Mon, Jun 27, 2022 at 12:02 PM Joseph Clark <jpclark186000 at gmail.com>
wrote:

> Thanks Dave, very helpful.
>
> For anyone who may want to know in the future (if such an individual
> exists), here is the relevant loop that I found in the source code:
>
> dplvl(1) = (plvl(1)+plvl(2))*0.5d0 - ptop do kl=2,klvl-1 dplvl(kl)= 0.5d0*
> (plvl(kl+1) - plvl(kl-1)) end do dplvl(klvl) = psfcmx -(plvl(klvl)+
> plvl(klvl-1))*0.5d0
>
> The dp's are centered differences except at the top and bottom. At the
> top, the dp's are a forward difference divided by 2, minus the pressure at
> the top (which is input into the function) and then at the bottom the dp's
> are equal to the surface pressure minus a forward difference divided by 2.
>
> Joe
>
> On Sun, Jun 26, 2022 at 7:52 PM Dave Allured - NOAA Affiliate <
> dave.allured at noaa.gov> wrote:
>
>> Download the 6.6.2 source code distribution from the NCL website.
>> Within, that function is in ni/src/lib/nfpfort/dpres_plevel_dp.f.  It is
>> called by a C wrapper function which can be found
>> in ni/src/lib/nfp/preshybW.c.
>>
>> You can also browse the NCL source code on the NCL github site.
>>
>>
>> On Sun, Jun 26, 2022 at 12:16 PM Joseph Clark via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>> Hi NCL users,
>>>
>>> Does anyone know where I can find the source code for the function
>>> dpres_plevel? I'm confused about the exact formula used to determine
>>> pressure weights, particularly near the surface and TOA. There clearly is a
>>> constraint that the sum of dp's adds to the surface pressure, but I'm not
>>> sure exactly what's going on in the function.
>>>
>>> Joseph\
>>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220627/583c34e4/attachment.htm>

From shea at ucar.edu  Mon Jun 27 18:24:02 2022
From: shea at ucar.edu (Dennis Shea)
Date: Mon, 27 Jun 2022 18:24:02 -0600
Subject: [ncl-talk] Troubles in vertical integral using NCL(vibeta,
 dpres) and Geocat (comp.dpres_plevel)
In-Reply-To: <CAHVpkHgX6cauobvxZcfPYfgx_hsALm95DFxERpRs6wj5BXXVSQ@mail.gmail.com>
References: <CAHVpkHgX6cauobvxZcfPYfgx_hsALm95DFxERpRs6wj5BXXVSQ@mail.gmail.com>
Message-ID: <CAOF1d_6PmkTGdRbz_QcpJSCq1Lg7J-=bn4WAQC_P4OXcU14tpA@mail.gmail.com>

Sorry: I have not looked at your codes or figures.

However, if you look at:

*   https://www.ncl.ucar.edu/Applications/wind.shtml*
<https://www.ncl.ucar.edu/Applications/wind.shtml>

Look at examples:

*mfc_div_1.ncl*
*mfc_div_2.ncl*

====*=====*

;---Calculate the Horizontal Moisture Flux Convergence [MFC]


Some extracted code lines

dp   = *dpres_plevel_Wrap*(plev, ps, ptop, 0) ; Pa; layar thickness
dpg  = dp/g
---
uq   = u*q

;---Integrated mass weighted moisture flux components

  uq_dpg = uq*dpg                ; mass weighted 'uq';
[m/s][g/kg][kg/m2]=>[m/s][g/kg]
  iuq    = *dim_sum_n*(uq_dpg, 1)  ; summation/integration

========

On Mon, Jun 27, 2022 at 3:52 AM Appo derbetini via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear NCL community,
>
> I want to integrate moisture flux using ncl or python.
> For this purpose i wrote ncl codes with VIBETA  and DPRES functions. I
> tried the same thing in python (GEOCAT.COMP).
>
> Surprisingly, each code is doing its own work !!!
>
> Attached, codes and figures produced.
>
> There is something I am doing wrong?
>
> Any suggestion is welcome.
>
> Best regards
> _______________________________________________
> 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/20220627/b1556346/attachment.htm>

From appopson4 at gmail.com  Tue Jun 28 08:37:24 2022
From: appopson4 at gmail.com (Appo derbetini)
Date: Tue, 28 Jun 2022 15:37:24 +0100
Subject: [ncl-talk] Troubles in vertical integral using NCL(vibeta,
 dpres) and Geocat (comp.dpres_plevel)
In-Reply-To: <CAOF1d_6PmkTGdRbz_QcpJSCq1Lg7J-=bn4WAQC_P4OXcU14tpA@mail.gmail.com>
References: <CAHVpkHgX6cauobvxZcfPYfgx_hsALm95DFxERpRs6wj5BXXVSQ@mail.gmail.com>
 <CAOF1d_6PmkTGdRbz_QcpJSCq1Lg7J-=bn4WAQC_P4OXcU14tpA@mail.gmail.com>
Message-ID: <CAHVpkHjC-tDFfqi394jH96pN1F2Vo+p2MDbQJjwbuJ1xv=ZASg@mail.gmail.com>

Dear Denis,

Thank you for your quick answer.
I double checked my code and examples you suggested, unfortunately there is
no change.
I am still in trouble !!!

Regards




Le mar. 28 juin 2022 ? 01:24, Dennis Shea <shea at ucar.edu> a ?crit :

> Sorry: I have not looked at your codes or figures.
>
> However, if you look at:
>
> *   https://www.ncl.ucar.edu/Applications/wind.shtml*
> <https://www.ncl.ucar.edu/Applications/wind.shtml>
>
> Look at examples:
>
> *mfc_div_1.ncl*
> *mfc_div_2.ncl*
>
> ====*=====*
>
> ;---Calculate the Horizontal Moisture Flux Convergence [MFC]
>
>
> Some extracted code lines
>
> dp   = *dpres_plevel_Wrap*(plev, ps, ptop, 0) ; Pa; layar thickness
> dpg  = dp/g
> ---
> uq   = u*q
>
> ;---Integrated mass weighted moisture flux components
>
>   uq_dpg = uq*dpg                ; mass weighted 'uq'; [m/s][g/kg][kg/m2]=>[m/s][g/kg]
>   iuq    = *dim_sum_n*(uq_dpg, 1)  ; summation/integration
>
> ========
>
> On Mon, Jun 27, 2022 at 3:52 AM Appo derbetini via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Dear NCL community,
>>
>> I want to integrate moisture flux using ncl or python.
>> For this purpose i wrote ncl codes with VIBETA  and DPRES functions. I
>> tried the same thing in python (GEOCAT.COMP).
>>
>> Surprisingly, each code is doing its own work !!!
>>
>> Attached, codes and figures produced.
>>
>> There is something I am doing wrong?
>>
>> Any suggestion is welcome.
>>
>> Best regards
>> _______________________________________________
>> 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/20220628/36a0bd1d/attachment.htm>

From setareh.rahimi at gmail.com  Tue Jun 28 10:25:04 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Tue, 28 Jun 2022 20:55:04 +0430
Subject: [ncl-talk] How to combine many separate nc files together and
 convert them to monthly?
Message-ID: <CAN+DFgzFiFekJo0yt2By47L00Q8Je0dm4vwiN3_DBsxuhRxRYg@mail.gmail.com>

Dear all,

I have downloaded the Normalized Difference Vegetation Index (NDVI), from
the Advanced Very High-Resolution Radiometer (AVHRR) sensor from the link
below:
https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html

As you can see this link provides daily data. And I have downloaded daily
data from 1986 to 2021. I put the daily data for each year in a separate
folder.
Now I tend to calculate the monthly anomaly of NDVI and do not know how to
combine all these files.

I used the below commands to call the files:
========================================
fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
   print(fils)

   f   = addfiles (fils, "r")

   x   = short2flt(f[:]->NDVI)

   printVarSummary(x)
=========================================
but got the attached error! Please have a look at the attached screenshot.
How can I combine all those files together, please?
I really appreciate any help you can provide.
Best wishes,

-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220628/4b9637de/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 1401-04-07 at 8.47.47 PM.png
Type: image/png
Size: 187042 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220628/4b9637de/attachment.png>

From shea at ucar.edu  Tue Jun 28 11:37:24 2022
From: shea at ucar.edu (Dennis Shea)
Date: Tue, 28 Jun 2022 11:37:24 -0600
Subject: [ncl-talk] How to combine many separate nc files together and
 convert them to monthly?
In-Reply-To: <CAN+DFgzFiFekJo0yt2By47L00Q8Je0dm4vwiN3_DBsxuhRxRYg@mail.gmail.com>
References: <CAN+DFgzFiFekJo0yt2By47L00Q8Je0dm4vwiN3_DBsxuhRxRYg@mail.gmail.com>
Message-ID: <CAOF1d_6X0NJrDC9gT93Xn0okxVGh=4y4S9LBfxof7n1+59O=FA@mail.gmail.com>

As indicated before, your path is wrong.
====
The attached works just fine.

%>  ncl NDVI_TimeSeries.ncl >&! NDVI.Daily_to_Monthly.txt

On Tue, Jun 28, 2022 at 10:25 AM Setareh Rahimi via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear all,
>
> I have downloaded the Normalized Difference Vegetation Index (NDVI), from
> the Advanced Very High-Resolution Radiometer (AVHRR) sensor from the link
> below:
>
> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
>
> As you can see this link provides daily data. And I have downloaded daily
> data from 1986 to 2021. I put the daily data for each year in a separate
> folder.
> Now I tend to calculate the monthly anomaly of NDVI and do not know how to
> combine all these files.
>
> I used the below commands to call the files:
> ========================================
> fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
>    print(fils)
>
>    f   = addfiles (fils, "r")
>
>    x   = short2flt(f[:]->NDVI)
>
>    printVarSummary(x)
> =========================================
> but got the attached error! Please have a look at the attached screenshot.
> How can I combine all those files together, please?
> I really appreciate any help you can provide.
> Best wishes,
>
> --
> S.Rahimi
>
> _______________________________________________
> 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/20220628/63573d06/attachment.htm>
-------------- next part --------------
 Copyright (C) 1995-2017 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 6.5.0 -29Mar2018_0148
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.


Variable: fils
Type: string
Total Size: 40 bytes
            5 values
Number of Dimensions: 1
Dimensions and sizes:	[5]
Coordinates: 
(0)	AVHRR-Land_v005_AVH13C1_NOAA-09_19860101_c20170612092103.nc
(1)	AVHRR-Land_v005_AVH13C1_NOAA-09_19860102_c20170612093403.nc
(2)	AVHRR-Land_v005_AVH13C1_NOAA-09_19861229_c20170612212717.nc
(3)	AVHRR-Land_v005_AVH13C1_NOAA-09_19861230_c20170612214313.nc
(4)	AVHRR-Land_v005_AVH13C1_NOAA-09_19861231_c20170612215936.nc

Variable: x
Type: float
Total Size: 518400000 bytes
            129600000 values
Number of Dimensions: 3
Dimensions and sizes:	[time | 5] x [latitude | 3600] x [longitude | 7200]
Coordinates: 
            time: [1826..2190]
            latitude: [89.975..-89.97499]
            longitude: [-179.975..179.975]
Number Of Attributes: 7
  long_name :	NOAA Climate Data Record of Normalized Difference Vegetation Index
  units :	1
  grid_mapping :	crs
  standard_name :	normalized_difference_vegetation_index
  _FillValue_original :	-9999
  _FillValue :	-9999
  valid_range :	( -0.09999999,  1 )

Variable: xM
Type: float
Total Size: 207360000 bytes
            51840000 values
Number of Dimensions: 3
Dimensions and sizes:	[time | 2] x [latitude | 3600] x [longitude | 7200]
Coordinates: 
            time: [1826..2188]
            latitude: [89.975..-89.97499]
            longitude: [-179.975..179.975]
Number Of Attributes: 9
  _FillValue :	-9999
  long_name :	NOAA Climate Data Record of Normalized Difference Vegetation Index
  units :	1
  grid_mapping :	crs
  standard_name :	normalized_difference_vegetation_index
  _FillValue_original :	-9999
  valid_range :	( -0.09999999,  1 )
  time :	1826
  NCL_tag :	calculate_monthly_values: arith=avg
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NDVI_TimeSeries.ncl
Type: application/octet-stream
Size: 533 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220628/63573d06/attachment.obj>

From setareh.rahimi at gmail.com  Tue Jun 28 11:46:03 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Tue, 28 Jun 2022 22:16:03 +0430
Subject: [ncl-talk] How to combine many separate nc files together and
 convert them to monthly?
In-Reply-To: <CAOF1d_6X0NJrDC9gT93Xn0okxVGh=4y4S9LBfxof7n1+59O=FA@mail.gmail.com>
References: <CAN+DFgzFiFekJo0yt2By47L00Q8Je0dm4vwiN3_DBsxuhRxRYg@mail.gmail.com>
 <CAOF1d_6X0NJrDC9gT93Xn0okxVGh=4y4S9LBfxof7n1+59O=FA@mail.gmail.com>
Message-ID: <CAN+DFgyYnJ5WRSiLkYv+i+CNvaeZMjf==miKvV0-pbPMoqA-Qw@mail.gmail.com>

Dear Dennis,
Many thanks for your email. I run the attached script and received the
attached error.
What is wrong, please?

On Tue, Jun 28, 2022 at 10:07 PM Dennis Shea <shea at ucar.edu> wrote:

> As indicated before, your path is wrong.
> ====
> The attached works just fine.
>
> %>  ncl NDVI_TimeSeries.ncl >&! NDVI.Daily_to_Monthly.txt
>
> On Tue, Jun 28, 2022 at 10:25 AM Setareh Rahimi via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Dear all,
>>
>> I have downloaded the Normalized Difference Vegetation Index (NDVI),
>> from the Advanced Very High-Resolution Radiometer (AVHRR) sensor from
>> the link below:
>>
>> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
>>
>> As you can see this link provides daily data. And I have downloaded daily
>> data from 1986 to 2021. I put the daily data for each year in a separate
>> folder.
>> Now I tend to calculate the monthly anomaly of NDVI and do not know how
>> to combine all these files.
>>
>> I used the below commands to call the files:
>> ========================================
>> fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
>>    print(fils)
>>
>>    f   = addfiles (fils, "r")
>>
>>    x   = short2flt(f[:]->NDVI)
>>
>>    printVarSummary(x)
>> =========================================
>> but got the attached error! Please have a look at the attached screenshot.
>> How can I combine all those files together, please?
>> I really appreciate any help you can provide.
>> Best wishes,
>>
>> --
>> S.Rahimi
>>
>> _______________________________________________
>> 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/20220628/ff0291b2/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 1401-04-07 at 10.14.13 PM.png
Type: image/png
Size: 97395 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220628/ff0291b2/attachment-0001.png>

From shea at ucar.edu  Tue Jun 28 12:53:05 2022
From: shea at ucar.edu (Dennis Shea)
Date: Tue, 28 Jun 2022 12:53:05 -0600
Subject: [ncl-talk] How to combine many separate nc files together and
 convert them to monthly?
In-Reply-To: <CAN+DFgyYnJ5WRSiLkYv+i+CNvaeZMjf==miKvV0-pbPMoqA-Qw@mail.gmail.com>
References: <CAN+DFgzFiFekJo0yt2By47L00Q8Je0dm4vwiN3_DBsxuhRxRYg@mail.gmail.com>
 <CAOF1d_6X0NJrDC9gT93Xn0okxVGh=4y4S9LBfxof7n1+59O=FA@mail.gmail.com>
 <CAN+DFgyYnJ5WRSiLkYv+i+CNvaeZMjf==miKvV0-pbPMoqA-Qw@mail.gmail.com>
Message-ID: <CAOF1d_66kRKOFQsKCGejxxCDxAy7LLpJPdsZL0+2UgR-X5KAcQ@mail.gmail.com>

NCL's *systemfunc* is invoking the Unix/Linux system 'ls'. It is
*independent* of NCL.

>From directory your NCL script is located, enter the following at the
command line:

%>  ls [1-2][8-9-0-1]*/*nc


On Tue, Jun 28, 2022 at 11:46 AM Setareh Rahimi <setareh.rahimi at gmail.com>
wrote:

> Dear Dennis,
> Many thanks for your email. I run the attached script and received the
> attached error.
> What is wrong, please?
>
> On Tue, Jun 28, 2022 at 10:07 PM Dennis Shea <shea at ucar.edu> wrote:
>
>> As indicated before, your path is wrong.
>> ====
>> The attached works just fine.
>>
>> %>  ncl NDVI_TimeSeries.ncl >&! NDVI.Daily_to_Monthly.txt
>>
>> On Tue, Jun 28, 2022 at 10:25 AM Setareh Rahimi via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>> Dear all,
>>>
>>> I have downloaded the Normalized Difference Vegetation Index (NDVI),
>>> from the Advanced Very High-Resolution Radiometer (AVHRR) sensor from
>>> the link below:
>>>
>>> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
>>>
>>> As you can see this link provides daily data. And I have downloaded
>>> daily data from 1986 to 2021. I put the daily data for each year in a
>>> separate folder.
>>> Now I tend to calculate the monthly anomaly of NDVI and do not know how
>>> to combine all these files.
>>>
>>> I used the below commands to call the files:
>>> ========================================
>>> fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
>>>    print(fils)
>>>
>>>    f   = addfiles (fils, "r")
>>>
>>>    x   = short2flt(f[:]->NDVI)
>>>
>>>    printVarSummary(x)
>>> =========================================
>>> but got the attached error! Please have a look at the attached
>>> screenshot.
>>> How can I combine all those files together, please?
>>> I really appreciate any help you can provide.
>>> Best wishes,
>>>
>>> --
>>> S.Rahimi
>>>
>>> _______________________________________________
>>> 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/20220628/997e46a1/attachment.htm>

From setareh.rahimi at gmail.com  Tue Jun 28 17:47:35 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Wed, 29 Jun 2022 04:17:35 +0430
Subject: [ncl-talk] How to combine many separate nc files together and
 convert them to monthly?
In-Reply-To: <CAOF1d_6nQvXTyP4jzAkE51jwVJWn=ya1ih5+D_HPC+A6LG4HEg@mail.gmail.com>
References: <CAN+DFgzFiFekJo0yt2By47L00Q8Je0dm4vwiN3_DBsxuhRxRYg@mail.gmail.com>
 <CAOF1d_6X0NJrDC9gT93Xn0okxVGh=4y4S9LBfxof7n1+59O=FA@mail.gmail.com>
 <CAN+DFgyYnJ5WRSiLkYv+i+CNvaeZMjf==miKvV0-pbPMoqA-Qw@mail.gmail.com>
 <CAOF1d_66kRKOFQsKCGejxxCDxAy7LLpJPdsZL0+2UgR-X5KAcQ@mail.gmail.com>
 <CAN+DFgyyXYEAXaTbQ6ZWnwCk2dg86nrNUYzF1dNc4-iQw_+rNg@mail.gmail.com>
 <CAOF1d_7Zp5kqWb30VryAs=OwEiJato-Hy-qH94EJE3HLjHC_pg@mail.gmail.com>
 <CAN+DFgwWtXMr6bt=vdfKL9ydWx7h4Br28PttuEC2o6p05StvgA@mail.gmail.com>
 <CAOF1d_6nQvXTyP4jzAkE51jwVJWn=ya1ih5+D_HPC+A6LG4HEg@mail.gmail.com>
Message-ID: <CAN+DFgx_Zb5JA4dfc52L9Wj8a13p-HSunG03YLrJOMd6zimdog@mail.gmail.com>

Dear Dennis,
So many thanks for your help. I could finally run the script successfully.
I attached the script and the plot. However, I would like to set the title
of the plot as "NDVI average from 2021-03-03 to 2021-03-12"
But this title overlaps with the main title of the file. How can I remove
the original title, please?
So many thanks for the time you allocate to check my request.
Best wishes,


On Wed, Jun 29, 2022 at 12:57 AM Dennis Shea <shea at ucar.edu> wrote:

> The file you sent,
>
> VIIRS-Land_v001-preliminary_NPP13C1_S-NPP_20210312_c20220419124706.nc
>
> does not have a problem.
>
> As you requested, I ran *https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl
> <https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl>*
> on the file you sent*, I*t is appended.
>
> %> ncl  setareh_VIIRS.ncl
>
> As you requested, I ran the script  on the VIIRS file. It ran
> successfully. The png is attached.
> I think it would be best to send help requests to ncl-talk.
> Maybe others can see issues that I am not seeing.
>
> Regards
> Dennis
>
>
>
> On Tue, Jun 28, 2022 at 1:42 PM Setareh Rahimi <setareh.rahimi at gmail.com>
> wrote:
>
>> Dear Dennis,
>> I tried to plot NDVI for a specific day, hence 12 March 2021 using the
>> file available at :
>> https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl
>> However, I got an error again. It seems that there is a problem with the
>> file. I attached the error message and the nc file. May I ask you please to
>> run the script for this data?
>> I need to make sure that the file is damaged.
>> Thanks in advance,
>>  VIIRS-Land_v001-preliminary_NPP13C1_S-NPP_20210...
>> <https://drive.google.com/file/d/12H5UHVPVoQisuBgo62hs4-66_B7OVXS_/view?usp=drive_web>
>>
>> On Tue, Jun 28, 2022 at 11:35 PM Dennis Shea <shea at ucar.edu> wrote:
>>
>>> These messages do not even look like NCL error messages.
>>>
>>> If you entered the following at the command line and got those messages,
>>> NCL in not involved.:
>>>
>>> %>  ls [1-2][8-9-0-1]*/*nc
>>>
>>> I'm not sure what to suggest.
>>>
>>>
>>>
>>>
>>>
>>>
>>> On Tue, Jun 28, 2022 at 12:59 PM Setareh Rahimi <
>>> setareh.rahimi at gmail.com> wrote:
>>>
>>>> Dear Dennis,
>>>> Thanks again for your advice. I did as you suggested and again got an
>>>> error. (attached)
>>>> Why should I receive this message, please?
>>>> ================================================
>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>> file: NclNetCDF4.c, line: 2493
>>>>
>>>> We Thought It Was Compound Data, But NOT.
>>>> ================================================
>>>>
>>>> On Tue, Jun 28, 2022 at 11:23 PM Dennis Shea <shea at ucar.edu> wrote:
>>>>
>>>>> NCL's *systemfunc* is invoking the Unix/Linux system 'ls'. It is
>>>>> *independent* of NCL.
>>>>>
>>>>> From directory your NCL script is located, enter the following at the
>>>>> command line:
>>>>>
>>>>> %>  ls [1-2][8-9-0-1]*/*nc
>>>>>
>>>>>
>>>>> On Tue, Jun 28, 2022 at 11:46 AM Setareh Rahimi <
>>>>> setareh.rahimi at gmail.com> wrote:
>>>>>
>>>>>> Dear Dennis,
>>>>>> Many thanks for your email. I run the attached script and received
>>>>>> the attached error.
>>>>>> What is wrong, please?
>>>>>>
>>>>>> On Tue, Jun 28, 2022 at 10:07 PM Dennis Shea <shea at ucar.edu> wrote:
>>>>>>
>>>>>>> As indicated before, your path is wrong.
>>>>>>> ====
>>>>>>> The attached works just fine.
>>>>>>>
>>>>>>> %>  ncl NDVI_TimeSeries.ncl >&! NDVI.Daily_to_Monthly.txt
>>>>>>>
>>>>>>> On Tue, Jun 28, 2022 at 10:25 AM Setareh Rahimi via ncl-talk <
>>>>>>> ncl-talk at mailman.ucar.edu> wrote:
>>>>>>>
>>>>>>>> Dear all,
>>>>>>>>
>>>>>>>> I have downloaded the Normalized Difference Vegetation Index
>>>>>>>> (NDVI), from the Advanced Very High-Resolution Radiometer (AVHRR) sensor
>>>>>>>> from the link below:
>>>>>>>>
>>>>>>>> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
>>>>>>>>
>>>>>>>> As you can see this link provides daily data. And I have downloaded
>>>>>>>> daily data from 1986 to 2021. I put the daily data for each year in a
>>>>>>>> separate folder.
>>>>>>>> Now I tend to calculate the monthly anomaly of NDVI and do not know
>>>>>>>> how to combine all these files.
>>>>>>>>
>>>>>>>> I used the below commands to call the files:
>>>>>>>> ========================================
>>>>>>>> fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
>>>>>>>>    print(fils)
>>>>>>>>
>>>>>>>>    f   = addfiles (fils, "r")
>>>>>>>>
>>>>>>>>    x   = short2flt(f[:]->NDVI)
>>>>>>>>
>>>>>>>>    printVarSummary(x)
>>>>>>>> =========================================
>>>>>>>> but got the attached error! Please have a look at the attached
>>>>>>>> screenshot.
>>>>>>>> How can I combine all those files together, please?
>>>>>>>> I really appreciate any help you can provide.
>>>>>>>> Best wishes,
>>>>>>>>
>>>>>>>> --
>>>>>>>> S.Rahimi
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> 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
>>>>>>
>>>>>>
>>>>
>>>> --
>>>> S.Rahimi
>>>>
>>>>
>>
>> --
>> S.Rahimi
>>
>>

-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220629/fa27a4fb/attachment-0001.htm>
-------------- next part --------------
fils = systemfunc("ls VIIRS*nc")
  print(fils)

   f   = addfiles (fils, "r")
  ;ListSetType (f, "cat")
   x   = short2flt(f[:]->NDVI)

   printVarSummary(x)
;************************************************
; Compute the climatology using a function in contributed.ncl
;************************************************
wks = gsn_open_wks("png","ndvi")           ; send graphics to PNG file

  res                  = True                ; plot mods desired
  res at gsnMaximize      = True

  res at cnFillOn         = True                ; color Fill
  res at cnFillMode       = "RasterFill"        ; Raster Mode
  res at cnLinesOn        = False               ; Turn off contour lines

  res at mpCenterLonF     = 0                   ; set map center at 180
  res at mpFillOn         = False


res at cnLevelSelectionMode = "ManualLevels"     ; set manual contour levels

  res at cnMinLevelValF   =  -0.4              ; set min contour level
  res at cnMaxLevelValF   =   0.5              ; one less than max
  res at cnLevelSpacingF  =   0.05              ; set contour spacing
  res at cnFillPalette    = "NEO_div_vegetation_a"   ; <<< VEGETATION ...place holder
setvalues NhlGetWorkspaceObjectId()
    "wsMaximumSize" : 500000000
  end setvalues

 res at gsnCenterString  = "NDVI average from 2021-03-03 to 2021-03-12"

  opt = True
  opt at nval_crit = 10   ; require at least 10 values for the "avg" is calculated.

  xM = dim_avg_n_Wrap( x, 0)
 printVarSummary(xM)

plot = gsn_csm_contour_map_ce(wks,xM, res) ; create plot
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ndvi.png
Type: image/png
Size: 311666 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220629/fa27a4fb/attachment-0001.png>

From asphilli at ucar.edu  Wed Jun 29 12:44:59 2022
From: asphilli at ucar.edu (Adam Phillips)
Date: Wed, 29 Jun 2022 12:44:59 -0600
Subject: [ncl-talk] How to combine many separate nc files together and
 convert them to monthly?
In-Reply-To: <CAN+DFgx_Zb5JA4dfc52L9Wj8a13p-HSunG03YLrJOMd6zimdog@mail.gmail.com>
References: <CAN+DFgzFiFekJo0yt2By47L00Q8Je0dm4vwiN3_DBsxuhRxRYg@mail.gmail.com>
 <CAOF1d_6X0NJrDC9gT93Xn0okxVGh=4y4S9LBfxof7n1+59O=FA@mail.gmail.com>
 <CAN+DFgyYnJ5WRSiLkYv+i+CNvaeZMjf==miKvV0-pbPMoqA-Qw@mail.gmail.com>
 <CAOF1d_66kRKOFQsKCGejxxCDxAy7LLpJPdsZL0+2UgR-X5KAcQ@mail.gmail.com>
 <CAN+DFgyyXYEAXaTbQ6ZWnwCk2dg86nrNUYzF1dNc4-iQw_+rNg@mail.gmail.com>
 <CAOF1d_7Zp5kqWb30VryAs=OwEiJato-Hy-qH94EJE3HLjHC_pg@mail.gmail.com>
 <CAN+DFgwWtXMr6bt=vdfKL9ydWx7h4Br28PttuEC2o6p05StvgA@mail.gmail.com>
 <CAOF1d_6nQvXTyP4jzAkE51jwVJWn=ya1ih5+D_HPC+A6LG4HEg@mail.gmail.com>
 <CAN+DFgx_Zb5JA4dfc52L9Wj8a13p-HSunG03YLrJOMd6zimdog@mail.gmail.com>
Message-ID: <CACvWuqDvYAvZempLRo0o-DqiwW3yt7bhF0z2gV6D6cnGwfac8g@mail.gmail.com>

Hi Setareh,
You can remove the left subtitle by setting gsnLeftString = "".
Alternatively, you can leave that subtitle alone and move the main title up
by using tiMainString instead of gsnCenterString.
See: https://www.ncl.ucar.edu/Applications/title.shtml
Adam

On Tue, Jun 28, 2022 at 5:47 PM Setareh Rahimi via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear Dennis,
> So many thanks for your help. I could finally run the script successfully.
> I attached the script and the plot. However, I would like to set the title
> of the plot as "NDVI average from 2021-03-03 to 2021-03-12"
> But this title overlaps with the main title of the file. How can I remove
> the original title, please?
> So many thanks for the time you allocate to check my request.
> Best wishes,
>
>
> On Wed, Jun 29, 2022 at 12:57 AM Dennis Shea <shea at ucar.edu> wrote:
>
>> The file you sent,
>>
>> VIIRS-Land_v001-preliminary_NPP13C1_S-NPP_20210312_c20220419124706.nc
>>
>> does not have a problem.
>>
>> As you requested, I ran *https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl
>> <https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl>*
>> on the file you sent*, I*t is appended.
>>
>> %> ncl  setareh_VIIRS.ncl
>>
>> As you requested, I ran the script  on the VIIRS file. It ran
>> successfully. The png is attached.
>> I think it would be best to send help requests to ncl-talk.
>> Maybe others can see issues that I am not seeing.
>>
>> Regards
>> Dennis
>>
>>
>>
>> On Tue, Jun 28, 2022 at 1:42 PM Setareh Rahimi <setareh.rahimi at gmail.com>
>> wrote:
>>
>>> Dear Dennis,
>>> I tried to plot NDVI for a specific day, hence 12 March 2021 using the
>>> file available at :
>>> https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl
>>> However, I got an error again. It seems that there is a problem with the
>>> file. I attached the error message and the nc file. May I ask you please to
>>> run the script for this data?
>>> I need to make sure that the file is damaged.
>>> Thanks in advance,
>>>  VIIRS-Land_v001-preliminary_NPP13C1_S-NPP_20210...
>>> <https://drive.google.com/file/d/12H5UHVPVoQisuBgo62hs4-66_B7OVXS_/view?usp=drive_web>
>>>
>>> On Tue, Jun 28, 2022 at 11:35 PM Dennis Shea <shea at ucar.edu> wrote:
>>>
>>>> These messages do not even look like NCL error messages.
>>>>
>>>> If you entered the following at the command line and got those
>>>> messages, NCL in not involved.:
>>>>
>>>> %>  ls [1-2][8-9-0-1]*/*nc
>>>>
>>>> I'm not sure what to suggest.
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> On Tue, Jun 28, 2022 at 12:59 PM Setareh Rahimi <
>>>> setareh.rahimi at gmail.com> wrote:
>>>>
>>>>> Dear Dennis,
>>>>> Thanks again for your advice. I did as you suggested and again got an
>>>>> error. (attached)
>>>>> Why should I receive this message, please?
>>>>> ================================================
>>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>>> file: NclNetCDF4.c, line: 2493
>>>>>
>>>>> We Thought It Was Compound Data, But NOT.
>>>>> ================================================
>>>>>
>>>>> On Tue, Jun 28, 2022 at 11:23 PM Dennis Shea <shea at ucar.edu> wrote:
>>>>>
>>>>>> NCL's *systemfunc* is invoking the Unix/Linux system 'ls'. It is
>>>>>> *independent* of NCL.
>>>>>>
>>>>>> From directory your NCL script is located, enter the following at the
>>>>>> command line:
>>>>>>
>>>>>> %>  ls [1-2][8-9-0-1]*/*nc
>>>>>>
>>>>>>
>>>>>> On Tue, Jun 28, 2022 at 11:46 AM Setareh Rahimi <
>>>>>> setareh.rahimi at gmail.com> wrote:
>>>>>>
>>>>>>> Dear Dennis,
>>>>>>> Many thanks for your email. I run the attached script and received
>>>>>>> the attached error.
>>>>>>> What is wrong, please?
>>>>>>>
>>>>>>> On Tue, Jun 28, 2022 at 10:07 PM Dennis Shea <shea at ucar.edu> wrote:
>>>>>>>
>>>>>>>> As indicated before, your path is wrong.
>>>>>>>> ====
>>>>>>>> The attached works just fine.
>>>>>>>>
>>>>>>>> %>  ncl NDVI_TimeSeries.ncl >&! NDVI.Daily_to_Monthly.txt
>>>>>>>>
>>>>>>>> On Tue, Jun 28, 2022 at 10:25 AM Setareh Rahimi via ncl-talk <
>>>>>>>> ncl-talk at mailman.ucar.edu> wrote:
>>>>>>>>
>>>>>>>>> Dear all,
>>>>>>>>>
>>>>>>>>> I have downloaded the Normalized Difference Vegetation Index
>>>>>>>>> (NDVI), from the Advanced Very High-Resolution Radiometer (AVHRR) sensor
>>>>>>>>> from the link below:
>>>>>>>>>
>>>>>>>>> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
>>>>>>>>>
>>>>>>>>> As you can see this link provides daily data. And I have
>>>>>>>>> downloaded daily data from 1986 to 2021. I put the daily data for each year
>>>>>>>>> in a separate folder.
>>>>>>>>> Now I tend to calculate the monthly anomaly of NDVI and do not
>>>>>>>>> know how to combine all these files.
>>>>>>>>>
>>>>>>>>> I used the below commands to call the files:
>>>>>>>>> ========================================
>>>>>>>>> fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
>>>>>>>>>    print(fils)
>>>>>>>>>
>>>>>>>>>    f   = addfiles (fils, "r")
>>>>>>>>>
>>>>>>>>>    x   = short2flt(f[:]->NDVI)
>>>>>>>>>
>>>>>>>>>    printVarSummary(x)
>>>>>>>>> =========================================
>>>>>>>>> but got the attached error! Please have a look at the attached
>>>>>>>>> screenshot.
>>>>>>>>> How can I combine all those files together, please?
>>>>>>>>> I really appreciate any help you can provide.
>>>>>>>>> Best wishes,
>>>>>>>>>
>>>>>>>>> --
>>>>>>>>> S.Rahimi
>>>>>>>>>
>>>>>>>>> _______________________________________________
>>>>>>>>> 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
>>>>>>>
>>>>>>>
>>>>>
>>>>> --
>>>>> S.Rahimi
>>>>>
>>>>>
>>>
>>> --
>>> S.Rahimi
>>>
>>>
>
> --
> S.Rahimi
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220629/29f8853c/attachment.htm>

From setareh.rahimi at gmail.com  Thu Jun 30 07:37:07 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Thu, 30 Jun 2022 18:07:07 +0430
Subject: [ncl-talk] How to combine many separate nc files together and
 convert them to monthly?
In-Reply-To: <CACvWuqDvYAvZempLRo0o-DqiwW3yt7bhF0z2gV6D6cnGwfac8g@mail.gmail.com>
References: <CAN+DFgzFiFekJo0yt2By47L00Q8Je0dm4vwiN3_DBsxuhRxRYg@mail.gmail.com>
 <CAOF1d_6X0NJrDC9gT93Xn0okxVGh=4y4S9LBfxof7n1+59O=FA@mail.gmail.com>
 <CAN+DFgyYnJ5WRSiLkYv+i+CNvaeZMjf==miKvV0-pbPMoqA-Qw@mail.gmail.com>
 <CAOF1d_66kRKOFQsKCGejxxCDxAy7LLpJPdsZL0+2UgR-X5KAcQ@mail.gmail.com>
 <CAN+DFgyyXYEAXaTbQ6ZWnwCk2dg86nrNUYzF1dNc4-iQw_+rNg@mail.gmail.com>
 <CAOF1d_7Zp5kqWb30VryAs=OwEiJato-Hy-qH94EJE3HLjHC_pg@mail.gmail.com>
 <CAN+DFgwWtXMr6bt=vdfKL9ydWx7h4Br28PttuEC2o6p05StvgA@mail.gmail.com>
 <CAOF1d_6nQvXTyP4jzAkE51jwVJWn=ya1ih5+D_HPC+A6LG4HEg@mail.gmail.com>
 <CAN+DFgx_Zb5JA4dfc52L9Wj8a13p-HSunG03YLrJOMd6zimdog@mail.gmail.com>
 <CACvWuqDvYAvZempLRo0o-DqiwW3yt7bhF0z2gV6D6cnGwfac8g@mail.gmail.com>
Message-ID: <CAN+DFgwJowR4gb5TBYAcqZh_e4JGGEHJR8QznJd7W5w9_NnvFA@mail.gmail.com>

Dear Adam,
So many thanks for your help. I could fix it now!
Best wishes,

On Wed, Jun 29, 2022 at 11:15 PM Adam Phillips <asphilli at ucar.edu> wrote:

> Hi Setareh,
> You can remove the left subtitle by setting gsnLeftString = "".
> Alternatively, you can leave that subtitle alone and move the main title up
> by using tiMainString instead of gsnCenterString.
> See: https://www.ncl.ucar.edu/Applications/title.shtml
> Adam
>
> On Tue, Jun 28, 2022 at 5:47 PM Setareh Rahimi via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Dear Dennis,
>> So many thanks for your help. I could finally run the script
>> successfully. I attached the script and the plot. However, I would like to
>> set the title of the plot as "NDVI average from 2021-03-03 to 2021-03-12"
>> But this title overlaps with the main title of the file. How can I remove
>> the original title, please?
>> So many thanks for the time you allocate to check my request.
>> Best wishes,
>>
>>
>> On Wed, Jun 29, 2022 at 12:57 AM Dennis Shea <shea at ucar.edu> wrote:
>>
>>> The file you sent,
>>>
>>> VIIRS-Land_v001-preliminary_NPP13C1_S-NPP_20210312_c20220419124706.nc
>>>
>>> does not have a problem.
>>>
>>> As you requested, I ran *https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl
>>> <https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl>*
>>> on the file you sent*, I*t is appended.
>>>
>>> %> ncl  setareh_VIIRS.ncl
>>>
>>> As you requested, I ran the script  on the VIIRS file. It ran
>>> successfully. The png is attached.
>>> I think it would be best to send help requests to ncl-talk.
>>> Maybe others can see issues that I am not seeing.
>>>
>>> Regards
>>> Dennis
>>>
>>>
>>>
>>> On Tue, Jun 28, 2022 at 1:42 PM Setareh Rahimi <setareh.rahimi at gmail.com>
>>> wrote:
>>>
>>>> Dear Dennis,
>>>> I tried to plot NDVI for a specific day, hence 12 March 2021 using the
>>>> file available at :
>>>> https://www.ncl.ucar.edu/Applications/Scripts/ndvi_3.ncl
>>>> However, I got an error again. It seems that there is a problem with
>>>> the file. I attached the error message and the nc file. May I ask you
>>>> please to run the script for this data?
>>>> I need to make sure that the file is damaged.
>>>> Thanks in advance,
>>>>  VIIRS-Land_v001-preliminary_NPP13C1_S-NPP_20210...
>>>> <https://drive.google.com/file/d/12H5UHVPVoQisuBgo62hs4-66_B7OVXS_/view?usp=drive_web>
>>>>
>>>> On Tue, Jun 28, 2022 at 11:35 PM Dennis Shea <shea at ucar.edu> wrote:
>>>>
>>>>> These messages do not even look like NCL error messages.
>>>>>
>>>>> If you entered the following at the command line and got those
>>>>> messages, NCL in not involved.:
>>>>>
>>>>> %>  ls [1-2][8-9-0-1]*/*nc
>>>>>
>>>>> I'm not sure what to suggest.
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> On Tue, Jun 28, 2022 at 12:59 PM Setareh Rahimi <
>>>>> setareh.rahimi at gmail.com> wrote:
>>>>>
>>>>>> Dear Dennis,
>>>>>> Thanks again for your advice. I did as you suggested and again got an
>>>>>> error. (attached)
>>>>>> Why should I receive this message, please?
>>>>>> ================================================
>>>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>>>> ncattname: ncid 23986176: NetCDF: Variable not found
>>>>>> file: NclNetCDF4.c, line: 2493
>>>>>>
>>>>>> We Thought It Was Compound Data, But NOT.
>>>>>> ================================================
>>>>>>
>>>>>> On Tue, Jun 28, 2022 at 11:23 PM Dennis Shea <shea at ucar.edu> wrote:
>>>>>>
>>>>>>> NCL's *systemfunc* is invoking the Unix/Linux system 'ls'. It is
>>>>>>> *independent* of NCL.
>>>>>>>
>>>>>>> From directory your NCL script is located, enter the following at
>>>>>>> the command line:
>>>>>>>
>>>>>>> %>  ls [1-2][8-9-0-1]*/*nc
>>>>>>>
>>>>>>>
>>>>>>> On Tue, Jun 28, 2022 at 11:46 AM Setareh Rahimi <
>>>>>>> setareh.rahimi at gmail.com> wrote:
>>>>>>>
>>>>>>>> Dear Dennis,
>>>>>>>> Many thanks for your email. I run the attached script and received
>>>>>>>> the attached error.
>>>>>>>> What is wrong, please?
>>>>>>>>
>>>>>>>> On Tue, Jun 28, 2022 at 10:07 PM Dennis Shea <shea at ucar.edu> wrote:
>>>>>>>>
>>>>>>>>> As indicated before, your path is wrong.
>>>>>>>>> ====
>>>>>>>>> The attached works just fine.
>>>>>>>>>
>>>>>>>>> %>  ncl NDVI_TimeSeries.ncl >&! NDVI.Daily_to_Monthly.txt
>>>>>>>>>
>>>>>>>>> On Tue, Jun 28, 2022 at 10:25 AM Setareh Rahimi via ncl-talk <
>>>>>>>>> ncl-talk at mailman.ucar.edu> wrote:
>>>>>>>>>
>>>>>>>>>> Dear all,
>>>>>>>>>>
>>>>>>>>>> I have downloaded the Normalized Difference Vegetation Index
>>>>>>>>>> (NDVI), from the Advanced Very High-Resolution Radiometer (AVHRR) sensor
>>>>>>>>>> from the link below:
>>>>>>>>>>
>>>>>>>>>> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
>>>>>>>>>>
>>>>>>>>>> As you can see this link provides daily data. And I have
>>>>>>>>>> downloaded daily data from 1986 to 2021. I put the daily data for each year
>>>>>>>>>> in a separate folder.
>>>>>>>>>> Now I tend to calculate the monthly anomaly of NDVI and do not
>>>>>>>>>> know how to combine all these files.
>>>>>>>>>>
>>>>>>>>>> I used the below commands to call the files:
>>>>>>>>>> ========================================
>>>>>>>>>> fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
>>>>>>>>>>    print(fils)
>>>>>>>>>>
>>>>>>>>>>    f   = addfiles (fils, "r")
>>>>>>>>>>
>>>>>>>>>>    x   = short2flt(f[:]->NDVI)
>>>>>>>>>>
>>>>>>>>>>    printVarSummary(x)
>>>>>>>>>> =========================================
>>>>>>>>>> but got the attached error! Please have a look at the attached
>>>>>>>>>> screenshot.
>>>>>>>>>> How can I combine all those files together, please?
>>>>>>>>>> I really appreciate any help you can provide.
>>>>>>>>>> Best wishes,
>>>>>>>>>>
>>>>>>>>>> --
>>>>>>>>>> S.Rahimi
>>>>>>>>>>
>>>>>>>>>> _______________________________________________
>>>>>>>>>> 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
>>>>>>>>
>>>>>>>>
>>>>>>
>>>>>> --
>>>>>> S.Rahimi
>>>>>>
>>>>>>
>>>>
>>>> --
>>>> S.Rahimi
>>>>
>>>>
>>
>> --
>> S.Rahimi
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>


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