diff --git a/docs/tech_guide/tecto.rst b/docs/tech_guide/tecto.rst
index 7d1c963..65de3e9 100644
--- a/docs/tech_guide/tecto.rst
+++ b/docs/tech_guide/tecto.rst
@@ -53,3 +53,7 @@ For regions of diverging or converging lateral movement, the elevation and strat
 .. note::
 
   Additional constrained for mesh refinement can be set if needed. As an example, one can use the horizontal velocity field magnitude as a parameter for refinement and create meshes with higher resolutions for regions with higher magnitudes in order to better represent the deformation in these areas.
+
+
+Flexural isostasy
+---------------------------------
\ No newline at end of file
diff --git a/gospl/tools/cgflex.py b/gospl/tools/cgflex.py
index 2baad19..c583990 100644
--- a/gospl/tools/cgflex.py
+++ b/gospl/tools/cgflex.py
@@ -3,16 +3,16 @@
 import sys
 import glob
 import pickle
-import pyproj
 import petsc4py
 import numpy as np
 
-import xarray as xr
 from mpi4py import MPI
 from pathlib import Path
 from time import process_time
 
 if "READTHEDOCS" not in os.environ:
+    import pyproj
+    import xarray as xr
     from gflex.f2d import F2D
 
 petsc4py.init(sys.argv)