FUNCTION cbrt(x)
  ! ! Return the cubic root for positive or negative x
  IMPLICIT NONE

  REAL :: x,cbrt

  cbrt=sign(1.,x)*(abs(x)**(1./3.))

  RETURN
END FUNCTION cbrt