Remove some commented-out code related to qsat calculation

Instead, refer to the relevant issue

Author: billsacks

cesm/flux_atmocn/flux_atmocn_Diurnal_mod.F90

@@ -26,7 +26,6 @@ module flux_atmocn_diurnal_mod
    use shr_const_mod,         only : shr_const_ocn_ref_sal, shr_const_zsrflyr, shr_const_rgas
    use shr_sys_mod,           only : shr_sys_abort
    use flux_atmocn_COARE_mod, only : cor30a
-   use shr_wv_sat_mod,        only : shr_wv_sat_qsat_liquid ! use saturation calculation consistent with CAM
 
    implicit none
    private
@@ -236,7 +235,8 @@ subroutine flux_atmOcn_diurnal(                      &
     real(R8)    :: tdiff(nMax)               ! tbot - ts
     real(R8)    :: vscl
 
-    ! NOTE: this should use the shr_wv_sat_qsat_liquid if this routine is ever used in production
+    ! NOTE: this should use the shr_wv_sat_qsat_liquid if this routine is ever used in
+    ! production (see https://github.com/ESCOMP/CMEPS/issues/624)
     qsat(Tk)   = 640380.0_R8 / exp(5107.4_R8/Tk)
     cdn(Umps)  =   0.0027_R8 / Umps + 0.000142_R8 + 0.0000764_R8 * Umps
     psimhu(xd) = log((1.0_R8+xd*(2.0_R8+xd))*(1.0_R8+xd*xd)/8.0_R8) - 2.0_R8*atan(xd) + 1.571_R8
@@ -354,10 +354,9 @@ subroutine flux_atmOcn_diurnal(                      &
              speed(n)   = 0.0_R8
           endif
 
-          ! This should be changed to use the subroutine below
+          ! This should be changed to use shr_wv_sat_qsat_liquid (see
+          ! https://github.com/ESCOMP/CMEPS/issues/624)
           ssq = 0.98_R8 * qsat(tBulk(n)) / rbot(n)   ! sea surf hum (kg/kg)
-          ! call shr_wv_sat_qsat_liquid(tBulk(n), pslv(n), qsat, ssq)
-          ! ssq    = 0.98_R8 * ssq                   ! sea surf hum (kg/kg)
 
           delt   = thbot(n) - tBulk(n)               ! pot temp diff (K)
           delq   = qbot(n) - ssq                     ! spec hum dif (kg/kg)
@@ -503,10 +502,9 @@ subroutine flux_atmOcn_diurnal(                      &
 
              !--need to update ssq,delt,delq as function of tBulk ----
 
-             ! This should be changed to use the subroutine below
+             ! This should be changed to use shr_wv_sat_qsat_liquid (see
+             ! https://github.com/ESCOMP/CMEPS/issues/624)
              ssq = 0.98_R8 * qsat(tBulk(n)) / rbot(n)      ! sea surf hum (kg/kg)
-             ! call shr_wv_sat_qsat_liquid(tBulk(n), pslv(n), qsat, ssq)
-             ! ssq    = 0.98_R8 * ssq                      ! sea surf hum (kg/kg)
 
              delt   = thbot(n) - tBulk(n)                  ! pot temp diff (K)
              delq   = qbot(n) - ssq                        ! spec hum dif (kg/kg)