From dfbc65829e9dd8e8067330152fa0e68a15834cca Mon Sep 17 00:00:00 2001
From: Michael Levy <mike.levy.work@gmail.com>
Date: Fri, 19 Apr 2024 17:08:39 -0600
Subject: [PATCH] Add chksum calls for MARBL forcings

Updated ice_ocn_bnd_type_chksum() in the NUOPC cap, though I don't think this
function is ever called
---
 .../nuopc_cap/mom_surface_forcing_nuopc.F90   | 55 +++++++++++++++++++
 1 file changed, 55 insertions(+)

diff --git a/config_src/drivers/nuopc_cap/mom_surface_forcing_nuopc.F90 b/config_src/drivers/nuopc_cap/mom_surface_forcing_nuopc.F90
index f9bd777a9b..ba8fbc750c 100644
--- a/config_src/drivers/nuopc_cap/mom_surface_forcing_nuopc.F90
+++ b/config_src/drivers/nuopc_cap/mom_surface_forcing_nuopc.F90
@@ -1533,6 +1533,60 @@ subroutine ice_ocn_bnd_type_chksum(id, timestep, iobt)
     chks = field_chksum( iobt%mass_berg  ) ; if (root) write(outunit,100) 'iobt%mass_berg      ', chks
   endif
 
+  ! MARBL forcing
+  if (associated(iobt%atm_fine_dust_flux)) then
+    chks = field_chksum(iobt%atm_fine_dust_flux)
+    if (root) write(outunit,110) 'iobt%atm_fine_dust_flux   ', chks
+  endif
+  if (associated(iobt%atm_coarse_dust_flux)) then
+    chks = field_chksum(iobt%atm_coarse_dust_flux)
+    if (root) write(outunit,110) 'iobt%atm_coarse_dust_flux ', chks
+  endif
+  if (associated(iobt%seaice_dust_flux)) then
+    chks = field_chksum(iobt%seaice_dust_flux)
+    if (root) write(outunit,110) 'iobt%seaice_dust_flux     ', chks
+  endif
+  if (associated(iobt%atm_bc_flux)) then
+    chks = field_chksum(iobt%atm_bc_flux)
+    if (root) write(outunit,110) 'iobt%atm_bc_flux          ', chks
+  endif
+  if (associated(iobt%seaice_bc_flux)) then
+    chks = field_chksum(iobt%seaice_bc_flux)
+    if (root) write(outunit,110) 'iobt%seaice_bc_flux       ', chks
+  endif
+  if (associated(iobt%nhx_dep)) then
+    chks = field_chksum(iobt%nhx_dep)
+    if (root) write(outunit,100) 'iobt%nhx_dep    ', chks
+  endif
+  if (associated(iobt%noy_dep)) then
+    chks = field_chksum(iobt%noy_dep)
+    if (root) write(outunit,100) 'iobt%noy_dep    ', chks
+  endif
+  if (associated(iobt%atm_co2_prog)) then
+    chks = field_chksum(iobt%atm_co2_prog)
+    if (root) write(outunit,110) 'iobt%atm_co2_prog         ', chks
+  endif
+  if (associated(iobt%atm_co2_diag)) then
+    chks = field_chksum(iobt%atm_co2_diag)
+    if (root) write(outunit,110) 'iobt%atm_co2_diag         ', chks
+  endif
+  if (associated(iobt%afracr)) then
+    chks = field_chksum(iobt%afracr)
+    if (root) write(outunit,100) 'iobt%afracr     ', chks
+  endif
+  if (associated(iobt%swnet_afracr)) then
+    chks = field_chksum(iobt%swnet_afracr)
+    if (root) write(outunit,110) 'iobt%swnet_afracr         ', chks
+  endif
+  if (associated(iobt%ifrac_n)) then
+    chks = field_chksum(iobt%ifrac_n)
+    if (root) write(outunit,100) 'iobt%ifrac_n    ', chks
+  endif
+  if (associated(iobt%swpen_ifrac_n)) then
+    chks = field_chksum(iobt%swpen_ifrac_n)
+    if (root) write(outunit,110) 'iobt%swpen_ifrac_n        ', chks
+  endif
+
   ! enthalpy
   if (associated(iobt%hrofl)) then
     chks = field_chksum( iobt%hrofl  ) ; if (root) write(outunit,100) 'iobt%hrofl      ', chks
@@ -1554,6 +1608,7 @@ subroutine ice_ocn_bnd_type_chksum(id, timestep, iobt)
   endif
 
 100 FORMAT("   CHECKSUM::",A20," = ",Z20)
+110 FORMAT("   CHECKSUM::",A30," = ",Z20)
 
   call coupler_type_write_chksums(iobt%fluxes, outunit, 'iobt%')