double precision kfac /2.4d0/
double precision x,x2,x3,x4,x5,licznik /1.12692801104297249644/
c facT=temp0/t_bath
double precision kfac /2.4d0/
double precision x,x2,x3,x4,x5,licznik /1.12692801104297249644/
c facT=temp0/t_bath
&(1.0d0/(-dsqrt(1.0d0-sinphi*sinthet))*(
& sinphi/sinthet*costhet*costhet_grad(j)
& +sinthet/sinphi*cosphi*cosphi_grad_long(j)))
&(1.0d0/(-dsqrt(1.0d0-sinphi*sinthet))*(
& sinphi/sinthet*costhet*costhet_grad(j)
& +sinthet/sinphi*cosphi*cosphi_grad_long(j)))
C grad_shield_side is Cbeta sidechain gradient
grad_shield_side(j,ishield_list(i),i)=
& (sh_frac_dist_grad(j)*-2.0d0
C grad_shield_side is Cbeta sidechain gradient
grad_shield_side(j,ishield_list(i),i)=
& (sh_frac_dist_grad(j)*-2.0d0
&(1.0d0/(-dsqrt(1.0d0-sinphi*sinthet))*(
& sinphi/sinthet*costhet*costhet_grad(j)
& +sinthet/sinphi*cosphi*cosphi_grad_long(j)))
&(1.0d0/(-dsqrt(1.0d0-sinphi*sinthet))*(
& sinphi/sinthet*costhet*costhet_grad(j)
& +sinthet/sinphi*cosphi*cosphi_grad_long(j)))
grad_shield_loc(j,ishield_list(i),i)=
& scale_fac_dist*VSolvSphere/VSolvSphere_div/2.0d0*
&(1.0d0/(dsqrt(1.0d0-sinphi*sinthet))*(
& sinthet/sinphi*cosphi*cosphi_grad_loc(j)
& ))
grad_shield_loc(j,ishield_list(i),i)=
& scale_fac_dist*VSolvSphere/VSolvSphere_div/2.0d0*
&(1.0d0/(dsqrt(1.0d0-sinphi*sinthet))*(
& sinthet/sinphi*cosphi*cosphi_grad_loc(j)
& ))
C write(2,*) "TOTAL VOLUME",i,VolumeTotal,fac_shield(i)
enddo
return
C write(2,*) "TOTAL VOLUME",i,VolumeTotal,fac_shield(i)
enddo
return