<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi Petroula,<br>
<br>
you do not have to create two plots, just use the min and max
function of NCL to get the minimum and maximum of lat and lon. See
attached script and plot.<br>
<br>
How to install the rangs/gshhs zip files read the page <br>
<br>
<a class="moz-txt-link-freetext" href="http://ncl.ucar.edu/Document/Graphics/rangs.shtml">http://ncl.ucar.edu/Document/Graphics/rangs.shtml</a><br>
<br>
Bye,<br>
Karin<br>
<br>
<div class="moz-cite-prefix">Am 10.03.17 um 13:16 schrieb Petroula
Louka:<br>
</div>
<blockquote
cite="mid:936383411.4728973.1489148208386@mail.yahoo.com"
type="cite">
<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>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
ncl-talk mailing list
<a class="moz-txt-link-abbreviated" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>
List instructions, subscriber options, unsubscribe:
<a class="moz-txt-link-freetext" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a>
</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
Dipl. Geophys. Karin Meier-Fleischer
Visualization, NCL
Application Support
Deutsches Klimarechenzentrum GmbH (DKRZ)
Bundesstrasse 45a - D20146 Hamburg - Germany
Phone: +49 (0)40 460094 126
Fax: +49 (0)40 460094 270
E-Mail: <a class="moz-txt-link-abbreviated" href="mailto:meier-fleischer@dkrz.de">meier-fleischer@dkrz.de</a>
URL: <a class="moz-txt-link-abbreviated" href="http://www.dkrz.de">www.dkrz.de</a>
Geschäftsführer: Prof. Dr. Thomas Ludwig
Sitz der Gesellschaft: Hamburg
Amtsgericht Hamburg HRB 39784
</pre>
</body>
</html>