School of GeoSciences

School of GeoSciences

Menu

PyGMT Examples

My main research interests are concerned with modelling paleo-icesheets. The model produces a large amount of output related to various aspects of the simulated ice sheet. Output typically includes scalar and vector fields.

Complex images are constructed from basic units. This example shows a plot of the surface topography used for the European experiments. Bedrock Topography
This image shows the ice surface at the simulated last glacial maximum. The area above sea level is displayed in grey. This image also shows two profile lines. Ice Sheet Surface
This image shows the internal temperature and velocity structure of the ice sheet along profile A, indicated in the previous image. Visualisation of the 3D fields is somewhat tricky because the vertical axis is scaled to the interval [0,1]. The scaled vertical axis has to be transformed to real vertical coordinates before visualisation. This was one of the main reasons why I embarked on writing a Python wrapper for GMT. Temperature and Velocity Profile
Another way to visualise simulated ice sheet behaviour is to extract data along a profile (Profile A, above) for each time slice and display it in a time-distance diagram. This image shows the basal temperature of the ice sheet. Red indicates ice at the melting point. The image also shows the location of the ice divide. Time-Distance Diagram
Complex diagrams can be generated by putting all the basic types together. This example shows surface run-off and melt water production for both sides of the ice divide. Discharge from Ice Sheet
This plot shows selected sea-level observations around Fennoscandinavia and the corresponding simulated sea-level curves. Discharge from Ice Sheet