Generalize Viscosity used to calculatate Prandtl Number in Thermohaline

Lines 116 to 119 in 1baa434

    
           nu_mol = 0.406d0*sqrt(amu)*pow(boltzm*T,2.5d0)/(qe4*loglambdah*rho)  
        
           ! From Spitzer "Physics of Fully Ionized Gases equation 5-54 
        
           ! Assumes pure H. Still trying to work out what it would be for a mixture.  
        
           vis = nu_mol + nu_rad   ! total viscosity

	nu_mol = 0.406d0sqrt(amu)pow(boltzmT,2.5d0)/(qe4loglambdah*rho)
	! From Spitzer "Physics of Fully Ionized Gases equation 5-54
	! Assumes pure H. Still trying to work out what it would be for a mixture.
	vis = nu_mol + nu_rad ! total viscosity