[ncl-talk] Wondering why I am getting this dimensions error

[Barry Lynn]

Hi:

The problem with the previous suggestion was that the sum_q became a 2
dimensional variable..  I needed it to be a 1D variables. I found this out
by using the "function" printVarSummary

Here is the correct solution for my needs, based on the very helpful
comments up to now.

Thanks

(sum_q was defined as a 2d variable -- note: sum_q is not really a sum, as
used here, but defined as 1 whereever qs > 0.1).

       sum_q = 0

       sum_q(0,0) = 1 ; just so I don't divide by 0!

       sum_q =  where(qs(k,:,:).ge.0.1,sum_q+1,sum_q)

       add_q = sum(sum_q); a 1d variable

 ;      printVarSummary(qs_1(k))

  ;     printVarSummary(sum_q)

  ;     printVarSummary(add_q)

       qs_1(k)  = sum(qs(k,:,:))/add_q+qs_1(k)

On Wed, Nov 27, 2019 at 8:18 PM Dennis Shea <shea at ucar.edu> wrote:

> re:   :=
> Yes.
>
> On Wed, Nov 27, 2019 at 10:49 AM Barry Lynn via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
>> Hi Dave:
>>
>> Thank you for your suggestions.
>>
>> I actually need to sum_q for all elements of the remaining (:,:) 2d array
>> to obtain a single number.
>>
>> I am dividing one number by it.
>>
>> I assume that it is okay to put the := and that I will get the sum over
>> all 2d elements?
>>
>> On Wed, Nov 27, 2019 at 7:40 PM Dave Allured - NOAA Affiliate <
>> dave.allured at noaa.gov> wrote:
>>
>>> "sum_q = 1." defines this variable as scalar.  Then the where statement
>>> attempts to write a 2-D array onto a scalar.  The assignment statement is
>>> not allowed to change dimensions.
>>>
>>> One way to fix this would be to use the reassignment operator ":=",
>>> which IS allowed to change dimensions on assignment.
>>>
>>> Another way to fix this would be to define sum_q as a 2-D array before
>>> the start of the loop.  I suspect that is what you really intended.
>>>
>>>
>>> On Wed, Nov 27, 2019 at 10:22 AM Barry Lynn via ncl-talk <
>>> ncl-talk at ucar.edu> wrote:
>>>
>>>> Hi:
>>>>
>>>> I have some simple code.
>>>>
>>>>       do k = 0,dims3d(0)-1
>>>>
>>>>        qs(k,:,:) = where(qs(k,:,:).lt.0.1,q_null(k,:,:),qs(k,:,:))
>>>>
>>>>        sum_q = 1.
>>>>
>>>>        sum_q =  where(qs(k,:,:).ge.0.1,sum_q+1,sum_q)
>>>>
>>>>        qs_1(k) = sum(qs(k,:,:))/sum_q+qs_1(k)
>>>>
>>>>       end do
>>>>
>>>>
>>>> Why can't I do the "where" sum?
>>>>
>>>>
>>>> fatal:Number of dimensions on right hand side do not match number of
>>>> dimension in left hand side
>>>>
>>>> --
>>>> Barry H. Lynn, Ph.D
>>>> Senior Associate Scientist, Lecturer,
>>>> The Institute of the Earth Science,
>>>> The Hebrew University of Jerusalem,
>>>> Givat Ram, Jerusalem 91904, Israel
>>>> Tel: 972 547 231 170
>>>> Fax: (972)-25662581
>>>>
>>>> C.E.O, Weather It Is, LTD
>>>> Weather and Climate Focus
>>>> 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 Associate Scientist, Lecturer,
>> The Institute of the Earth Science,
>> The Hebrew University of Jerusalem,
>> Givat Ram, Jerusalem 91904, Israel
>> Tel: 972 547 231 170
>> Fax: (972)-25662581
>>
>> C.E.O, Weather It Is, LTD
>> Weather and Climate Focus
>> http://weather-it-is.com
>> Jerusalem, Israel
>> Local: 02 930 9525
>> Cell: 054 7 231 170
>> Int-IS: x972 2 930 9525
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

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

C.E.O, Weather It Is, LTD
Weather and Climate Focus
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: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191127/18480ddd/attachment.html>