[ncl-talk] Help with copying array into a larger array

[Laura Fowler]

Hi Mary:

Thanks. This makes complete sense. I thought that since the latitude
resolution of precip was 400 whereas that of data was 720, they could
not have the same name "lat". This is the reason why I used "lat2" for
data instead. It works much better now.

Have a good week-end.
Laura


On Fri, Jun 8, 2018 at 12:08 PM, Mary Haley <haley at ucar.edu> wrote:
> Hi Laura,
>
> For reasons I don't understand right now, it doesn't like that one variable
> has a "lat" dimension, while the other one is "lat2".
>
> If you rename the data "lat2" dimension to "lat" before you do the copy, I
> think it should work (I tried it on a small test program).
>
> data!1 = "lat"
>
> -Mary
>
>
>
>
>
> On Fri, Jun 8, 2018 at 10:43 AM, Laura Fowler <laura at ucar.edu> wrote:
>>
>> Good morning:
>>
>> I would like to copy TRMM data into a global array:
>>
>> printVarSummary of my TRMM data shows:
>>
>> Variable: precip
>> Type: float
>> Total Size: 2304000 bytes
>>             576000 values
>> Number of Dimensions: 2
>> Dimensions and sizes: [lon | 1440] x [lat | 400]
>> Coordinates:
>>             lon: [-179.875..179.875]
>>             lat: [-49.875..49.875]
>> Number Of Attributes: 5
>>   _FillValue : -9999.9
>>   units : mm
>>   long_name : Daily accumulated precipitation (combined microwave-IR)
>> estimate with gauge calibration over land
>>   coordinates : lat lon
>>   average_op_ncl : dim_avg_n over dimension(s): ncl_join
>>
>>
>> I would like to put the TRMM data (limited between the latitudes
>> -49.875 to 49.875 into an array between the latitudes  -89.875 to
>> 89.875. printVarSummary of my output TRMM data shows:
>>
>> Variable: data
>> Type: float
>> Total Size: 4147200 bytes
>>             1036800 values
>> Number of Dimensions: 2
>> Dimensions and sizes: [lon | 1440] x [lat2 | 720]
>> Coordinates:
>>             lon: [-179.875..179.875]
>>             lat2: [-89.875..89.875]
>> Number Of Attributes: 1
>>   _FillValue : 9.96921e+36
>>
>>
>> I tried the following. To put precip into data, I tried:
>> data(:,{-49.875:49.875}) = precip(:,{-49.875:49.875})
>> or
>> data(:,{lat2|-49.875:49.875}) = precip(:,{lat|-49.875:49.875})
>>
>> but it is not working. The last error I get is:
>>
>> fatal:syntax error: line 59 in file
>> global_mean.precipitation_minus_TRMM-3B42.ncl before or near lat2
>> data(:,{lat2
>>
>>
>> After looking at the documentation, I am not sure what I am doing wrong.
>>
>> 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

!-------------------------------------------------------------------------------------------------------------