def METHOD_NAME():
    sc = SpillContainer()
    new_subs = GnomeOil('oil_crude')
    new_subs.densities = [1004.0]
    new_subs.density_ref_temps = [288.15]
    sc.substance= new_subs
    sc.prepare_for_model_run(weathering_array_types)
    sc._append_data_arrays(100)
    w = Water()
    w.set('temperature', 288, 'K')
    w.set('salinity', 0, 'psu')
    init_density(sc, 100, water=w)
    assert np.all(sc['density'] == w.get('density'))