anchor.diff:       &         luminosity(numberOfMassesWithColors,2),c_2)
             call extrap1(mass,c_1,c_2,mtrk(numberOfMassesWithColors-1),
        &         mtrk(numberOfMassesWithColors),c3)
  -          if(c3.lt.0.0) c3=0.0
  +          if (c3 < 0.0) then 
  +            c3 = 0.0
  +          end if
             call extrap1(lumi,color_R(numberOfMassesWithColors-1,1),
        &         color_R(numberOfMassesWithColors-1,2),
        &         luminosity(numberOfMassesWithColors-1,1),
  
positive.code:      &         luminosity(numberOfMassesWithColors,2),c_2)
            call extrap1(mass,c_1,c_2,mtrk(numberOfMassesWithColors-1),
       &         mtrk(numberOfMassesWithColors),c4)
            if(c4.lt.0.0) c4=0.0
            call extrap1(lumi,color_I(numberOfMassesWithColors-1,1),
       &         color_I(numberOfMassesWithColors-1,2),
       &         luminosity(numberOfMassesWithColors-1,1),
  
negative.code:      &         luminosity(numberOfMassesWithColors,ns2-1),
       &         luminosity(numberOfMassesWithColors,ns2),c_2)
            call extrap1(mass,c_1,c_2,mtrk(1),mtrk(2),c5)      
            check3 = 1
            GOTO 45
  C       Luminosity between known values
          else
            do i = 1, ns1 - 1
              if (lumi >= luminosity(numberOfMassesWithColors - 1, i + 1)
       &          .and. lumi 
       &                <= luminosity(numberOfMassesWithColors - 1, i)) 
       &      then
                i1 = i
                check1 = 1
                GOTO 15
              end if
            end do
  15        continue
            do i = 1, ns2 - 1
              if (lumi >= luminosity(numberOfMassesWithColors, i + 1) 
       &          .and. lumi <= luminosity(numberOfMassesWithColors, i)) 
       &      then
                i2 = i
                check2 = 1
                GOTO 20
              end if
            end do
  20        continue
            if (check1 == 1 .and. check2 == 1) then
              check3 = 1
              a1 = lumi - luminosity(numberOfMassesWithColors - 1, i1)
              a2 = lumi - luminosity(numberOfMassesWithColors, i2)
              b1 = luminosity(numberOfMassesWithColors - 1, i1 + 1) 
       &           - luminosity(numberOfMassesWithColors - 1, i1)
              b2 = luminosity(numberOfMassesWithColors, i2 + 1)
       &           - luminosity(numberOfMassesWithColors, i2)
              c_1 = color_U(numberOfMassesWithColors - 1, i1)
       &            + (color_U(numberOfMassesWithColors - 1, i1 + 1)
       &               - color_U(numberOfMassesWithColors - 1, i1)) * a1 
       &              / b1
              c_2 = color_U(numberOfMassesWithColors, i2)
       &            + (color_U(numberOfMassesWithColors, i2 + 1)
       &               - color_U(numberOfMassesWithColors, i2)) * a2 / b2
              call extrap1(mass,c_1,c_2,mtrk(1),
       &           mtrk(numberOfMassesWithColors),c1)
              c_1 = color_B(numberOfMassesWithColors - 1, i1)
       &            + (color_B(numberOfMassesWithColors - 1, i1 + 1)
       &               - color_B(numberOfMassesWithColors - 1, i1)) * a1 
       &              / b1
              c_2 = color_B(numberOfMassesWithColors, i2)
       &            + (color_B(numberOfMassesWithColors, i2 + 1)
       &               - color_B(numberOfMassesWithColors, i2)) * a2 / b2
              call extrap1(mass,c_1,c_2,mtrk(numberOfMassesWithColors-1),
       &           mtrk(numberOfMassesWithColors),c2)
              c_1 = color_V(numberOfMassesWithColors - 1, i1)
       &            + (color_V(numberOfMassesWithColors - 1, i1 + 1)
       &               - color_V(numberOfMassesWithColors - 1, i1)) * a1 
       &              / b1
              c_2 = color_V(numberOfMassesWithColors, i2) 
       &            + (color_V(numberOfMassesWithColors, i2 + 1)
       &               - color_V(numberOfMassesWithColors, i2)) * a2 / b2
              call extrap1(mass,c_1,c_2,mtrk(1),
       &           mtrk(numberOfMassesWithColors),c3)
              c_1 = color_R(numberOfMassesWithColors - 1, i1)
       &            + (color_R(numberOfMassesWithColors - 1, i1 + 1)
       &               - color_R(numberOfMassesWithColors - 1, i1)) * a1 
       &              / b1
              c_2 = color_R(numberOfMassesWithColors, i2)
       &            + (color_R(numberOfMassesWithColors, i2 + 1)
       &               - color_R(numberOfMassesWithColors, i2)) * a2 / b2
              call extrap1(mass,c_1,c_2,mtrk(numberOfMassesWithColors-1),
       &           mtrk(numberOfMassesWithColors),c4)
              c_1 = color_I(numberOfMassesWithColors - 1, i1)
       &            + (color_I(numberOfMassesWithColors - 1, i1 + 1)
       &               - color_I(numberOfMassesWithColors - 1, i1)) * a1 
       &              / b1
              c_2 = color_I(numberOfMassesWithColors, i2)
       &            + (color_I(numberOfMassesWithColors, i2 + 1)
       &               - color_I(numberOfMassesWithColors, i2)) * a2 / b2
              call extrap1(mass,c_1,c_2,mtrk(numberOfMassesWithColors-1),
       &           mtrk(numberOfMassesWithColors),c5)
              GOTO 45
  