c c program difjk version 2 D. C. Griffin Last modified: 08/08/08 c c program to calculate the angular distribution of the differential c excitation cross section between individual levels for electron- c atom and electron-ion scattering in the jk coupling scheme. c implicit real*8(a-h,o-z) c parameter dimensions c c nlvm--maximum number of target levels c nlvfm--maximum number of final target levels to which differential c cross sections are to be calculated c ncm--maximum number of channels in k-matrices c nclm--maximum number of channels for a given level c nlm--maximum number of angular momenta for continuum electron c njm--maximum number of momentum transfer values c nwm--maximum number of partial waves actually used in calculation c nftm--maximum number of different factorials used in coupling c coefficient calculations - 800 should be enough c ntem-maximum number of different angles c parameter(nlvm=120,nlvfm=10,ncm=620,nclm=20,nlm=46,njm=2*nlm, 1 nlbm=2*nlm,nwm=102,nftm=800,ntem=181) character name(nlvm)*10,pltunt(5)*12 dimension title(10),fj(nlvm),enj(nlvm),lpt(nlvm) common/blok1/thmn,thmx,einc,theta(ntem),fji,fjf(nlvfm), 1 sigbri(nlm),sigbrf(nlm,nlvfm),nthta,iprint,icharg,j2mx,nlvf, 2 ilvi,ilvf(nlvfm),nthmx,limn,limx,lfmx(nlm,nlvfm), 3 lfmn(nlm,nlvfm),ielas,iphr,nipc,nfpc common/blok2/gamma(nftm) common/blok3/fjtot(nwm),tr(nwm,ncm,ncm),ti(nwm,ncm,ncm), 1 flcs(nwm,ncm),fkvs(nwm,ncm),ncci(nwm),nccf(nwm,nlvfm), 2 icci(nwm,nclm),iccf(nwm,nclm,nlvfm),iwu,imrd common/blok4/fmr(nlm,nlvfm,nlm,njm,4),fmi(nlm,nlvfm,nlm,njm,4), 1 qr(nlm),qi(nlm),alm(nlvfm,nlbm),dsg(ntem,nlvfm),sig(nlvfm), 2 sigic(nlvfm) common/blok5/pl(nlbm),imbmx namelist/difjk/iphr,iprint,icharg,j2mx,nlv,ilvi,nlvf,imrd,nthta, 1 thmn,thmx,ein,iunit data zero,threehf,two,onehndrd,one80/0.0d0,1.5d0,2.0d0,1.00d2, 1 1.80d2/ data convev,convsg,convsg1,convsg2,convsg3,convsg4/13.6058d0, 1 2.80028d-17,1.0d-19,1.0d-18,1.0d-17,1.0d-16/ data (pltunt(i),i=1,5) /' a0^2/sr ','10-19 cm2/sr', 1 '10-18 cm2/sr','10-17 cm2/sr','10-16 cm2/sr'/ open(unit=5,file='difjin') open(unit=6,file='difjout') open(unit=7,file='difjplt.dat') c c input from difjin c c title(10)--title for run - 80 characters c read(5,100) title 100 format(10a8) write(6,200) title 200 format(/80('*')//10a8//80('*')//) c c iphr--if input r-matrices were generated with r-matrix code, set c iphr equal to 1. this corrects the differential cross c section for the fano phase convention in which the spherical c harmonic includes the factor i**l. c (default:1) c iprint--if equal to 1, print k-matrices; if equal to 2, print c t-matrices; if equal to 3, print ql and a-lambda c coefficents; if equal to 4, print m array. c (default:0) c icharg--charge of target c j2mx--maximum value of 2xjtot to be included in the cross section c calculation c nlv--number of target levels included in input r-matrices c ilvi--level number of the initial level from list of target levels c must be in the range 1<=ilvi<=nlv c nlvf--number of final levels considered c imrd--set equal to zero to read k-matrices from unitarized c distorted-wave or close-coupling calculation c set equal to 1 to read t-matrices from non-unitarized c distorted-wave run from file tmatnu c (default:0) c nthta--number of angles for 0