[ncl-talk] Computing max of monthly NO2 data each year

[Jack Ogaja]

Yes, I oversaw it.

Jack

> On 29 Sep 2015, at 18:22, Doumbia <thiernodoumbia at yahoo.fr> wrote:
> 
> I add some thing in your code and it gives me the good results.
> yr = i would be yr = i+1.
> 
> The new script is:
> do n=0,7
>    yr=0
>    k =0
>    do i = 11,275,12
>       EUmax(n,k) = max(EU_plot(n,yr:i))
>       k=k+1
>       yr=i+1 
>    end do
> end do
> 
> Jack, thank again.
> Thierno
> 
>> Le 29 Sep 2015 à 18:02, Doumbia <thiernodoumbia at yahoo.fr> a écrit :
>> 
>> Hi Jack,
>> I really appreciate your teaching skills.
>> 
>> I recognize that there’s weird things in my script but I apologize for that, I am new in NCL.
>> 
>> Unfortunately, your suggestions don’t give the good results i.e. the maximum of monthly NO2 data for each year. I compared these results with those obtained from excel.
>> 
>> You are right, the "EU(120:395,0)” is for 1990 to 2012 (23 years).
>> 
>> Thank you for your help.
>> 
>> Thierno
>> 
>> 
>>> 
>>> ...
>>> “The results I’ve obtained aren’t what I expected.”
>>> ...
>>> Well, even Aristotle would struggle to figure out what you expected if you didn't tell him…
>>> 
>>> Your script still got visible mistakes though. For example the loop is quite extravagant:
>>> …
>>> k = ispan(0,11,1)
>>> do n=0,7
>>> do i = 0,275,12
>>>  do j=0,11
>>>    EUmax(n,:) = max(EU_plot(n,i+k(j))) <== “j” will always be the same type as, and equal to “k(j)”?
>>>  end do
>>> end do
>>> end do
>>> …
>>> and since the left side of your equation has a free dimension, the variable will always be overwritten up to the end of the loop.
>>> 
>>> You can clean the loop and do something like,
>>> do n=0,7
>>>   yr=0
>>>   k =0
>>>   do i = 11,275,12
>>>      EUmax(n,k) = max(EU_plot(n,yr:i))
>>>      k=k+1
>>>      yr=i
>>>   end do
>>> end do
>>> 
>>> I was also wondering whether from your file “EU(120:395,0)” will give you information from 2000 to 2012(12 years)? Probably that is the unexpected result you mentioned. Who knows?
>>> 
>>> I hope that helps.
>>> 
>>> Jack
>>> 
>>>> On 29 Sep 2015, at 15:02, Doumbia <thiernodoumbia at yahoo.fr> wrote:
>>>> 
>>>> Dear ncl users,
>>>> 
>>>> I am trying to compute the monthly maximum NO2 in a yearly basis.
>>>> 
>>>> My input file contains 7 columns representing different variables and 396 lines (monthly mean NO2 data from 1980 to 2012). I am only interested by the 2000-2012 period .
>>>> 
>>>> The results I’ve obtained aren’t what I expected.
>>>> <max_of_monthly_data.ncl><Europe_bis_ncl.txt>
>>>> I will appreciate any help to resolve this.
>>>> 
>>>> Please find attached my script as well as my input file.
>>>> 
>>>> Thank you in advance,
>>>> Thierno
>>>> 
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>> 
>> 
>