rhoIcefront = rhoShelfice
& rhoShelfice, SHELFICEkappa,
_RL rhoShelfice
CALL WRITE_0D_RL( rhoShelfice, INDEX_NONE,