From 8bfccd22293fbddde2878e6000059d3d353989ac Mon Sep 17 00:00:00 2001
From: alperaltuntas <alperaltuntas@gmail.com>
Date: Tue, 31 Oct 2023 12:51:31 -0600
Subject: [PATCH] set %label in register_netcdf_field and register_netcdf_axis

---
 src/framework/MOM_netcdf.F90 | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/src/framework/MOM_netcdf.F90 b/src/framework/MOM_netcdf.F90
index 95e6aa7bb7..73f276aba9 100644
--- a/src/framework/MOM_netcdf.F90
+++ b/src/framework/MOM_netcdf.F90
@@ -217,6 +217,8 @@ function register_netcdf_field(handle, label, axes, longname, units) &
   allocate(dimids(size(axes)))
   dimids(:) = [(axes(i)%dimid, i = 1, size(axes))]
 
+  field%label = label
+
   ! Determine the corresponding netCDF data type
   ! TODO: Support a `pack`-like argument
   select case (kind(1.0))
@@ -225,7 +227,7 @@ function register_netcdf_field(handle, label, axes, longname, units) &
     case (real64)
       xtype = NF90_DOUBLE
     case default
-      call MOM_error(FATAL, "register_netcdf_axis: Unknown kind(real).")
+      call MOM_error(FATAL, "register_netcdf_field: Unknown kind(real).")
   end select
 
   ! Register the field variable
@@ -293,6 +295,8 @@ function register_netcdf_axis(handle, label, units, longname, points, &
         "Axis must either have explicit points or be a time axis ('T').")
   endif
 
+  axis%label = label
+
   if (present(points)) then
     axis_size = size(points)
     allocate(axis%points(axis_size))