corrected initial conditions in ic_gb

2013-04-11 13:43:23 +00:00 · 2013-04-11 13:43:23 +00:00 · 912e9e1ec3
commit 912e9e1ec3
parent b83b6f0f5b
1 changed files with 5 additions and 2 deletions
--- 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