Re: Having difficulty in generating a 3D plot from the input (Oceanic Temperature) Netcdf file:

Hi John,

I am not knowledgeable with the use of the tdez2d command, but in
looking through the underlying code (TDSTRI in particular), it looks
like a regular grid gets broken up into triangles, 2 per rectangular
grid cell. It also appears there's a hard limit on the maximum number
of triangles, at 110,000 triangles. If I understand everything, your
input with NX=318, NY=191 produces a grid of 60738 cells, which
times-2, exceeds the 110000 limit.

If you built NCL from source, I can tell you where/how to bump up the
number of triangles. If not, and if you have the flexibility in your
data, perhaps see if you can find settings for NX, NY such that
NX*NY*2 < 110,000.

I hope that helps...
Rick

On Feb 17, 2010, at 12:13 PM, McCormick.JohnW@epamail.epa.gov wrote:

>
> ######################################################################
> ### After running this script on a Linux 2.6.18-164.11.1.el5 OS,
> (Dell Precision Hardware Station)
> ### I got the following message below:
> ###
> #####################################################################
> iris1.rtpnc.epa.gov:/home/mccormik/GulfBreeze/Software02/Iso-Plots%
> ncl 20100212_3d_Test00.ncl
> Copyright (C) 1995-2009 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 5.1.1
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
> Triangle list overflow in TDSTRI
> iris1.rtpnc.epa.gov:/home/mccormik/GulfBreeze/Software02/Iso-Plots%
>
>
>
> However, one observation that I had noticed about this program is
> that when I reduced the indices (sizes)
> to about 80 or 90 (for all three dimensions), it ran okay. The
> above problem occurred when I set the following:
> NX = 318
> NY = 191
>
> Also, where the system error message says: "Triangle list overflow
> in TDSTRI", could that mean that I'm not using the
> correct function ?? I have been assuming that the function would
> accept values from unit-cubes (as grids) and not triangles.
>
> Please let me know if you (or anyone who sees this message) have
> come across this condition and if there is any way to correct or
> reconstruct the NCL script in "fixing" this overflow problem.
>
>
> John W. McCormick, Systems Programmer / DBA
> High Performance Computing & Scientific Visualization
> Lockheed Martin Information Technology, Supporting the EPA
> Research Triangle Park, NC
> 919-541-0890
>
> Federal Infrastructure Contact - Ravi Nair
> 919-541-5467 - nair.ravi@epa.gov
> Federal Infrastructure Contact - Joe Retzer, Ph.D.
> 919-541-4190 - retzer.joseph@epa.gov
>
> ############################################
> ### Work
> Information ###
> ############################################
> John W. McCormick, Systems Programmer, MD: N127-01
> Lockheed Martin, Contractor to US-EPA
> US-EPA, Bldg NCC, Rm N135B
> Durham, NC 27711
> Office: 919-541-0890
> Fax : 919-541-0056
> Work Email: mccormick.johnw@epa.gov
>
>
> John W. McCormick, Systems Programmer / DBA
> High Performance Computing & Scientific Visualization
> Lockheed Martin Information Technology, Supporting the EPA
> Research Triangle Park, NC
> 919-541-0890
>
> Federal Infrastructure Contact - Ravi Nair
> 919-541-5467 - nair.ravi@epa.gov
> Federal Infrastructure Contact - Joe Retzer, Ph.D.
> 919-541-4190 - retzer.joseph@epa.gov
>
> ############################################
> ### Work
> Information ###
> ############################################
> John W. McCormick, Systems Programmer, MD: N127-01
> Lockheed Martin, Contractor to US-EPA
> US-EPA, Bldg NCC, Rm N135B
> Durham, NC 27711
> Office: 919-541-0890
> Fax : 919-541-0056
> Work Email: mccormick.johnw@epa.gov
>
> -----Forwarded by JohnW McCormick/RTP/USEPA/US on 2010.02.17 13:59
> -----
>
>
>
>
> cc: herchenroder.barry@epa.gov
> Subject: I have attempted to run the following script which would
> generate a 3D plot from the input (Oceanic Temperature) Netcdf file:
>
> To: William Hobbs:
> From: John McCormick
> Date: 02/16/2010
> Subject: 3D plot from the input (Oceanic Temperature) Netcdf file
> gives overflow system message:
>
> ######################################################################
> ### Script to run in attempting to get a 3-D Plot of Temperature in
> the
> ### Northern Gulf (off the coast of Louisiana)....
> ###
> #####################################################################
> ;*************************************************
> ;
> ;************************************************
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"
> ;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/ccm_func.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/popRemap.ncl"
> 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"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"
> begin
> NX = 318 ; Longitude Dimension (x) -- 318
> NY = 191 ; Latitude Dimension (y) -- 191
> x = new((/NX/),float)
> y = new((/NY/),float)
> t2 = new((/NX,NY/),float)
> ; x = (/0.00, 1.00, 2.00/)
> ; y = (/0.00, 1.00, 2.00/)
>
> ;*************************************************
> ; Read in Temperature Values at depth of 50 meters
> ;************************************************
>
> nc_file_in ="/home/mccormik/GulfBreeze/Data/TemperatureOceanFiles/
> 20090916_20090918/t3d_2009091600.nc"
> ; a = addfile("/home/mccormik/GulfBreeze/
> t3d_2009061100_nc_ascii.nc","r")
> ; a = addfile("/amber/home/mccormik/GulfBreeze/
> t3d_2009061100_nc_ascii.nc","r")
> a = addfile(nc_file_in,"r")
> idepth = 1
> T = a->Temperature(0,idepth,:,:) ; (time,lev,lat,lon)
> lon = a->Longitude
> lat = a->Latitude
> dep = a->Depth
> ;print(T)
> ;print(lon)
> ;print(lat)
>
> do i = 0,317 ; Longitude (x from 0 thru 317)
> x(i) = (i)*1.0
> do j = 0,190 ; Latitude
> y(j) = (j)*1.0
> t2(i,j)=T(j,i)*1.0
> end do
> end do
>
> ; print(x)
> ; print(y)
> ; print(t2)
> ; Create the output grid.
> xmin = -2.0
> ymin = -2.0
> zmin = -2.0
> xmax = 2.0
> ymax = 2.0
> zmax = 2.0
> ; ii = fspan(0,20.,21)
> ; xo = xmin + (ii/(NX-1)) * (xmax-xmin)
> ; yo = ymin + (ii/(NY-1)) * (ymax-ymin)
> ; zo = zmin + (ii/(NZ-1)) * (zmax-zmin)
> ;print(x)
> ;print(y)
> ;print(z)
> wks = gsn_open_wks("x11","test")
> rho = 2.5
> theta = -136.
> phi = -127.
> ;
> ; Interpolate. "output" will be a 3-dimensional grid dimensioned
> 21x21x21.
> ;
> tdez2d(wks, x, y, t2, rho, theta, phi, -6)
> frame(wks)
> end
>
>
>
>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Feb 17 14:02:28 2010