Add some comments

Author: billsacks

cesm/flux_atmocn/flux_atmocn_COARE_mod.F90

@@ -153,9 +153,12 @@ subroutine flux_atmOcn_COARE(                       &
           endif
 
           if (aofluxes_use_shr_wv_sat) then
+             ! This version uses a qsat calculation method consistent with what's used in CAM
              call shr_wv_sat_qsat_liquid(ts(n), pslv(n), esat_val, qsat_val)
              ssq = 0.98_R8 * qsat_val   ! sea surf hum (kg/kg)
           else
+             ! This version uses the qsat calculation method that was used for many years,
+             ! prior to Aug 2025, and which is still being used by default in NorESM
              ssq = 0.98_R8 * qsat(ts(n)) / rbot(n)   ! sea surf hum (kg/kg)
           end if
 

cesm/flux_atmocn/flux_atmocn_Large.F90

@@ -215,9 +215,12 @@ subroutine flux_atmOcn_large(              &
           endif
 
           if (aofluxes_use_shr_wv_sat) then
+             ! This version uses a qsat calculation method consistent with what's used in CAM
              call shr_wv_sat_qsat_liquid(ts(n), pslv(n), esat_val, qsat_val)
              ssq    = 0.98_R8 * qsat_val                ! sea surf hum (kg/kg)
           else
+             ! This version uses the qsat calculation method that was used for many years,
+             ! prior to Aug 2025, and which is still being used by default in NorESM
              ssq    = 0.98_R8 * qsat(ts(n)) / rbot(n)   ! sea surf hum (kg/kg)
           end if
           delt   = thbot(n) - ts(n)                  ! pot temp diff (K)