Evolution equation for

We break this equation up into 5 parts(5 modules)

C evolve_k1
      call diff(nx,ny,nz,h,alpha,temp1,'x',ml,dl,nl)
      call diff(nx,ny,nz,h,alpha,temp2,'y',ml,dl,nl)
      call diff(nx,ny,nz,h,alpha,temp3,'z',ml,dl,nl)
C evolve_k2
      call diff(nx,ny,nz,h,temp1,t11,'x',ml,dl,nl)
      call diff(nx,ny,nz,h,temp1,t12,'y',ml,dl,nl)
      call diff(nx,ny,nz,h,temp1,t13,'z',ml,dl,nl)
      call diff(nx,ny,nz,h,temp2,t22,'y',ml,dl,nl)
      call diff(nx,ny,nz,h,temp2,t23,'z',ml,dl,nl)
      call diff(nx,ny,nz,h,temp3,t33,'z',ml,dl,nl)
      t21 = t12
      t31 = t13
      t32 = t23
      FORALL (j=1:3,i=1:3,i<=j)
         t(i,j) = t(i,j) - c(1,i,j) *dalpha(1) -
     .              c(2,i,j)*dalpha(2) - c(3,i,j) * dalpha(3)
         t(j,i) = t(i,j)
      END FORALL
      call rais02i1() 
      FORALL (i=1:3,j=1:3)
         t(i,j) = -t(i,j)
      END FORALL