Added missing factor in spin_doctor

repepo · Oct 16, 2023 · 1799d47 · 1799d47
1 parent 2f4a379
commit 1799d47
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/bin/spin_doctor.py b/bin/spin_doctor.py
@@ -180,13 +180,13 @@ def main(ncpus):
         if par.thermal:
 
             [ TE[i], Dthm0, Wadv_thm[i] ] = np.sum( tdgn, 0) 
-            Dthm[i] = Dthm0 * par.Etherm
+            Dthm[i] = Dthm0 * par.OmgTau * par.Etherm
 
 
         if par.compositional:
 
             [ CE[i], Dcmp0, Wadv_cmp[i] ] = np.sum( cdgn, 0)
-            Dcmp[i] = Dcmp0 * par.Ecomp
+            Dcmp[i] = Dcmp0 * par.OmgTau * par.Ecomp
 
 
         # --------------------------------------------------------- Computing residuals to check the power balance: