From 65ec9214a909f3407152a27a5ebeb99c8cf09e0e Mon Sep 17 00:00:00 2001
From: David Bailey <dbailey@ucar.edu>
Date: Fri, 25 Oct 2024 16:34:24 -0600
Subject: [PATCH] Additional bug fix to floe_area_c

---
 columnphysics/icepack_fsd.F90 | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/columnphysics/icepack_fsd.F90 b/columnphysics/icepack_fsd.F90
index 9735090f..291ee6fb 100644
--- a/columnphysics/icepack_fsd.F90
+++ b/columnphysics/icepack_fsd.F90
@@ -194,8 +194,11 @@ subroutine icepack_init_fsd_bounds( &
       floe_rad_c = (floe_rad_h+floe_rad_l)/c2
 
       floe_area_l = c4*floeshape*floe_rad_l**2
-      floe_area_c = c4*floeshape*floe_rad_c**2
       floe_area_h = c4*floeshape*floe_rad_h**2
+!     floe_area_c = c4*floeshape*floe_rad_c**2
+!     This is exactly in the middle of floe_area_h and floe_area_l
+!     Whereas the above calculation is closer to floe_area_l.
+      floe_area_c = (floe_area_h+floe_area_l)/c2
 
       floe_binwidth = floe_rad_h - floe_rad_l