Adds meltw (water flux from seaice and snow)

The corresponding variable from CIME is *Fioi_meltw
NOAA-GFDL · Oct 12, 2018 · 918f218 · 918f218
1 parent 67e70c5
commit 918f218
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 0 deletions.
diff --git a/config_src/mct_driver/MOM_surface_forcing.F90 b/config_src/mct_driver/MOM_surface_forcing.F90
@@ -158,6 +158,7 @@ module MOM_surface_forcing
   real, pointer, dimension(:,:) :: v_flux          =>NULL() !< j-direction wind stress (Pa)
   real, pointer, dimension(:,:) :: t_flux          =>NULL() !< sensible heat flux (W/m2)
   real, pointer, dimension(:,:) :: melth           =>NULL() !< sea ice and snow melt heat flux (W/m2)
+  real, pointer, dimension(:,:) :: meltw           =>NULL() !< water flux due to sea ice and snow melting (kg/m2/s)
   real, pointer, dimension(:,:) :: q_flux          =>NULL() !< specific humidity flux (kg/m2/s)
   real, pointer, dimension(:,:) :: salt_flux       =>NULL() !< salt flux (kg/m2/s)
   real, pointer, dimension(:,:) :: lw_flux         =>NULL() !< long wave radiation (W/m2)
@@ -446,6 +447,10 @@ subroutine convert_IOB_to_fluxes(IOB, fluxes, Time, G, CS, &
     if (associated(fluxes%melth)) &
       fluxes%melth(i,j) = G%mask2dT(i,j) * IOB%melth(i-i0,j-j0)
 
+    ! water flux due to sea ice and snow melt (kg/m2/s)
+    if (associated(fluxes%meltw)) &
+      fluxes%meltw(i,j) = G%mask2dT(i,j) * IOB%meltw(i-i0,j-j0)
+
     ! latent heat flux (W/m^2)
     if (associated(fluxes%latent)) &
       fluxes%latent(i,j) = G%mask2dT(i,j) * IOB%latent_flux(i-i0,j-j0)
@@ -776,6 +781,7 @@ subroutine IOB_allocate(IOB, isc, iec, jsc, jec)
              IOB% v_flux (isc:iec,jsc:jec),          &
              IOB% t_flux (isc:iec,jsc:jec),          &
              IOB% melth (isc:iec,jsc:jec),          &
+             IOB% meltw (isc:iec,jsc:jec),          &
              IOB% q_flux (isc:iec,jsc:jec),          &
              IOB% salt_flux (isc:iec,jsc:jec),       &
              IOB% lw_flux (isc:iec,jsc:jec),         &
@@ -801,6 +807,7 @@ subroutine IOB_allocate(IOB, isc, iec, jsc, jec)
   IOB%v_flux          = 0.0
   IOB%t_flux          = 0.0
   IOB%melth           = 0.0
+  IOB%meltw           = 0.0
   IOB%q_flux          = 0.0
   IOB%salt_flux       = 0.0
   IOB%lw_flux         = 0.0
@@ -1316,6 +1323,7 @@ subroutine ice_ocn_bnd_type_chksum(id, timestep, iobt)
   write(outunit,100) 'iobt%v_flux         ', mpp_chksum( iobt%v_flux         )
   write(outunit,100) 'iobt%t_flux         ', mpp_chksum( iobt%t_flux         )
   write(outunit,100) 'iobt%melth          ', mpp_chksum( iobt%melth          )
+  write(outunit,100) 'iobt%meltw          ', mpp_chksum( iobt%meltw          )
   write(outunit,100) 'iobt%q_flux         ', mpp_chksum( iobt%q_flux         )
   write(outunit,100) 'iobt%rofl_flux      ', mpp_chksum( iobt%rofl_flux      )
   write(outunit,100) 'iobt%rofi_flux      ', mpp_chksum( iobt%rofi_flux      )

diff --git a/config_src/mct_driver/ocn_cap_methods.F90 b/config_src/mct_driver/ocn_cap_methods.F90
@@ -76,6 +76,9 @@ subroutine ocn_import(x2o, ind, grid, ice_ocean_boundary, ocean_public, logunit,
       ! snow&ice melt heat flux  (W/m^2)
       ice_ocean_boundary%melth(i,j) = x2o(ind%x2o_Fioi_melth,k)
 
+      ! water flux from snow&ice melt (kg/m2/s)
+      ice_ocean_boundary%meltw(i,j) = x2o(ind%x2o_Fioi_meltw,k)
+
       ! liquid runoff
       ice_ocean_boundary%rofl_flux(i,j) = x2o(ind%x2o_Foxx_rofl,k) * GRID%mask2dT(ig,jg)
 
@@ -120,6 +123,7 @@ subroutine ocn_import(x2o, ind, grid, ice_ocean_boundary, ocean_public, logunit,
         write(logunit,F01)'import: day, secs, j, i, q_flux          = ',day,secs,j,i,ice_ocean_boundary%q_flux(i,j)
         write(logunit,F01)'import: day, secs, j, i, t_flux          = ',day,secs,j,i,ice_ocean_boundary%t_flux(i,j)
         write(logunit,F01)'import: day, secs, j, i, melth           = ',day,secs,j,i,ice_ocean_boundary%melth(i,j)
+        write(logunit,F01)'import: day, secs, j, i, meltw           = ',day,secs,j,i,ice_ocean_boundary%meltw(i,j)
         write(logunit,F01)'import: day, secs, j, i, latent_flux     = ',&
                           day,secs,j,i,ice_ocean_boundary%latent_flux(i,j)
         write(logunit,F01)'import: day, secs, j, i, runoff          = ',&