Ticket #84: SIGMA2.FOR
| File SIGMA2.FOR, 485 bytes (added by jhjcho, 12 years ago) |
|---|
| Line | |
|---|---|
| 1 | subroutine new_i0_term |
| 2 | |
| 3 | include 'frills_sources:function.inc' |
| 4 | |
| 5 | real*8 cal,dens,rho,rp,phi,gamma,moment,bkgd |
| 6 | real*8 a,b,c,d,e,f,g,h,pi |
| 7 | |
| 8 | cal=p(1) |
| 9 | dens=p(2) |
| 10 | rho=p(3) |
| 11 | rp=p(4) |
| 12 | phi=p(5) |
| 13 | gamma=p(6) |
| 14 | moment=p(7) |
| 15 | bkgd=p(8) |
| 16 | |
| 17 | pi=3.141592654 |
| 18 | |
| 19 | do j=1,nv |
| 20 | |
| 21 | a=x(j)*x(j) |
| 22 | b=moment*moment |
| 23 | c=rho*rho |
| 24 | d=gamma*gamma |
| 25 | e=dens*dens |
| 26 | f=cal*6.0E-22*pi*c*phi*d |
| 27 | g=rp*a*e |
| 28 | h=dexp(-1.*a*b) |
| 29 | ycal(j)=((f/g)*h)+bkgd |
| 30 | end do |
| 31 | |
| 32 | return |
| 33 | end |
