subroutine dload(f,kstep,kinc,time,noel,npt,layer,kspt,coords,
     $     jltyp,sname)
C     
      include 'ABA_PARAM.INC'
      dimension time(2), coords(3)
      character*80 sname
      parameter (pi=3.14159265358979d0,two=2.0d0,zero=0.d0)
      parameter(onetmp=1.d0, abig=1.d36)
      parameter(asmall=onetmp/abig)
      parameter(pten=5.7297)
      angnp1=0.0
      dangBeg=7*pi/9
      dangEnd1=dangBeg+ pi/18
      dangEnd2=dangBeg+pi/6
      dangEnd3=dangBeg+pi/4.5
C     
      c=100.0d0
      ca=70.0d0
      f=0.0d0
C
      y = coords(1)
      x = coords(3)
      dydx = abs(y/x)
      angnp1 = atan(dydx)
      if (y .gt. zero) then
        if (x .lt. zero) then
          angnp1 = pi-angnp1
        end if
      else
        if (x .lt. zero) then
          angnp1 = pi+angnp1
        else
          angnp1 = two*pi-angnp1
        end if
      end if
      if(angnp1.ge.dangBeg.and.angnp1.le.dangEnd1)Then
         f=c-pten*(c-ca)*(dangEnd1-angnp1)
      else if (angnp1.ge.dangEnd1.and.angnp1.le.dangEnd2)Then
         f=c
      else if (angnp1.ge.dangEnd2.and.angnp1.le.dangEnd3)Then 
         f=ca+pten*(c-ca)*(dangEnd3-angnp1)
      else
         f=0.0
      endif

      rad=sqrt(coords(1)**2+coords(3)**2)
      if(f.gt.asmall)Then
         f=f+8.0*(rad-1.0)
      endif

      return
      end