About PIG sensitivity tutorial

in-woo

Hello, ISSM users.

I have a question with 'PIG sensitivity tutorial' that shows floating ice figure.

According to https://issm.jpl.nasa.gov/documentation/tutorials/pigsensitivity/, transient results were showed as map with changed floating ice.

And how can I show this figure?

As a result of transient solution, struct of 'md.results.TransientSolution' has following :

StressbalanceConvergenceNumSteps

step

time

Vel

Pressure

Thickness

Surface

Base

IceVolume

IceVolumeAboveFloatation

SmbMassBalance

MaskGroundediceLevelset

SolutionType

errlog

outlog

Thank you

heleneseroussi

Dear In-Woo,

The extent of the floating and grounded ice are stored in md.results.TransientSolution(XX).MaskGroundediceLevelset, so you can simply plot this field at the beginning and the end of the simulation to see how the floating ice extent evolves in the simulation.

Best,

Helene

in-woo

Dear Helene.

Thank you for your reply.

But when I plot figure with using. 'md.results.TransientSolution(end).MaskGroundediceLevelset' and 'md.results.TransientSolution(1).MaskGroundediceLevelset', the figure shows like this (Figure 1).

run1.ini&en.floatingice.jpg

In tutorial, figure seems to be normalized MaskGroundediceLevelset.

Screenshot from 2016-12-20 16-46-25.png

I want to know how to show like this.

Do I have to convert MaskGroundediceLevelset with following code?


% grounded ice
pos0 = find(md.results.TransientSolution(xx).MaskGroundediceLevelset > 0)
md.results.TransientSolution(1).MaskGroundediceLevelset(pos0) = 1;
% floating ice
pos0 = find(md.results.TransientSolution(xx).MaskGroundediceLevelset < 0)
md.results.TransientSolution(1).MaskGroundediceLevelset(pos0) = -1;

Thank you.

heleneseroussi

Hello In-Woo,

This is the same field and you don't have to convert anything. The only difference is that I used a different 'caxis' in order to better highlight the grounded and floating parts. So just add: 'caxis', in plotmodel when you plot these fields. The grounding line is where the 0 values are.

Best,

Helene