[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