update

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
+---------------------------------

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)