Rgas=8.314d7
mu=.62
unit_rho=1d-24
unit_pp=1d-14
unit_TT=29.98
rhoT=ppm*unit_pp*mu/(Rgas*unit_rho*unit_TT)
mu=rhoT/(rhom*TTm)-1.
;plot,tt,mu,yr=[-1,0.]
;
plot,tt,rhom*TTm*unit_TT,yr=[0,500]
oplot,tt,rhoT*unit_TT
;
rhoC=14.3
rhoW=0.19
f=.1
;
rhom0=f*rhoC+(1.-f)*rhoW
Rgas=8.314e7
mu=0.62
P=24.2d-14
rhoT0=P*mu/(Rgas*unit_rho)
;
TC=P*mu/(Rgas*rhoC*unit_rho)
TW=P*mu/(Rgas*rhoW*unit_rho)
TTm0=f*TC+(1.-f)*TW
;
print,rhoT0/(TTm0*rhom0)
END