How to read the Bootstrap SIC data from NSIDC correctly?

From: gmail <leo.aries.g_at_nyahnyahspammersnyahnyah>
Date: Wed Aug 21 2013 - 08:49:14 MDT

Hi, all

      I want to read the bootstrap sea ice concentration from NSIDC (http://nsidc.org/data/docs/daac/nsidc0079_bootstrap_seaice.gd.html). Here below

is my scripts :

;*************************************************
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"

begin

data_path = "/Data/NSIDC/nsidc0079_gsfc_bootstrap_seaice/final-gsfc/north/monthly/"
data_filename = "bt_201212_f17_v02_n.bin"
plot_filename = data_filename

nrows = 448
ncols = 304

setfileoption("bin","ReadByteOrder","LittleEndian")
lat = fbindirread("psn25lats_v3.dat",0,(/nrows,ncols/), "float") ; about the .dat , please refer to http://nsidc.org/data/polar_stereo/tools_geo_pixel.html#geolocate
lon = fbindirread("psn25lons_v3.dat",0,(/nrows,ncols/), "float") ; I also try the other function such as fbinread, fbinrecread,

lat = lat/100000.
lon = lon/100000.
lat@units = "degrees_north"
lon@units = "degrees_east"
printVarSummary(lon)
printVarSummary(lat)
print(max(lon)+" " +min(lon))
print(max(lat)+" " +min(lat))

 data = fbindirread(data_path + data_filename, 0, (/nrows,ncols/), "float")
 printVarSummary(data)
 print(max(data) + " " + min(data))

¡­¡­..
; ****************************************

However, I got the error information as listed in the following :

Variable: lon
Type: float
Total Size: 544768 bytes
            136192 values
Number of Dimensions: 2
Dimensions and sizes: [448] x [304]
Coordinates:
Number Of Attributes: 1
  units : degrees_east

Variable: lat
Type: float
Total Size: 544768 bytes
            136192 values
Number of Dimensions: 2
Dimensions and sizes: [448] x [304]
Coordinates:
Number Of Attributes: 1
  units : degrees_north
(0) 2.69049e-43 -3.40272e+33
(0) 1.26117e-43 4.34403e-44
fatal:fbindirread: The size implied by the dimension array and record number is greater that the size of the file, can't continue
fatal:["Execute.c":7743]:Execute: Error occurred at or near line 53 in file read.ncl

I've read through the instruction in the related url, and also try the other function such as fbinread, fbinrecread, cbinread, but got no correct results.
Anyone can give me some advise£¿

Thanks!

Leo

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Wed Aug 21 12:29:55 2013

This archive was generated by hypermail 2.1.8 : Wed Aug 21 2013 - 13:54:15 MDT