diff --git a/src/LIB/SURCOUCHE/src/mode_io_write_nc4.f90 b/src/LIB/SURCOUCHE/src/mode_io_write_nc4.f90
index 727fde505c6a7c8de9f5babf75f421c8215bf7a3..deb8dfa689174893e9ebeb04917b8bbb22e39b45 100644
--- a/src/LIB/SURCOUCHE/src/mode_io_write_nc4.f90
+++ b/src/LIB/SURCOUCHE/src/mode_io_write_nc4.f90
@@ -302,8 +302,14 @@ IF (.NOT.GISCOORD) THEN
     IF (TPFIELD%CDIR=='XY') THEN
       if ( iboxid == 0 ) then
         ! Main domain (or no box)
-        gcoordprint =       kshape(1) == tpfile%tncdims%tdims( NMNHDIM_ARAKAWA(tpfield%ngrid,1) )%nlen &
-                      .and. kshape(2) == tpfile%tncdims%tdims( NMNHDIM_ARAKAWA(tpfield%ngrid,2) )%nlen
+        if ( tpfile%ctype == 'MNHOUTPUT' ) then
+          ! Do the hypothesis that if the field is has CDIR= 'XY', it is in the coordinate system
+          ! Checking that the size corresponds is a bit complex (need to take into account if some points were removed)
+          gcoordprint = .true.
+        else
+          gcoordprint =       kshape(1) == tpfile%tncdims%tdims( NMNHDIM_ARAKAWA(tpfield%ngrid,1) )%nlen &
+                        .and. kshape(2) == tpfile%tncdims%tdims( NMNHDIM_ARAKAWA(tpfield%ngrid,2) )%nlen
+        end if
       else
         ! Box (subdomain)
         ishift = NMNHDIM_BOX_FIRST_ENTRY - NMNHDIM_NI
@@ -340,8 +346,8 @@ IF (.NOT.GISCOORD) THEN
           DEALLOCATE(YCOORDS)
         end if
       ELSE
-        CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_attr_write_nc4','coordinates not implemented for variable ' &
-                                                                    //TRIM(TPFIELD%CMNHNAME))
+        CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Field_attr_write_nc4', TPFILE%CNAME //  ': coordinates not implemented for variable ' &
+                                                                   //TRIM(TPFIELD%CMNHNAME))
       END IF
     ELSE
       !No YCOORDS for CDIR/='XY'