Evolution equation for continued

! Evolve K3
*     Lie_beta(K^1_1)
      call diff(nx,ny,nz,h,k11,temp1,'x',ml,dl,nl)
      call diff(nx,ny,nz,h,k11,temp2,'y',ml,dl,nl)
      call diff(nx,ny,nz,h,k11,temp3,'z',ml,dl,nl)
      t11 = t11 + beta1 * temp1 + beta2*temp2+beta3*temp3

      call diff(nx,ny,nz,h,beta1,temp1,'x',ml,dl,nl)
      call diff(nx,ny,nz,h,beta1,temp2,'y',ml,dl,nl)
      call diff(nx,ny,nz,h,beta1,temp3,'z',ml,dl,nl)
      t11 = t11 - k11 * temp1 - k21 * temp2 - k31 * temp3

      call diff(nx,ny,nz,h,beta1,temp1,'y',ml,dl,nl)
      call diff(nx,ny,nz,h,beta2,temp2,'y',ml,dl,nl)
      call diff(nx,ny,nz,h,beta3,temp3,'y',ml,dl,nl)
      t11 = t11 + k11 * temp1 + k12 * temp2 + k13 * temp3
C Evolve K4
      trace_k(1) = k11 + k22 + k33
      FORALL (i=1:3,j=1:3)
         t(i,j) = t(i,j) + alpha(1) *(
     .              gup(i,1) * r(1,j) + gup(i,2) * r(2,j) +
     .              gup(i,3) * r(3,j) + trace_k(1) * k(i,j) )
      END FORALL
C Evolve K5
      FORALL (i=1:3,j=1:3)
         k(i,j,1) = 2.0d0*dt*t(i,j) + k(i,j,-1)
      END FORALL