1 subroutine getFld (T, cp, rmu, rlm, rlm2mu, 2 & con) 3c 4c---------------------------------------------------------------------- 5c 6c This routine calculates the fluid material properties. 7c 8c input: 9c T (npro) : temperature 10c cp (npro) : specific heat at constant pressure 11c 12c output: 13c rmu (npro) : Mu 14c rlm (npro) : Lambda 15c rlm2mu (npro) : Lambda + 2 Mu 16c con (npro) : Conductivity 17c 18c Note: material type flags 19c matflg(2): 20c eq. 0, constant viscosity 21c eq. 1, generalized Sutherland viscosity 22c matflg(3): 23c eq. 0, Stokes approximation 24c eq. 1, shear proportional bulk viscosity 25c 26c 27c Farzin Shakib, Winter 1987. 28c Zdenek Johan, Winter 1991. (Fortran 90) 29c---------------------------------------------------------------------- 30c 31 include "common.h" 32c 33 dimension T(npro), cp(npro), 34 & rmu(npro), rlm(npro), 35 & rlm2mu(npro), con(npro) 36c 37c 38c.... constant viscosity 39c 40 if (matflg(2,1) .eq. 0) then 41c 42 rmu = datmat(1,2,1) 43c 44 else 45c 46c.... generalized Sutherland viscosity 47c 48 rmu = datmat(1,2,1) * (T/datmat(2,2,1))*sqrt(T/datmat(2,2,1)) 49 & * ( datmat(2,2,1) + datmat(3,2,1) ) / (T + datmat(3,2,1)) 50c 51 endif 52c 53c.... calculate the second viscosity coefficient 54c 55 if (matflg(3,1) .eq. 0) then 56c 57 rlm = -pt66 * rmu 58c 59 else 60c 61 rlm = (datmat(1,3,1) - pt66) * rmu 62c 63 endif 64c 65c.... calculate the remaining quantities 66c 67 cp = datmat(1,3,1) 68 rlm2mu = rlm + two * rmu 69 con = datmat(1,4,1) 70c 71c.... return 72c 73 return 74 end 75