We could optionally write the model references to a NetCDF file just before the solving process and reassign afterwards. This could reduce the total memory usage by 30%. (Example: for a 112 GB peak memory example solving, 31 GB was held by linopy and 81 GB by gurobi.)