Fwd: Re: error at mjoclivar_15.ncl

I inadvertently did not include ncl-talk on my response.

On Mar 18, 2013, at 5:56 PM, Dennis Shea wrote:

> mjoclivar_14, the precursor to mjoclivar_15, applies a Lanczos
> filter to the data. This results in data at the beginning and
> end being set to _FillValue (missing)
>
> ;************************************************
> ; create BandPass Filter
> ;************************************************
> ihp = 2 ; bpf=>band pass filter
> nWgt = 201
> sigma = 1.0 ; Lanczos sigma
> fca = 1./100.
> fcb = 1./20.
> wgt = filwgts_lanczos (nWgt, ihp, fca, fcb, sigma )
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/filwgts_lanczos.shtml
>
> ====
>
> mjoclivar_15 wants a specific MJO event. This is a subset of your entire
> daily time series. It looks like you used the entire series
>
> ymdStrt = 19750101 ; start yyyymmdd
> ymdLast = 19821231 ; last
>
> You should select a subset period.
>
> Note that the example used the the 6-month period
>
> ymdStrt = 19961016 ; start yyyymmdd
> ymdLast = 19970415 ; last
>
>
>
>
>
> On 3/18/13 12:54 PM, Mary Haley wrote:
>> Marilla,
>>
>> There's something odd going on with the mjoclivar_14.ncl script, because it appears to be writing values equal to 32767 to the file. These values are apparently the missing value, which is odd, because the variable on the file should be float. 32767 is not a "normal" missing value for floats.
>>
>> Did you get errors when you ran the mjoclivar_14.ncl file?
>>
>> --Mary
>>
>> On Mar 18, 2013, at 12:02 PM, marilia.gregorio@cptec.inpe.br wrote:
>>
>>> Mary
>>>
>>> now I just have this error.
>>> The line 108 is the same as yours.
>>>
>>> Adding
>>> printMinMax (pc1, 0)
>>> printMinMax (pc2, 0)
>>> print ("pc1 (0) =" + pc1 (0))
>>> print ("pc2 (0) =" + pc2 (0))
>>>
>>> I got:
>>>
>>> (0) PC1: min=-3.14385 max=3.43989
>>> (0) PC2: min=-2.67514 max=2.90698
>>> (0) pc1 (0) =32766
>>> (0) pc2 (0) =32766
>>> (0) xBegin=32766
>>> (0) yBegin=32766
>>> (0) nt=1
>>> (0) pc1(nt)=32766
>>> (0) pc2(nt)=32766
>>> fatal:The result of the conditional expression yields a missing value. NCL can not determine branch, see ismissing function
>>> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 4195 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl
>>>
>>> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 4502 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl
>>>
>>> Thanks
>>>
>>> Marilia
>>>
>>>
>>>
>>>
>>>
>>> Quoting Mary Haley <haley@ucar.edu>:
>>>
>>>> Marilla,
>>>>
>>>> The error is coming from when mjoclivar_15.ncl script is trying to attach polylines to your plot.
>>>>
>>>> Are you getting any other errors?
>>>>
>>>> You need to do some debugging, because we don't have enough information to give you a fix.
>>>> The problem may be with the mjoclivar_14.ncl code (which writes the file you need to run this
>>>> script), or with mjoclivar_15.ncl itself.
>>>>
>>>> I think it is getting some missing values for whatever is being passed to gsn_add_polyline.
>>>>
>>>> On line 108 of the mjoclivar_15.ncl code I have:
>>>>
>>>> plot@$unique_string("dum")$ = gsn_add_polyline(wks, plot, (/xBegin,pc1(nt)/), (/yBegin,pc2(nt)/), plLine)
>>>>
>>>> Do you have this as well?
>>>>
>>>> If so, then you need to look at "xBegin", "pc1(nt)", "yBegin", and "pc2(nt)".
>>>>
>>>> These are all coming out of the "pc1" and "pc2" arrays, which are coming from the "MJO_PC_INDEX.nc" file which was written by mjoclivar_14.ncl.
>>>>
>>>> To debug this, after these two lines:
>>>>
>>>> pc1 = f->PC1(iStrt:iLast)
>>>> pc2 = f->PC2(iStrt:iLast)
>>>>
>>>> add:
>>>>
>>>> printMinMax(pc1,0)
>>>> printMinMax(pc2,0)
>>>> print("pc1(0) = " + pc1(0))
>>>> print("pc2(0) = " + pc2(0))
>>>>
>>>>
>>>> Inside the "do nt=1,ntim-1" do loop, does it fail on the first iteration, or later?
>>>> You may want to print "pc2(nt)" and "pc1(nt)" to see if those values are missing.
>>>>
>>>> --Mary
>>>>
>>>> On Mar 18, 2013, at 11:26 AM, marilia.gregorio@cptec.inpe.br wrote:
>>>>
>>>>> Hy...
>>>>>
>>>>> I'm trying to run the mjoclivar_15.ncl script. I have already ran the
>>>>> mjoclivar_14.ncl
>>>>> and its work very well.
>>>>>
>>>>> At mjoclivar_15.ncl I got the message:
>>>>>
>>>>> fatal:The result of the conditional expression yields a missing value.
>>>>> NCL can not
>>>>> determine branch, see ismissing function
>>>>> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 4195 in file
>>>>> $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl
>>>>>
>>>>> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 4502 in file
>>>>> $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl
>>>>>
>>>>> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 108 in file
>>>>> mjoclivar_15.ncl
>>>>>
>>>>>
>>>>> I really don't know what these means.
>>>>> I'm using Linux with NCL 6.1.2
>>>>>
>>>>> Thanks very much
>>>>>
>>>>> Marilia
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>
>>>
>>>
>>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Mar 18 18:18:17 2013