
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr">

<div id="divtagdefaultwrapper" style="font-size:10pt;color:#2672EC;background-color:#FFFFFF;font-family:Arial,Helvetica,sans-serif;">

<p></p>

<div>Hi all,<br>

<br>

I'm trying to plot a text file containing three columns lon , lat, value using ncl.<br>

<br>

First I converted my text file to nc using CDO. The nc file is like:<br>

<br>

##################################################################################<br>

<br>

<div>Variable: f<br>

Type: file<br>

filename:&nbsp;&nbsp; &nbsp;tir1<br>

path:&nbsp;&nbsp; &nbsp;tir1.nc<br>

&nbsp;&nbsp; file global attributes:<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; CDI : Climate Data Interface version ?? (http://mpimet.mpg.de/cdi)<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Conventions : CF-1.4<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; history : Tue Sep 06 11:56:42 2016: cdo -f nc input,my_gridfile tir1.nc<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; CDO : Climate Data Operators version 1.7.0 (http://mpimet.mpg.de/cdo)<br>

&nbsp;&nbsp; dimensions:<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ncells = 53248<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time = 1&nbsp; // unlimited<br>

&nbsp;&nbsp; variables:<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; double lon ( ncells )<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; standard_name :&nbsp;&nbsp; &nbsp;longitude<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; units :&nbsp;&nbsp; &nbsp;degrees_east<br>

<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; double lat ( ncells )<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; standard_name :&nbsp;&nbsp; &nbsp;latitude<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; units :&nbsp;&nbsp; &nbsp;degrees_north<br>

<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; double time ( time )<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; standard_name :&nbsp;&nbsp; &nbsp;time<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; units :&nbsp;&nbsp; &nbsp;hours since 0-0-0 00:00:00<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calendar :&nbsp;&nbsp; &nbsp;proleptic_gregorian<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; axis :&nbsp;&nbsp; &nbsp;T<br>

<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; float var1 ( time, ncells )<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; grid_type :&nbsp;&nbsp; &nbsp;unstructured<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; coordinates :&nbsp;&nbsp; &nbsp;lat lon<br>

###################################################################################<br>

<br>

</div>

but I'm getting following errors:<br>

<br>

<span style="color: rgb(255, 0, 0);">fatal:ScalarFieldInitialize: Insufficient number of elements in sfDataArray</span><br>

<span style="color: rgb(255, 0, 0);">fatal:Unable to initialize layer-Can't Create</span><br>

<span style="color: rgb(255, 0, 0);">fatal:Unable to access object with id:-4</span><br>

<span style="color: rgb(255, 0, 0);">fatal:PID #-4 can't be found in NhlSetValues</span><br>

<span style="color: rgb(255, 0, 0);">fatal:Unable to access object with id:-4</span><br>

<span style="color: rgb(255, 0, 0);">warning:Unable to add DataItem &quot;(null)&quot; to DataList &quot;cnScalarFieldData&quot;</span><br>

<span style="color: rgb(255, 0, 0);">fatal:CompileDataList:DataList has no valid members</span><br>

</div>

<p></p>

<p><br>

</p>

<p><br>

</p>

<p>######################################################################################</p>

<p><br>

</p>

<p>my script is :</p>

<p><br>

</p>

<p></p>

<div>load &quot;/usr/share/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>

load &quot;/usr/share/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br>

load &quot;/usr/share/ncarg/nclscripts/csm/contributed.ncl&quot;<br>

<br>

<br>

<br>

begin<br>

;---read in zonal winds<br>

&nbsp; a &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= addfile(&quot;/home/amit/Desktop/tir1.nc&quot;,&quot;r&quot;)<br>

<br>

&nbsp; filename = &quot;/media/amit/E/India_shapefile/State_bound(Satmet).shp&quot;<br>

<br>

&nbsp; ;filename = &quot;/media/amit/E/India_shapefile/Districtbound(Satmet).shp&quot;<br>

<br>

<br>

&nbsp; var &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= doubletofloat(a-&gt;var1)<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br>

lat &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;=&nbsp; a-&gt;lat<br>

lon &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;=&nbsp; a-&gt;lon<br>

<br>

<br>

&nbsp; wks = gsn_open_wks(&quot;png&quot;,&quot;TIR1&quot;)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; send graphics to PNG file<br>

<br>

&nbsp; res&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; plot mods desired<br>

&nbsp; res@cnFillOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; turn on color fill<br>

&nbsp; res@cnLinesOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>

&nbsp; res@gsnMaximize&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp; ; maximize plot in frame<br>

&nbsp; res@gsnDraw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>

&nbsp; res@gsnFrame&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>

<br>

&nbsp; res@mpOutlineOn&nbsp;&nbsp; = False&nbsp;&nbsp; ; Use outlines from shapefile<br>

&nbsp; res@mpFillOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False &nbsp;<br>

&nbsp; res@cnFillMode = &quot;RasterFill&quot;<br>

&nbsp; res@trGridType = &quot;TriangularMesh&quot;<br>

&nbsp; res@cnRasterSmoothingOn = False <br>

<br>

res@cnLevelSelectionMode = &quot;ExplicitLevels&quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br>

&nbsp; res@cnLevels&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = (/190,220,240,250,260,265,270,275,280,285,290,300,320,360/) ;&nbsp;&nbsp; ; 14 contour values<br>

&nbsp; res@cnFillPalette&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = (/&quot;Grey&quot;,&quot;BlueViolet&quot;,&quot;Blue4&quot;,&quot;Blue2&quot;,&quot;Cornflowerblue&quot;,&quot;darkslategray1&quot;,&quot;Cyan&quot;\<br>

&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;,&quot;chartreuse&quot;,&quot;khaki1&quot;,&quot;Yellow&quot;\ &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;; 15 contour colors<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ,&quot;Orange&quot;,&quot;darkorange4&quot;, &quot;hotpink&quot;&nbsp;&nbsp;&nbsp; \<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ,&quot;Red&quot;,&quot;Black&quot; /)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br>

<br>

res@tiMainString&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;TIR1 Plot&quot;<br>

&nbsp; res@tiMainFontHeightF = 0.016 <br>

&nbsp; res@gsnAddCyclic&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>

&nbsp; res@lbOrientation&nbsp;&nbsp; = &quot;Vertical&quot;<br>

<br>

&nbsp; res@sfXArray &nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = lon<br>

&nbsp; res@sfYArray&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;=&nbsp; lat<br>

<br>

<br>

res@mpMinLonF &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= 68.0<br>

res@mpMaxLonF &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= 98.0<br>

res@mpMinLatF &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= 0.0<br>

res@mpMaxLatF &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= 38.0<br>

<br>

res@pmTickMarkDisplayMode &nbsp;&nbsp; &nbsp;= &quot;Always&quot;<br>

<br>

res@tiXAxisString&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= &quot;Longitude&quot;<br>

res@tiYAxisString&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= &quot;Lattitude&quot;<br>

<br>

res@tiXAxisFontHeightF &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= 0.02<br>

res@tiYAxisFontHeightF &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= 0.02 <br>

<br>

res@gsnStringFontHeightF &nbsp;&nbsp; &nbsp;= 0.005 &nbsp;<br>

<br>

&nbsp; plot = gsn_csm_contour_map(wks,var, res)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; create plot<br>

<br>

;---Resource list for shapefile outlines<br>

&nbsp;plres&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<br>

&nbsp;plres@gsLineColor&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;black&quot;<br>

&nbsp;plres@gsLineThicknessF&nbsp; = &quot;1.5&quot;<br>

&nbsp;dumstr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; = unique_string(&quot;poly&quot;)<br>

&nbsp; <br>

&nbsp;id = gsn_add_shapefile_polylines(wks,plot,filename,plres)<br>

<br>

&nbsp;draw(plot)<br>

&nbsp; frame(wks)<br>

<br>

end</div>

####################################################################################################################<br>

<p></p>

<p><br>

</p>

</div>

</body>

</html>