Commits

John Wise committed c93b567

A similar fix to the metal cooling rate (MetalCooling = 1) in a
routine that calculates separate cooling rates for projections.

Comments (0)

Files changed (1)

src/enzo/cool1d_sep.F

                metalc(i) = metalc(i) - metalc_cmb
 
                xi = min(3.d0, metal(i,j,k)/(d(i,j,k)*ZSOLAR))
-               edot(i,18) = -metalc(i) * d(i,j,k) * 
-     $              (xi/JHW_METALS_NORMZ) / dom
+               edot(i,18) = -metalc(i) * d(i,j,k) * d(i,j,k)
+     $              (xi/JHW_METALS_NORMZ)
 #ifdef UNUSED
                if (xi .gt. 1e-06) then
                   write(6,777) i,j,k, exp(log_xe), exp(logtem(i)),