reading gtopo30 data

Hi,
I am having some problem trying to read the gtopo30 binary data, could
anyone help me out?
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
 setfileoption("bin","ReadByteOrder","BigEndian")
 nlat=6000-1
 nlon=4800-1
 z=fbindirread ("E060N40.DEM",0,(/nlon,nlat/),"byte")
 res=0.00833333333333
 ulat=39.99583333333333
 ulon=60.00416666666667
 lat = ulat + ispan(-nlat,0,1)*res
 lon = ispan(0,nlon,1)*res + ulon

 ;=============================
 ; create lat and long coordinate variables
 ;============================
 lon!0 = "lon"
 lon@long_name = "lon"
 lon@units = "degrees-east"
 lon&lon = lon

 lat!0 = "lat"
 lat@long_name = "lat"
 lat@units = "degrees_north"
 lat&lat = lat
 print(lat)
 print(lon)
 ;=============================
 ; name dimensions of z and assign coordinate variables
 ;============================
  z!0 = "lat"
  z!1 = "lon"
  z&lat = lat
  z&lon = lon
  z@long_name = "Elevation"
  z@units = "meters"

Its giving a error message saying
fatal:Coordinate variables must be the same dimension as their dimension
fatal:No coordinate variable exists for dimension (lat) in variable (z)
Segmentation fault

Also I am not sure if the format is "byte" if I use double or float while
reading the binary file, it gives error message.
Regards,
Prabhakar

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Nov 26 19:52:45 2009