Moments_varendo estimation option

Hi,
I am in the process of getting familiar with Dynare’s tools. I have downloaded the lastest Dynare version: 4.1.0 and I have installed Matlab 7.0 in my computer. I have estimated several versions of the New Keynesian model and parameter estimation procedures always run perfectly. However, I have problems when I introduce the option “moments_varendo” in the estimation procedure.

The estimation command used is standard and it reads as follows:

estimation(datafile=dci2,nobs=98,mh_replic=10000,mh_jscale=0.49,moments_varendo);

The error messages I am getting read as follows:

Warning: Divide by zero.

In dsge_simulated_theoretical_conditional_variance_decomposition at 71
In posterior_analysis>job at 75
In posterior_analysis at 34
In compute_moments_varendo at 110
In dynare_estimation_1 at 1070
In dynare_estimation at 62
In icasuno at 195
In dynare at 132
??? Reference to non-existent field ‘sigma_e_is_diagonal’.

Error in ==> dsge_simulated_theoretical_conditional_variance_decomposition at 120
StateSpaceModel.sigma_e_is_diagonal = M_.sigma_e_is_diagonal;

Error in ==> posterior_analysis>job at 75
[nvar,vartan,NumberOfFiles] = …

Error in ==> posterior_analysis at 34
oo_ = job(type,SampleSize,arg1,arg2,arg3,options_,M_,oo_);

Error in ==> compute_moments_varendo at 110
oo_ = posterior_analysis(‘conditional decomposition’,var_list_(i,:),M_.exo_names(j,:),Steps,options_,M_,oo_);

Error in ==> dynare_estimation_1 at 1070
oo_ = compute_moments_varendo(‘posterior’,options_,M_,oo_,var_list_);

Error in ==> dynare_estimation at 62
dynare_estimation_1(var_list,varargin{:});

Error in ==> icasuno at 195
dynare_estimation(var_list_);

Error in ==> dynare at 132
evalin(‘base’,fname) ;

Anyone knows what is going on? Thanks for your help.