if(rad.ge.5.0001.and.rad.lt.600.)then
      if(vmaxt.lt.65)then
        if(shamp.lt.9.7)then
          a_1 = -17.5*exp(-(rad-100)**2/((rad+50)*0.5)**2)+7.5
          b_1 = 1.1*(rad-5)**0.6-exp((rad-125)/100)
          a_2 = 4-10*exp(-(rad-245)**2/250**2)
          b_2 = 1-10*exp(-(rad-335)**2/60**2)
        elseif(shamp.ge.9.7.and.shamp.lt.19.4)then       
          a_1 = -20*exp(-(rad-65)**2/((rad+50)*0.65)**2)+2.5
          b_1 = 1.1*(rad-5)**0.75-((rad-25)/47.5)**2
          a_2 = -2.5-12.5*exp(0-(rad-245)**2/100**2)
          b_2 = -1.5*(rad-5)**0.25
        else
          a_1 = -35*exp(-(rad-85)**2/((rad+150)*0.35)**2)-5
          b_1 = 1.1*(rad-5)**0.75-((rad-25)/50)**2+12.5
          a_2 = -1.5-20*exp(-(rad-205)**2/100**2)
          b_2 = -2.5*(rad-5)**0.32
        endif
      elseif(vmaxt.ge.65.and.vmaxt.lt.95)then
        if(shamp.lt.9.7)then
          a_1 = -15*exp(-(rad-100)**2/((rad+300)*0.25)**2)+(-8+rad/20)-
     1             exp((rad-175)/85)
          b_1 = 1.1*(rad-5)**0.35
          a_2 = 10-8*exp((rad)**2/600**2)
          b_2 = 6.5-8*exp(-(rad-245)**2/250**2)-exp((rad-250)/100) 
        elseif(shamp.ge.9.7.and.shamp.lt.19.4)then       
          a_1 = -15*exp(-(rad-100)**2/((rad+300)*0.2)**2)-
     1             9-exp((rad-250)/100)
          b_1 = 1.1*(rad-5)**0.75-2.5*(rad/50)**1.5+2.5
          a_2 = -2.5-10*exp(-(rad-225)**2/(100-(rad-100)/10)**2)
          b_2 = -1.25*(rad-5)**0.2-exp((rad-400)/100)
        else
          a_1 = -25*exp(-(rad-100)**2/((rad+100)*0.75)**2)-
     1             9+exp((rad-200)/75)
          b_1 = 1.1*(rad-5)**0.7-((rad-25)/75)**2+5
          a_2 = -1.5-16*exp(-(rad-375)**2/200**2)
          b_2 = -2.5*(rad-5)**0.23
        endif
      else
        if(shamp.lt.9.7)then
          a_1 = -15*exp(-(rad-100)**2/((rad+300)*0.5)**2)+(5+rad/20)-
     1             2*exp((rad-75)/150)
          b_1 = 1.1*(rad-5)**0.35
          a_2 = 1.5-11.5*exp(-(rad-375)**2/125**2)
          b_2 = 6.5-8*exp(-(rad-245)**2/250**2)+
     1            22.5*exp(-(rad-410)**2/50**2)
        elseif(shamp.ge.9.7.and.shamp.lt.19.4)then       
          a_1 = -15*exp(-(rad-150)**2/((rad+500)*0.15)**2)-11
          b_1 = 1.1*((rad-5)/10)**1.3-7*((rad-50)/100)**2
          a_2 = 1-10*exp(-(rad-350)**2/125**2)
          b_2 = 4.5-20*exp(-(rad-300)**2/100**2)
        else
          a_1 = -15*exp(-(rad-100)**2/((rad+300)*0.35)**2)-15
          b_1 = 1.1*(rad-5)**0.72-((rad-25)/50)**2+5
          a_2 = -1.5-20*exp(-(rad-245)**2/50**2)+
     1            15*exp(-(rad-340)**2/50**2)
          b_2 = -5-10*exp(-(rad-245)**2/100**2)
        endif
      endif
      else
          a_1=0.
          b_1=0.
          a_2=0.
          b_2=0.
      endif