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

adriana.bossolasco at cnrs-orleans.fr adriana.bossolasco at cnrs-orleans.fr
Thu Aug 23 04:39:29 MDT 2018


Thanks for the correction

Now I get

abossola at olivia:~$ gcc --version
gcc (Ubuntu 7.3.0-16ubuntu3) 7.3.0
Copyright (C) 2017 Free Software Foundation, Inc.
This is free software; see the source for copying conditions.  There is 
NO
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR 
PURPOSE.

abossola at olivia:~$ uname -a
Linux olivia 4.15.0-32-generic #35-Ubuntu SMP Fri Aug 10 17:58:07 UTC 
2018 x86_64 x86_64 x86_64 GNU/Linux

I tryed the script that Denis send me and is the same I have the contour 
lines but not the color (see the attached graphic)
Le 23/08/2018 12:09, Guido Cioni a écrit :
> Be careful!
> 
>> abossola at olivia:~/Adriana/CESM/MAM7/MAM7_tes4$  gcc--version
>> gcc--version: command not found
> 
> It should be gcc --version (NOTE THE SPACE!!) and not gcc--version
> 
> You do the same error afterwards
> 
>> ncl 1> ncl-V
> 
> instead than ncl -V, although also this command should be executed
> oustide of NCL (only in the bash). The commands uname -a and gcc
> --version as well need to be executed outside of NCL.
> 
> I'd suggest you to have a general introduction on Linux and NCL before
> attempting to fix your plot.
> 
> However, did you try the script that Dennis sent?
> 
>> On 23. Aug 2018, at 11:12, adriana.bossolasco at cnrs-orleans.fr wrote:
>> 
>> Hi
>> 
>> I try the commands and this is what I obtain:
>> 
>> abossola at olivia:~/Adriana/CESM/MAM7/MAM7_tes4$  ncl -V
>> 6.4.0
>> abossola at olivia:~/Adriana/CESM/MAM7/MAM7_tes4$ uname -a
>> Linux olivia 4.15.0-32-generic #35-Ubuntu SMP Fri Aug 10 17:58:07
>> UTC 2018 x86_64 x86_64 x86_64 GNU/Linux
>> abossola at olivia:~/Adriana/CESM/MAM7/MAM7_tes4$  gcc--version
>> gcc--version: command not found
>> 
>> Other wise if I type this command in the ncl shell, then says:
>> 
>> ncl 1> ncl-V
>> fatal:syntax error: line 1 before or near -
>> ncl-
>> 
>> ncl 3> uname -a
>> fatal:syntax error: line 3 before or near -
>> uname -
>> ------^
>> 
>> fatal:error in statement
>> ncl 4>  gcc--version
>> fatal:syntax error: line 4 before or near -
>> gcc-
>> ----^
>> 
>> fatal:error in statement
>> 
>> I know that it is very strange because in my old computer everything
>> working, it could be due to the version of Linux?
>> 
>> This is my version now:
>> NAME="Ubuntu"
>> VERSION="18.04.1 LTS (Bionic Beaver)"
>> ID=ubuntu
>> ID_LIKE=debian
>> PRETTY_NAME="Ubuntu 18.04.1 LTS"
>> VERSION_ID="18.04"
>> HOME_URL="https://www.ubuntu.com/"
>> SUPPORT_URL="https://help.ubuntu.com/"
>> BUG_REPORT_URL="https://bugs.launchpad.net/ubuntu/"
>> 
> PRIVACY_POLICY_URL="https://www.ubuntu.com/legal/terms-and-policies/privacy-policy"
>> VERSION_CODENAME=bionic
>> UBUNTU_CODENAME=bionic
>> 
>> I really appreciate any help!
>> 
>> Le 22/08/2018 17:58, Dennis Shea a écrit :
>> Sorry. I accidentally hit "Send".
>> Attached is the 2nd plot.
>> On Wed, Aug 22, 2018 at 9:56 AM, Dennis Shea <shea at ucar.edu> wrote:
>> 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 [1]","r")
>> ;u = a->U(0,35,:,:)
>> a = addfile("uwnd.2008.nc [2]","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 [1] [3] 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 [1]
>> 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/ [4] [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 [5] [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 [1] [3] [8]
>> anonymous
>> your_email
>> cd incoming
>> put MAM_GEOS_1anee.cam.h0.2000-12-07-00000.nc [1] [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] [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 [2] [6] [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]
>> [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 [7] [2]
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk [7] [2]
>  Links:
> ------
> [1] http://MAM_GEOS_1anee.cam.h0.2000-12-07-00000.nc [3] [1]
> [2] http://mailman.ucar.edu/mailman/listinfo/ncl-talk [7]
> [3] http://www.eccforum.org
> [4] http://www.toni-klemm.de
> [5] http://twitter.com/toniklemm
> [6] http://www.ncl.ucar.edu/Document/Manuals/ [8]
> [7] http://www.ncl.ucar.edu/Applications/cylineq.shtml [9]
> [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
> Links:
> ------
> [1] http://MAM_GEOS_1anee.cam.h0.2000-12-07-00000.nc
> [2] http://uwnd.2008.nc
> [3] http://ftp.cgd.ucar.edu
> [4] HTTP://WWW.NCL.UCAR.EDU/DOCUMENT/MANUALS/
> [5] HTTP://WWW.NCL.UCAR.EDU/APPLICATIONS/CYLINEQ.SHTML
> [6] http://www.toni-klemm.de
> [7] http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> [8] http://www.ncl.ucar.edu/Document/Manuals/
> [9] http://www.ncl.ucar.edu/Applications/cylineq.shtml
> 
> --
> 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
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> 
> Guido Cioni
> http://guidocioni.altervista.org
> 
> 
> Links:
> ------
> [1] http://ftp.cgd.ucar.edu
> [2] http://www.toni-klemm.de
> [3] http://mam_geos_1anee.cam.h0.2000-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 --------------
A non-text attachment was scrubbed...
Name: adriana.000001.png
Type: image/png
Size: 388561 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180823/2ebddc13/attachment-0001.png>


More information about the ncl-talk mailing list