@@ -70,6 +70,13 @@ function ice_terminal_velocity_number_weighted(
7070 bnds = integral_bounds (state, logλ; p)
7171 return integrate (number_weighted_integrand, bnds... ; ∫kwargs... ) / N_ice
7272end
73+ function ice_terminal_velocity_number_weighted (
74+ velocity_params:: CMP.Chen2022VelType , ρₐ, params:: CMP.ParametersP3 , L_ice, N_ice, F_rim, ρ_rim, logλ;
75+ use_aspect_ratio = true , ∫kwargs... ,
76+ )
77+ state = get_state (params; L_ice, N_ice, F_rim, ρ_rim)
78+ return ice_terminal_velocity_number_weighted (velocity_params, ρₐ, state, logλ; use_aspect_ratio, ∫kwargs... )
79+ end
7380
7481"""
7582 ice_terminal_velocity_mass_weighted(velocity_params::CMP.Chen2022VelType, ρₐ, state::P3State, logλ; [use_aspect_ratio], [∫kwargs...])
@@ -107,3 +114,10 @@ function ice_terminal_velocity_mass_weighted(
107114 bnds = integral_bounds (state, logλ; p)
108115 return integrate (mass_weighted_integrand, bnds... ; ∫kwargs... ) / L_ice
109116end
117+ function ice_terminal_velocity_mass_weighted (
118+ velocity_params:: CMP.Chen2022VelType , ρₐ, params:: CMP.ParametersP3 , L_ice, N_ice, F_rim, ρ_rim, logλ;
119+ use_aspect_ratio = true , ∫kwargs... ,
120+ )
121+ state = get_state (params; L_ice, N_ice, F_rim, ρ_rim)
122+ return ice_terminal_velocity_mass_weighted (velocity_params, ρₐ, state, logλ; use_aspect_ratio, ∫kwargs... )
123+ end
0 commit comments