diff --git a/src/gray.f b/src/gray.f index 20aa448..46505e2 100644 --- a/src/gray.f +++ b/src/gray.f @@ -3895,8 +3895,14 @@ c c pppx=x0t*rcixx+y0t*rcixy pppy=x0t*rcixy+y0t*rciyy - ppx=-pppx*gzt/(pppx*gxt+pppy*gyt) - ppy=-pppy*gzt/(pppx*gxt+pppy*gyt) + denpp=pppx*gxt+pppy*gyt + if (denpp.ne.0.0d0) then + ppx=-pppx*gzt/denpp + ppy=-pppy*gzt/denpp + else + ppx=0.0d0 + ppy=0.0d0 + end if c anzt=sqrt((1.0d0+gr2)/(1.0d0+ppx**2+ppy**2)) anxt=ppx*anzt