[ncl-talk] Missing color fill when using gsn_csm_contour_map (ncl V 6.4.0)

Dennis Shea shea at ucar.edu
Wed Aug 22 09:56:06 MDT 2018 

Not sure why you are not getting color.

If the following does not create color then please send the output from

%> ncl -V
%> uname -a
%> gcc--version

=========

;---read in zonal winds
 ;a = addfile("MAM_GEOS_1anee.cam.h0.2000-12-07-00000.nc","r")
 ;u = a->U(0,35,:,:)
  a = addfile("uwnd.2008.nc","r")
  u = a->uwnd(0,10,:,:)
  printVarSummary (u)
 ;print(u)

  wks  = gsn_open_wks("png","adriana")  ; send graphics to PNG file

  res                 = True
  res at gsnMaximize     = True       ; maximize in frame
  res at cnLinesOn       = True
  res at cnFillOn        = True
  res at cnFillPalette   = "BlueYellowRed"   ; change the color palette
  res at tiMainString    = "Default"
  plot = gsn_csm_contour_map(wks,u,res)


  res at mpFillDrawOrder = "PostDraw" ; Make sure land fill is drawn
                                   ; on top of filled contours

  res at tiMainString    = "add gray land fill after contouring"

  plot = gsn_csm_contour_map(wks,u,res)

On Wed, Aug 22, 2018 at 2:22 AM, <adriana.bossolasco at cnrs-orleans.fr> wrote:

