source: LMDZ6/branches/Optimisation_LMDZ/libf/misc/cbrt.F @ 3703

Last change on this file since 3703 was 2291, checked in by Ehouarn Millour, 9 years ago

Added utility routines cbrt to compute cubic roots (used by some planetary physics packages) and getin_p (to read from def files in parallel).
EM
File size: 185 bytes
Line 
1      FUNCTION cbrt(x)
2      ! Return the cubic root for positive or negative x
3      IMPLICIT NONE
4
5      REAL x,cbrt
6
7      cbrt=sign(1.,x)*(abs(x)**(1./3.))
8
9      RETURN
10      END
11
Note: See TracBrowser for help on using the repository browser.