function [sig,beta]=molec(lambda,T,p);
% Calculate the extinction coefficien and 
% back scatter based on temperature and pressure


R=287;
g=9.81;
P=3/(8*pi);
sig=g./(R*T).*rayleigh(lambda,p);
beta=sig*P;