Lines Matching refs:k

1045       do k = 1, ipord
1049          Tau(:,1,1,k) = rho*T*(STau(:,1,k)+u1*STau(:,2,k)+
1050      &        u2*STau(:,4,k)+u3*STau(:,7,k)+betaT*STau(:,11,k))
1051          Tau(:,1,2,k) = rho*T*(STau(:,2,k)+u1*STau(:,3,k)+
1052      &        u2*STau(:,5,k)+u3*STau(:,8,k)+betaT*STau(:,12,k))
1053          Tau(:,1,3,k) = rho*T*(STau(:,4,k)+u1*STau(:,5,k)+
1054      &        u2*STau(:,6,k)+u3*STau(:,9,k)+betaT*STau(:,13,k))
1055          Tau(:,1,4,k) = rho*T*(STau(:,7,k)+u1*STau(:,8,k)+
1056      &        u2*STau(:,9,k)+u3*STau(:,10,k)+betaT*STau(:,14,k))
1057          Tau(:,1,5,k) = rho*T*(STau(:,11,k)+u1*STau(:,12,k)+
1058      &        u2*STau(:,13,k)+u3*STau(:,14,k)+betaT*STau(:,15,k))
1061          Tau(:,2,1,k) = T(:)*(STau(:,2,k) + u1(:)*STau(:,11,k))
1062          Tau(:,2,2,k) = T(:)*(STau(:,3,k) + u1(:)*STau(:,12,k))
1063          Tau(:,2,3,k) = T(:)*(STau(:,5,k) + u1(:)*STau(:,13,k))
1064          Tau(:,2,4,k) = T(:)*(STau(:,8,k) + u1(:)*STau(:,14,k))
1065          Tau(:,2,5,k) = T(:)*(STau(:,12,k) + u1(:)*STau(:,15,k))
1067          Tau(:,3,1,k) = T(:)*(STau(:,4,k) + u2(:)*STau(:,11,k))
1068          Tau(:,3,2,k) = T(:)*(STau(:,5,k) + u2(:)*STau(:,12,k))
1069          Tau(:,3,3,k) = T(:)*(STau(:,6,k) + u2(:)*STau(:,13,k))
1070          Tau(:,3,4,k) = T(:)*(STau(:,9,k) + u2(:)*STau(:,14,k))
1071          Tau(:,3,5,k) = T(:)*(STau(:,13,k) + u2(:)*STau(:,15,k))
1073          Tau(:,4,1,k) = T(:)*(STau(:,7,k) + u3(:)*STau(:,11,k))
1074          Tau(:,4,2,k) = T(:)*(STau(:,8,k) + u3(:)*STau(:,12,k))
1075          Tau(:,4,3,k) = T(:)*(STau(:,9,k) + u3(:)*STau(:,13,k))
1076          Tau(:,4,4,k) = T(:)*(STau(:,10,k) + u3(:)*STau(:,14,k))
1077          Tau(:,4,5,k) = T(:)*(STau(:,14,k) + u3(:)*STau(:,15,k))
1081          Tau(:,5,1,k) = betaT(:)*STau(:,11,k)
1082          Tau(:,5,2,k) = betaT(:)*STau(:,12,k)
1083          Tau(:,5,3,k) = betaT(:)*STau(:,13,k)
1084          Tau(:,5,4,k) = betaT(:)*STau(:,14,k)
1085          Tau(:,5,5,k) = betaT(:)*STau(:,15,k)
1106            do k = 1, ipord
1109                     rLyi(:,i,k) = rLyi(:,i,k)+Tau(:,i,j,k)*eigmax(:,j) 
1191            do k = 1, ipord
1194          rLymi(:,i,k) = rLymi(:,i,k)+Tau(:,i,j,k)*eigmax(:,j)