disp('Enter Freq [GHz], T [K], Pres [mb], h2omix [kg/kg]:')
prompt=input('','s');
dat=textscan(prompt,'%f','delimiter',',');

atm_ext=absorb(dat{1}(1),dat{1}(2),dat{1}(3),dat{1}(4),0);
disp(['Atm_ext [km^{-1}] = ',num2str(atm_ext,'%10.5e\n')])