Philippe 06/05/2024: outputs: do not write 2D coordinates for main domain if not needed

src/LIB/SURCOUCHE/src/mode_io_write_nc4.f90

@@ -1863,7 +1863,7 @@ if ( .not. lcartesian ) then
     call Gather_hor_coord2d( zxhat,  zyhat,  zlatf_glob, zlonf_glob )
   end if
 
-  if ( tzfile%ctype /= 'MNHOUTPUT' .or. lout_bigbox_write ) then !Do not write coordinates of main domain if it is not written
+  if ( tzfile%ctype /= 'MNHOUTPUT' ) then !Do not write coordinates of main domain if it is not written
     ! Mass point
     call Write_hor_coord2d( zlatm_glob, zlonm_glob,  'latitude',  'longitude',  iim_min, iim_max, ijm_min, ijm_max )
     ! u point
@@ -1872,6 +1872,22 @@ if ( .not. lcartesian ) then
     call Write_hor_coord2d( zlatv_glob, zlonv_glob, 'latitude_v','longitude_v', iiv_min, iiv_max, ijv_min, ijv_max )
     ! xi vorticity point (=f point =uv point)
     call Write_hor_coord2d( zlatf_glob, zlonf_glob, 'latitude_f','longitude_f', iiu_min, iiu_max, ijv_min, ijv_max )
+  else if ( tzfile%ctype == 'MNHOUTPUT' .and. lout_bigbox_write ) then !Do not write coordinates of main domain if it is not written
+    ! Write both coordinates (lat and lon even if only one of them is needed).
+    ! In practice, this simplification should not have real impact (most of the time, the 2 coordinates appear simultaneously)
+
+    ! Mass point
+    if ( tout_boxes(0)%lcoord_latm .or. tout_boxes(0)%lcoord_latm ) &
+      call Write_hor_coord2d( zlatm_glob, zlonm_glob,  'latitude',  'longitude',  iim_min, iim_max, ijm_min, ijm_max )
+    ! u point
+    if ( tout_boxes(0)%lcoord_latu .or. tout_boxes(0)%lcoord_latu ) &
+      call Write_hor_coord2d( zlatu_glob, zlonu_glob, 'latitude_u','longitude_u', iiu_min, iiu_max, iju_min, iju_max )
+    ! v point
+    if ( tout_boxes(0)%lcoord_latv .or. tout_boxes(0)%lcoord_latv ) &
+      call Write_hor_coord2d( zlatv_glob, zlonv_glob, 'latitude_v','longitude_v', iiv_min, iiv_max, ijv_min, ijv_max )
+    ! xi vorticity point (=f point =uv point)
+    if ( tout_boxes(0)%lcoord_latf .or. tout_boxes(0)%lcoord_latf ) &
+      call Write_hor_coord2d( zlatf_glob, zlonf_glob, 'latitude_f','longitude_f', iiu_min, iiu_max, ijv_min, ijv_max )
   end if
 
   Deallocate( zlat, zlon )