c..for silly compilers that don't have the circular functions in degrees.

      double precision function cosd(angle)
      implicit none
      save
      double precision angle
      double precision pi,a2rad,rad2a
      parameter        (pi    = 3.1415926535897932384d0, 
     1                  a2rad = pi/180.0d0,  
     2                  rad2a = 180.0d0/pi)
      cosd = cos(angle*a2rad)
      return
      end

      double precision function sind(angle)
      implicit none
      save
      double precision angle
      double precision pi,a2rad,rad2a
      parameter        (pi    = 3.1415926535897932384d0, 
     1                  a2rad = pi/180.0d0,  
     2                  rad2a = 180.0d0/pi)
      sind = sin(angle*a2rad)
      return
      end

      double precision function tand(angle)
      implicit none
      save
      double precision angle
      double precision pi,a2rad,rad2a
      parameter        (pi    = 3.1415926535897932384d0, 
     1                  a2rad = pi/180.0d0,  
     2                  rad2a = 180.0d0/pi)
      tand = tan(angle*a2rad)
      return
      end





      double precision function acosd(rads)
      implicit none
      save
      double precision rads
      double precision pi,a2rad,rad2a
      parameter        (pi    = 3.1415926535897932384d0, 
     1                  a2rad = pi/180.0d0,  
     2                  rad2a = 180.0d0/pi)
      acosd = acos(rads) * rad2a
      return
      end

      double precision function asind(rads)
      implicit none
      save
      double precision rads
      double precision pi,a2rad,rad2a
      parameter        (pi    = 3.1415926535897932384d0, 
     1                  a2rad = pi/180.0d0,  
     2                  rad2a = 180.0d0/pi)
      asind = asin(rads) * rad2a
      return
      end

      double precision function atan2d(denom,numer)
      implicit none
      save
      double precision denom,numer
      double precision pi,a2rad,rad2a
      parameter        (pi    = 3.1415926535897932384d0, 
     1                  a2rad = pi/180.0d0,  
     2                  rad2a = 180.0d0/pi)
      atan2d = atan2(denom,numer)*rad2a
      return
      end