<html><head></head><body><div style="color:#000; background-color:#fff; font-family:verdana, helvetica, sans-serif;font-size:16px"><div id="yui_3_16_0_ym19_1_1489139220873_82512"><span><br></span></div><div id="yui_3_16_0_ym19_1_1489139220873_82628">Dear sir/madam</div><div id="yui_3_16_0_ym19_1_1489139220873_82629"><br id="yui_3_16_0_ym19_1_1489139220873_82630"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82631">I
am new to the ncl graphics and I'm trying to plot a field from the
ascii file attached. As you may see it is a 3-column file
(lon,lat,indicator). I have used the file plotAscii.ncl given below to
create the plot and receive the following error:</div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82632"><br id="yui_3_16_0_ym19_1_1489139220873_82633"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82634">-----------------------</div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82635">fatal:Cannot create triangular mesh: supply additional resources or build with Triangle package<br id="yui_3_16_0_ym19_1_1489139220873_82636">fatal:ContourPlotDraw: draw error<br id="yui_3_16_0_ym19_1_1489139220873_82637">fatal:ContourPlotDraw: draw error<br id="yui_3_16_0_ym19_1_1489139220873_82638">fatal:PlotManagerDraw: error in plot draw<br id="yui_3_16_0_ym19_1_1489139220873_82639">fatal:_NhlPlotManagerDraw: Draw error<br id="yui_3_16_0_ym19_1_1489139220873_82640"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82641">-----------------------<br id="yui_3_16_0_ym19_1_1489139220873_82642"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82643"><br id="yui_3_16_0_ym19_1_1489139220873_82644"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82645">I only get a png image with the map of the corresponding area in low res (as it's the default).<br id="yui_3_16_0_ym19_1_1489139220873_82646"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82647">Could you please help me with the script?</div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82648">Also, if I want to use highres for the coeastlines in which directory should I add the GSHHS file?</div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82649"><br id="yui_3_16_0_ym19_1_1489139220873_82650"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82651">The version of ncl installed (already built from ncl homepage) is 6.3.0 .<br id="yui_3_16_0_ym19_1_1489139220873_82652"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82653"><br id="yui_3_16_0_ym19_1_1489139220873_82654"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82655">Thank you for any help.</div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82656"><br id="yui_3_16_0_ym19_1_1489139220873_82657"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82658">Best regards,<br id="yui_3_16_0_ym19_1_1489139220873_82659"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82660">Petroula<br id="yui_3_16_0_ym19_1_1489139220873_82661"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82662"><br id="yui_3_16_0_ym19_1_1489139220873_82663"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82664"><br id="yui_3_16_0_ym19_1_1489139220873_82665"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82666">-----------------------</div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82667">plotAscii.ncl <br id="yui_3_16_0_ym19_1_1489139220873_82668"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82669">-----------------------</div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82670"><br id="yui_3_16_0_ym19_1_1489139220873_82671"></div>load "~/WRF/LIBRARIES/ncl_ncarg-6.3.0/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br id="yui_3_16_0_ym19_1_1489139220873_82672">load "~/WRF/LIBRARIES/ncl_ncarg-6.3.0/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br id="yui_3_16_0_ym19_1_1489139220873_82673"><br id="yui_3_16_0_ym19_1_1489139220873_82674">begin<br id="yui_3_16_0_ym19_1_1489139220873_82675"><br id="yui_3_16_0_ym19_1_1489139220873_82676"> fname = "~/FlyingRoute/OUTPUT/storm_2014071412_12.dat"<br id="yui_3_16_0_ym19_1_1489139220873_82677"> n_col = numAsciiCol(fname)<br id="yui_3_16_0_ym19_1_1489139220873_82678"> n_var = n_col<br id="yui_3_16_0_ym19_1_1489139220873_82679"> n_pt = numAsciiRow(fname)<br id="yui_3_16_0_ym19_1_1489139220873_82680"> dataf = asciiread(fname,(/n_pt,n_var/),"float")<br id="yui_3_16_0_ym19_1_1489139220873_82681">;---read longitude<br id="yui_3_16_0_ym19_1_1489139220873_82682"> lon = dataf (:,0)<br id="yui_3_16_0_ym19_1_1489139220873_82683">;---read latitude<br id="yui_3_16_0_ym19_1_1489139220873_82684"> lat = dataf (:,1)<br id="yui_3_16_0_ym19_1_1489139220873_82685">;---read data for presence of convective cloud<br id="yui_3_16_0_ym19_1_1489139220873_82686"> storm = dataf (:,2)<br id="yui_3_16_0_ym19_1_1489139220873_82687">;---define the workstation (plot type and name)<br id="yui_3_16_0_ym19_1_1489139220873_82688"> wks = gsn_open_wks("png","plot_field")<br id="yui_3_16_0_ym19_1_1489139220873_82689">;--set resources<br id="yui_3_16_0_ym19_1_1489139220873_82690"> res = True<br id="yui_3_16_0_ym19_1_1489139220873_82691"> res@sfXArray = lon<br id="yui_3_16_0_ym19_1_1489139220873_82692"> res@sfYArray = lat<br id="yui_3_16_0_ym19_1_1489139220873_82693"> res@cnFillOn = True ; Turn on contour fill.<br id="yui_3_16_0_ym19_1_1489139220873_82694"> res@cnLinesOn = False ; Turn off contour lines<br id="yui_3_16_0_ym19_1_1489139220873_82695">;-- First create just a contour plot in order to retrieve<br id="yui_3_16_0_ym19_1_1489139220873_82696">; --the start/end limits of the actual location of the data.<br id="yui_3_16_0_ym19_1_1489139220873_82697">; --Don't draw the contour plot.<br id="yui_3_16_0_ym19_1_1489139220873_82698"> res@gsnDraw = False<br id="yui_3_16_0_ym19_1_1489139220873_82699"> res@gsnFrame = False<br id="yui_3_16_0_ym19_1_1489139220873_82700"> contour = gsn_csm_contour(wks,storm,res)<br id="yui_3_16_0_ym19_1_1489139220873_82701">;-- gsnDraw and gsnFrame are on for the next plot.<br id="yui_3_16_0_ym19_1_1489139220873_82702"> res@gsnDraw = True<br id="yui_3_16_0_ym19_1_1489139220873_82703"> res@gsnFrame = True<br id="yui_3_16_0_ym19_1_1489139220873_82704">;-- Retrieve the actual start and end limits of the data<br id="yui_3_16_0_ym19_1_1489139220873_82705"> getvalues contour@data<br id="yui_3_16_0_ym19_1_1489139220873_82706"> "sfXCActualStartF" : xs<br id="yui_3_16_0_ym19_1_1489139220873_82707"> "sfXCActualEndF" : xe<br id="yui_3_16_0_ym19_1_1489139220873_82708"> "sfYCActualStartF" : ys<br id="yui_3_16_0_ym19_1_1489139220873_82709"> "sfYCActualEndF" : ye<br id="yui_3_16_0_ym19_1_1489139220873_82710"> end getvalues<br id="yui_3_16_0_ym19_1_1489139220873_82711"> <br id="yui_3_16_0_ym19_1_1489139220873_82712"> res@mpLimitMode = "LatLon"<br id="yui_3_16_0_ym19_1_1489139220873_82713"> res@mpMinLonF = xs<br id="yui_3_16_0_ym19_1_1489139220873_82714"> res@mpMaxLonF = xe<br id="yui_3_16_0_ym19_1_1489139220873_82715"> res@mpMinLatF = ys<br id="yui_3_16_0_ym19_1_1489139220873_82716"> res@mpMaxLatF = ye<br id="yui_3_16_0_ym19_1_1489139220873_82717"><br id="yui_3_16_0_ym19_1_1489139220873_82718"> res@tiMainString = "Areas of convective cloud"<br id="yui_3_16_0_ym19_1_1489139220873_82719"> res@pmTickMarkDisplayMode = "Always" ; Turn on map tickmarks.<br id="yui_3_16_0_ym19_1_1489139220873_82720"><br id="yui_3_16_0_ym19_1_1489139220873_82721"> plot = gsn_csm_contour_map(wks,storm,res)<br id="yui_3_16_0_ym19_1_1489139220873_82722"><br id="yui_3_16_0_ym19_1_1489139220873_82723"><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82724">end</div><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82725"><div dir="ltr" id="yui_3_16_0_ym19_1_1489139220873_82800"><br id="yui_3_16_0_ym19_1_1489139220873_82726"></div></div></div></body></html>