ESCOMP · whlipscomb · Feb 20, 2026 · Feb 20, 2026 · Feb 20, 2026 · Mar 10, 2026
diff --git a/libglide/glide.F90 b/libglide/glide.F90
@@ -281,7 +281,7 @@ subroutine glide_initialise(model)
     ! Initialise isostasy first
     call init_isostasy(model)
 
-    select case(model%isostasy%whichrelaxed)
+    select case(model%isostasy%which_relaxed)
 
     case(RELAXED_TOPO_INPUT)   ! Supplied input topography is relaxed
 

diff --git a/libglide/glide_diagnostics.F90 b/libglide/glide_diagnostics.F90
diff --git a/libglide/glide_nc_custom.F90 b/libglide/glide_nc_custom.F90
@@ -199,6 +199,14 @@ subroutine glide_nc_filldvars(outfile,    model)
        call nc_errorhandle(__FILE__,__LINE__,status)
     end if
 
+    ! axis dimension (used for CalvingMIP output)
+
+    if (model%options%which_ho_calvingmip_domain /= HO_CALVINGMIP_DOMAIN_NONE) then
+       status = parallel_inq_varid(NCO%id,'axis',varid)
+       status= parallel_put_var(NCO%id,varid,model%calving%axis)
+       call nc_errorhandle(__FILE__,__LINE__,status)
+    end if
+
   end subroutine glide_nc_filldvars
 
 end module glide_nc_custom
diff --git a/libglide/glide_setup.F90 b/libglide/glide_setup.F90
diff --git a/libglide/glide_types.F90 b/libglide/glide_types.F90