diff --git a/libglissade/glissade.F90 b/libglissade/glissade.F90
index ef0dc5dd..4093a137 100644
--- a/libglissade/glissade.F90
+++ b/libglissade/glissade.F90
@@ -457,7 +457,11 @@ subroutine glissade_initialise(model, evolve_ice)
     ! handle relaxed/equilibrium topo
     ! Initialise isostasy first
 
-    call init_isostasy(model)
+    if (model%options%isostasy == ISOSTASY_COMPUTE) then
+
+       call init_isostasy(model)
+
+    endif
 
     select case(model%isostasy%whichrelaxed)