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