Hi, Fan
Thank you for the reply! I've modified my scripts as below:
setfileoption("bin","ReadByteOrder","LittleEndian")
lat = fbindirread("psn25lats_v3.dat", 0, (/nrows,ncols/) , "integer") / 100000. ; lat and lon stored in 4-byte integer, scaled by 100,000
lon = fbindirread("psn25lons_v3.dat", 0, (/nrows,ncols/) , "integer") / 100000. ; http://nsidc.org/data/polar_stereo/tools_geo_pixel.html#psn25_pss25_lats
data = fbindirread(data_path + data_filename, 0, (/nrows,ncols/), "short") / 10. ; 2-byte integer, scaled by 10, http://nsidc.org/data/docs/daac/nsidc0079_bootstrap_seaice.gd.html
and it works.
If "LittleEndian" was changed to "Native", the results will not be changed. I don't know why there are no difference between them.
Thank
Leo
ÔÚ 2013-8-22£¬ÉÏÎç2:59£¬Tingting Fan <tingting@ucar.edu> дµÀ£º
> Hi Leo,
> I used to use SIC data based on Bootstrap algorithm, I can shear the script.
> I think the problem is that you did not describe data correctly, while you
> should re-read the document. If you modify your script as:
> ********************************
> ...
> setfileoption("bin","ReadByteOrder","Native")
> y = fbindirread("./NSIDC/south/bt_201112_f17_v02_s.bin", 0, -1, "short")
> s = onedtond(y,(/nrows,ncols/))
> ...
> ********************************
> It may work!
> If you have any questions, feel free to ask!
> Tingting
>
> --
> Tingting Fan
> grad student visitor
> NCAR/CGD/CAS
> tel:303-497-1745<binread.ncl>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Aug 21 19:04:05 2013
This archive was generated by hypermail 2.1.8 : Fri Aug 30 2013 - 14:04:57 MDT