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
This archive was generated by hypermail 2.1.8 : Wed Aug 21 2013 - 13:54:15 MDT