getind_latlon2d

The following script, using getind_latlon2d, results in this error:

fatal:Conditional statements (if and do while) require SCALAR logical
values, see all and any functions
fatal:Execute: Error occurred at or near line 8538 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl
fatal:Execute: Error occurred at or near line 53

Any ideas why getind_latlon2d would produce this error?

Thanks, Mike

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"

begin

nlat=1801
nlon=3601

a = fbindirread("/Users/notaro/ndvi/namer_albers.lat",0,(/
1024,1280/),"short")
lat2d=short2flt(a)
lat2d=lat2d/100.
b = fbindirread("/Users/notaro/ndvi/namer_albers.lon",0,(/
1024,1280/),"short")
lon2d=short2flt(b)
lon2d=lon2d/100.
lat2d!0="lat"
lat2d!1="lon"
lon2d!0="lat"
lon2d!1="lon"
lat2d_at_units="degrees_north"
lon2d_at_units="degrees_east"

lats=ndtooned(lat2d)
lons=ndtooned(lon2d)

lat=fspan(-90.,90.,nlat)
lon=fspan(-180.,180.,nlon)
lat!0="lat"
lat_at_units="degrees_north"
lat&lat=lat
lon!0="lon"
lon_at_units="degrees_east"
lon&lon=lon

latit=new((/nlat,nlon/),float)
longit=new((/nlat,nlon/),float)
do i=0,nlat-1
   latit(i,:)=lat(i)
end do
do j=0,nlon-1
   longit(:,j)=lon(j)
end do

nm=getind_latlon2d(latit,longit,lats,lons)

print(nm)

end

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Mar 14 2008 - 08:49:57 MDT