[ncl-talk] warning:cnConstFEnableFill to enable fill
Kunal Bali
kunal.bali9 at gmail.com
Thu Feb 22 07:29:33 MST 2018
please ignore my previous emails.
I have sorted out this issue. Now getting values between 0 to 1.
---
Kunal Bali
On Thu, Feb 22, 2018 at 6:06 PM, Kunal Bali <kunal.bali9 at gmail.com> wrote:
> One more thing I would like to add.
>
> As I mentioned in the earlier mail that I run the same script for both AOD
> and PM, which having same grid structure. In the script, I add the
> coordinated so that I can get the CoD values around that grid point.
> such as latz = (/28.545/)
> lonz = (/77.192/)
>
> But if you see the attached figure, right side showing CoD summary of AOD
> which extracting the required coordinated
>
> BUT
>
> the left side showing the CoD of PM which is extracting some other
> coordinated, not as mentioned in the script and like AOD.
> SO
>
> I got confused because both files have the same grid structure then the
> script should run the same way. Why the script is not running the same way
> in both cases.
>
> please provide some information on that.
>
>
>
> ---
> Kunal Bali
>
>
>
>
> On Thu, Feb 22, 2018 at 12:03 PM, Kunal Bali <kunal.bali9 at gmail.com>
> wrote:
>
>> Hello Ma'am,
>>
>> I have one doubt regarding your question.
>>
>> I run this script for both AOD and PM2.5 variable. The script is giving
>> CoD values between 0 to 1 for AOD but not for PM2.5. If I assume that there
>> is something wrong with the code then it should give some error in AOD data
>> as well.
>>
>> Both files have the same grid structure
>>
>> for AOD
>>
>> #
>> # gridID 1
>> #
>> gridtype = lonlat
>> gridsize = 5273664
>> xsize = 2497
>> ysize = 2112
>> xname = lon
>> xlongname = "longitude"
>> xunits = "degrees_east"
>> yname = lat
>> ylongname = "latitude"
>> yunits = "degrees_north"
>> xfirst = 59
>> xinc = 0.015625
>> yfirst = 5.0078125
>> yinc = 0.015625
>>
>>
>> for PM
>>
>>
>> #
>> # gridID 1
>> #
>> gridtype = lonlat
>> gridsize = 5273664
>> xsize = 2497
>> ysize = 2112
>> xname = lon
>> xlongname = "longitude"
>> xunits = "degrees_east"
>> yname = lat
>> ylongname = "latitude"
>> yunits = "degrees_north"
>> xfirst = 59
>> xinc = 0.015625
>> yfirst = 5.0078125
>> yinc = 0.015625
>>
>>
>>
>>
>>
>> ---
>> Kunal Bali
>>
>>
>>
>>
>>
>> On Wed, Feb 21, 2018 at 10:26 PM, Mary Haley <haley at ucar.edu> wrote:
>>
>>> Kunal,
>>>
>>> The issue is obviously coming somewhere from the calculations, and you
>>> are going to have to debug this code to figure out where the problem is.
>>>
>>> As you can see from the printMinMax output, The CoD1 variable has a
>>> min/max of 1, so taking the average of this is also going to give you a
>>> min/max of 1 (CoD2), and then doing a sqrt on this is also going to give
>>> you a min/max of 1 (CoD_final).
>>>
>>> You need to go back to the data_single, data_subset, and DATA_SINGLE_3D
>>> and make sure these are being calculated correctly. You should also print
>>> the min/max of these arrays to further pinpoint the problem.
>>>
>>> --Mary
>>>
>>>
>>> On Wed, Feb 21, 2018 at 4:49 AM, Kunal Bali <kunal.bali9 at gmail.com>
>>> wrote:
>>>
>>>> Dear Barry Sir,
>>>>
>>>> I set the levels via
>>>>
>>>> res at cnLevelSelectionMode = "ManualLevels"
>>>> res at cnMinLevelValF = 0.0
>>>> res at cnMaxLevelValF = 1
>>>> res at cnLevelSpacingF = 0.005
>>>>
>>>>
>>>> ---
>>>> Kunal Bali
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> On Wed, Feb 21, 2018 at 5:16 PM, Barry Lynn <barry.h.lynn at gmail.com>
>>>> wrote:
>>>>
>>>>> Hi:
>>>>>
>>>>> My apologies for asking simple question, but where do you set the
>>>>> levels you want to plot?
>>>>>
>>>>> I see that these two lines are commented out.
>>>>>
>>>>> ; res at cnLevelSelectionMode = "ExplicitLevels" ; set explicit
>>>>> contour levels
>>>>>
>>>>> ; res at cnLevels = (/0.04,0.06,0.08,0.1,0.12,0.14
>>>>> ,0.16,0.18,0.2,0.24,0.26,0.28,0.3,0.32,0.34,0.36,0.38,0.4/)
>>>>>
>>>>>
>>>>> On Wed, Feb 21, 2018 at 1:23 PM, Kunal Bali <kunal.bali9 at gmail.com>
>>>>> wrote:
>>>>>
>>>>>> Dear NCL users,
>>>>>>
>>>>>> I am trying to calculate the coefficient of divergence using the
>>>>>> attached script. The script is running successfully but I am getting only
>>>>>> one value i.e. 1 as a minimum and maximum value which should come in
>>>>>> between 0 to 1. And because of that I am also receiving the warning such as
>>>>>>
>>>>>> warning:ContourPlotInitialize: scalar field is constant; no contour
>>>>>> lines will appear; use cnConstFEnableFill to enable fill
>>>>>> if I set
>>>>>>
>>>>>> res at cnConstFEnableFill = True
>>>>>>
>>>>>> then the whole field is drawn in a single color as expected it will
>>>>>> give only one color.
>>>>>>
>>>>>>
>>>>>>
>>>>>> Regarding data
>>>>>> The data is in different 5 grid structure covering different Indian
>>>>>> region so I merged them and got one single file which cover whole Indian
>>>>>> region and then runs the attached script.
>>>>>>
>>>>>> Before merging the data, I used the same script for different grid
>>>>>> structures separately and it was giving the correct values between 0 to 1.
>>>>>>
>>>>>> But this time I am getting only one value.
>>>>>>
>>>>>> So could anyone please provide some information on that.
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> *Copyright (C) 1995-2017 - All Rights Reserved University Corporation
>>>>>> for Atmospheric Research NCAR Command Language Version 6.4.0 The use of
>>>>>> this software is governed by a License Agreement. See
>>>>>> http://www.ncl.ucar.edu/ <http://www.ncl.ucar.edu/> for more
>>>>>> details.Variable: totcType: doubleTotal Size: 4050173952
>>>>>> <040%205017%203952> bytes 506271744 valuesNumber of Dimensions:
>>>>>> 3Dimensions and sizes: [time | 96] x [lat | 2112] x [lon |
>>>>>> 2497]Coordinates: time: [35430.125..104813.7916666667
>>>>>> <079%201666%206667>] lat: [5.0078125..37.9921875]
>>>>>> lon: [ 59.. 98]Number Of Attributes: 6 _FillValue : -28672 long_name
>>>>>> : AOT at 0.55 micron hdfeos_name : Optical_Depth_055 projection
>>>>>> : Albers Conical Equal_Area unit : None missing_value :
>>>>>> -28672Variable: totcType: doubleTotal Size: 4050173952 <040%205017%203952>
>>>>>> bytes 506271744 valuesNumber of Dimensions: 3Dimensions and
>>>>>> sizes: [time | 96] x [lat | 2112] x [lon | 2497]Coordinates:
>>>>>> time: [35430.125..104813.7916666667 <079%201666%206667>] lat:
>>>>>> [5.0078125..37.9921875] lon: [ 59.. 98]Number Of Attributes:
>>>>>> 6 _FillValue : -28672 long_name : AOT at 0.55 micron hdfeos_name
>>>>>> : Optical_Depth_055 projection : Albers Conical Equal_Area unit
>>>>>> : None missing_value : -28672Variable: data_subsetType: doubleTotal
>>>>>> Size: 295196160 bytes 36899520 valuesNumber of Dimensions:
>>>>>> 3Dimensions and sizes: [time | 96] x [lat | 646] x [lon |
>>>>>> 595]Coordinates: time: [35430.125..104813.7916666667
>>>>>> <079%201666%206667>] lat: [29.5234375..19.4453125]
>>>>>> lon: [91.390625..82.109375]Number Of Attributes: 6 missing_value :
>>>>>> -28672 unit : None projection : Albers Conical Equal_Area
>>>>>> hdfeos_name : Optical_Depth_055 long_name : AOT at 0.55 micron
>>>>>> _FillValue : -28672(0) nearest location: k=0 n=481 m=882 156.836
>>>>>> 5297.11(0) -----(0) -----Variable: CoD1Type: doubleTotal Size:
>>>>>> 295196160 bytes 36899520 valuesNumber of Dimensions: 3Dimensions
>>>>>> and sizes: [96] x [646] x [595]Coordinates: Number Of Attributes: 1
>>>>>> _FillValue : -28672(0) min=1 max=1(0) ---Variable: CoD2Type:
>>>>>> doubleTotal Size: 3074960 bytes 384370 valuesNumber of
>>>>>> Dimensions: 2Dimensions and sizes: [646] x [595]Coordinates: Number Of
>>>>>> Attributes: 3 long_name : Coef of Divergence _FillValue : -28672
>>>>>> average_op_ncl : dim_avg_n function was applied(0) Coef of Divergence
>>>>>> : min=1 max=1Variable: CoD_finalType: doubleTotal Size: 3074960
>>>>>> bytes 384370 valuesNumber of Dimensions: 2Dimensions and
>>>>>> sizes: [646] x [595]Coordinates: Number Of Attributes: 1 _FillValue
>>>>>> : -28672(0) min=1 max=1warning:ContourPlotInitialize: scalar field
>>>>>> is constant; no contour lines will appear; use cnConstFEnableFill to enable
>>>>>> fill*
>>>>>>
>>>>>> ---
>>>>>> Kunal
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> _______________________________________________
>>>>>> ncl-talk mailing list
>>>>>> ncl-talk at ucar.edu
>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> Barry H. Lynn, Ph.D
>>>>> Senior Lecturer,
>>>>> The Institute of the Earth Science,
>>>>> The Hebrew University of Jerusalem,
>>>>> Givat Ram, Jerusalem 91904, Israel
>>>>> Tel: 972 547 231 170
>>>>> Fax: (972)-25662581
>>>>>
>>>>> C.E.O, Weather It Is, LTD
>>>>> Weather and Climate Focus
>>>>> http://weather-it-is.com
>>>>> Jerusalem, Israel
>>>>> Local: 02 930 9525
>>>>> Cell: 054 7 231 170
>>>>> Int-IS: x972 2 930 9525
>>>>> US 914 432 3108 <091443%2023108>
>>>>>
>>>>
>>>>
>>>> _______________________________________________
>>>> 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/20180222/01ac750f/attachment.html>
More information about the ncl-talk
mailing list