Implement Floe Size Distribution (FSD) in Icepack

columnphysics/icepack_intfc.F90

@@ -53,9 +53,16 @@ module icepack_intfc
       use icepack_itd, only: icepack_init_itd_hist
       use icepack_itd, only: icepack_aggregate
 
+      use icepack_fsd, only: icepack_init_fsd_bounds
+      use icepack_fsd, only: icepack_init_fsd
+      use icepack_fsd, only: icepack_cleanup_fsd
+
       use icepack_mechred, only: icepack_step_ridge
       use icepack_mechred, only: icepack_ice_strength
 
+      use icepack_wavefracspec, only: icepack_init_wave
+      use icepack_wavefracspec, only: icepack_step_wavefracture
+
       use icepack_shortwave, only: icepack_prep_radiation
       use icepack_shortwave, only: icepack_step_radiation
 

columnphysics/icepack_parameters.F90

@@ -266,6 +266,22 @@ module icepack_parameters
                          !   2 = WMO standard
                          !   3 = asymptotic formula
 
+!-----------------------------------------------------------------------
+! Parameters for the floe size distribution
+!-----------------------------------------------------------------------
+
+      integer (kind=int_kind), public :: &
+         nfreq = 25                   ! number of frequencies
+
+      real (kind=dbl_kind), public :: &
+         floeshape = 0.666_dbl_kind   ! constant from Steele (unitless)
+
+      logical (kind=log_kind), public :: &
+         wave_spec = .false.          ! if true, use wave forcing
+
+      character (len=char_len), public :: &
+         wave_spec_type = 'constant'  ! 'none', 'constant', or 'random'
+
 !-----------------------------------------------------------------------
 ! Parameters for melt ponds
 !-----------------------------------------------------------------------
@@ -369,6 +385,7 @@ subroutine icepack_init_parameters(   &
          kalg_in, kstrength_in, krdg_partic_in, krdg_redist_in, mu_rdg_in, &
          atmbndy_in, calc_strair_in, formdrag_in, highfreq_in, natmiter_in, &
          tfrz_option_in, kitd_in, kcatbound_in, hs0_in, frzpnd_in, &
+         floeshape_in, wave_spec_in, wave_spec_type_in, nfreq_in, &
          dpscale_in, rfracmin_in, rfracmax_in, pndaspect_in, hs1_in, hp1_in, &
          bgc_flux_type_in, z_tracers_in, scale_bgc_in, solve_zbgc_in, &
          modal_aero_in, skl_bgc_in, solve_zsal_in, grid_o_in, l_sk_in, &
@@ -567,6 +584,22 @@ subroutine icepack_init_parameters(   &
                             !   2 = WMO standard
                             !   3 = asymptotic formula
 
+!-----------------------------------------------------------------------
+! Parameters for the floe size distribution
+!-----------------------------------------------------------------------
+
+      integer (kind=int_kind), intent(in), optional :: &
+         nfreq_in           ! number of frequencies
+
+      real (kind=dbl_kind), intent(in), optional :: &
+         floeshape_in       ! constant from Steele (unitless)
+
+      logical (kind=log_kind), intent(in), optional :: &
+         wave_spec_in       ! if true, use wave forcing
+
+      character (len=char_len), intent(in), optional :: &
+         wave_spec_type_in  ! type of wave spectrum forcing 
+
 !-----------------------------------------------------------------------
 ! Parameters for biogeochemistry
 !-----------------------------------------------------------------------
@@ -745,6 +778,10 @@ subroutine icepack_init_parameters(   &
       if (present(tfrz_option_in)       ) tfrz_option      = tfrz_option_in
       if (present(kitd_in)              ) kitd             = kitd_in
       if (present(kcatbound_in)         ) kcatbound        = kcatbound_in
+      if (present(floeshape_in)         ) floeshape        = floeshape_in
+      if (present(wave_spec_in)         ) wave_spec        = wave_spec_in
+      if (present(wave_spec_type_in)    ) wave_spec_type   = wave_spec_type_in
+      if (present(nfreq_in)             ) nfreq            = nfreq_in
       if (present(hs0_in)               ) hs0              = hs0_in
       if (present(frzpnd_in)            ) frzpnd           = frzpnd_in
       if (present(dpscale_in)           ) dpscale          = dpscale_in
@@ -828,6 +865,7 @@ subroutine icepack_query_parameters(   &
          kalg_out, kstrength_out, krdg_partic_out, krdg_redist_out, mu_rdg_out, &
          atmbndy_out, calc_strair_out, formdrag_out, highfreq_out, natmiter_out, &
          tfrz_option_out, kitd_out, kcatbound_out, hs0_out, frzpnd_out, &
+         floeshape_out, wave_spec_out, wave_spec_type_out, nfreq_out, &
          dpscale_out, rfracmin_out, rfracmax_out, pndaspect_out, hs1_out, hp1_out, &
          bgc_flux_type_out, z_tracers_out, scale_bgc_out, solve_zbgc_out, &
          modal_aero_out, skl_bgc_out, solve_zsal_out, grid_o_out, l_sk_out, &
@@ -1035,6 +1073,22 @@ subroutine icepack_query_parameters(   &
                              !   2 = WMO standard
                              !   3 = asymptotic formula
 
+!-----------------------------------------------------------------------
+! Parameters for the floe size distribution
+!-----------------------------------------------------------------------
+
+      integer (kind=int_kind), intent(out), optional :: &
+         nfreq_out          ! number of frequencies
+
+      real (kind=dbl_kind), intent(out), optional :: &
+         floeshape_out      ! constant from Steele (unitless)
+
+      logical (kind=log_kind), intent(out), optional :: &
+         wave_spec_out      ! if true, use wave forcing
+
+      character (len=char_len), intent(out), optional :: &
+         wave_spec_type_out ! type of wave spectrum forcing
+
 !-----------------------------------------------------------------------
 ! Parameters for biogeochemistry
 !-----------------------------------------------------------------------
@@ -1254,6 +1308,10 @@ subroutine icepack_query_parameters(   &
       if (present(tfrz_option_out)       ) tfrz_option_out  = tfrz_option
       if (present(kitd_out)              ) kitd_out         = kitd
       if (present(kcatbound_out)         ) kcatbound_out    = kcatbound
+      if (present(floeshape_out)         ) floeshape_out    = floeshape
+      if (present(wave_spec_out)         ) wave_spec_out    = wave_spec
+      if (present(wave_spec_type_out)    ) wave_spec_type_out = wave_spec_type
+      if (present(nfreq_out)             ) nfreq_out        = nfreq
       if (present(hs0_out)               ) hs0_out          = hs0
       if (present(frzpnd_out)            ) frzpnd_out       = frzpnd
       if (present(dpscale_out)           ) dpscale_out      = dpscale
@@ -1423,6 +1481,10 @@ subroutine icepack_write_parameters(iounit)
         write(iounit,*) "  tfrz_option   = ", tfrz_option
         write(iounit,*) "  kitd          = ", kitd
         write(iounit,*) "  kcatbound     = ", kcatbound
+        write(iounit,*) "  floeshape     = ", floeshape
+        write(iounit,*) "  wave_spec     = ", wave_spec
+        write(iounit,*) "  wave_spec_type= ", wave_spec_type
+        write(iounit,*) "  nfreq         = ", nfreq
         write(iounit,*) "  hs0           = ", hs0
         write(iounit,*) "  frzpnd        = ", frzpnd
         write(iounit,*) "  dpscale       = ", dpscale