def add_Inorg_and_Org_totals2array(ds, InOrgVar='Inorg_Total', OrgVar='Org_Total'):
    """
    Add inorganic and organic sub totals to dataset
    """
    # Add aggregated values to ds
    OrgVars = [
        'EmisCH2IBr_Ocean', 'EmisCH2ICl_Ocean', 'EmisCH2I2_Ocean', 'EmisCH3I_Ocean',
    ]
    InOrgVars = ['EmisI2_Ocean', 'EmisHOI_Ocean', ]
    # - Inorganic
    # template off the first species
    ds[InOrgVar] = ds[InOrgVars[0]].copy()
    # sum values to this
    arr = ds[InOrgVar].values
    for var_ in InOrgVars[1:]:
        print(var_)
        arr = arr + ds[var_].values
    ds[InOrgVar].values = arr
    attrs = ds[InOrgVar].attrs
    attrs['long_name'] = InOrgVar
    ds[InOrgVar].attrs = attrs
    # - Organic
    # template off the first species
    ds[OrgVar] = ds[OrgVars[0]].copy()
    # sum values to this
    arr = ds[OrgVar].values
    for var_ in OrgVars[1:]:
        print(var_)
        arr = arr + ds[var_].values
    ds[OrgVar].values = arr
    attrs = ds[OrgVar].attrs
    attrs['long_name'] = OrgVar
    ds[OrgVar].attrs = attrs
    return ds