disp('Which parameters?');
disp('  1 = Original ');
disp('  2 = theta=.216, corresponding shocks');
icase = input('? ');

if icase ==1;
  cd origparam
else
  cd altparam
end;   
lwedge
cd ..
