Merge pull request #614 from FESOM/fix/ib_fw_flux

Fix/ib fw flux
FESOM · Jul 22, 2024 · 7b801dc · 7b801dc
2 parents eecc017 + e98f4bc
commit 7b801dc
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/icb_coupling.F90 b/src/icb_coupling.F90
@@ -96,10 +96,10 @@ subroutine icb2fesom(mesh, partit, ice)
 
     do n=1, myDim_nod2d+eDim_nod2D
         if (.not.turn_off_hf) then
-                net_heat_flux(n)   = net_heat_flux(n) + ibhf(n) 
+                net_heat_flux(n)   = net_heat_flux(n) + ibhf(n) * steps_per_ib_step 
         end if
         if (.not.turn_off_fw) then
-                fresh_wa_flux(n)   = fresh_wa_flux(n) + (ibfwb(n)+ibfwl(n)+ibfwe(n)+ibfwbv(n))
+                fresh_wa_flux(n)   = fresh_wa_flux(n) + (ibfwb(n)+ibfwl(n)+ibfwe(n)+ibfwbv(n)) * steps_per_ib_step
         end if
     end do
 !---wiso-code-begin

diff --git a/src/io_meandata.F90 b/src/io_meandata.F90
@@ -622,7 +622,7 @@ subroutine ini_mean_io(ice, dynamics, tracers, partit, mesh)
     call def_stream(nod2D, myDim_nod2D, 'ibfwbv',  'basal iceberg melting',            'm/s',    ibfwbv(:),        1, 'm', i_real4, partit, mesh)
     call def_stream(nod2D, myDim_nod2D, 'ibfwl',   'lateral iceberg melting',          'm/s',    ibfwl(:),         1, 'm', i_real4, partit, mesh)
     call def_stream(nod2D, myDim_nod2D, 'ibfwe',   'iceberg erosion',                  'm/s',    ibfwe(:),         1, 'm', i_real4, partit, mesh)
-    call def_stream(nod2D, myDim_nod2D, 'ibhf',    'heat flux from iceberg melting',   'm/s',    ibhf(:),          1, 'm', i_real4, partit, mesh)
+    call def_stream(nod2D, myDim_nod2D, 'ibhf',    'heat flux from iceberg melting',   'W/m2',   ibhf(:),          1, 'm', i_real4, partit, mesh)
   end if
 !------------------------------------------
 !_______________________________________________________________________________