From 38e1ec6b1e9295a328f5925aa034f82689f56f4d Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?St=C3=A9phane=20Adjemian=20=28Hermes=29?= Date: Tue, 23 Aug 2016 16:38:26 +0200 Subject: [PATCH] Fixed bug. Routine forcst was crashing if at least one of the observed variables, with measurement error, was not in varlist. --- matlab/forcst.m | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/matlab/forcst.m b/matlab/forcst.m index c2e46ddf3..9c7654ae6 100644 --- a/matlab/forcst.m +++ b/matlab/forcst.m @@ -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