Do not call init_isostasy unless isostasy = 1

The code was calling subroutine init_isostasy with isostasy = 0 = ISOSTASY_NONE. This subroutine is now called only if isostasy = 1 = ISOSTASY_COMPUTE.
ESCOMP · Sep 18, 2021 · 44f771c · 44f771c
1 parent 165fb8a
commit 44f771c
Showing 1 changed file with 5 additions and 1 deletion.
diff --git 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)