		load ./highfinput   
        X=highf;
		%
		% Options for regem
		%
		OPTIONS.regress = 'ttls';
		OPTIONS.regpar = 15
		OPTIONS.stagtol = 5e-2;
		OPTIONS.maxit = 40;
		OPTIONS.inflation = 1; % this works pretty well all around. One could do a short 
		%%%%%%%%%%%%%%%%%%%%%%%%%%calibration (1900-1980 in our case) with validation to adjust for a given 
		%%%%%%%%%%%%%%%%%%%%%%%%%application if necessary
		OPTIONS.disp = 1;
		%
		[X, M, C, Xerr] = regemhigh(X, OPTIONS);
         save highf.filled

		
