Philippe 11/04/2024: Fill_tfieldmetadata: populate ndimlist if CDIR=XY for 4D...

Philippe 11/04/2024: Fill_tfieldmetadata: populate ndimlist if CDIR=XY for 4D fields (4th dimension set to NMNHDIM_NOTLISTED)

src/LIB/SURCOUCHE/src/modd_field.f90

@@ -425,6 +425,9 @@ type(tfieldmetadata) function Fill_tfieldmetadata( cmnhname, cstdname, clongname
           tpfield%ndimlist(1:2) = NMNHDIM_ARAKAWA(tpfield%ngrid,1:2)
         else if ( tpfield%ndims == 3 .and. tpfield%ngrid /= NGRIDUNKNOWN ) then
           tpfield%ndimlist(1:3) = NMNHDIM_ARAKAWA(tpfield%ngrid,1:3)
+        else if ( tpfield%ndims == 4 .and. tpfield%ngrid /= NGRIDUNKNOWN ) then
+          tpfield%ndimlist(1:3) = NMNHDIM_ARAKAWA(tpfield%ngrid,1:3)
+          tpfield%ndimlist(4) = NMNHDIM_NOTLISTED
         else
           call Print_msg( NVERB_DEBUG, 'GEN', 'Fill_tfielddata', 'ndimlist not filled for variable ' // Trim( ymnhname ) )
         end if