<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=NL-BE link="#0563C1" vlink="#954F72"><div class=WordSection1><p class=MsoNormal>Hello all,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><span lang=EN-GB>Lately I am experiencing some issues with plotting GFS data on charts. I used the cbarn.gs from the old opengrads version and the newer version as well as cbar cbarb cbarc etc… It all comes down to the same error: Cannot plot colorbar, no shading information<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>Once I start plotting a timeseries some of them go well and get the colorbar, others don’t. I have been searching for an explanation but couldn’t come up with one. I have it on the 0.25 0.50 and 1 degree GFS data. I have it for some charts, others don’t seem to have the issue. For example the Cape and Cinh charts do struggle with it, the temperature, dewpoint and low level shear don’t.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>An example is below. It is the script that is started through a .bat file that plots the cape 180mb data. Maybe someone of you guys might find the error? <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>**<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>function main(args)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>* Parse the arguments: date, hour<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB> if (args = '')<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB> prompt 'Enter forecast date (example, 20110717) --> ' <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB> pull date<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB> prompt 'Enter forecast hour (example, 00 or 06 or 12 or 18) --> ' <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB> pull hour<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB> else <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB> date = subwrd(args,1)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB> hour = subwrd(args,2)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB> endif<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>'reinit'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'sdfopen https://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs'date'/gfs_0p25_'hour'z'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>SAY 180Cape<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>* Get info from the descriptor file<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'q ctlinfo'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>_ctl = result<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>_undef = getctl(undef)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>_tdef = getctl(tdef)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>_zdef = getctl(zdef)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>maps = 66<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>* Get the Time axis info<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>tsize = subwrd(_tdef,2)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>_t1 = 1 ;* 2nd half of timeseries<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>_t2 = 45<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>tsize = _t2 - _t1 + 1<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set t '_t1' '_t2<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'q dims'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>times = sublin(result,5)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>_time1 = subwrd(times,6) <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>_time2 = subwrd(times,8)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>_tdim = _time1' '_time2<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>tincr = subwrd(_tdef,5)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>_tdef = 'tdef 'tsize' linear '_time1' 'tincr<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>huh = subwrd(_tdef,4)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>****************************************<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>* 180-0mb Caoe<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>****************************************<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB> i = 3<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB> while ( i<maps )<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>* Set up a few preliminary characteristics<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>setmap(1)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 16 0 100 254'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 17 0 150 254'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 18 0 200 254'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 19 0 230 240'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 20 0 230 160'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 21 0 230 120'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 21 0 230 80'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 22 0 240 40'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 23 0 250 0'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 24 254 254 0'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 25 254 225 0'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 26 254 200 0'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 27 254 175 0'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 28 254 150 0'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 29 230 125 0'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 30 230 100 0'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 31 220 75 30'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 32 200 50 30'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 33 180 25 30'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 34 170 0 30'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 35 180 0 50'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 36 200 0 100'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 37 254 0 150'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 38 254 0 200'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 39 254 0 254'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 40 254 90 254'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set rgb 41 255 150 255'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set t 'i<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set gxout shaded'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set csmooth on'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set clevs 0 200 400 600 800 1000 1200 1400 1600 1800 2000 2200 2400 2600 2800 3000 3200 3400 3600 3800 4000'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set ccols 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'d cape180_0mb'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>'cbarn'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>'query time'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>time = subwrd(result,3)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>day = subwrd(result,6)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>fcsttime = substr(time,1,13)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>fcstday = substr(day,1,3)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>'q dims'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>times = sublin(result,5)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>hub = subwrd(times,6) <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set strsiz 0.12'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set string 0 r 3 90' ; 'draw string 9.3 6 Cape in Joules/kg'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set string 0 r 3 0' ; 'draw string 8.7 8.4 Valid: 'fcstday' 'fcsttime<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set string 0 l 3 0' ; 'draw string 0.1 0.1 Data: NOAA NCEP GFS 0.25 - WWW.STORMCELL.BE'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set string 2 r 3 0' ; 'draw string 8.7 0.1 GFS Run: 'huh<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>* Draw Labels at the top of the page<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set string 0 r 11 0'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set strsiz 0.12'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'set string 4 l 3 0' ; 'draw string 0.1 8.4 180-0mb above ground cape'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>'gxprint bel_180cape'i'.png x1000 y844 -f overlay.png'<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>i = i+1<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>endwhile<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>* END OF MAIN SCRIPT<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB>**<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB style='mso-fareast-language:NL-BE'><o:p> </o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:NL-BE'>Met vriendelijke groeten,<o:p></o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:NL-BE'><o:p> </o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:NL-BE'>Tim Van Damme<o:p></o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:NL-BE'>+32 499 61 25 61<o:p></o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:NL-BE'><a href="mailto:tim@stormcell.be"><span style='color:#0563C1'>tim@stormcell.be</span></a><o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>