1 SUBROUTINE grad(gradbs1, gradbs2, gradbs3, gradbs4, irho)
19 REAL(rprec) :: gradbs1, gradbs2, gradbs3, gradbs4
23 REAL(rprec) :: gradbs, de, de1, di, di1, p
26 gradbs =-2.5_dp*gpsi(irho)*qsafety(irho)*betar(irho)*sign_jacobian
29 di = dense(irho)/zeff1
31 p = de*tempe1(irho) + di*tempi1(irho) + 1.e-36_dp
32 gradbs1 = gradbs*de1*tempe1(irho)/p
33 gradbs2 = gradbs*di1*tempi1(irho)/p
35 gradbs3 = alphae*gradbs*temperho1*de/p
37 gradbs4 = alphai*gradbs*tempirho1*di/p