
<html><head></head><body><div style="font-family: Verdana;font-size: 12.0px;"><div>Hello !!!</div>


<div>&nbsp;</div>


<div>I am plotting precipitation over Africa, I am using the world mask but I have the netcdf mask of the Congo and I would like to plot the area of this country shaded or just the contour in a different color. How could I do that ? Thank you !!!</div>


<div>&nbsp;</div>


<div>

<div>&#39;sdfopen H:&#92;precip.mon.mean.nc&#39;<br/>

&#39;sdfopen H:&#92;lsmask.nc&#39;<br/>

* Estes s&Atilde;&pound;o tons de AZUL<br/>

&#39;set rgb 16&nbsp;&nbsp; 0&nbsp;&nbsp; 0 255&#39;<br/>

&#39;set rgb 17&nbsp; 20&nbsp; 20 255&#39;<br/>

&#39;set rgb 18&nbsp; 55&nbsp; 55 255&#39;<br/>

&#39;set rgb 19 110 110 255&#39;<br/>

&#39;set rgb 20 165 165 255&#39;<br/>

&#39;set rgb 21 220 220 255&#39;<br/>

*<br/>

&#39;set mpdset lowres&#39;<br/>

&#39;set lat -30 30&#39;<br/>

&#39;set lon 320 410&#39;<br/>

&#39;set map 1 1 6&#39;<br/>

&#39;set clopts 1 1 0.12&#39;<br/>

&#39;set xlopts 1 2 0.14&#39;<br/>

&#39;set ylopts 1 2 0.14&#39;<br/>

&#39;set grid off&#39;<br/>

*<br/>

tempo=1<br/>

while(tempo&lt;13)<br/>

* clima 79 a 13 (jan79 t=1) (jan13 t=240)<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;define mediau=ave(precip.1,t=&#39;tempo&#39;,t=&#39;tempo+239&#39;,12)&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set grads off&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set gxout shaded&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set clevs 1 3 6 9 12&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set ccols 0 20 19 18 17 16&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;d mediau&#39;&nbsp; &nbsp;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;run cbarn.gs&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set gxout contour&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set grads off&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set cmax 12&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set cint 3&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set ccolor 1&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set cstyle 1&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set cthick 7&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set clab off&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;d mediau&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;sdfopen H:&#92;co_datos.nc&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; *If I difine here the new lat and long of co_datos then then print just the area selected and no the whole africa !!<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;d band1&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set clevs 54029&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;set gxout shaded&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;run&nbsp; G:&#92;cbarn.gs.gs&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;draw title Precipitation Rate (mm/day)&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;printim precip&#39;tempo&#39;.gif gif white&#39;<br/>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &#39;c&#39;<br/>

&nbsp;&nbsp;&nbsp; tempo=tempo+1<br/>

endwhile<br/>

&#39;quit&#39;</div>


<div>&nbsp;</div>

</div></div></body></html>