[ncl-talk] The result of the conditional expression yields a missing value (mjoclivar_15)

reza tisa rezatisa1 at gmail.com
Tue Dec 31 17:15:15 MST 2019


Thank You again Mr.Dave, It works.
I am following your suggesstion to took 12-months periods, and it looks
great.
Best Regards,
Reza


Pada tanggal Rab, 1 Jan 2020 pukul 02.16 Dave Allured - NOAA Affiliate <
dave.allured at noaa.gov> menulis:

> This error is more subtle.  I took a quick look at the code for
> mjoclivar_15.  It looks like it is currently set up for a maximum of 11 or
> 12 months, but you are trying to plot 24 months.  To begin, try reducing
> the time range.
>
> Later you could experiment with increasing certain arrays to accommodate a
> longer time range.  However, notice that there will be a trade-off.  The
> MJO plot will start to become overcrowded and hard to read.  What do you
> think would be better for your project?  A single dense plot, or separate
> plots for 12-month periods?
>
>
> On Tue, Dec 31, 2019 at 4:10 AM reza tisa <rezatisa1 at gmail.com> wrote:
>
>> Now i already adding at least 100 days before and after the date i want
>> to get (20160101-20171231),
>> i set the mjoclivar_14 like this
>> mjoclivar_14:
>>
>>
>> *   ymdStrt = 20150910                         ; start yyyymmdd   ymdLast
>> = 20180415                         ; last  *
>>
>> Result file :
>> MJO_PC_INDEX.nc *:   *
>> https://drive.google.com/open?id=1pr6wTdRXE0AWaPbh-Q8wWHca8NZLTLQw
>>
>> Then i set the time for mjoclivar_15
>> mjoclivar_15.ncl  (
>> https://drive.google.com/open?id=1TZO_iumtPSLuXITY8J4hRCKvfzBlzMrh):
>>
>> *  ymdStrt = 20160101                        ; start yyyymmdd  ymdLast =
>> 20171231                         ; last  *
>>
>> and i got this error for mjoclivar_15:
>> .........
>> fatal:Subscript out of range, error in subscript #0
>> fatal:["Execute.c":8635]:Execute: Error occurred at or near line 149
>> ..........
>> ncl 147>       if (iday(nt).eq.1) then
>> ncl 148>           nMon = nMon+1
>> *ncl 149>           monInfo(nMon,0) = namMonth(imon(nt)-1) *
>> ncl 150>           monInfo(nMon,1) = colMonth(imon(nt)-1)
>> ncl 151>       end if
>> ...........
>>
>> *This is my printVarSummary:*
>> *ncl 148> printVarSummary(namMonth(imon(nt)-1))*
>> Variable: namMonth (subsection)
>> Type: string
>> Total Size: 8 bytes
>>             1 values
>> Number of Dimensions: 1
>> Dimensions and sizes: [1]
>> Coordinates:
>>
>> *I also printVarSummary for pc1 and pc2 variable:*
>> *Variable: pc1*
>> Type: float
>> Total Size: 2924 bytes
>>             731 values
>> Number of Dimensions: 1
>> Dimensions and sizes: [time | 731]
>> Coordinates:
>>             time: [1893408..1910928]
>> Number Of Attributes: 5
>>   matrix : covariance
>>   ts_mean : ( 3.883529e-09, -1.027128e-09 )
>>   long_name : PC1
>>   info : PC1/stddev(PC1)
>>   _FillValue : 32766
>>
>> *Variable: pc2*
>> Type: float
>> Total Size: 2924 bytes
>>             731 values
>> Number of Dimensions: 1
>> Dimensions and sizes: [time | 731]
>> Coordinates:
>>             time: [1893408..1910928]
>> Number Of Attributes: 5
>>   matrix : covariance
>>   ts_mean : ( 3.883529e-09, -1.027128e-09 )
>>   long_name : PC2
>>   info : PC2/stddev(PC2)
>>   _FillValue : 32766
>> (0) PC1 : min=-2.94831   max=3.18312
>> (0) PC2 : min=-2.46549   max=2.53068
>> (0) Missing value count = 0
>> (0) Missing value count = 0
>>
>> Any advice for this error?
>>
>
> <snip>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200101/cbf7c170/attachment.html>


More information about the ncl-talk mailing list