=== Header part === L1 Title. Not read in. L2- You can put the information on the external file(s). Format: "A1,I3,1X,A8,2X,A50" The 1st row is usually a blank, then you can put the unit # (I3), status (A8), and the file path (A50). If you put any character/number on the 1st row, the line is disregarded. If the unit # is larger than 99, the reading procedure for external file(s) information is terminated. ---- INPUT ---- This line has no meaning in the code. === Main part === L1 : COMMENT [A50] COMMENT User's comment. Appears in outlist as it is. L2 : LIMFS IONS IFRM IMIR ICAL [5I5] LIMFS Limit of size of output file (unit # of KIBTBL) in unit of MB If 0, replaced by 1000000 (roughly 1TB) One line is estimated to be 128 bytes. IONS=0 No output for the configuration that is kinematically not allowed <>0 Output for the configuration that is kinematically not allowed IFRM=0 Kinematics/Observable output given in the Laboratory frame (L) =1 the c.m. frame (G) =2 the Projectile-rest frame (V) IMIR<>0 +z -> -z conversion in kinematics output ICAL=0 Kinematics survey mode; observables not calculated <>0 Observables calculated L3 : ZP AP ZA AA [F5.0,F10.0,F5.0,F10.0] ZP AP Atomic and mass #'s of the probe (particle 0) ZA AA Atomic and mass #'s of the nucleus A to be broken L4 : IKIN ELAB ICTREIN [I5,F10.0,I5] IKIN=0 Forward kinematics (probe is the projectile) <>0 Inverse kinematics (nucleus to be broken is the projectile) ELAB Kinetic energy per nucleon of the projectile ICTREIN=0 The total kinetic energy Ein of the projectile is calculated with Ein = ELAB*nint(AP) or ELAB*nint(AA). <>0 Ein = ELAB*AP or ELAB*AA is used. L5 : ISH EBIND ZSP ASP BETASP ICTRM [I5,F10.0,F5.0,2F10.0,I5] ISH=0 Depth of s.p. pot. is specified. =1 Depth of s.p. pot. is changed to reproduce the separation energy given by EBIND (>0). >9 S.p. wave function is read from unit ISH. EBIND Binding energy of the struck nucleon (when ISH=1) When ISH=0 and IBMC=0, regarded as the depth of the central s.p. pot. When ISH=0 and IBMC<>0, it has no meaning. ZSP Atomic # of the struck nucleon ASP Mass # of the struck nucleon BETASP Range of non-locality for the nucleon in A (Perey-Buck correction) If negative. the correction function is read from unit ISH (if ISH<10, process is terminated). ICTRM Mass parameter in the Sch. Eq. for the bound-state problem: =1 Mass of the particle to be struck 2 Reduced mass L6 : FJ FL SFAC NOD [2F5.0,F10.0,I5] FJ FL Total spin and orbital angular momentum of the s.p. state SFAC Spectroscopic factor (maximum=2j+1 for nucleon-knockout) NOD node (0 for the lowest state) L7 : IBMC RC ICTRC A0C RCL ICTRCL [I5,F10.0,I5,2F10.0,I5] IBMC=1 Bohr-Mottelson (BM) pot. used for central part <>1 Parameters read-in for central part RC Reduced radial parameter ICTRC Control of RC =0: fac=AA^1/3 1: fac=(AA-ASP)^1/3 2: fac=(AA-ASP)^1/3 + ASP^1/3 3: fac=1 the radial parameter is given by RC*fac. A0C Diffuseness parameter RCL Reduced Coulomb radius ICTRCL Control of RCL (similar to ICTRC) L8 : IBMS V0LS RS ICTRS AS [I5,2F10.0,I5,F10.0] IBMS=1 BM pot. used for spin-orbit part <>1 Parameters read-in for spin-orbit part V0LS Depth RS Reduced radial parameter ICTRS Control of RS (similar to ICTRC) AS Diffuseness parameter L9 : LMAX0 LMAX1 LMAX2 [3I5] LMAXi Maximum orbital angular momentum for particle i (=0,1,2) If negative, LMAXi=min( nint(K_i*R_max), |LMAXi| ) (automatic setting) L10: IVAR IEX FKNCUT IXUNT KUNT [2I5,F5.0,2I5] IVAR=1 T1 and Omega1 are controlled with IVVAR. =9 K_Bz in the A-rest frame is controlled with IVVAR and K_Bb in the A-rest frame is controlled with IVTHX. IVPHX, IVTH2, and IVPH2 have no meaning. >29 Kinematics profile (9 quantities) in L-frame read from unit IVAR (9F10.0); IVVAR, IVTHX, IVPHX,IVTH2, IVPH2 have no meaning. IEX<>1 Particle 1 is particle 0 =1 Particle 1 is the particle inside A in the initial state (applicable only when IPOT(2)=IPOT(3)=0 or IPOT(2)=IPOT(3)=1) FKNCUT Cutoff for the missing momentum (K_N) (in fm^-1) IXUNT Control for the unit of X-sec. in output =1 : micro barn <>1: milli barn KUNT Control for the unit of K in output =0: fm^-1 =1: MeV/c =2: GeV/c L11: IVVAR VARMIN VARMAX DVAR [I5,3F10.0] L12: IVTHX THXMIN THXMAX DTHX [I5,3F10.0] L13: IVPHX PHXMIN PHXMAX DPHX [I5,3F10.0] L14: IVTH2 TH2MIN TH2MAX DTH2 [I5,3F10.0] L15: IVPH2 PH2MIN PH2MAX DPH2 [I5,3F10.0] IVx=0 x is fixed at xMIN <>0 x varies from xMIN to xMAX with the increment Dx [NOTE: Units of energy (T), angle (TH and PH), and wave number (K) are MeV, degree, and fm^-1, respectively. These quantities must be given in the L-frame.] L16: KIB: TBL OUT TMD LG PX TR TL [7I5] KIBTBL Unit # for TDX table when IVAR<>9 Unit # for DDX when IVAR=9 KIBOUT Unit # for outlist file KIBTMD Unit # for transition matrix density (when > 0) Allowed only when single or no degree of freedom in kinematics is varied. KIBLG Unit # for longitudinal momentum distribution (when > 0) KIBPX Unit # for p_x momentum distribution (when > 0) if > 0, IVTHX=1 and nint(THXMAX-THXMIN/DTHX)+1 > 3 must be satisfied KIBTR Unit # for transverse momentum distribution (when > 0) KIBTL Unit # for total momentum distribution (when > 0) if > 0, IVTHX=1, nint(THXMAX-THXMIN/DTHX)+1 > 3, IVVAR=1, nint(VARMAX-VARMIN/DVAR)+1 > 3 must be satisfied [NOTE: KIBTBL and KIBOUT must be positive. KIBLG, KIBPX, KIBTR, KIBTL are effective only when IVAR=9. Trivially, all unit #s must be different from each other.] L17: IELM KIBELM IONSH KINELM IELMEDG [5I5] IELM=0 Isotropic free NN c.s. at ELAB (in unit of mb) =3 Free differential NN c.s. read from unit KIBELM (as a function of E and theta, see below) =4 Free (on-shell) NN t-matrix read from unit KIBELM Available only when reaction is in coplanar (see below) KIBELM Unit # of the external file for elementary process (needed unless IELM=0) IONSH Choice of on-shell approximation (needed when IELM=3,4) =1 Final-state prescription =2 Initial-state prescription =3 Energy-average prescription =4 momentum-average prescription KINELM=1 Output for the kinematics of the elementary process <>1 No output for the kinematics of the elementary process IELMEDG=0 If the scattering energy is out of the range prepared, the process is terminated. <>0 If the scattering energy is out of the range prepared, the value at the nearest energy (on an edge) is adopted. L18: RMAX DR NGR NGTH NGPH NGK1 NGPH1Q [2F10.0,5I5] RMAX DR Maximum and increment of the integration variable NG-x Mesh number of Gauss integration for x=R, theta_R, phi_R, K1, or phi_1Q NG for K1 and phi1Q effective only when IVAR=9. L19: IPOT0 FACV0 FACW0 FACVS0 FACWS0 BETA0 IMS0 IEDG0 [I5,4F5.0,F10.0,2I5] L20: IPOT1 FACV1 FACW1 FACVS1 FACWS1 BETA1 IMS1 IEDG1 [I5,4F5.0,F10.0,2I5] L21: IPOT2 FACV2 FACW2 FACVS2 FACWS2 BETA2 IMS2 IEDG2 [I5,4F5.0,F10.0,2I5] IPOTi=0 Plane wave is used for particle i (i=0-2). 1 Koning-Delaroche pot. (w/ Coulomb) is adopted for particle i. >9 Opt. pot. for particle i is read from unit IPOTi (must cover the energy range appropriately) FACVi V_i is multiplied by FACVi If negative, the Coulomb potential is turned off and FACVi = |FACVi| FACWi W_i is multiplied by FACWi FACVSi VS_i is multiplied by FACVSi FACWSi WS_i is multiplied by FACWSi BETAi Range of non-locality for particle i. If negative. the correction function is read from ext. file (if IPOTi<10, process is terminated). IMSi=0 Reduced energy is used as the kinetic factor. <>0 Reduced mass is used as the kinetic factor. IEDGi=0 If the scattering energy is out of the range prepared, the process is terminated. <>0 If the scattering energy is out of the range prepared, the value at the nearest energy (on an edge) is adopted. === External file for bound state wave function === L1 : RSPMAX DRSP [free] RSPMAX Maximum radius (must be not smaller than RMAX) DRSP Increment of radius (can be different from DR) L2-: RSP WFSP FNLSPW (for RSPMAX/DRSP+1 lines) [F10.0,2E20.12] RSP Radius WFSP Radial W.Fn. (not multiplied by R) FNLSPW Nonlocal correction factor (effective only when BETASP<0) === External file for optical potential === L1 : NEPMAXW RPMAXW DRPW R0CLW ICTRCLW [I5,3F10.0,I5] NEPMAXW # of energy points RPMAXW Maximum radius (must be not smaller than RMAX) DRPW Increment of radius (can be different from DR) R0CLW Reduced Coulomb radius If 0 and the Coulomb potential (see below) is not given, Koning-Delaroche's Coulomb radius is adopted. ICTRCLW Control of R0CLW =1: fac=(AX)^1/3 2: fac=(ASC)^1/3 + AX^1/3 3: fac=1 AX : AA (for particle 0) or AB (for particle 1 and 2) ASC: Mass # of particle 0, 1, or 2 the radial parameter is given by R0CLW*fac. L2 : EPW [free] L3-: VPOTW WPOTW VLSW WLSW VCOUW FNLWR FNLWI (for RPMAXW/DRPW+1 lines) [7E20.12] VPOTW Real central WPOTW Imaginary central VLSW Real spin-orbit WLSW Imaginary spin-orbit VCOUW Coulomb potential FNLWR Real nonlocal correction factor (effective only when BETA<0) FNLWI Imaginary nonlocal correction factor (effective only when BETA<0) Repeat "L2 L3-" for NEPMAXW times. EPW must be sorted in ascending order. === External file for elementary process (IELM=3) === L1 : NENNMAX THNNMAX DTHNN [free] NENNMAX # of energy points THNNMAX Maximum NN scattering angle (in degrees) DTHNN Increment L2 : EELM [free] EELM NN scattering energy (in NN laboratory frame) L3-: DSIGPP DSIGPN (for THNNMAX/DTHNN+1 lines) [free] DSIGPP pp differential X-sec. (mb/sr) DSIGPN pn differential X-sec. (mb/sr) Repeat "L2 L3-" for NENNMAX times. EELM must be sorted in ascending order. === External file for NN scattering amplitude from KIBELM (IELM=4) === L1 : q0mn q0mx dq0 [free] q0mn Minimum value of the relative wave number q0mx Maximum value of the relative wave number dq0 increment of the wave number L2 : th0mn th0mx thmn thmx dth [free] th0mn Minimum polar angle of the relative wave number in the initial state th0mx Maximum polar angle of the relative wave number in the initial state thmn Minimum polar angle of the relative wave number in the final state thmx Maximum polar angle of the relative wave number in the final state dth increment of the polar angle (for both the initial and the final states) [NOTE: IELM=4 is allowed only for reaction kinematics in coplanar. NN amplitudes with phi = 0, 180, and 360 deg. (azimuthal angle of the relative momentum in the initial and the final states) have to be given. ] L3- : Re[t_{NN}] Im[t_{NN}] [free] Data are read as follows. They should be prepared accordingly. do iq = 1,nqmx (q0mn -- q0mx) do itht0 = 1,ntht0 (th0mn -- th0mx) do iphi0 = 1,nphi0 (For 0, 180, 360 deg.) do itht = 1,ntht (thmn -- thmx) do iphi = 1,nphi (For 0, 180, 360 deg.) do ispin_1=1,2 do ispin_sp=1,2 do ispin_2=1,2 do ispin_3=1,2 read(kibelm,*) Re[t_{NN}] Im[t_{NN}] end do end do end do end do end do end do end do end do end do spin_x=1 when particle x has spin up, =2 when down x = 1: incident particle =sp: struck particle (initial state) = 2: emitted particle in the final state (that used to be incident particle) = 3: emitted particle in the final state (that used to be struck particle) "L3" is repeated for nqmx*ntht0*3*ntht*3*16 times.