diff --git a/src/cvmix_driver/gen_modules_cvmix_kpp.F90 b/src/cvmix_driver/gen_modules_cvmix_kpp.F90
index 62d0ff1a8..5769230f2 100644
--- a/src/cvmix_driver/gen_modules_cvmix_kpp.F90
+++ b/src/cvmix_driver/gen_modules_cvmix_kpp.F90
@@ -931,8 +931,8 @@ subroutine calc_cvmix_kpp(ice, dynamics, tracers, partit, mesh)
                 ! applied and where surface averages are taken — not something that 
                 ! varies with depth.
                 sldepth = kpp_surf_layer_ext*max(kpp_obldepth(node), kpp_minOBLdepth)
-                nzsfc = kpp_obldepth(node)
-                do nztmp = nun, int(kpp_obldepth(node))
+                nzsfc = int(kpp_nzobldepth(node))
+                do nztmp = nun, int(kpp_nzobldepth(node))
                     if (-zbar_3d_n(nztmp+1,node)>=sldepth) then
                         nzsfc = nztmp
                         exit