Fix expression merge with outer join for disjoint coordinates (#514)

* fix: explicitly set outer join for expression merge

Set join="outer" when merging expressions without common coordinate
structure to ensure consistent behavior across xarray versions and
properly handle expressions with disjoint coordinates.

* doc: add release note for outer join fix in expression merge

Author: FabianHofmann

doc/release_notes.rst

@@ -9,6 +9,7 @@ Upcoming Version
 * Improved constraint equality check in `linopy.testing.assert_conequal` to less strict optionally
 * Minor bugfix for multiplying variables with numpy type constants
 * Harmonize dtypes before concatenation in lp file writing to avoid dtype mismatch errors. This error occurred when creating and storing models in netcdf format using windows machines and loading and solving them on linux machines.
+* Fix expression merge to explicitly use outer join when combining expressions with disjoint coordinates for consistent behavior across xarray versions
 
 Version 0.5.6
 --------------

linopy/expressions.py

@@ -1999,6 +1999,8 @@ def merge(
 
         if override:
             kwargs["join"] = "override"
+        else:
+            kwargs.setdefault("join", "outer")
 
     if dim == TERM_DIM:
         ds = xr.concat([d[["coeffs", "vars"]] for d in data], dim, **kwargs)