def roll_da_to_0360(xr_da, lon_name="longitude"):
    """Roll DataArray's data with longitudes from (-180, 180) to (0, 360)."""
    # Roll longitudes and corresponding data
    out = xr_da.roll(**{lon_name: xr_da[lon_name].shape[0] // 2}, roll_coords=True)
    # Reset western (negative) longitudes to values within (180, 360) range
    out[lon_name] = out[lon_name] - (out[lon_name] // 360) * 360
    return out