phSolver/compressible/e3tau.f

2      &                     u1,     u2,     u3,
41      &            u2(npro),                  u3(npro),
69      &         + u3*u3*(gijd(:,1)*gijd(:,3)-gijd(:,2)*gijd(:,2))
71      &         + u1*u3*(gijd(:,2)*gijd(:,5)-gijd(:,4)*gijd(:,3))
72      &         + u2*u3*(gijd(:,4)*gijd(:,2)-gijd(:,1)*gijd(:,5)))
84      &     + u3*u3*gijd(:,6)
86      &     + u1*u3*gijd(:,4)
87      &     + u2*u3*gijd(:,5))*two  + 1.0e-15 !FIX FOR INVALID MESHES
163      &        + (u1*(u1*gijd(:,1) + two*(u2*gijd(:,2)+u3*gijd(:,4)))
164      &        + u2*(u2*gijd(:,3) + two*u3*gijd(:,5))
165      &        + u3*u3*gijd(:,6)))
285      &     u1,              u2,             u3,
326      &            u2(npro),                  u3(npro),
352       u = sqrt(u1**2+u2**2+u3**2)
353       eb1 = cp*T - pt5*(u1**2+u2**2+u3**2)
372      &               u2,              u3,              a1,
585         betaT = cp*T + pt5*(u1**2+u2**2+u3**2) !reuse betaT
588      &         u2*STau(:,4)+u3*STau(:,7)+betaT*STau(:,11))
590      &         u2*STau(:,5)+u3*STau(:,8)+betaT*STau(:,12))
592      &         u2*STau(:,6)+u3*STau(:,9)+betaT*STau(:,13))
594      &         u2*STau(:,9)+u3*STau(:,10)+betaT*STau(:,14))
596      &         u2*STau(:,13)+u3*STau(:,14)+betaT*STau(:,15))
611         Tau(:,4,1) = T(:)*(STau(:,7) + u3(:)*STau(:,11))
612         Tau(:,4,2) = T(:)*(STau(:,8) + u3(:)*STau(:,12))
613         Tau(:,4,3) = T(:)*(STau(:,9) + u3(:)*STau(:,13))
614         Tau(:,4,4) = T(:)*(STau(:,10) + u3(:)*STau(:,14))
615         Tau(:,4,5) = T(:)*(STau(:,14) + u3(:)*STau(:,15))
746      &     u1,              u2,             u3,
788      &            u2(npro),                  u3(npro),
814       u = sqrt(u1**2+u2**2+u3**2)
815       eb1 = cp*T - pt5*(u1**2+u2**2+u3**2)
834      &               u2,              u3,              a1,
1047          betaT = cp*T + pt5*(u1**2+u2**2+u3**2) !reuse betaT
1050      &        u2*STau(:,4,k)+u3*STau(:,7,k)+betaT*STau(:,11,k))
1052      &        u2*STau(:,5,k)+u3*STau(:,8,k)+betaT*STau(:,12,k))
1054      &        u2*STau(:,6,k)+u3*STau(:,9,k)+betaT*STau(:,13,k))
1056      &        u2*STau(:,9,k)+u3*STau(:,10,k)+betaT*STau(:,14,k))
1058      &        u2*STau(:,13,k)+u3*STau(:,14,k)+betaT*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))
1212      &                       u1,     u2,     u3,
1252      &            u2(npro),                  u3(npro),
1280      &     + u3*u3*(gijd(:,1)*gijd(:,3)-gijd(:,2)*gijd(:,2))
1282      &     + u1*u3*(gijd(:,2)*gijd(:,5)-gijd(:,4)*gijd(:,3))
1283      &     + u2*u3*(gijd(:,4)*gijd(:,2)-gijd(:,1)*gijd(:,5)))
1297      &     + u3*u3*gijd(:,6)
1299      &     + u1*u3*gijd(:,4)
1300      &     + u2*u3*gijd(:,5))*two  + 1.0e-15 !FIX FOR INVALID MESHES