* *=== ncdnvp ===========================================================* * *----------------------------------------------------------------------* * * * Copyright (C) 2003-2019: CERN & INFN * * All Rights Reserved. * * * * Used only by PReequilibrium INitialization routines * * * * Created on 22 June 1994 by Alfredo Ferrari & Paola Sala * * Infn - Milan * * * * Rhapsf (i,l) = a (skin) parameter for the density distribution, * * or s-shell armonic oscillator constant (a_s) for * * nuclei using shell model density distributions * * already including (possible) folding distributions* * Rhcpsf (i,l) = c (R1/2) parameter for the density distribution, * * or p-shell armonic oscillator constant (a_p) for * * nuclei using shell model density distributions * * already including (possible) folding distributions* * Rh0psf (i,l) = Rh0 (central density) parameter for the density * * distribution, or p-shell armonic oscillator * * constant (a_p) for nuclei using shell model dens- * * ity distributions NOT including (possible) fold- * * ing distributions * * Rhcesf (i,l) = Rhocen, central density for the p/n density * * distribution * * Rhrdsf (i,l) = radius at which the i,l density distribution drops* * at Rhfrsf (i,l) times the central value * * Rrmssf (i,l) = r-rms () parameter for the density distri- * * bution * * Rreqsf (i,l) = equivalent sharp sphere radius for the density * * distributions * * ( i = 1 : proton, i = 2 : neutron * * l = 1 charge (matter) proton (neutron) distribu- * * tion, integral = Z (N) * * l = 2 charge (matter) proton (neutron) distribu- * * tion, folded with folrms, integral = Z (N) * * l = 3 charge (matter) proton (neutron) distribu- * * tion, folded with folrms, integral = Z (N) * * l = 4 charge (matter) proton (neutron) distribu- * * tion, folded with folrms, integral = Z (N) * * l = 5 point/centre proton (neutron) distribution, * * integral = Z (N) * * l = 6 proton (neutron) distribution, folded for * * pion potential, integral = Z (N) ) * * Sshlnc (i) = Number of nucleons of type i in the s shell (for * * shell model like calculations, i = 1 : proton, * * i = 2 : neutron ) * * Pshlnc (i) = Number of nucleons of type i in the p shell (for * * shell model like calculations, i = 1 : proton, * * i = 2 : neutron ) * * Vpapsf (i,k) = a (skin) parameter for the potential distribution * * or s-shell armonic oscillator constant (a_s) for * * nuclei using shell model density like potential * * distributions already including (possible) fold- * * ing distributions * * Vpcpsf (i,k) = c (R1/2) parameter for the potential distribution * * or p-shell armonic oscillator constant (a_p) for * * nuclei using shell model density like potential * * distributions already including (possible) fold- * * ing distributions * * Vp0psf (i,k) = Vp0 (central potential) parameter for the poten- * * tial distribution distribution, or p-shell armo- * * nic oscillator constant (a_p) for nuclei using * * shell model density distributions NOT including * * (possible) folding distributions * * Vpdpsf (i,k) = Vp (central potential) depth for the potential * * distribution * * Vpedsf (i,k) = Vp "edge" depth at the nucleus p/n radius * * Vprdsf (i,k) = Radius defining the potential "edge" * * Rcvpsf (i,k) = correlation r_0 for saturation correction for po- * * tential distributions * * Vpbpsf (i,k) = b (skin-like) parameter for the potential distri- * * bution * * Vrmssf (i,k) = r-rms () parameter for the potential distri- * * bution * * Vreqsf (i,k) = sharp radius parameter for the potential distri- * * bution * * Pcvpsf (k) = density exponent for saturation correction for po-* * tential distributions * * Vpexsf (k) = volume integral of the contact like (Pauli) part * * of the exchange interaction * * Akvppp = relative strength of pp(nn) interaction * * Akvppn = relative strength of pn(np) interaction * * Bkvppp = relative strength of pp(nn) exchange interaction * * Bkvppn = relative strength of pn(np) exchange interaction * * Vbarnn (k) = volume integral of the NN interaction (GeV x fm^3)* * Jnrhvp = present nucleon index (1=proton, 2=neutron) * * Jbouvp = present bound or unbound nucleon index (1=bound, * * 2=unbound) * * Jrhodn = present density index 1-6 * * Jrhofl = density index (1-6) to be (un)folded * * Lshmdu = flag for shell model like density calculations * * * *----------------------------------------------------------------------* * * R_0 for the equivalent sharp sphere radius for the proton density * (or charge) distribution PARAMETER ( R0PRCH = 1.128 D+00 ) * R_0 fundamental one according to Myers PARAMETER ( R0DRPM = 1.18 D+00 ) * Various parameters of the Droplet model * PARAMETER ( CKDRPM = 240. D+00 ) PARAMETER ( CJDRPM = 36.8 D+00 ) PARAMETER ( CLDRPM = 100. D+00 ) PARAMETER ( CMDRPM = 0. D+00 ) PARAMETER ( CQDRPM = 17. D+00 ) PARAMETER ( CA1DRP = 15.96 D+00 ) PARAMETER ( CA2DRP = 20.69 D+00 ) PARAMETER ( CA3DRP = 0. D+00 ) PARAMETER ( CC1DRP = THRTHR / FIVFIV * COUGFM / R0DRPM * GEVMEV ) PARAMETER ( CC2DRP = CC1DRP / 336.D+00 * ( ONEONE / CJDRPM & + 18.D+00 / CKDRPM ) ) PARAMETER ( HC3DRP = CC1DRP * HLFHLF * FIVFIV / R0DRPM / R0DRPM ) * ( 3 / 2 pi )^2/3 PARAMETER ( HC4DRP = 0.610887057710857 D+00 ) PARAMETER ( CC4DRP = FIVFIV / FOUFOU * HC4DRP * CC1DRP ) PARAMETER ( CC5DRP = CC1DRP * CC1DRP / CQDRPM / 64.D+00 ) PARAMETER ( R0SSHM = HLFHLF / PIPIPI / ERFA00 ) PARAMETER ( R0PSHM = ONETHI / PIPIPI / ERFA00 ) LOGICAL LSHMDU, LUNFRU COMMON / NCDNVP / RHCPSF (2,8), RHAPSF (2,8), RH0PSF (2,8), & RHCESF (2,8), RHRDSF (2,8), RHFRSF (2,8), & RRMSSF (2,8), RREQSF (2,8), VPCPSF (2,2), & VPAPSF (2,2), VP0PSF (2,2), VPDPSF (2,2), & VPEDSF (2,2), VPRDSF (2,2), VRMSSF (2,2), & VREQSF (2,2), VPBPSF (2,2), RCVPSF (2,2), & PCVPSF (2) , VBARNN (2) , VPEXSF (2), & SSHLNC (2) , PSHLNC (2) , AKVPPP, AKVPPN, & BKVPPP, BKVPPN, & JRHDNF (2) , JNRHVP, JBOUVP, JRHODN, JRHOFL, & LSHMDU, LUNFRU SAVE / NCDNVP /