The model is a c-shell script, run at linux command line:
>csh flowline_model.csh
Required input data sets, all in same projected coordinate system. Examples given were used to calculate results for
Bevan, S.L. et al., 'Centuries of melt on Larsen C Ice Shelf', 2017. Submitted to The Cryosphere Discussions.
Par file with projection information
- e.g. larsenc_250.dem_par
Surface trajectory vector
- e.g. Cabinet_Inlet_forward_flowline.ps.txt
- e.g. Whirlwind_Inlet_forward_flowline.ps.txt
Gridded surface mass balance in metres weq per year
Gridded ice thickness and surface elevation data
- e.g. Bedmap2 dataset from British Antarctic Survey (https://secure.antarctica.ac.uk/data/bedmap2/)
Gridded surface velocity magnitude, and +/- 10% values.
- e.g. MEaSURES Antarctic velocity dataset, available from the NSIDC (https:\/nsidc.org/data/docs/measures/nsidc0484_rignot/)
Gridded vertical strain rate
Depth Density data (Depth, density, mean density from surface to depth)
- e.g. lcis_density_summary_CI_120.mean.txt
- e.g. lcis_density_summary_WI_70.mean.txt
Software
Generic Mapping Tools (http://www.soest.hawaii.edu/gmt/) and R (https://www.r-project.org/).
Output data columns are:
#$1 = easting
#$2 = northing
#$3 = segment length
#$4 = along flow velocity magnitude
#$5 = smb rate
#$6 = vstrain rate
#$7 = ice shelf thickness
#$8 = ice shelf surface elevation
#$9 = cell residence time
#$10 = smb per cell
#$11 = accumulated time
#$12 = accumulated distance
#$13 = accumulated smb
#$14 = accumulated z1 from start point 1 etc.
#$15 = accumulated z2
#$16 = accumulated z3
#$17 = accumulated z4
#$18 = accumulated z5
#$19 = accumulated z6