[ncl-talk] how to fill an array with an other

[Laura Fowler]

Hi Rashed:
What your proposed only works if I were to loop over cells. This would take
quite some time if cells were a large number? There must be some kind of an
array syntax to do this.
Thanks,
Laura

On Tue, May 5, 2020 at 12:08 PM Rashed Mahmood <rashidcomsis at gmail.com>
wrote:

> Hi Laura,
> The error message is telling what is wrong, you defining a one dimensional
> and then trying to assign it two dimensional data:
> *ppbl = new((/cells/),typeof(p))*
>
> *ppbl(0:cells-1) = p(0:cells-1,kpbl(0:cells-1)) ; this is 2D and even
> wrong dimension sizes note the nVertLevels*
>
>
> I am not sure if the following is what you want?
>
> vlevel = 0 ; 0 or any other VertLevels
> ppbl = new((/cells/),typeof(p))
> ppbl(0:cells-1) = p(0:cells-1,vlevel)
>
>
> Cheers,
> Rashed
>
>
>
> On Tue, May 5, 2020 at 8:37 AM Laura Fowler via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
>> Hello:
>> I cannot figure out how to do the following and would love to get some
>> help:
>>
>> I have an array kpbl which contains the vertical index of the PBL top:
>> printVarSummary gives:
>>
>> Variable: kpbl
>>
>> Type: integer
>>
>> Total Size: 15728832 bytes
>>
>>             3932208 values
>>
>> Number of Dimensions: 2
>>
>> Dimensions and sizes: [Time | 13] x [cells | 163842]
>>
>> Coordinates:
>>
>> Number Of Attributes: 2
>>
>>   units : unitless
>>
>>   long_name : index level of PBL top
>>
>>
>>
>> I have a second array p which contains the 3D pressure: print VarSummary
>> gives:
>>
>>
>> Variable: p
>>
>> Type: float
>>
>> Total Size: 865085760 bytes
>>
>>             216271440 values
>>
>> Number of Dimensions: 3
>>
>> Dimensions and sizes: [Time | 13] x [cells | 163842] x [nVertLevels | 55]
>>
>> Coordinates:
>>
>> Number Of Attributes: 2
>>
>>   units : hPa
>>
>>   long_name : Pressure
>>
>>
>>
>> What I would like to get is the pressure at the index level kpbl, i.e.
>> pressure at the PBL top. In fortran and for one time-level (itime), I would
>> have a loop like:
>>
>> do icell = 0, cells-1
>>      ppbl(itime,icell) = p(itime,icell,kpbl(iCell))
>> end do
>>
>> In NCL, I tried the following (for one time-level):
>> ppbl = new((/cells/),typeof(p))
>> ppbl(0:cells-1) = p(0:cells-1,kpbl(0:cells-1))
>>
>> but this does not work. I get an error message that the left and right
>> side do not have the same dimensions. I am definitely running out of ideas
>> on how to solve this and would appreciate a hint.
>>
>> Many thanks,
>> Laura
>>
>> --
>>
>> !-------------------------------------------------------------------------------------------------------------
>> Laura D. Fowler
>>
>> Mesoscale and Microscale Meteorology Division (MMM)
>> National Center for Atmospheric Research
>> P.O. Box 3000, Boulder CO 80307-3000
>>
>> e-mail: laura at ucar.edu
>> phone: 303-497-1628
>>
>>
>> !-------------------------------------------------------------------------------------------------------------
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

-- 
!-------------------------------------------------------------------------------------------------------------
Laura D. Fowler

Mesoscale and Microscale Meteorology Division (MMM)
National Center for Atmospheric Research
P.O. Box 3000, Boulder CO 80307-3000

e-mail: laura at ucar.edu
phone: 303-497-1628

!-------------------------------------------------------------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200506/566413c2/attachment.html>