subroutine prepot c subroutine inifit c c..... this subroutine simply reads the fit parameters and should be c..... referenced before the potential is called fot the first time c..... the fitting parameters are stored in the common blocks common/cparm/ c......and common/cint/ c implicit real*8(a-h,o-z) parameter(nmx=198,npmx=12) common/cparm/ a(nmx),p(npmx) common/cint/ ix(nmx),iy(nmx),iz(nmx),mmax data ix / * 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, * 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, * 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, * 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, * 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, * 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, * 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, * 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, * 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, * 7, 8, 8, 0, 0, 0, 0, 0, 0/ data iy / * 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, * 4, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 7, 7, 8, 0, 0, 0, 0, 0, 0, * 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, * 3, 3, 4, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 7, 7, 8, 0, 0, 0, 0, * 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, * 4, 4, 4, 4, 5, 5, 5, 6, 6, 7, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, * 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 5, 5, 6, 0, 0, 0, 0, 0, * 1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 4, 4, 5, 0, 0, 0, 0, 1, 1, * 1, 1, 2, 2, 2, 3, 3, 4, 0, 0, 0, 1, 1, 1, 2, 2, 3, 0, 0, 1, 1, * 2, 0, 1, 0, 0, 0, 0, 0, 0/ data iz / * 1, 2, 3, 4, 5, 6, 7, 8, 1, 2, 3, 4, 5, 6, 7, 1, 2, 3, 4, 5, 6, * 1, 2, 3, 4, 5, 1, 2, 3, 4, 1, 2, 3, 1, 2, 1, 1, 2, 3, 4, 5, 6, * 7, 8, 0, 1, 2, 3, 4, 5, 6, 7, 0, 1, 2, 3, 4, 5, 6, 0, 1, 2, 3, * 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 3, 0, 1, 2, 0, 1, 0, 1, 2, 3, 4, * 5, 6, 7, 0, 1, 2, 3, 4, 5, 6, 0, 1, 2, 3, 4, 5, 0, 1, 2, 3, 4, * 0, 1, 2, 3, 0, 1, 2, 0, 1, 0, 1, 2, 3, 4, 5, 6, 0, 1, 2, 3, 4, * 5, 0, 1, 2, 3, 4, 0, 1, 2, 3, 0, 1, 2, 0, 1, 0, 1, 2, 3, 4, 5, * 0, 1, 2, 3, 4, 0, 1, 2, 3, 0, 1, 2, 0, 1, 0, 1, 2, 3, 4, 0, 1, * 2, 3, 0, 1, 2, 0, 1, 0, 1, 2, 3, 0, 1, 2, 0, 1, 0, 1, 2, 0, 1, * 0, 1, 0, 0, 0, 0, 0, 0, 0/ c The data has been splitted since the old sun compiler does not c accept more than 19 continuation card. data (a(i),i=1,72) / * 3.446353e-01, -1.030216e+00, 8.746902e-01, -1.756380e-01, * -1.283421e-01, 4.046147e-01, -3.708786e-01, 9.815560e-02, * -1.042546e+00, 2.706368e+00, -7.156916e-01, 1.898295e-01, * -7.962587e-01, 3.527891e-01, 2.730624e-03, 1.174452e+00, * -1.918118e+00, -1.360439e+00, 6.168840e-01, 8.927372e-01, * -3.493863e-01, -4.248271e-01, 4.809767e-01, 2.485677e+00, * -1.145523e+00, 3.007968e-02, -6.998328e-01, -3.865584e-01, * -1.587307e+00, 1.513058e-01, 1.218593e+00, 4.061863e-01, * 6.164074e-01, -7.124100e-01, -1.912015e-01, 1.394301e-01, * -1.786764e-01, 2.603647e-01, -2.889760e-02, 1.128801e+00, * -2.424831e+00, 1.562723e+00, -2.747457e-01, -2.459826e-02, * 3.446353e-01, 1.423155e+00, -4.763302e-01, -2.017482e+00, * 2.356091e+00, -2.048582e+00, 1.250254e+00, -2.760210e-01, * -1.042546e+00, -2.125961e+00, 3.722102e-01, -2.426763e+00, * 3.657578e+00, -1.606100e+00, 2.173684e-02, 1.174452e+00, * 2.087784e-01, 6.280354e-01, 1.981022e+00, -3.933349e+00, * 1.073628e+00, -4.248271e-01, 6.207242e-01, -1.662749e+00, * 2.109088e+00, 1.107771e+00, -6.998328e-01, 7.194553e-01/ data (a(i),i=73,144) / * 1.685128e+00, -1.472575e+00, 1.218593e+00, -1.434811e+00, * -7.761216e-01, -7.124100e-01, 7.146970e-01, 1.394301e-01, * 2.603647e-01, -1.120778e+00, 2.014198e-01, 1.639369e+00, * -2.333689e-01, -1.000990e+00, 3.675041e-01, -1.030216e+00, * -4.763302e-01, 2.652506e+00, -1.811591e+00, 9.287747e-01, * -4.655039e-01, 1.799556e-02, 2.706368e+00, 3.722102e-01, * 5.783930e+00, -4.246126e+00, 1.993008e+00, 4.117272e-01, * -1.918118e+00, 6.280354e-01, -2.702635e+00, 3.510761e+00, * -2.223585e+00, 4.809767e-01, -1.662749e+00, -4.833982e+00, * -8.589023e-01, -3.865584e-01, 1.685128e+00, 3.171338e+00, * 4.061863e-01, -7.761216e-01, -1.912015e-01, -2.889760e-02, * 2.014198e-01, -1.384754e+00, -5.094533e-01, 1.715063e+00, * -4.619253e-01, 8.746902e-01, -2.017482e+00, -1.811591e+00, * 2.138828e+00, -6.286188e-01, 1.121651e-01, -7.156916e-01, * -2.426763e+00, -4.246126e+00, -7.668529e-01, -3.120883e-01, * -1.360439e+00, 1.981022e+00, 3.510761e+00, 2.237935e+00, * 2.485677e+00, 2.109088e+00, -8.589023e-01, -1.587307e+00, * -1.472575e+00, 6.164074e-01, 1.128801e+00, 1.639369e+00/ data (a(i),i=145,198) / * -5.094533e-01, -1.355390e+00, 1.525463e-01, -1.756380e-01, * 2.356091e+00, 9.287747e-01, -6.286188e-01, 7.764178e-02, * 1.898295e-01, 3.657578e+00, 1.993008e+00, -3.120883e-01, * 6.168840e-01, -3.933349e+00, -2.223585e+00, -1.145523e+00, * 1.107771e+00, 1.513058e-01, -2.424831e+00, -2.333689e-01, * 1.715063e+00, 1.525463e-01, -1.283421e-01, -2.048582e+00, * -4.655039e-01, 1.121651e-01, -7.962587e-01, -1.606100e+00, * 4.117272e-01, 8.927372e-01, 1.073628e+00, 3.007968e-02, * 1.562723e+00, -1.000990e+00, -4.619253e-01, 4.046147e-01, * 1.250254e+00, 1.799556e-02, 3.527891e-01, 2.173684e-02, * -3.493863e-01, -2.747457e-01, 3.675041e-01, -3.708786e-01, * -2.760210e-01, 2.730624e-03, -2.459826e-02, 9.815560e-02, * 8.970869e-01, 1.333088e+00, 8.970869e-01, 1.628670e+00, * 1.171451e+00, 1.628670e+00/ data p /.2236136690, 2.3412233303, 1.7347116549, 1.3636202850, * 0.7400188359, 0.0, .1736921184, 2.19220027263, 1.4001303861, * 1.3601150624, .6960248781, .1736921184/ m=ix(nmx-6) mmax=m+1 return end