Fixed bug.

Routine forcst was crashing if at least one of the observed variables,
with measurement error, was not in varlist.
time-shift
Stéphane Adjemian (Hermes) 2016-08-23 16:38:26 +02:00
parent 5df056b38f
commit 38e1ec6b1e
1 changed files with 4 additions and 2 deletions

View File

@ -86,8 +86,10 @@ end
if nargout==3
var_yf_ME=var_yf;
[loc_H,loc_varlist]=ismember(options_.varobs',options_.varlist);
loc_varlist(loc_varlist==0)=[];
var_yf_ME(:,loc_varlist)=var_yf(:,loc_varlist)+repmat(diag(M_.H(loc_H,loc_H))',horizon,1);
loc_varlist(loc_varlist==0)=[];
if ~isempty(loc_varlist)
var_yf_ME(:,loc_varlist)=var_yf(:,loc_varlist)+ repmat(diag(M_.H(loc_H,loc_H))',horizon,1);
end
int_width_ME = zeros(horizon,nvar);
end