> Hi thank you very much for the answer
> I corrected U and print the variable for see if this have values
>
> now is U(0,34,:,:)
>
> Variable: u
> Type: float
> Total Size: 55296 bytes
>             13824 values
> Number of Dimensions: 2
> Dimensions and sizes:    [lat | 96] x [lon | 144]
> Coordinates:
>             lat: [ -90..  90]
>             lon: [   0..357.5]
> Number Of Attributes: 6
>   lev :    581.2501460313797
>   time :     327
>   mdims :    1
>   units :    m/s
>   long_name :    Zonal wind
>   cell_methods :    time: mean
>
> And also have values, I checked it with print(u)
>
> But still the graphic is empty (white, no contour fill)
> On the other hand I can not access to ftp.cgd.ucar.edu for attach my
> file. Say:
> ERROR
> The requested URL could not be retrieved
> Any other idea how can I pass the file MAM_GEOS_1anee.cam.h0.2000-12-
> 07-00000.nc
>
> Thank you very much
>
> Le 21/08/2018 20:36, Dennis Shea a écrit :
>
>> [0]
>> If you are just learning NCL, I suggest you look at the USER MANUAL
>> at:
>>       HTTP://WWW.NCL.UCAR.EDU/DOCUMENT/MANUALS/ [6]
>>
>> [1]
>> As noted by Toni,
>>
>> u = U(0,,:,:)
>>
>> is not correct. In fact, it should have generated a FATAL subscript
>> syntax error.
>>
>> Note: although "fatal", NCL will continue. In fact, it will still try
>> to create the plot. Since it did not have valid data, NCL generated a
>> background and a default color lable bar. That is what you attached.
>>
>> [2]
>> Your following resource setting needs some clarification:
>>
>>   res at tiMainString    = "NCL version 6.1.x will give you gray land"
>>
>> For historical reasons, NCL's default graphic mode is black-and-white
>> (B/W). This is true of *all* NCL versions. Not just 6.1.x.
>>
>> When, 'gsn_csm_contour_map' is used, the default background is
>> continental outlines with light gray used to fill land areas. See:
>>
>>    HTTP://WWW.NCL.UCAR.EDU/APPLICATIONS/CYLINEQ.SHTML [7]
>>
>>    Examples 1 and 2
>>
>> Color is turned on via:
>>    res at cnFillOn = True     ; make color graphic
>>
>> When color is activated, (a) the background map including the light
>> gray land fil is generated; then, (b) the color contour are drawn over
>> that background.  The light gray land fill is overwritten (no big
>> deal). The following resource setting tells the 'gsn_csm_contour_map'
>> do the colr contours and then force the gray land fill last. [There
>> are some reasons for doing that.]
>>
>> res at mpFillDrawOrder = "PostDraw" ; Make sure land fill is drawn
>>                                    ; on top of filled contours
>>
>> [3]
>>
>> If you want you can ftp a sample file to:
>>
>>   ftp ftp.cgd.ucar.edu [8]
>>
>>   anonymous
>>
>>   your_email
>>
>>   cd incoming
>>
>>   put MAM_GEOS_1anee.cam.h0.2000-12-07-00000.nc [9]
>>
>>   quit
>>
>> On Tue, Aug 21, 2018 at 9:45 AM, Toni Klemm <toni-klemm at tamu.edu>
>> wrote:
>>
>> Hi Adriana,
>>>
>>> I noticed there are two commas in the command to extract the sind
>>> vectors: u = a->U(0,,:,:)  Could that be the problem?
>>>
>>> Could you attach the netCDF file from the example that you’re
>>> using? (MAM_GEOS_1anee.cam.h0.2000-12-07-00000.nc [1])
>>>
>>> Toni
>>>
>>> Toni Klemm, Ph.D.
>>> Postdoctoral Research Associate
>>> Department of Ecosystem Science and Management
>>> College of Agriculture and Life Sciences
>>> Texas A&M University, College Station, TX
>>> Contributor to the Early Career Climate Forum [3]
>>> www.toni-klemm.de [4] | @toniklemm [5]
>>>
>>> On Aug 21, 2018, at 10:01 AM, adriana.bossolasco at cnrs-orleans.fr
>>>> wrote:
>>>>
>>>> Hi all
>>>>
>>>> I have a problem with ncl when I try to plot a contour graphic. I
>>>> thought it was a mistake I was making, but then when I try with a
>>>> example very simple still appears blank. And ncl does not tell me
>>>> any error.
>>>> Here is the script (have my file but is exactly the same that
>>>> newcolor_1.ncl)
>>>>
>>>> ;************************************************
>>>> ; zonal_2.ncl
>>>> ;
>>>> ; Concepts illustrated:
>>>> ;   - Drawing a zonal means plot
>>>> ;   - Using dim_avg_Wrap to calculate a zonal average
>>>> ;************************************************
>>>> ;load "/usr/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>>> ;load "/usr/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>>> ;load "/usr/lib/ncarg/nclscripts/csm/contributed.ncl"
>>>> ;************************************************
>>>> begin
>>>> ;************************************************
>>>> ; variable and file handling
>>>> ;************************************************
>>>> ;---read in zonal winds
>>>> a = addfile("MAM_GEOS_1anee.cam.h0.2000-12-07-00000.nc [1]","r")
>>>>
>>>> u = a->U(0,,:,:)                                ; July zonal
>>>> winds
>>>>
>>>> wks  = gsn_open_wks("png","newcolor")  ; send graphics to PNG
>>>> file
>>>>
>>>> res                 = True
>>>> res at gsnMaximize     = True       ; maximize in frame
>>>> res at cnLinesOn       = False
>>>> res at cnFillOn        = True
>>>> res at cnFillPalette   = "BlueYellowRed"   ; change the color
>>>> palette
>>>>
>>>> res at mpFillDrawOrder = "PostDraw" ; Make sure land fill is drawn
>>>> ; on top of filled contours
>>>>
>>>> res at tiMainString    = "NCL version 6.1.x will give you gray land"
>>>>
>>>> plot = gsn_csm_contour_map(wks,u,res)
>>>> end
>>>>
>>>> I obtained a map without color fill, I was looking in some
>>>> previous talk but nothing works. Can Anyone help me on this?
>>>> thanks
>>>> I attached the maps
>>>>
>>>> Thank you
>>>> --
>>>> Adriana Bossolasco
>>>> LPC2E UMR 7328 (CNRS, Université d'Orléans)
>>>> 3A avenue de la Recherche Scientifique
>>>> 45071 Orléans Cedex 2 - FRANCE
>>>> Tel  +33 (0)2 38 25 52 85
>>>> e-mail:
>>>>
>>>>
>>> adriana.bossolasco at cnrs-orleans.fr<newcolor.png>____________
>> ___________________________________
>>
>>> 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://MAM_GEOS_1anee.cam.h0.2000-12-07-00000.nc
>> [2] http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>> [3] http://www.eccforum.org
>> [4] http://www.toni-klemm.de
>> [5] http://twitter.com/toniklemm
>> [6] http://www.ncl.ucar.edu/Document/Manuals/
>> [7] http://www.ncl.ucar.edu/Applications/cylineq.shtml
>> [8] http://ftp.cgd.ucar.edu
>> [9] http://12-07-00000.nc
>>
>
> --
> Adriana Bossolasco
> LPC2E UMR 7328 (CNRS, Université d'Orléans)
> 3A avenue de la Recherche Scientifique
> 45071 Orléans Cedex 2 - FRANCE
> Tel  +33 (0)2 38 25 52 85
> e-mail: adriana.bossolasco at cnrs-orleans.fr
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180822/ab33f81e/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: adriana.ncl
Type: application/octet-stream
Size: 797 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180822/ab33f81e/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: adriana.000001.png
Type: image/png
Size: 274140 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180822/ab33f81e/attachment.png>

More information about the ncl-talk mailing list