[ncl-talk] Segmentation fault

Karin Meier-Fleischer meier-fleischer at dkrz.de
Mon May 25 01:51:44 MDT 2015


Please take a look at the following .hlures web page: 

http://www.ncl.ucar.edu/Document/Graphics/hlures.shtml

Bye,
Karin


> Am 25.05.2015 um 07:23 schrieb Amadou Coulibaly <mpapin24 at gmail.com>:
> 
> Hi Dennis and Rick,
> 
> I have been able to plot one my wavelet spectrum using 5 years of hourly data, instead of all the 30 years, which were bringing the following error message: Segmentation fault (core dumped).
> 
> But I till have some warring and errors messages such as: 
> 
> warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for X axis failed: consider adjusting trXTensionF value
> warning:IrTransInitialize: error creating spline approximation for trXCoordPoints; defaulting to linear
> warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for X axis failed: consider adjusting trXTensionF value
> warning:IrTransInitialize: error creating spline approximation for trXCoordPoints; defaulting to linear
> warning:XyPlotSetValues:xyYStyle is NhlLOG:trYMinF can't be <= 0.0:Setting xyComputeYMin to True
> fatal:ContourPlotDraw: Workspace reallocation would exceed maximum size 100000000
> fatal:ContourPlotDraw: draw error
> fatal:ContourPlotDraw: draw error
> fatal:PlotManagerDraw: error in plot draw
> fatal:_NhlPlotManagerDraw: Draw error
> 
> I saw somewhere that to solve the first fatal (fatal:ContourPlotDraw: Workspace reallocation would exceed maximum size 100000000), I must change that maximum size 100000000 to  *wsMaximumSize : 500000000 in the file:.hluresfile and put it in my home directory. I don't understand how to put it in my home directory?
> 
> How to solve the issues of those errors?
> 
> Best
> 
>> On 23 May 2015 at 19:09, Dennis Shea <shea at ucar.edu> wrote:
>> [1] If your data has missing values, you can not use either specx_anal or the wavelet tool.
>> 
>> [2] http://www.ncl.ucar.edu/Document/Functions/Built-in/dtrend.shtml
>> 
>>      states: "Missing values (_FillValue) are not allowed."
>> 
>> [3] If there are **isolated** missing values you could likely do a simple linear interpolation
>>      to fill them. However, **any** interpolation will affect the spectra.
>> 
>>      If you have many _FillValue .... I would not even try to fill them in.
>>      The resulting spectra will reflect the effects of the interpolation as much as the 'rea;' data.
>> 
>> ---
>> Punch line: I don't think you can get what you want ... given your data.
>> 
>> --
>> 
>> Your statement:
>> 
>> "My script is from this link: 
>> http://paos.colorado.edu/research/wavelets/
>> 
>> There is no NCL script at that web site.
>> 
>> In fact, NCL call the fortran code presented at the website.
>> 
>> ---
>> The interactive code there is IDL ... not NCL
>> 
>> Regards
>>      
>>      
>> 
>>> On Sat, May 23, 2015 at 9:54 AM, Amadou Coulibaly <mpapin24 at gmail.com> wrote:
>>> Hi Dennis,
>>> 
>>> My script is from this link: 
>>> http://paos.colorado.edu/research/wavelets/
>>> 
>>> My Temp look like this:
>>> 
>>> Variable: Temp
>>> Type: float
>>> Total Size: 1051968 bytes
>>>             262992 values
>>> Number of Dimensions: 1
>>> Dimensions and sizes:    [time | 262992]
>>> Coordinates: 
>>>             time: [410227200..1356994800]
>>> Number Of Attributes: 7
>>>   coordinates :    lon lat
>>>   grid_type :    unstructured
>>>   comment :    Air temperature is the bulk temperature of the air, not the surf
>>> ace (skin) temperature.
>>>   shortname :    ta
>>>   code :      11
>>>   units :    K
>>>   _FillValue :    9.96921e+36
>>> (0)    299.15
>>> (1)    9.96921e+36
>>> (2)    9.96921e+36
>>> (3)    296.35
>>> 
>>> It contains a lot of missing data. That's why I used the following command:
>>> ; Remove mean and linear trend (recommended, not required)
>>> ; Missing values are not allowed. 
>>> ;************************************  
>>>  Temp = dtrend(Temp,False)
>>> 
>>> After that I printed Temp, but look stange for me with negative values:
>>> 
>>> Variable: Temp
>>> Type: float
>>> Total Size: 1051968 bytes
>>>             262992 values
>>> Number of Dimensions: 1
>>> Dimensions and sizes:    [time | 262992]
>>> Coordinates: 
>>>             time: [410227200..1356994800]
>>> Number Of Attributes: 7
>>>   coordinates :    lon lat
>>>   grid_type :    unstructured
>>>   comment :    Air temperature is the bulk temperature of the air, not the surf
>>> ace (skin) temperature.
>>>   shortname :    ta
>>>   code :      11
>>>   units :    K
>>>   _FillValue :    9.96921e+36
>>> (0)    -8.256236e+36
>>> (1)    1.713005e+36
>>> (2)    1.713035e+36
>>> (3)    -8.256145e+36
>>> 
>>> I think the problem is from that command trying to remove mean and linear trend. Because without that I can get the wavelet function (w), but it contains only zero as values.
>>> 
>>> How to handle this issue?
>>> 
>>> Best
>>> 
>>> 
>>>> On 23 May 2015 at 14:17, Dennis Shea <shea at ucar.edu> wrote:
>>>> [1] I have no idea why the script should seg fault. Maybe your parameter settings are not appropriate for the input data.
>>>> 
>>>> [2] What does "Temp" look like? printVarSummary(Temp)
>>>> 
>>>> [3] Are there missing values?
>>>> 
>>>> [4] If you want the "power spectrum", why not look at 
>>>> 
>>>> http://www.ncl.ucar.edu/Document/Functions/Built-in/specx_anal.shtml
>>>> http://www.ncl.ucar.edu/Applications/spec.shtml
>>>> 
>>>> Good luck
>>>> 
>>>>> On Sat, May 23, 2015 at 8:00 AM, Amadou Coulibaly <mpapin24 at gmail.com> wrote:
>>>>> Dear ncl users,
>>>>> 
>>>>> When I tried ti run my script (wavelet1.ncl), I am getting a strange error which is: 
>>>>> Segmentation fault (core dumped)
>>>>> 
>>>>> I am using ncl.6.2.1 and my system is: Linux amadou-Satellite-C855 3.19.0-16-generic #16-Ubuntu SMP Thu Apr 30 16:09:58 UTC 2015 x86_64 x86_64 x86_64 GNU/Linux
>>>>> 
>>>>> My dataset contains hourly data (1983-2012) of wind(speed/direction), air_temperature, ... for 46 sites. I just wanted to plot the power spectrum of only one site and for the all 30 years. I think the error comes from line 45 when I tried to compute the wavelet function (w). Please find attached my script.
>>>>> 
>>>>> How can I solve this issue?
>>>>> 
>>>>> 
>>>>> -- 
>>>>>  COULIBALY   AMADOU   
>>>>> PhD  Student  on  West  African  Climate  System (WACS)
>>>>> FUTA - Federal  University  of  Technology of Akure, Nigeria
>>>>> 
>>>>> Visiting Student - University of Cologne, Germany
>>>>> Institute of Geophysics and Meteorology
>>>>> Pohligstr. 3 / Office 3.102
>>>>> D-50969 Köln
>>>>> 
>>>>> Project: WASCAL (West African Science Service Centre on Climate Change and Adapted Land Use) 
>>>>> Phone:(+234) 810 795 2836 /(+223) 65 67 27 57 /(+226) 64 57 37 27 /+49 15218352574
>>>>> 
>>>>> E-mail: mpapin24 at gmail.com / coulibalya68 at yahoo.com 
>>>>> "The time is always right to do right": Nelson Mandela
>>>>> 
>>>>> "Character is like a tree and reputation like a shadow.  The shadow is what we think of it; the tree is the real thing" : Abraham Lincoln
>>>>> 
>>>>> "Do what you can, with what you have, where you are" Theodore Roosevelt
>>>>> 
>>>>>  
>>>>> 
>>>>>  
>>>>> 
>>>>>  
>>>>> 
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>> 
>>>> 
>>> 
>>> 
>>> 
>>> -- 
>>>  COULIBALY   AMADOU   
>>> PhD  Student  on  West  African  Climate  System (WACS)
>>> FUTA - Federal  University  of  Technology of Akure, Nigeria
>>> 
>>> Visiting Student - University of Cologne, Germany
>>> Institute of Geophysics and Meteorology
>>> Pohligstr. 3 / Office 3.102
>>> D-50969 Köln
>>> 
>>> Project: WASCAL (West African Science Service Centre on Climate Change and Adapted Land Use) 
>>> Phone:(+234) 810 795 2836 /(+223) 65 67 27 57 /(+226) 64 57 37 27 /+49 15218352574
>>> 
>>> E-mail: mpapin24 at gmail.com / coulibalya68 at yahoo.com 
>>> "The time is always right to do right": Nelson Mandela
>>> 
>>> "Character is like a tree and reputation like a shadow.  The shadow is what we think of it; the tree is the real thing" : Abraham Lincoln
>>> 
>>> "Do what you can, with what you have, where you are" Theodore Roosevelt
>>> 
>>>  
>>> 
>>>  
>>> 
>>>  
>> 
> 
> 
> 
> -- 
>  COULIBALY   AMADOU   
> PhD  Student  on  West  African  Climate  System (WACS)
> FUTA - Federal  University  of  Technology of Akure, Nigeria
> 
> Visiting Student - University of Cologne, Germany
> Institute of Geophysics and Meteorology
> Pohligstr. 3 / Office 3.102
> D-50969 Köln
> 
> Project: WASCAL (West African Science Service Centre on Climate Change and Adapted Land Use) 
> Phone:(+234) 810 795 2836 /(+223) 65 67 27 57 /(+226) 64 57 37 27 /+49 15218352574
> 
> E-mail: mpapin24 at gmail.com / coulibalya68 at yahoo.com 
> "The time is always right to do right": Nelson Mandela
> 
> "Character is like a tree and reputation like a shadow.  The shadow is what we think of it; the tree is the real thing" : Abraham Lincoln
> 
> "Do what you can, with what you have, where you are" Theodore Roosevelt
> 
>  
> 
>  
> 
>  
> <wavelet2.ncl>
> _______________________________________________
> ncl-talk mailing list
> 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/20150525/19ccbb8e/attachment.html 


More information about the ncl-talk mailing list