From 912e9e1ec347872013f00c3bf5997190d3945124 Mon Sep 17 00:00:00 2001
From: Lorenzo Figini <lorenzo.figini@istp.cnr.it>
Date: Thu, 11 Apr 2013 13:43:23 +0000
Subject: [PATCH] corrected initial conditions in ic_gb

---
 src/gray.f | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/src/gray.f b/src/gray.f
index 16e3505..02db2b3 100644
--- a/src/gray.f
+++ b/src/gray.f
@@ -3893,8 +3893,11 @@ c
             du1tz=0.5d0*u*(dx0t**2*dwwxx+dy0t**2*dwwyy+
      .            2.0d0*dx0t*dy0t*dwwxy)/ddfu
 c
-          ppx=x0t*rcixx+y0t*rcixy
-          ppy=x0t*rcixy+y0t*rciyy
+          pppx=x0t*rcixx+y0t*rcixy
+          pppy=x0t*rcixy+y0t*rciyy
+          ppx=-pppx*gzt/(pppx*gxt+pppy*gyt)
+          ppy=-pppy*gzt/(pppx*gxt+pppy*gyt)
+c
           anzt=sqrt((1.0d0+gr2)/(1.0d0+ppx**2+ppy**2))
           anxt=ppx*anzt
           anyt=ppy*anzt