Re: How to read the Bootstrap SIC data from NSIDC correctly

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

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-222:59Tingting 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