The tdez2d routine is supposed to be an "easy" interface to the
Tdpack functions, but in a case like this, it might be better to not
use it, because as Rick pointed out, it is hard-coding the value for
NTRI.
Instead of using tdez2d, you can call the individual routines that
tdez2d uses, which will allow you to bump up "NTRI" yourself.
Example 4 on the Tdpack applications page shows how to call the
individual routines.
This example additionally does an isosurface (see call to "tditri"),
but you can leave this part out:
http://www.ncl.ucar.edu/Applications/tdpack.shtml
--Mary
On Feb 17, 2010, at 3:21 PM, Dennis Shea wrote:
> For crude testing purposes, you could thin the arrays
> using NCL's array syntax
>
> tdez2d(wks, x(::2), y(::2), t2(::2,::2), rho, theta, phi, -6)
>
> This decimates the input array sizes. The decimation need not be
> the same in x and y, eg
>
> tdez2d(wks, x, y(::2), t2(:,::2), rho, theta, phi, -6)
>
> Good luck
>
> Rick Brownrigg wrote:
>> 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
>> <mailto: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 <mailto:nair.ravi@epa.gov>
>>> Federal Infrastructure Contact - Joe Retzer, Ph.D.
>>> 919-541-4190 - retzer.joseph@epa.gov
>>> <mailto: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
>>> <mailto: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 <mailto:nair.ravi@epa.gov>
>>> Federal Infrastructure Contact - Joe Retzer, Ph.D.
>>> 919-541-4190 - retzer.joseph@epa.gov
>>> <mailto: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
>>> <mailto:mccormick.johnw@epa.gov>
>>>
>>> -----Forwarded by JohnW McCormick/RTP/USEPA/US on 2010.02.17
>>> 13:59 -----
>>>
>>>
>>>
>>>
>>> cc: herchenroder.barry@epa.gov <mailto: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
> _______________________________________________
> 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 15:34:36 2010
This archive was generated by hypermail 2.1.8 : Thu Feb 18 2010 - 10:33:29 MST