diff --git a/earthdiagnostics/ocean/interpolatecdo.py b/earthdiagnostics/ocean/interpolatecdo.py
index 42dd17a99020e021b2296c43bae56a7ca06254cd..dcf8dd53991126666d621939ab0522136cd596bf 100644
--- a/earthdiagnostics/ocean/interpolatecdo.py
+++ b/earthdiagnostics/ocean/interpolatecdo.py
@@ -126,6 +126,8 @@ class InterpolateCDO(Diagnostic):
         :type options: list[str]
         :return:
         """
+        # Initialize CDO on the main thread
+        Utils.cdo()
         options_available = (
             DiagnosticDomainOption(default_value=ModelingRealms.ocean),
             DiagnosticVariableListOption(