def H_fr_frbd(freeboard, rho_sw, rho_i):
    """
    Compute iceberg total height based on freeboard and seawater and ice density
    """
    berg_H = freeboard * (rho_sw/(rho_sw-rho_i))
    return berg_H