Additional bug fix to floe_area_c

CICE-Consortium · Oct 25, 2024 · 65ec921 · 65ec921
1 parent 1bff0fb
commit 65ec921
Showing 1 changed file with 4 additions and 1 deletion.
diff --git 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