fixed crash in RtoEuler for PHI=180°

2014-01-15 08:32:43 +00:00 · 2014-01-15 08:32:43 +00:00 · 98f677916c
parent cf318f4392
commit 98f677916c
1 changed files with 1 additions and 1 deletions
--- a/code/math.f90
+++ b/code/math.f90
@ -1418,7 +1418,7 @@ pure function math_RtoEuler(R)

 math_RtoEuler(2) = acos(myVal)

- if(math_RtoEuler(2) < 1.0e-8_pReal) then
+ if((math_RtoEuler(2) < 1.0e-8_pReal) .or. (pi-math_RtoEuler(2) < 1.0e-8_pReal)) then
 ! calculate phi2
     math_RtoEuler(3) = 0.0_pReal
 ! calculate phi1