Update README.md

Author: gerlero

README.md

@@ -150,20 +150,21 @@ with my_pitz.fv_schemes as f:
     f["snGradSchemes"]["default"] = "uncorrected"
 ```
 
-### ⏳ Run multiple cases in parallel
-
-Outside an [asyncio](https://docs.python.org/3/library/asyncio.html) context:
+### 🔢 Direct field file access without `FoamCase`
 
 ```python
-from foamlib import AsyncFoamCase
-
-case1 = AsyncFoamCase("path/to/case1")
-case2 = AsyncFoamCase("path/to/case2")
+import numpy as np
+from foamlib import FoamFieldFile
 
-AsyncFoamCase.run_all_wait([case1, case2])
+# Read field data directly
+U = FoamFieldFile("0/U")
+print(f"Velocity field shape: {np.shape(U.internal_field)}")
+print(f"Boundaries: {list(U.boundary_field)}")
 ```
 
-Within an asyncio context (e.g. in a [Jupyter](https://jupyter.org/) notebook):
+### ⏳ Run multiple cases in parallel
+
+In an `asyncio` context (e.g. [`asyncio`](https://docs.python.org/3/library/asyncio.html) REPL or [Jupyter](https://jupyter.org/) notebook):
 
 ```python
 from foamlib import AsyncFoamCase
@@ -174,19 +175,9 @@ case2 = AsyncFoamCase("path/to/case2")
 await AsyncFoamCase.run_all([case1, case2])
 ```
 
-### 🔢 Direct field file access
-
-```python
-import numpy as np
-from foamlib import FoamFieldFile
-
-# Read field data directly
-U = FoamFieldFile("0/U")
-print(f"Velocity field shape: {np.shape(U.internal_field)}")
-print(f"Boundaries: {list(U.boundary_field)}")
-```
+Note: outside of an `asyncio` context, you can use [`asyncio.run()`](https://docs.python.org/3/library/asyncio-runner.html#asyncio.run).
 
-### 🎯 Optimization with HPC clusters
+### 🎯 Full optimization run on a Slurm-based HPC cluster
 
 ```python
 import os