The fopprofile function creates a vertical profile of surface elevation, bed elevation, and aircraft elevation along a flight line given waypoints. This function does NOT require Matlab's Mapping Toolbox but it does require Antarctic Mapping Tools and the Bedmap2 Toolbox.
fopprofile(lat,lon,z) fopprofile(x,y,z) fopprofile(lat,lon,'height',h) fopprofile(x,y,'height',h) fopprofile(...,'turns',turns) fopprofile(...,'turns',turns,'radius',TurningRadius_m) fopprofile(...,'speedkn',Speed_in_knots) fopprofile(...,'science')
fopprofile(lat,lon,z) plots a Bedmap2 profile of ice surface and bed elevations along the transect given by lat,lon. Geo coordinates can be in the form of decimal degrees or at two-column matrices of degrees and decimal minutes (i.e., [latdeg latmin],[londeg lonmin] ). The third input z is altitude in meters above the GL04c geoid. z can be a scalar to set flying height to a constant, or can be set at each waypoint. Starting and ending z values are always set to the Bedmap2 surface elevation.
fopprofile(x,y,z) specifies input waypoints in ps71 coordinates. Polar stereo- graphic coordinates are automatically determined via the islatlon function.
fopprofile(lat,lon,'height',h) sets flying height relative to the surface rathern than relative to the ellipsoid. If h is a scalar, flying altitude is the surface elevation + h at each waypoint. If h is a vector heights are set for each waypoint. Height is automatically set to zero at the start and end points. Units are meters.
fopprofile(x,y,'height',h) as above, but with ps71 input coordinates.
fopprofile(...,'turns',turns) specifies a vector of turn types for calculation by the fopturns function. Use zeros for each waypoint except 270 where at waypoints where you want big loopy turns and set turns to 180 at the at the beginning of a 180 degree turn. For more info see the fopturns documentation. If turns are not specified, straight-line distances and times are calculated.
fopprofile(...,'turns',turns,'radius',TurningRadius_m) specifies an aircraft turning radius in meters. Default turning radius is 5000 m.
fopprofile(...,'speedkn',Speed_in_knots) specifies flying speed in knots. Must be constant (scalar). Default speed is 165 knots.
fopprofile(...,'science') generates a profile tailored for science. By default, if the science option is not selected, the profile is tailored for pilots.
Specify waypoints in decimal degrees and a constant 3003 meter flying elevation:
lats = [-79.4677 -79.5170 -75.2492 -75.3733 -79.6945 -79.8210 -78.0170 -79.4677]; lons = [-112.0860 -121.7898 -131.0195 -132.3613 -123.5588 -124.9253 -129.2335 -112.0860]; z = 3003; modismoa(mean(lats),mean(lons),900,'inset','nw') plotm(lats,lons,'g-','linewidth',2)
Get an elevation profile along the green straight-line path:
Specify waypoints as ps71 coordinates, set flying height above surface at each waypoint, and estimate flight time based on 140 knot flying speed:
[x,y] = ll2ps(lats,lons); h = [0 1000 1600 1300 1100 1100 800 0]; figure fopprofile(x,y,'height',h,'speedkn',140)
Same as Example 1, but specifying 270 degree turns at waypoints 2, 3, and 7 and start a 180 degree turn at waypoint 5:
turns = [0 0 270 270 180 0 270 0]; figure fopprofile(lats,lons,z,'turns',turns);
Note that waypoints 2, 3, and 7 are marked twice in the profile because they are passed twice--once before the loop and once after the loop.
Same as Example 3, but specify a 12 km turning radius and plot a science profile:
The fopprofile function and supporting documentation were written by Chad A. Greene of the University of Texas at Austin's Institute for Geophysics (UTIG), November 2015.