WebSVN – lagranto.ecmwf – Diff – /trunk/density/density.f

           klo=k
         endif
       goto 1
       endif
       h=xa(khi)-xa(klo)
+      if (h.eq.0.) then
-      if (h.eq.0.) pause 'bad xa input in splint'
+         print*,'bad xa input in splint'
+         stop
+      endif
       a=(xa(khi)-x)/h
       b=(x-xa(klo))/h
       y=a*ya(klo)+b*ya(khi)+((a**3-a)*y2a(klo)+(b**3-b)*y2a(khi))*(h**
      *2)/6.
       return
       PHTOPHS = ZRPI18*ASIN(ZARG)
       RETURN
       END
+c     ------------------------------------------------------------------
+c     Compute Cos/Sin of an argument in Degree instead of Radian
+c     ------------------------------------------------------------------
+      real function cosd(arg)
+      real,intent(IN) :: arg
+      real,parameter :: grad2rad=3.1415926/360.
+      cosd=cos(arg*grad2rad)
+      return
+      end
+      real function sind(arg)
+      real,intent(IN) :: arg
+      real,parameter :: grad2rad=3.1415926/360.
+      sind=sin(arg*grad2rad)
+      return
+      end

Subversion Repositories lagranto.ecmwf