diff --git a/environment.yml b/environment.yml
index e4f8f1a74..c566bd845 100644
--- a/environment.yml
+++ b/environment.yml
@@ -16,7 +16,7 @@ dependencies:
   - pyopenfast>=4.1.0
   # - openraft>=2.0.2
   - osqp  
-  - pcrunch>=2.1.4
+  - pcrunch>=2.1.5
   - pip
   - pyhams>=1.3
   #- pyoptsparse
diff --git a/pyproject.toml b/pyproject.toml
index 9fe2625c5..f2ccc8f9c 100644
--- a/pyproject.toml
+++ b/pyproject.toml
@@ -51,12 +51,12 @@ dependencies = [
   "openpyxl",
   "openraft",
   "osqp",
-  "pcrunch",
+  "pcrunch>=2.1.5",
   "pyhams",
-  "rosco",
+  "rosco>=2.10.1",
   "scipy",
   "smt",
-  "wisdem",
+  "wisdem>4.0.4",
 ]
 
 # List additional groups of dependencies here (e.g. development
diff --git a/weis/aeroelasticse/openmdao_openfast.py b/weis/aeroelasticse/openmdao_openfast.py
index ad4283e58..85abd4b62 100644
--- a/weis/aeroelasticse/openmdao_openfast.py
+++ b/weis/aeroelasticse/openmdao_openfast.py
@@ -2861,6 +2861,7 @@ def run_FAST(self, inputs, discrete_inputs, fst_vt):
             lss_fatigue = FatigueParams(load2stress=1.0,
                                         dnv_name='B1',
                                         dnv_type='air',
+                                        units='kPa',
                                         ultimate_stress=1e-3*inputs['lss_ultstress'],
                                         S_intercept=1e-3*inputs['lss_wohlerA'])
             for s in ['Ax','Sh']:
@@ -2886,6 +2887,7 @@ def run_FAST(self, inputs, discrete_inputs, fst_vt):
             tower_fatigue_base = FatigueParams(load2stress=1.0,
                                                dnv_name='D',
                                                dnv_type='air',
+                                               units='kPa',
                                                ultimate_stress=1e-3*inputs['tower_ultstress'][0],
                                                S_intercept=1e-3*inputs['tower_wohlerA'][0])
             for s in ['Ax','Sh']:
@@ -2905,6 +2907,7 @@ def run_FAST(self, inputs, discrete_inputs, fst_vt):
                 monopile_fatigue_base = FatigueParams(load2stress=1.0,
                                                       dnv_name='D',
                                                       dnv_type='sea',
+                                                      units='Pa',
                                                       ultimate_stress=inputs['monopile_ultstress'][0],
                                                       S_intercept=inputs['monopile_wohlerA'][0])
                 for s in ['Ax','Sh']: