diff --git a/codes/coarse_k_point_implementation.ipynb b/analysis/coarse_k_point_implementation.ipynb
similarity index 100%
rename from codes/coarse_k_point_implementation.ipynb
rename to analysis/coarse_k_point_implementation.ipynb
diff --git a/analysis/codes b/analysis/codes
new file mode 120000
index 0000000000000000000000000000000000000000..f496e1c1ce26b6f6cb6504053420ebf88f11ea50
--- /dev/null
+++ b/analysis/codes
@@ -0,0 +1 @@
+../codes/
\ No newline at end of file
diff --git a/codes/data/graphene_example.nc b/analysis/data/graphene_example.nc
similarity index 100%
rename from codes/data/graphene_example.nc
rename to analysis/data/graphene_example.nc
diff --git a/codes/graphene_test_extended_hubbard.ipynb b/analysis/graphene_test_extended_hubbard.ipynb
similarity index 100%
rename from codes/graphene_test_extended_hubbard.ipynb
rename to analysis/graphene_test_extended_hubbard.ipynb
diff --git a/codes/mean_field_demo.ipynb b/analysis/mean_field_demo.ipynb
similarity index 100%
rename from codes/mean_field_demo.ipynb
rename to analysis/mean_field_demo.ipynb
diff --git a/codes/scipy_optimizer.ipynb b/analysis/scipy_optimizer.ipynb
similarity index 99%
rename from codes/scipy_optimizer.ipynb
rename to analysis/scipy_optimizer.ipynb
index fdf7394d66b86bffadb9efcd75e04e350c5b3468..83e465e86d2bff9835146550ecaca9e910ec0f1e 100644
--- a/codes/scipy_optimizer.ipynb
+++ b/analysis/scipy_optimizer.ipynb
@@ -308,11 +308,24 @@
   },
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 32,
    "id": "e17fc96c-c463-4e1f-8250-c254d761b92a",
    "metadata": {},
    "outputs": [],
-   "source": []
+   "source": [
+    "import xarray as xr\n",
+    "gap_da = xr.DataArray(data=gap, coords=dict(Us=Us, Vs=Vs))"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 33,
+   "id": "0cb395cd-84d1-49b4-89dd-da7a2d09c8d0",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "gap_da.to_netcdf('./data/graphene_example.nc')"
+   ]
   }
  ],
  "metadata": {