diff --git a/src/gray.f b/src/gray.f
index 47dd34f..9883f15 100644
--- a/src/gray.f
+++ b/src/gray.f
@@ -376,7 +376,7 @@ c
       common/parwv/ak0,akinv,fhz
       common/dxgyg/derxg,deryg
       common/toxmxx/toxxmx,xgmax,ygxmx,anplxmx,anprxmx,thkbxmx,ixmx
-      common/toxmxt/toxmax,xgtmx,ygtmx,anpltmx,anprtmx,thkbtmx,itmx
+!      common/toxmxt/toxmax,xgtmx,ygtmx,anpltmx,anprtmx,thkbtmx,itmx
 
       anplopt=sqrt(yg/(yg+1.0d0))
       thkb=1.8d2/pi*acos(anpl/sqrt(an2))
@@ -398,16 +398,16 @@ c
           anprxmx=anpr
           thkbxmx=thkb
         end if
-        if (tox.gt.toxmax) then
-          itmx=i
-          toxmax=tox
-          xgtmx=xg
-          ygtmx=yg
-          anpltmx=anpl
-          anprtmx=anpr
-          thkbtmx=thkb
-        end if
-        write(46,98) i,tox,xg,yg,anpl,anpr,thkb,xv/1.d2,anv,derxg,deryg
+!        if (tox.gt.toxmax) then
+!          itmx=i
+!          toxmax=tox
+!          xgtmx=xg
+!          ygtmx=yg
+!          anpltmx=anpl
+!          anprtmx=anpr
+!          thkbtmx=thkb
+!        end if
+!        write(46,98) i,tox,xg,yg,anpl,anpr,thkb,xv/1.d2,anv,derxg,deryg
       end if
 c
       return
@@ -419,9 +419,9 @@ c
       subroutine print_tox
       implicit real*8 (a-h,o-z)
       common/toxmxx/toxxmx,xgmax,ygxmx,anplxmx,anprxmx,thkbxmx,ixmx
-      common/toxmxt/toxmax,xgtmx,ygtmx,anpltmx,anprtmx,thkbtmx,itmx
+!      common/toxmxt/toxmax,xgtmx,ygtmx,anpltmx,anprtmx,thkbtmx,itmx
       write(44,98) ixmx,toxxmx,xgmax,ygxmx,anplxmx,anprxmx,thkbxmx
-      write(45,98) itmx,toxmax,xgtmx,ygtmx,anpltmx,anprtmx,thkbtmx
+!      write(45,98) itmx,toxmax,xgtmx,ygtmx,anpltmx,anprtmx,thkbtmx
 c
       return
  98   format(1x,i8,30(1x,e16.8e3))
@@ -586,6 +586,7 @@ c
       write(7,*)'#Icd Pa Jphimx dPdVmx '//
      .'rhotj rhotjava rhotp rhotpav '//
      .'drhotjava drhotpav ratjbmx stmx psipol chipol index_rt'
+      write(44,*) '#i Tox X Y Npl Nperp thetakB'
       write(48,*) '#psi rhot Jphi Jcdb dPdV Icdins Pins P% index_rt'
 
       else
@@ -2698,7 +2699,7 @@ c
       common/ipol/ipolc
 c used only for computation of O-X transmission function
       common/toxmxx/toxxmx,xgmax,ygxmx,anplxmx,anprxmx,thkbxmx,ixmx
-      common/toxmxt/toxmax,xgtmx,ygtmx,anpltmx,anprtmx,thkbtmx,itmx
+!      common/toxmxt/toxmax,xgtmx,ygtmx,anpltmx,anprtmx,thkbtmx,itmx
 c
       istpr=0
       istpl=1
@@ -2708,9 +2709,9 @@ c
       ipolc=0
 c
       ixmx=0
-      itmx=0
+!      itmx=0
       xgmax=0.0d0
-      toxmax=0.0d0
+!      toxmax=0.0d0
 c
       return
       end
@@ -3555,7 +3556,7 @@ c
       parameter(lw20=nnw*2+nnh*4+nnw*nnh,lw02=nnw*4+nnh*2+nnw*nnh)
       parameter(lw11=nnw*3+nnh*3+nnw*nnh)
       parameter(nrs=1,nzs=1)
-      dimension cc01(lw10),cc10(lw01),cc02(lw02),cc20(lw20),cc11(lw11)
+      dimension cc01(lw01),cc10(lw10),cc02(lw02),cc20(lw20),cc11(lw11)
       dimension tfp(nrest),cfp(nrest),wrkfd(nrest)
 c
       common/eqnn/nr,nz,npp,nintp