load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"
begin
;
; Create the input arrays.
;
xmin = -1.4
xmax = 1.4
ymin = -1.2
ymax = 1.2
ndata = 500
xi = new(ndata,float)
yi = new(ndata,float)
zi = new(ndata,float)
do i=0,ndata-1
xi(i) = xmin + (xmax-xmin)*rand()/32767.
yi(i) = ymin + (ymax-ymin)*rand()/32767.
zi(i) = xi(i) + yi(i)
t1 = 1.0/((fabs(xi(i)-0.1))^2.75 + fabs(yi(i))^2.75 + 0.09)
t2 = 1.0/((fabs(xi(i)+0.1))^2.75 + fabs(yi(i))^2.75 + 0.09)
zi(i) = 0.3*(zi(i)+t1-t2)
end do
;
; Set up a list of output points that corresponds to a 2D grid.
;
nx = 29
ny = 25
no = nx*ny
xo = new(no,float)
yo = new(no,float)
scalex = 1./(nx-1)*(xmax-xmin)
do i=0,nx-1
istart = i*ny
iend = istart+ny-1
xo(istart:iend) = xmin + i*scalex
yo(istart:iend) = fspan(ymin,ymax,ny)
end do
;
; Find the approximation surface.
;
knots = (/10,10/)
zo = csa2ls(xi,yi,zi,knots,xo,yo)
;
; Convert to a 2D array for plotting.
;
zo2 = onedtond(zo,(/nx,ny/))
;
; Create workstation object.
;
NCGM=1
X11=0
PS=0
if (NCGM .eq. 1) then
wid = gsn_open_wks("ncgm","nm19n")
end if
if (X11 .eq. 1) then
wid = gsn_open_wks("x11","nm19n")
end if
if (PS .eq. 1) then
wid = gsn_open_wks("ps","nm19n")
end if
;
; Draw plot.
;
rho = 2.5
theta = -154.
phi = 80.
xo2 = fspan(xmin,xmax,nx)
yo2 = fspan(ymin,ymax,ny)
tdez2d(wid, xo2, yo2, zo2, rho, theta, phi, 6)
frame(wid)
end
home |
contents |
defs |
procedures |
examples |
errors