[ncl-talk] warning:cnConstFEnableFill to enable fill
Kunal Bali
kunal.bali9 at gmail.com
Thu Feb 22 05:36:13 MST 2018
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/00f699d9/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2018-02-22 at 5.36.57 PM.png
Type: image/png
Size: 622851 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180222/00f699d9/attachment-0001.png>
More information about the ncl-talk
mailing list