[ncl-talk] Wondering why I am getting this dimensions error
Barry Lynn
barry.h.lynn at gmail.com
Wed Nov 27 10:49:18 MST 2019
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191127/44298271/attachment-0001.html>
More information about the ncl-talk
mailing list