[ncl-talk] Cannot create triangular mesh
Adam Phillips
asphilli at ucar.edu
Tue Apr 2 11:01:59 MDT 2019
Hello,
I ran your script using version 6.6.2. I didn't get any error messages, and
it produced the attached plot. NCL v6.0.0 is almost 8 years old, I highly
recommend updating to the newest version (6.6.2).
Adam
On Tue, Apr 2, 2019 at 5:18 AM MM Ali <mmali110 at gmail.com> wrote:
> Hi,
> I took a program of Dr. Petroula modified by Prof. Thomas from ncl-talk
> and slightly modified to suit my requirements (attached below). The program
> plots the continent and the legent without any contour. The error is:
>
> fatal:Cannot create triangular mesh: supply additional resources or build
> with Triangle package
> fatal:ContourPlotDraw: draw error
> warning:WorkstationDeactivate: workstation not active or not opened
>
> The data set I used is attached. The ncl version I use is 6.0.0.0
>
> It would be great if someone can help me in solving this problem.
>
> Thanks
> Ali
> The Code (also attached):
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> begin
>
> fname = "tst.txt"
>
> n_col = numAsciiCol(fname)
> n_var = n_col
> n_pt = numAsciiRow(fname)
> dataf = asciiread(fname,(/n_pt,n_var/),"float")
>
> ;---read longitude
> lon = dataf (:,0)
>
> ;---read latitude
> lat = dataf (:,1)
>
> ;---read data for presence of convective cloud
> storm = dataf (:,2)
>
> ;---define the workstation (plot type and name)
> wks = gsn_open_wks("png","plot_field")
>
> ;--set resources
> wks = gsn_open_wks("x11","ce")
> res = True
> res at gsnMaximize = True
>
> res at mpDataBaseVersion = "MediumRes" ;-- map resolution
> res at mpMinLonF = min(lon)
> res at mpMaxLonF = max(lon)
> res at mpMinLatF = min(lat)
> res at mpMaxLatF = max(lat)
>
> res at sfXArray = lon
> res at sfYArray = lat
>
> res at cnFillOn = True ; Turn on contour fill.
> res at cnLinesOn = False ; Turn off contour lines
>
> ;-- create the contour map plot
> ; plot = gsn_csm_contour_map(wks,storm,res) ; create a default
> plot
> contour = gsn_csm_contour_map(wks,storm,res)
>
> end
>
>
>
>
> --
> *Dr. MM Ali*
> Formerly:
> Scientist G & Group Director
> NRSC/ISRO
> email: mmali110 at gmail.com
> mmali at coaps.fsu.edu
> Mob: +91 9949045110
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
--
Adam Phillips
Associate Scientist, Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/ 303-497-1726
<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190402/cf8d0bdd/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ce.png
Type: image/png
Size: 95373 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190402/cf8d0bdd/attachment-0001.png>
More information about the ncl-talk
mailing list