ICESat data import

The Fortran 90/95 algorithm reads data from a set of binary ICESat files, filters them according to user defined geographic, temporal and quality settings and writes the resulting records into ascii files (for GIS import as point layers) and binary files (for analysis with the algorithm for the detection of height differences between neighbouring elevations). The original version of the algorithm is described in Alberti & Biscaro (2010).


Language: Fortran 90/95
Compiler: g95 - Windows
Version: 2016-04-30

Since data in GLA12 files are in big-endian notation (the Unix default), when working in Windows you must compile the program using the big-endian data option.
For instance, with the free g95 compiler, the command line should include the "-fendian=BIG" option, e.g.:
  prompt> g95 ReadGLA12_v5.0.f95 -fendian=BIG -o readgla12.exe

Program input

Two ascii file, the former with the list of GLA12 binary files to read (one record for line), e.g.:


the latter with the analysis parameters, e.g.:

-90 90
-180 180
2005 1 1 2009 1 1
300 100


-90 90 #minimum and maxixum latitude
-180 180 #minimum and maxixum longitude
2005 1 1 2009 1 1 #time window: initial and final year month day; if not required: -1 as first value
300 100 #filters on saturation elevation correction and gain value; when not desired substitute value with -1 (e.g. 300 -1 or -1 100 or -1 -1 )

Program output

A set of ascii files for GIS import, each one corresponding to an input Glas binary file, and a binary file for height variation analysis with the HeightVarGLA_v1.0 program. You can find further methodology details in Alberti & Biscaro (2010).


Mauro Alberti