[ncl-talk] Error subscript mjoclivar_14 : Dimension sizes of left hand side and right hand side of assignment do not match
reza tisa
rezatisa1 at gmail.com
Mon Dec 30 08:46:25 MST 2019
Thanks Mr.Dennis
Pada tanggal Sen, 30 Des 2019 pukul 21.45 Dennis Shea <shea at ucar.edu>
menulis:
> Whenever NCL tells you
>
> f*atal:Dimension sizes of left hand side and right hand side of
> assignment do not match*
> fatal:["Execute.c":8635]:Execute: Error occurred at or near *line 82*
>
> You should use printVarSummary() before the offending line
> *.*
>
> As noted by DaveA:
>
>
> Dimensions and sizes: *[time | 731] x [lat | 73] x [lon | 145]*
> Coordinates:
> *lon: [ 0.. 360]*
> Number Of Attributes: 2
> _FillValue : -9.99e+08
> long_name : Anomalies from Daily Climatology
>
> Dimensions and sizes: *[time | 731] x [lat | 73] x [lon | 144]*
> Coordinates:
> * lon: [ 0..357.5]*
> Number Of Attributes: 3
> _FillValue : 32766
> long_name : Anomalies: Daily OLR
> units : W/m^2
>
> On Mon, Dec 30, 2019 at 4:12 AM reza tisa via ncl-talk <ncl-talk at ucar.edu>
> wrote:
>
>> It works! Thank you so much Mr.Dave. May you blessed with good health
>> and happiness
>>
>>
>> Pada tanggal Sen, 30 Des 2019 pukul 16.43 Dave Allured - NOAA Affiliate <
>> dave.allured at noaa.gov> menulis:
>>
>>> Those array summaries are helpful, thanks. It looks like the original
>>> program assumed that all three input files are on exactly the same grid.
>>> However, you can see that your U and V files have one extra longitude,
>>> compared to OLR. This caused the dimension error.
>>>
>>> You will need all three variables on the same grid because they are
>>> involved in common calculations, later in the program. Fortunately in this
>>> case, it is unusually easy to make U and V on (0...360) align to OLR on
>>> (0...357.5), because all the longitudes are exactly the same, except for
>>> the extra OLR longitude at 360. Simply remove the extra longitude by
>>> subsetting at the right places in the program.
>>>
>>> Here is one way to do it (untested). Use the reassignment operator
>>> *(:=)* to avoid conflicts with the preceding use of "work", then add a
>>> new line to subset the array to a conforming grid alignment. Repeat this
>>> for V as well as U.
>>>
>>> ncl 75> work = read_rename(f,"OLR_anom",iStrt,iLast,latS,latN) ;
>>> (time,lat,lon)
>>> ncl 76> OLR = dim_avg_n_Wrap(work, 1) ;
>>> (time,lon)
>>> ncl 77>
>>> dncl 78> f = addfile (diri+filu850 , "r")
>>>
>>> ncl 79> * work := read_rename(f,"U_anom",iStrt,iLast,latS,latN) ;
>>> (time,lat,lon)*
>>> ++++++ *work := work(:,:,{0:357.5})*
>>> ncl 80> U850 = dim_avg_n_Wrap(work, 1) ; (time,lon)
>>>
>>>
>>> On Mon, Dec 30, 2019 at 12:16 AM reza tisa <rezatisa1 at gmail.com> wrote:
>>>
>>>> This is what i got, but sorry i still dont get what should i do
>>>> Script mjoclivar_14:
>>>> https://drive.google.com/open?id=1UdrtSvNioUmbkL3xuwF_QZbRrkNlEfCr
>>>>
>>>> printVarSummary (work)
>>>> printVarSummary (read_rename(f,"U_anom",iStrt,iLast,latS,latN))
>>>>
>>>> Variable: work
>>>> Type: float
>>>> Total Size: 10939968 bytes
>>>> 2734992 values
>>>> Number of Dimensions: 3
>>>> Dimensions and sizes: [time | 1461] x [lat | 13] x [lon | 144]
>>>> Coordinates:
>>>> time: [1884648..1919688]
>>>> lat: [-15..15]
>>>> lon: [ 0..357.5]
>>>> Number Of Attributes: 3
>>>> units : W/m^2
>>>> long_name : Anomalies: Daily OLR
>>>> _FillValue : 32766
>>>>
>>>> Variable: unnamed (return)
>>>> Type: float
>>>> Total Size: 11015940 bytes
>>>> 2753985 values
>>>> Number of Dimensions: 3
>>>> Dimensions and sizes: [time | 1461] x [lat | 13] x [lon | 145]
>>>> Coordinates:
>>>> time: [ 0..2102400]
>>>> lat: [ -15.. 15]
>>>> lon: [ 0.. 360]
>>>> Number Of Attributes: 2
>>>> long_name : Anomalies from Daily Climatology
>>>> _FillValue : -9.99e+08
>>>> fatal:Dimension sizes of left hand side and right hand side of
>>>> assignment do not match
>>>> fatal:["Execute.c":8635]:Execute: Error occurred at or near line 82
>>>>
>>>>
>>>> ncl 74>
>>>> ncl 75> work = read_rename(f,"OLR_anom",iStrt,iLast,latS,latN)
>>>> ;(time,lat,lon)
>>>> ncl 76> OLR = dim_avg_n_Wrap(work, 1) ;
>>>> (time,lon)
>>>> ncl 77>
>>>> ncl 78> f = addfile (diri+filu850 , "r")
>>>> ncl 79> printVarSummary (work)
>>>> ncl 80> printVarSummary (read_rename(f,"U_anom",iStrt,iLast,latS,latN))
>>>> ncl 81>
>>>> ncl 82> work = read_rename(f,"U_anom",iStrt,iLast,latS,latN) ;
>>>> (time,lat,lon)
>>>> ncl 83> U850 = dim_avg_n_Wrap(work, 1) ; (time,lon)
>>>>
>>>>
>>>>
>>>> Pada tanggal Sen, 30 Des 2019 pukul 12.55 Dave Allured - NOAA Affiliate
>>>> <dave.allured at noaa.gov> menulis:
>>>>
>>>>> We need to show you how to debug, instead of just giving you the
>>>>> answer. Whenever you get this ...
>>>>>
>>>>> *"Dimension sizes of left hand side and right hand side of assignment
>>>>> do not match"*
>>>>>
>>>>> ... then you need to show printVarSummary of the left hand and right
>>>>> hand sides, NOT the original data. Here you would add two statements
>>>>> between lines 78 and 79, because the program breaks on line 79.
>>>>>
>>>>> printVarSummary (work)
>>>>> printVarSummary (read_rename(f,"U_anom",iStrt,iLast,latS,latN))
>>>>>
>>>>> Line 79 will work only if all dimensions are exactly the same size,
>>>>> between left and right. Now you will need to explain why one or more of
>>>>> the dimension sizes are different here.
>>>>>
>>>>> Note to others: There is more going on here than just "use the
>>>>> reassignment operator".
>>>>>
>>>>>
>>>>> On Sun, Dec 29, 2019 at 10:20 PM reza tisa via ncl-talk <
>>>>> ncl-talk at ucar.edu> wrote:
>>>>>
>>>>>> Hello NCL User, I really need help for doing mjoclivar_14.
>>>>>> I already success on doing mjoclivar_2 to get data anomalies for *OLR
>>>>>> and zonal wind* (uwnd 200 hPa & uwnd 850 hPa) for 2015-2018.
>>>>>> This is my data :
>>>>>> My mjoclivar_14 script :
>>>>>>
>>>>>> https://drive.google.com/open?id=0B6v1GPWVNcu7eGo0V0J0VTlpNmstYTFFVmJiVGpRTGdXdlBjOLR
>>>>>> anomalies
>>>>>> Original mjoclivar_14 script :
>>>>>> https://www.ncl.ucar.edu/Applications/Scripts/mjoclivar_14.ncl
>>>>>> OLR anomalies :
>>>>>> https://drive.google.com/open?id=1lgzfmmadfmXUWkUT1n8DAsTbvp0oQ-Fc
>>>>>> uwnd 200 anomalies :
>>>>>> https://drive.google.com/open?id=1Eb9a4VD96n8Wg9iDoRg7pOHr60tl9OWw
>>>>>> uwnd 850 anomalies :
>>>>>> https://drive.google.com/open?id=1PKpTYxGT5VAE6F5ydu0j89VHlRN7gC9p
>>>>>>
>>>>>> That anomalies data i use for running mjoclivar_14 . Then i got this
>>>>>> error.
>>>>>>
>>>>>> fatal:Dimension sizes of left hand side and right hand side of
>>>>>> assignment do not match
>>>>>> fatal:["Execute.c":8635]:Execute: Error occurred at or near *line 79*
>>>>>>
>>>>>> ....................
>>>>>> ;***********************************************************
>>>>>> ; Read anomalies
>>>>>> ;***********************************************************
>>>>>> ncl 75> work = read_rename(f,"OLR_anom",iStrt,iLast,latS,latN)
>>>>>> ; (time,lat,lon)
>>>>>> ncl 76> OLR = dim_avg_n_Wrap(work, 1)
>>>>>> ; (time,lon)
>>>>>> ncl 77>
>>>>>> dncl 78> f = addfile (diri+filu850 , "r")
>>>>>>
>>>>>> ncl 79> * work = read_rename(f,"U_anom",iStrt,iLast,latS,latN)
>>>>>> ; (time,lat,lon)*
>>>>>> ncl 80> U850 = dim_avg_n_Wrap(work, 1) ; (time,lon)
>>>>>> .....................
>>>>>>
>>>>>> This is my printVarSummary :
>>>>>> Variable *OLR_anom* :
>>>>>> ncl 5> printVarSummary(T)
>>>>>> Variable: T
>>>>>> Type: float
>>>>>> Total Size: 61432128 bytes
>>>>>> 15358032 values
>>>>>> Number of Dimensions: 3
>>>>>> Dimensions and sizes: [time | 1461] x [lat | 73] x [lon | 144]
>>>>>> Coordinates:
>>>>>> time: [1884648..1919688]
>>>>>> lat: [90..-90]
>>>>>> lon: [ 0..357.5]
>>>>>> Number Of Attributes: 3
>>>>>> _FillValue : 32766
>>>>>> long_name : Anomalies: Daily OLR
>>>>>> units : W/m^2
>>>>>>
>>>>>>
>>>>>> Variable *U_anom *(200 hPa) :
>>>>>> ncl 11> printVarSummary(T2)
>>>>>>
>>>>>> Variable: T2
>>>>>> Type: float
>>>>>> Total Size: 61858740 bytes
>>>>>> 15464685 values
>>>>>> Number of Dimensions: 3
>>>>>> Dimensions and sizes: [time | 1461] x [lat | 73] x [lon | 145]
>>>>>> Coordinates:
>>>>>> time: [ 0..2102400]
>>>>>> lat: [ -90.. 90]
>>>>>> lon: [ 0.. 360]
>>>>>> Number Of Attributes: 2
>>>>>> _FillValue : -9.99e+08
>>>>>> long_name : Anomalies from Daily Climatology
>>>>>>
>>>>>> Variable *U_anom* (850 hPa) :
>>>>>> ncl 15> printVarSummary(T3)
>>>>>>
>>>>>> Variable: T3
>>>>>> Type: float
>>>>>> Total Size: 61858740 bytes
>>>>>> 15464685 values
>>>>>> Number of Dimensions: 3
>>>>>> Dimensions and sizes: [time | 1461] x [lat | 73] x [lon | 145]
>>>>>> Coordinates:
>>>>>> time: [ 0..2102400]
>>>>>> lat: [ -90.. 90]
>>>>>> lon: [ 0.. 360]
>>>>>> Number Of Attributes: 2
>>>>>> _FillValue : -9.99e+08
>>>>>> long_name : Anomalies from Daily Climatology
>>>>>>
>>>>>> Huge thanks for anybody can make this work, really need help this for
>>>>>> my mini thesis. Thank you
>>>>>>
>>>>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191230/7fcc6b0e/attachment.html>
More information about the ncl-talk
mailing list