From rubin Wed Aug 14 10:22:51 1991 Received: by thyphy.physics.orst.edu (AIX 3.1/UCB 5.61/4.03) id AA16728; Wed, 14 Aug 91 10:22:51 -0700 Date: Wed, 14 Aug 91 10:22:51 -0700 From: rubin Message-Id: <9108141722.AA16728@thyphy.physics.orst.edu> To: lu Status: R Received: by thyphy.PHYSICS.ORST.EDU (AIX 3.1/UCB 5.61/4.03) id AA16046; Thu, 6 Dec 90 01:59:53 -0800 Date: Thu, 6 Dec 90 01:59:53 -0800 From: tim Message-Id: <9012060959.AA16046@thyphy.PHYSICS.ORST.EDU> To: rubin Differences between old and new. -------------------------------------- getins 61,62c61,62 < if (xnang .eq. 0) then < xnang = 3 --- > if (xnang .eq. 0.) then > xnang = 3. 90,94c90,94 < write(10,741) nr,lxmax < write(10,770) Tlab < write(10,750) ngp, kode, b, xnang, ymin1, ymin2 < write(10,880) achp, acmp, wsp, achn, acmn, wsn < write(10,830) nz, na, (nifty(n),n=1,20) --- > c write(10,741) nr,lxmax > c write(10,770) Tlab > c write(10,750) ngp, kode, b, xnang, ymin1, ymin2 > c write(10,880) achp, acmp, wsp, achn, acmn, wsn > c write(10,830) nz, na, (nifty(n),n=1,20) 101c101 < write(10,910) nes, nwaves --- > c write(10,910) nes, nwaves 135,137c135,137 < 741 format (1h0,5i4) < 750 format (1h ,i4,i3,e10.2,i4,2f5.2) < 770 format (1h ,f10.4) --- > c741 format (1h0,5i4) > c750 format (1h ,i4,i3,e10.2,f8.2,2f5.2) > c770 format (1h ,f10.4) 139,141c139,141 < 830 format (1h ,2i5,1x,i2,i2,8i1,10i3) < 880 format (1h ,8f10.4) < 910 format (' nes, nwaves = ',2i5) --- > c830 format (1h ,2i5,1x,i2,i2,8i1,10i3) > c880 format (1h ,8f10.4) > c910 format (' nes, nwaves = ',2i5) --------------------------------------- legpol ---------------------------------- lpotp 73c73 < c nifty(16) = 0, modified aussian rho and rho**2 & for p* z=n --- > c nifty(16) = 0, modified Gaussian rho and rho**2 & for p* z=n 162c162 < data hbarc, pi /197.3289e+0, 3.141593e+0/ --- > c data hbarc, pi /197.3289e+0, 3.141593e+0/ 165c165,166 < --- > hbarc = 197.3289e+0 > pi = 3.141593e+0 229c230 < call GAUSS2 (ngp, kode, 0.0d0, a, gp, wt) --- > call GAUSS2 (ngp, kode, a, b, gp, wt) 578c579 < 227 format(i3, 8e15.6 / 200(3x, 8e15.6 / )) --- > c227 format(i3, 8e15.6 / 200(3x, 8e15.6 / )) 586c587 < 820 format (47h R matrix calculated with ms series up to order,i3) --- > c820 format (47h R matrix calculated with ms series up to order,i3) 588,589c589,590 < 928 format(1h0,33h re/im f matrix before inversion ) < 930 format (1h ,26h normalized r(k,k..k)= ,24x,2e14.7) --- > c928 format(1h0,33h re/im f matrix before inversion ) > c930 format (1h ,26h normalized r(k,k..k)= ,24x,2e14.7) 597,599c598,600 < 1010 format (1h+,85x,4heta=,f10.4,8h delta=,f10.4) < 1020 format (41h0********** eta greater than 1 **********,//) < 1030 format (1h ,15h t(k,k,k)= ,35x,4e14.7) --- > c1010 format (1h+,85x,4heta=,f10.4,8h delta=,f10.4) > c1020 format (41h0********** eta greater than 1 **********,//) > c1030 format (1h ,15h t(k,k,k)= ,35x,4e14.7) 601,603c602,604 < 1040 format (38h now with exact coulomb after matching) < 1050 format (20x,40hcoulomb phase modified nuclear amplitude/16h t(k, < 1k.k) = ,4e14.7) --- > c1040 format (38h now with exact coulomb after matching) > c1050 format (20x,40hcoulomb phase modified nuclear amplitude/16h t(k, > c 1k.k) = ,4e14.7) 605c606 < 1080 format (1h ,8x,23h born approximation ) --- > c1080 format (1h ,8x,23h born approximation ) 617,619c618,620 < 1093 format(1h ,f7.0,2f8.3,3e11.3,f8.3,2e11.3,3f9.3,e11.3) < 1100 format(1h ,f7.0,f8.3,2e11.3,2f8.1,f10.3,2e11.3,2f9.3,e12.4 < 1,e10.2) --- > c1093 format(1h ,f7.0,2f8.3,3e11.3,f8.3,2e11.3,3f9.3,e11.3) > c1100 format(1h ,f7.0,f8.3,2e11.3,2f8.1,f10.3,2e11.3,2f9.3,e12.4 > c 1,e10.2) 622,626c623,627 < 1120 format (1h ,18x,10e11.3) < 1498 format(' ???main ur(lid,lid,1)= ',e10.4,' ui(same)= 'e10.4) < 1921 format(1h ,7h sigot=,e13.5,7h sigit=,e13.5, < 27h sig2t=,e13.5,2he=,f8.2,5hxgam=,e13.5) < 1979 format(1h ,17h ur ui(k0,k0,k0)= ,2e13.5) --- > c1120 format (1h ,18x,10e11.3) > c1498 format(' ???main ur(lid,lid,1)= ',e10.4,' ui(same)= 'e10.4) > c1921 format(1h ,7h sigot=,e13.5,7h sigit=,e13.5, > c 27h sig2t=,e13.5,2he=,f8.2,5hxgam=,e13.5) > c1979 format(1h ,17h ur ui(k0,k0,k0)= ,2e13.5) 628,630c629,631 < 3335 format(1h ,19h tr ti (nspin=3-6)=,8e13.5) < 4000 format(' ', f7.2, 3e11.3) < 4010 format(1h ,f7.2,2e11.3) --- > c3335 format(1h ,19h tr ti (nspin=3-6)=,8e13.5) > c4000 format(' ', f7.2, 3e11.3) > c4010 format(1h ,f7.2,2e11.3) ----------------------------------- matchc 54c54 < return --- > c return 58,59c58,60 < 1159 format (1h ,6hldum =,I10.4,6hlcmax=,I10.4,6h eta=,f10.4) < 1211 format (1h ,20h%%%%%%%%%%%%%%rhoc=,f10.4) --- > c1159 format (1h ,6hldum =,I10,6hlcmax=,I10,6h eta=,f10.4) > c1211 format (1h ,20h% %%%%%%%%%%%%%rhoc=,f10.4) > return -------------------------------- optp 456,458c456,458 < 350 format(' C13 neutron spin ff for k = kp = ',e15.5/10(10x,6e15.6/)) < 400 format (3x, 2e19.4, 22x, e19.4) < 401 format(1x, i2, 19x, e19.4, 41x, e19.4) --- > c350 format(' C13 neutron spin ff for k = kp = ',e15.5/10(10x,6e15.6/)) > c400 format (3x, 2e19.4, 22x, e19.4) > c401 format(1x, i2, 19x, e19.4, 41x, e19.4) 464,470c464,470 < 1001 format(' Vll(ldum,1,1)=',e12.6,' Vll(ldum,2,1)='e12.6) < 1002 format(' ',e10.4,' ',e10.4,' 'e10.4) < 1005 format(' Vll(ldum,1,6)=',e12.6,' Vll(1,2,6)='e12.6) < 1010 format(' ???2 lid=',i3,' lpia=',i3) < 1015 format(' f(1,1)= ',e10.4,' f(2,1)= ',e10.4) < 1200 format(' ???3 lpia= ',i3,'lborn= ',i3) < 1210 format(' Ur(lid,lid,1)=',e13.7,' Ui(lid,lid,1)=',e13.7) --- > c1001 format(' Vll(ldum,1,1)=',e12.6,' Vll(ldum,2,1)='e12.6) > 1002 format(' ',e10.4,' ',e10.4,' ',e10.4) > 1005 format(' Vll(ldum,1,6)=',e12.6,' Vll(1,2,6)=',e12.6) > c1010 format(' ???2 lid=',i3,' lpia=',i3) > c1015 format(' f(1,1)= ',e10.4,' f(2,1)= ',e10.4) > c1200 format(' ???3 lpia= ',i3,'lborn= ',i3) > c1210 format(' Ur(lid,lid,1)=',e13.7,' Ui(lid,lid,1)=',e13.7) ----------------------------------- vcoul 7c7 < c all units are mev-1 --- > c all units are mev-1???the unit of rcut is fm 18c18 < if (q.gt.1.) go to 30 --- > c if (q.gt.1.) go to 30 20,22c20,22 < vcoul = rcut*rcut*zalpha/4./pi/pi < go to 50 < 30 continue --- > c vcoul = rcut*rcut*zalpha/4./pi/pi > c go to 50 > c30 continue 36c36,41 < 45 vcoul = (0.5/pi/pi)*(zalpha/q/q)*(ff(1)-cos(q*rcut)) --- > CCCCCCCCCCCCCCCC > if (q.gt.1.) go to 45 > vcoul = (0.25/pi/pi)*zalpha*ff(1)*(rcut/hbarc)**2 > go to 50 > CCCCCCCCCCCCCCCC > 45 vcoul = (0.5/pi/pi)*(zalpha/q/q)*ff(1)*(1.-cos(q*rcut/hbarc)) ------------------------------------- voptth.f 13c13 < dimension v10(2), v1m1(2), v01(2), v11(2), v00(2), vss(2) --- > c dimension v10(2), v1m1(2), v01(2), v11(2), v00(2), vss(2) 103c103 < v11a = apb(ja) * nz * ff(1) + apb(je) * nn * ff(2) --- > v11a = apb(ja) * nz * ff(1) + apb(je) * nn * ff(2)+vc 106c106 < v11(j) = aovera * (v11a + s2th2 * v11b + c2th2 * v11c)+vc --- > v11(j) = aovera * (v11a + s2th2 * v11b + c2th2 * v11c) 110,111c110,111 < v00(j) = aovera * (v11a + v1m1a + x * (v11b - v11c))+vc < vss(j) = aovera * (v11a - v1m1a - v11b - v11c)+vc --- > v00(j) = aovera * (v11a + v1m1a + x * (v11b - v11c)) > vss(j) = aovera * (v11a - v1m1a - v11b - v11c) -------------------------------------- xsects