1 subroutine electra(ro,t,ye,ee,pe,se 2 a ,eer,eet,per,pet,ser,set,keyps) 3 implicit real*8 (a-h,o-z) 4 common /nunu/ nu,dnudr,dnudb,eta,detadnu,nup 5 data facen,facpr,facs,rg /2.037300d+24,1.358200d+24,1.686304d-10 6 1,8.314339d+07/ 7 data a1,a2,a3,a4 /2.059815d-03,-7.027778d-03 8 1,4.219747d-02,-1.132427d+00/ 9 beta=facs*t 10 b32=b12*beta 11 u=(f62/f52)**2 12 dudnu=2.0d0*u*(df62/f62-df52/f52) 13 x=beta*u 14 f=1.0d0+x*(2.5d0+x*(2.0d0+0.5d0*x)) 15 df=2.5d0+x*(4.0d0+1.5d0*x) 16 dfdb=u*df 17 fi32=f32+(f-1.0d0)*f52/u 18 dfidnu=dfidu*dudnu+df32+(f-1.0d0)*df52/u 19 dfidb=dfdb*f52/u 20 dfidbet=dfidb+dfidnu*dnudb 21 gs=sqrt(g) 22 dg=0.75d0*gs 23 dgdb=u*dg 24 dgdu=beta*dg 25 gi32=f32+(g-1.0d0)*f52/u 26 dgidu=(u*dgdu-g+1.0d0)*f52/us 27 dgidnu=dgidu*dudnu+df32+(g-1.0d0)*df52/u 28 dgidb=dgdb*f52/u 29 dgidbet=dgidb+dgidnu*dnudb 30 dgidroe=dgidnu*dnudr 31 em=facen*b52*fi32 32 demdbet=facen*b32*(2.5d0*fi32+beta*dfidbet) 33 dpmdbet=facpr*b32*(2.5d0*gi32+beta*dgidbet) 34 demdroe=facen*b52*dfidroe 35 dpmdroe=facpr*b52*dgidroe 36 call divine(nup,fp12,dfp12,s12) 37 s42=2.0d0 38 call divine(nup,fp42,dfp42,s42) 39 eer=(ye*(demdroe+depdroe)-(em+ep)/ro)/ro 40 eet=facs*(demdbet+depdbet)/ro 41 per=ye*(dpmdroe+dppdroe) 42 pet=facs*(dpmdbet+dppdbet) 43 end 44