diff --git a/src/LIB/SURCOUCHE/src/mode_io_tools.f90 b/src/LIB/SURCOUCHE/src/mode_io_tools.f90
index 0bc64e81353760e79b944df399dfbbbae5ab8d88..baa85c944d1ba14e02d3911e74215f7f1e4409da 100644
--- a/src/LIB/SURCOUCHE/src/mode_io_tools.f90
+++ b/src/LIB/SURCOUCHE/src/mode_io_tools.f90
@@ -226,7 +226,7 @@ SUBROUTINE IO_Dim_box_bounds( TPBOX, TPFIELD, KIINF, KISUP, KJINF, KJSUP, KKINF,
   ! Subroutine to return the box boundaries in global coordinates
   USE MODD_BAKOUT,     ONLY: LOUT_PHYSICAL_SIMPLIFIED
   USE MODD_DIM_n,      ONLY: NIMAX_ll, NJMAX_ll, NKMAX
-  USE MODD_FIELD,      ONLY: NMNHDIM_LEVEL, NMNHDIM_LEVEL_W, TFIELDMETADATA
+  USE MODD_FIELD,      ONLY: NMNHDIM_LEVEL, NMNHDIM_LEVEL_W, NMNHDIM_BOX_LEVEL, NMNHDIM_BOX_LEVEL_W, TFIELDMETADATA
   USE MODD_OUT_n,      ONLY: CMAINDOMAINNAME, TOUTBOXMETADATA
   USE MODD_PARAMETERS, ONLY: JPHEXT, JPVEXT
 
@@ -267,7 +267,8 @@ SUBROUTINE IO_Dim_box_bounds( TPBOX, TPFIELD, KIINF, KISUP, KJINF, KJSUP, KKINF,
         IF ( .NOT.PRESENT(KKINF) .OR. .NOT.PRESENT(KKSUP) ) &
           call Print_msg( NVERB_FATAL, 'IO', 'IO_Dim_box_bounds', trim(tpfield%cmnhname) // ': missing needed dummy arguments' )
 
-        IF ( TPFIELD%NDIMLIST(3) == NMNHDIM_LEVEL .OR. TPFIELD%NDIMLIST(3) == NMNHDIM_LEVEL_W ) THEN
+        IF ( TPFIELD%NDIMLIST(3) == NMNHDIM_LEVEL     .OR. TPFIELD%NDIMLIST(3) == NMNHDIM_LEVEL_W .OR. &
+             TPFIELD%NDIMLIST(3) == NMNHDIM_BOX_LEVEL .OR. TPFIELD%NDIMLIST(3) == NMNHDIM_BOX_LEVEL_W  ) THEN
           KKINF = JPVEXT + TPBOX%NKINF
           IF ( TPBOX%CNAME == CMAINDOMAINNAME .AND. TPFIELD%NGRID == 4 .AND. .NOT.LOUT_PHYSICAL_SIMPLIFIED ) THEN
             ! There is one more PHYSICAL vertical layer for w points (but the same number for physical + unphysical borders)