[ncl-talk] Errors while trying to plot IAGOS data. Only the contour map plot, but the data is missing.
najib.yusuf at carnasrda.com
najib.yusuf at carnasrda.com
Wed Sep 27 12:22:54 MDT 2017
On 2017-09-27 13:41, Dennis Shea wrote:
> You stated: "I tried solving but couldn't get it"
>
> Please look at the output from your commands.
>
> The output is for *you* to look at.
>
> ----
>
> More on David's comment ...
>
> re: GOLDEN RULE OF DATA PROCESSING... LOOK AT YOUR DATA
>
> You used:
>
> var = f->O3_PM(0) <==== the reads only the 1st value [index (0) ]
>
> Variable: var
> Type: double
> Total Size: 8 bytes
> 1 values <=== ONE VALUE ... how do you contour
> one value?
>
> ---
>
> The min & max are _FillValue. Look at the output from your script.
>
> (0) min/max var = -9999/-9999
> ---
>
> Even if you read the entire *time series* ... all values are
> _FillValue
>
> var = f->O3_PM ; entire series
>
> printVarSummary(var)
>
> Variable: var
> Type: double
> Total Size: 63024 bytes
> 7878 values
> Number of Dimensions: 1
> Dimensions and sizes: [UTC_time | 7878] <=========
> Coordinates:
> UTC_time: [70841..106359]
> Number Of Attributes: 8
> standard_name : mole_fraction_of_ozone_in_air
> long_name : Ozone mixing ratio
> comment : Measured by MOZAIC package
> units : 0.000000001
> missing_value : -9999
> _FillValue : -9999
> scale_factor : 1
> availability : 0
>
> print ("min/max var = " + min(var) + "/" + max(var))
> (0) MIN/MAX VAR = -9999/-9999
>
> ---
>
> It should be plotted as a time series not as a contour plot.
>
> On Wed, Sep 27, 2017 at 10:34 AM, David Brown <dbrown at ucar.edu> wrote:
>
>> You will never be able to contour the O3_PM variable from this file
>> because it contains only fill values:
>>
>> ncl 0> f = addfile("IAGOS_timeseries_2013101819404151.nc [1]","r")
>> ncl 1> printMinMax(f->O3_PM,0)
>> (0) Ozone mixing ratio (0.000000001) : min=-9999 max=-9999
>>
>> Presumably this is what the attribute "availability" means: no valid
>> values for this variable.
>>
>> Also for some reason your variable "var" contains only a single
>> value.
>> There is no way to make a contour plot with only one point even if
>> it
>> were a valid value. Since the only dimension in the file is a time
>> dimension, maybe this data set would be better displayed as a time
>> series x-y plot.
>> -dave
>>
>> On Wed, Sep 27, 2017 at 8:16 AM, <najib.yusuf at carnasrda.com> wrote:
>>>
>>> Dear Sirs,
>>>
>>> I try to plot contour map of O3_PM variable from .nc IAGOS file
>> but I am
>>> getting strange errors, I tried solving but couldn't get it.
>> Kindly assist
>>> to guide. Attached is the file and the IAGOS.ncl as script.Thank
>> you.
>>>
>>> The summary of the print var is;
>>>
>>> Variable: var
>>> Type: double
>>> Total Size: 8 bytes
>>> 1 values
>>> Number of Dimensions: 1
>>> Dimensions and sizes: [UTC_time | 1]
>>> Coordinates:
>>> Number Of Attributes: 9
>>> UTC_time : 70841
>>> standard_name : mole_fraction_of_ozone_in_air
>>> long_name : Ozone mixing ratio
>>> comment : Measured by MOZAIC package
>>> units : 0.000000001
>>> missing_value : -9999
>>> _FillValue : -9999
>>> scale_factor : 1
>>> availability : 0
>>>
>>>
>>> Variable: var
>>> Type: double
>>> Total Size: 8 bytes
>>> 1 values
>>> Number of Dimensions: 1
>>> Dimensions and sizes: [UTC_time | 1]
>>> Coordinates:
>>> Number Of Attributes: 9
>>> UTC_time : 70841
>>> standard_name : mole_fraction_of_ozone_in_air
>>> long_name : Ozone mixing ratio
>>> comment : Measured by MOZAIC package
>>> units : 0.000000001
>>> missing_value : -9999
>>> _FillValue : -9999
>>> scale_factor : 1
>>> availability : 0
>>> (0) -9999
>>> (0) min/max var = -9999/-9999
>>> fatal:MeshScalarFieldInitialize: Insufficient number of elements
>> in
>>> sfDataArray
>>> fatal:Unable to initialize layer-Can't Create
>>> fatal:Unable to access object with id:-4
>>> fatal:Unable to access object with id:-4
>>> warning:Unable to add DataItem "(null)" to DataList
>> "cnScalarFieldData"
>>> fatal:CompileDataList:DataList has no valid members
>>> warning:ContourPlotSetValues: cnLevelSpacingF value causes level
>> count to
>>> exceed maximum: using AUTOMATICLEVELS mode
>>> ncl 48>
>>>
>>> Thank you
>>>
>>> Najib
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk [2]
>>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk [2]
>
>
>
> Links:
> ------
> [1] http://2013101819404151.nc
> [2] http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Thank you sir David and Dennis.
Got it. Thank you
Najib
More information about the ncl-talk
mailing list