@@ -70,6 +70,13 @@ function ice_terminal_velocity_number_weighted(
70
70
bnds = integral_bounds (state, logλ; p)
71
71
return integrate (number_weighted_integrand, bnds... ; ∫kwargs... ) / N_ice
72
72
end
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
73
80
74
81
"""
75
82
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(
107
114
bnds = integral_bounds (state, logλ; p)
108
115
return integrate (mass_weighted_integrand, bnds... ; ∫kwargs... ) / L_ice
109
116
end
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