options nocenter ps = 80 ls = 76;
data thawall;
   infile 'thawall.dat' firstobs = 2;
   input ex p a t r ilc ilt chlc chlt ruc rut rchlc rchlt puc put;
   pfold = put / puc;
   pfold2 = sqrt(pfold + 0.5);
   ifold = ilt / ilc;
   ifold2 = sqrt(ifold + 0.5);
   rfold = rut / ruc;
   rfold2 = sqrt(rfold + 0.5);
   chlac = chlc * 1000 / 32.0762143;
   chlat = chlt * 1000 / 32.0762143;
   pchlc = 0.9936 * puc / chlc;
   pchlt = 0.9936 * put / chlt;
   pareac = puc * 0.9936 / 32.0762143;
   pareat = put * 0.9936 / 32.0762143;
   rareac = ruc * 0.9936 / 32.0762143;
   rareat = rut * 0.9936 / 32.0762143;
   parel = (pareac - pareat) / pareac;
   rarel = (rareac - rareat) / rareac;
   idiff = ilt - ilc;
   pdiff = pchlc - pchlt;
   rdiff = rchlc - rchlt;
   irel = idiff / ilc;
   prel = pdiff / pchlc;
   rrel = rdiff / rchlc;
   
proc sort; by a;
proc print;
   var a t idiff pdiff;