add d0 to integer quantities to ensure dp.

MESAHub · Sep 6, 2024 · a49ed93 · a49ed93
1 parent a7d6526
commit a49ed93
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/rates/private/rates_support.f90 b/rates/private/rates_support.f90
@@ -174,7 +174,7 @@ subroutine get_rates_from_table(r1, r2)
             rate_raw_dRho(i) = rate_raw(i) * (ddtab(i) / dtab(i))
 
             ! Derivative with respect to temperature (dT)
-            rate_raw_dT(i) = rate_raw(i) * (rattab_f(2,k,i) + 2*dt*rattab_f(3,k,i) + 3*dt**2*rattab_f(4,k,i)) / (btemp)
+            rate_raw_dT(i) = rate_raw(i) * (rattab_f(2,k,i) + 2*dt*rattab_f(3,k,i) + 3d0*dt**2*rattab_f(4,k,i)) / (btemp)
 
             end do
 
@@ -208,7 +208,8 @@ subroutine do_make_rate_tables( &
          include 'formats'
 
          ierr = 0
-         rate_logR = 0
+         rate_logR = 0d0
+
          rattab_f(1:4,1:nrattab,1:num_reactions) =>  &
                rattab_f1(1:4*nrattab*num_reactions)
 
@@ -284,7 +285,7 @@ subroutine do_make_rate_tables( &
                       if (rattab(j, i) > 0.0_dp) then  ! Only take log of positive values
                           rate_logR = log10(rattab(j, i))
                       else if (rattab(j, i) == 0.0_dp) then
-                          rate_logR = -999
+                          rate_logR = -999d0
                       else
                           ! warning for rates set to -1
                           ! like rni56ec_to_co56,rco56ec_to_fe56