diff --git a/src/gray.f b/src/gray.f
index 1257aaf..e4a2271 100644
--- a/src/gray.f
+++ b/src/gray.f
@@ -186,7 +186,7 @@ c
       dimension iiv(jmx,kmx),tauv(jmx,kmx,nmx),alphav(jmx,kmx,nmx)
       dimension iop(jmx,kmx),iow(jmx,kmx),ihcd(jmx,kmx),istore(jmx,kmx)
       dimension tau1v(jmx,kmx),yyrfl(jmx,kmx,6),y(6),dery(6)
-      dimension ext(jmx,kmx,0:3),eyt(jmx,kmx,0:3)
+      dimension ext(jmx,kmx,0:4),eyt(jmx,kmx,0:4)
       dimension xv(3),anv(3),xvrfl(3),anvrfl(3),xvvac(3),anw(3),anwcl(3)
       dimension xwcl(3),xvjk(3,jmx,kmx),anvjk(3,jmx,kmx)
       dimension ywrk(6,jmx,kmx),ypwrk(6,jmx,kmx)
@@ -3969,7 +3969,7 @@ c
       dimension dffiu(jmx),ddffiu(jmx)
       dimension grad2(jmx,kmx),dgrad2v(ndimm,jmx,kmx)
       dimension gri(3,jmx,kmx),ggri(3,3,jmx,kmx)
-      complex*16 ext(jmx,kmx,0:3),eyt(jmx,kmx,0:3)
+      complex*16 ext(jmx,kmx,0:4),eyt(jmx,kmx,0:4)
       complex*16 ui,sss,ddd,phic,qi1,qi2,tc,ts,qqxx,qqxy,qqyy
       complex*16 dqi1,dqi2,dqqxx,dqqyy,dqqxy
       complex*16 d2qi1,d2qi2,d2qqxx,d2qqyy,d2qqxy
@@ -4285,7 +4285,7 @@ c
       dimension dffiu(jmx),ddffiu(jmx)
       dimension grad2(jmx,kmx),dgrad2v(ndimm,jmx,kmx)
       dimension gri(3,jmx,kmx),ggri(3,3,jmx,kmx)
-      complex*16 ext(jmx,kmx,0:3),eyt(jmx,kmx,0:3)
+      complex*16 ext(jmx,kmx,0:4),eyt(jmx,kmx,0:4)
 c
       common/nray/nrayr,nrayth
       common/rwmax/rwmax
@@ -4481,7 +4481,7 @@ c
       dimension dffiu(jmx),ddffiu(jmx)
       dimension grad2(jmx,kmx),dgrad2v(ndimm,jmx,kmx)
       dimension gri(3,jmx,kmx),ggri(3,3,jmx,kmx)
-      complex*16 ext(jmx,kmx,0:3),eyt(jmx,kmx,0:3)
+      complex*16 ext(jmx,kmx,0:4),eyt(jmx,kmx,0:4)
 c
       common/nray/nrayr,nrayth
       common/wrk/ywrk0,ypwrk0