1 MODULE visc_g 2 3 4 Use param 5 Use param1 6 7 8 ! 9 ! trace of D_g at i, j, k 10 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: trD_g 11 ! 12 ! Turbulent viscosity of fluid phase: sum of molecular 13 ! and eddy viscosities 14 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: MU_gt 15 ! 16 ! Granular second coefficient of viscosity 17 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: LAMBDA_gt 18 ! 19 ! Characteristic length for turbulence 20 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: L_scale 21 22 ! 23 ! gas viscosity needs to be updated during iterations? 24 LOGICAL Recalc_visc_g 25 26 27 28 !!!HPF$ align trD_g(:) with TT(:) 29 !!!HPF$ align MU_gt(:) with TT(:) 30 !!!HPF$ align LAMBDA_gt(:) with TT(:) 31 !!!HPF$ align L_scale(:) with TT(:) 32 33 END MODULE visc_g 34