--- imach/src/imach.c	2003/06/27 13:00:02	1.94
+++ imach/src/imach.c	2004/02/20 13:25:42	1.97
@@ -1,6 +1,19 @@
-/* $Id: imach.c,v 1.94 2003/06/27 13:00:02 brouard Exp $
+/* $Id: imach.c,v 1.97 2004/02/20 13:25:42 lievre Exp $
   $State: Exp $
   $Log: imach.c,v $
+  Revision 1.97  2004/02/20 13:25:42  lievre
+  Version 0.96d. Population forecasting command line is (temporarily)
+  suppressed.
+
+  Revision 1.96  2003/07/15 15:38:55  brouard
+  * imach.c (Repository): Errors in subdirf, 2, 3 while printing tmpout is
+  rewritten within the same printf. Workaround: many printfs.
+
+  Revision 1.95  2003/07/08 07:54:34  brouard
+  * imach.c (Repository):
+  (Repository): Using imachwizard code to output a more meaningful covariance
+  matrix (cov(a12,c31) instead of numbers.
+
   Revision 1.94  2003/06/27 13:00:02  brouard
   Just cleaning
 
@@ -165,6 +178,9 @@
 #include <time.h>
 #include "timeval.h"
 
+/* #include <libintl.h> */
+/* #define _(String) gettext (String) */
+
 #define MAXLINE 256
 #define GNUPLOTPROGRAM "gnuplot"
 /*#define GNUPLOTPROGRAM "..\\gp37mgw\\wgnuplot"*/
@@ -193,11 +209,11 @@
 #define ODIRSEPARATOR '/'
 #endif
 
-/* $Id: imach.c,v 1.94 2003/06/27 13:00:02 brouard Exp $ */
+/* $Id: imach.c,v 1.97 2004/02/20 13:25:42 lievre Exp $ */
 /* $State: Exp $ */
 
-char version[]="Imach version 0.96b, June 2003, INED-EUROREVES ";
-char fullversion[]="$Revision: 1.94 $ $Date: 2003/06/27 13:00:02 $"; 
+char version[]="Imach version 0.96d, February 2004, INED-EUROREVES ";
+char fullversion[]="$Revision: 1.97 $ $Date: 2004/02/20 13:25:42 $"; 
 int erreur, nberr=0, nbwarn=0; /* Error number, number of errors number of warnings  */
 int nvar;
 int cptcovn=0, cptcovage=0, cptcoveff=0,cptcov;
@@ -240,7 +256,7 @@ char fileresvpl[FILENAMELENGTH];
 char title[MAXLINE];
 char optionfile[FILENAMELENGTH], datafile[FILENAMELENGTH],  filerespl[FILENAMELENGTH];
 char optionfilext[10], optionfilefiname[FILENAMELENGTH], plotcmd[FILENAMELENGTH];
-char tmpout[FILENAMELENGTH]; 
+char tmpout[FILENAMELENGTH],  tmpout2[FILENAMELENGTH]; 
 char command[FILENAMELENGTH];
 int  outcmd=0;
 
@@ -2291,7 +2307,7 @@ void evsij(char fileres[], double ***eij
 
     hf=hstepm*stepm/YEARM;  /* Duration of hstepm expressed in year unit. */
 
-    /* Computing Variances of health expectancies */
+    /* Computing  Variances of health expectancies */
 
      for(theta=1; theta <=npar; theta++){
       for(i=1; i<=npar; i++){ 
@@ -2845,8 +2861,8 @@ void varprob(char optionfilefiname[], do
   fprintf(fichtm,"\n<li><h4> Computing and drawing one step probabilities with their confidence intervals</h4></li>\n");
   fprintf(fichtm,"\n");
 
-  fprintf(fichtm,"\n<li><h4> <a href=\"%s\">Computing matrix of variance-covariance of step probabilities</a></h4></li>\n",optionfilehtmcov);
-  fprintf(fichtmcov,"\n<h4>Computing matrix of variance-covariance of step probabilities</h4>\n\
+  fprintf(fichtm,"\n<li><h4> <a href=\"%s\">Matrix of variance-covariance of pairs of step probabilities (drawings)</a></h4></li>\n",optionfilehtmcov);
+  fprintf(fichtmcov,"\n<h4>Matrix of variance-covariance of pairs of step probabilities</h4>\n\
   file %s<br>\n",optionfilehtmcov);
   fprintf(fichtmcov,"\nEllipsoids of confidence centered on point (p<inf>ij</inf>, p<inf>kl</inf>) are estimated\
 and drawn. It helps understanding how is the covariance between two incidences.\
@@ -2879,9 +2895,9 @@ To be simple, these graphs help to under
 	fprintf(ficgp, "**********\n#\n");
 	
 	
-	fprintf(fichtm, "\n<hr  size=\"2\" color=\"#EC5E5E\">********** Variable "); 
+	fprintf(fichtmcov, "\n<hr  size=\"2\" color=\"#EC5E5E\">********** Variable "); 
 	for (z1=1; z1<=cptcoveff; z1++) fprintf(fichtm, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtab[j1][z1]]);
-	fprintf(fichtm, "**********\n<hr size=\"2\" color=\"#EC5E5E\">");
+	fprintf(fichtmcov, "**********\n<hr size=\"2\" color=\"#EC5E5E\">");
 	
 	fprintf(ficresprobcor, "\n#********** Variable ");    
 	for (z1=1; z1<=cptcoveff; z1++) fprintf(ficresprobcor, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtab[j1][z1]]);
@@ -3090,21 +3106,19 @@ void printinghtml(char fileres[], char t
 		  double jprev1, double mprev1,double anprev1, \
 		  double jprev2, double mprev2,double anprev2){
   int jj1, k1, i1, cpt;
-  /*char optionfilehtm[FILENAMELENGTH];*/
-/*   if((fichtm=fopen(optionfilehtm,"a"))==NULL)    { */
-/*     printf("Problem with %s \n",optionfilehtm), exit(0); */
-/*     fprintf(ficlog,"Problem with %s \n",optionfilehtm), exit(0); */
-/*   } */
 
    fprintf(fichtm,"<ul><li><h4>Result files (first order: no variance)</h4>\n \
- - Observed prevalence in each state (during the period defined between %.lf/%.lf/%.lf and %.lf/%.lf/%.lf): <a href=\"%s\">%s</a> <br>\n \
- - Estimated transition probabilities over %d (stepm) months: <a href=\"%s\">%s</a><br>\n \
- - Stable prevalence in each health state: <a href=\"%s\">%s</a> <br>\n \
+ - Observed prevalence in each state (during the period defined between %.lf/%.lf/%.lf and %.lf/%.lf/%.lf): <a href=\"%s\">%s</a> <br>\n ",
+	   jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,subdirf2(fileres,"p"),subdirf2(fileres,"p"));
+   fprintf(fichtm,"\
+ - Estimated transition probabilities over %d (stepm) months: <a href=\"%s\">%s</a><br>\n ",
+	   stepm,subdirf2(fileres,"pij"),subdirf2(fileres,"pij"));
+   fprintf(fichtm,"\
+ - Stable prevalence in each health state: <a href=\"%s\">%s</a> <br>\n",
+	   subdirf2(fileres,"pl"),subdirf2(fileres,"pl"));
+   fprintf(fichtm,"\
  - Life expectancies by age and initial health status (estepm=%2d months): \
-   <a href=\"%s\">%s</a> <br>\n</li>", \
-	   jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,subdirf2(fileres,"p"),subdirf2(fileres,"p"),\
-	   stepm,subdirf2(fileres,"pij"),subdirf2(fileres,"pij"),\
-	   subdirf2(fileres,"pl"),subdirf2(fileres,"pl"),\
+   <a href=\"%s\">%s</a> <br>\n</li>",
 	   estepm,subdirf2(fileres,"e"),subdirf2(fileres,"e"));
 
 fprintf(fichtm," \n<ul><li><b>Graphs</b></li><p>");
@@ -3138,28 +3152,32 @@ fprintf(fichtm," \n<ul><li><b>Graphs</b>
         fprintf(fichtm,"\n<br>- Health life expectancies by age and initial health state (%d): %s%d%d.png <br> \
 <img src=\"%s%d%d.png\">",cpt,subdirf2(optionfilefiname,"exp"),cpt,jj1,subdirf2(optionfilefiname,"exp"),cpt,jj1);
      }
-     fprintf(fichtm,"\n<br>- Total life expectancy by age and \
-health expectancies in states (1) and (2): %s%d.png<br>\
-<img src=\"%s%d.png\">",subdirf2(optionfilefiname,"e"),jj1,subdirf2(optionfilefiname,"e"),jj1);
    } /* end i1 */
  }/* End k1 */
  fprintf(fichtm,"</ul>");
 
 
- fprintf(fichtm,"\n<br><li><h4> Result files (second order: variances)</h4>\n\
- - Parameter file with estimated parameters and covariance matrix: <a href=\"%s\">%s</a> <br>\n\
- - Variance of one-step probabilities: <a href=\"%s\">%s</a> <br>\n\
- - Variance-covariance of one-step probabilities: <a href=\"%s\">%s</a> <br>\n\
- - Correlation matrix of one-step probabilities: <a href=\"%s\">%s</a> <br>\n\
- - Variances and covariances of life expectancies by age and initial health status (estepm=%d months): <a href=\"%s\">%s</a><br>\n\
- - Health expectancies with their variances (no covariance): <a href=\"%s\">%s</a> <br>\n\
+ fprintf(fichtm,"\
+\n<br><li><h4> Result files (second order: variances)</h4>\n\
+ - Parameter file with estimated parameters and covariance matrix: <a href=\"%s\">%s</a> <br>\n", rfileres,rfileres);
+
+ fprintf(fichtm," - Variance of one-step probabilities: <a href=\"%s\">%s</a> <br>\n",
+	 subdirf2(fileres,"prob"),subdirf2(fileres,"prob"));
+ fprintf(fichtm,"\
+ - Variance-covariance of one-step probabilities: <a href=\"%s\">%s</a> <br>\n",
+	 subdirf2(fileres,"probcov"),subdirf2(fileres,"probcov"));
+
+ fprintf(fichtm,"\
+ - Correlation matrix of one-step probabilities: <a href=\"%s\">%s</a> <br>\n",
+	 subdirf2(fileres,"probcor"),subdirf2(fileres,"probcor"));
+ fprintf(fichtm,"\
+ - Variances and covariances of life expectancies by age and initial health status (estepm=%d months): <a href=\"%s\">%s</a><br>\n",
+	 estepm, subdirf2(fileres,"v"),subdirf2(fileres,"v"));
+ fprintf(fichtm,"\
+ - Health expectancies with their variances (no covariance): <a href=\"%s\">%s</a> <br>\n",
+	 subdirf2(fileres,"t"),subdirf2(fileres,"t"));
+ fprintf(fichtm,"\
  - Standard deviation of stable prevalences: <a href=\"%s\">%s</a> <br>\n",\
-	 rfileres,rfileres,\
-	 subdirf2(fileres,"prob"),subdirf2(fileres,"prob"),\
-	 subdirf2(fileres,"probcov"),subdirf2(fileres,"probcov"),\
-	 subdirf2(fileres,"probcor"),subdirf2(fileres,"probcor"),\
-	 estepm, subdirf2(fileres,"v"),subdirf2(fileres,"v"),\
-	 subdirf2(fileres,"t"),subdirf2(fileres,"t"),\
 	 subdirf2(fileres,"vpl"),subdirf2(fileres,"vpl"));
 
 /*  if(popforecast==1) fprintf(fichtm,"\n */
@@ -3168,7 +3186,8 @@ health expectancies in states (1) and (2
 /* 	<br>",fileres,fileres,fileres,fileres); */
 /*  else  */
 /*    fprintf(fichtm,"\n No population forecast: popforecast = %d (instead of 1) or stepm = %d (instead of 1) or model=%s (instead of .)<br><br></li>\n",popforecast, stepm, model); */
-fprintf(fichtm," <ul><li><b>Graphs</b></li><p>");
+ fflush(fichtm);
+ fprintf(fichtm," <ul><li><b>Graphs</b></li><p>");
 
  m=cptcoveff;
  if (cptcovn < 1) {m=1;ncodemax[1]=1;}
@@ -3184,10 +3203,13 @@ fprintf(fichtm," <ul><li><b>Graphs</b></
        fprintf(fichtm," ************\n<hr size=\"2\" color=\"#EC5E5E\">");
      }
      for(cpt=1; cpt<=nlstate;cpt++) {
-       fprintf(fichtm,"<br>- Observed and period prevalence (with confident\
-interval) in state (%d): %s%d%d.png <br>\
+       fprintf(fichtm,"<br>- Observed (cross-sectional) and period (incidence based) \
+prevalence (with 95%% confidence interval) in state (%d): %s%d%d.png <br>\
 <img src=\"%s%d%d.png\">",cpt,subdirf2(optionfilefiname,"v"),cpt,jj1,subdirf2(optionfilefiname,"v"),cpt,jj1);  
      }
+     fprintf(fichtm,"\n<br>- Total life expectancy by age and \
+health expectancies in states (1) and (2): %s%d.png<br>\
+<img src=\"%s%d.png\">",subdirf2(optionfilefiname,"e"),jj1,subdirf2(optionfilefiname,"e"),jj1);
    } /* end i1 */
  }/* End k1 */
  fprintf(fichtm,"</ul>");
@@ -3842,7 +3864,7 @@ void prwizard(int ncovmodel, int nlstate
 	} /* end k*/
       } /*end j */
     } /* end i */
-  }
+  } /* end itimes */
 
 } /* end of prwizard */
 
@@ -3855,8 +3877,11 @@ int main(int argc, char *argv[])
 {
   int movingaverage(double ***probs, double bage,double fage, double ***mobaverage, int mobilav);
   int i,j, k, n=MAXN,iter,m,size=100,cptcode, cptcod;
-  int jj, imk;
+  int jj, ll, li, lj, lk, imk;
   int numlinepar=0; /* Current linenumber of parameter file */
+  int itimes;
+
+  char ca[32], cb[32], cc[32];
   /*  FILE *fichtm; *//* Html File */
   /* FILE *ficgp;*/ /*Gnuplot File */
   double agedeb, agefin,hf;
@@ -3911,6 +3936,12 @@ int main(int argc, char *argv[])
 
   long total_usecs;
  
+/*   setlocale (LC_ALL, ""); */
+/*   bindtextdomain (PACKAGE, LOCALEDIR); */
+/*   textdomain (PACKAGE); */
+/*   setlocale (LC_CTYPE, ""); */
+/*   setlocale (LC_MESSAGES, ""); */
+
   /*   gettimeofday(&start_time, (struct timezone*)0); */ /* at first time */
   (void) gettimeofday(&start_time,&tzp);
   curr_time=start_time;
@@ -4457,6 +4488,7 @@ int main(int argc, char *argv[])
 
 }*/
 
+
   printf("Total number of individuals= %d, Agemin = %.2f, Agemax= %.2f\n\n", imx, agemin, agemax);
   fprintf(ficlog,"Total number of individuals= %d, Agemin = %.2f, Agemax= %.2f\n\n", imx, agemin, agemax); 
 
@@ -4602,27 +4634,25 @@ Interval (in months) between two waves:
   fprintf(ficres,"title=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%e stepm=%d ncovcol=%d nlstate=%d ndeath=%d maxwav=%d mle= 0 weight=%d\nmodel=%s\n", title, datafile, lastobs, firstpass,lastpass,ftol, stepm, ncovcol, nlstate, ndeath, maxwav, weightopt,model);
   
 
-  jk=1;
   fprintf(ficres,"# Parameters nlstate*nlstate*ncov a12*1 + b12 * age + ...\n");
   printf("# Parameters nlstate*nlstate*ncov a12*1 + b12 * age + ...\n");
   fprintf(ficlog,"# Parameters nlstate*nlstate*ncov a12*1 + b12 * age + ...\n");
   for(i=1,jk=1; i <=nlstate; i++){
     for(k=1; k <=(nlstate+ndeath); k++){
-      if (k != i) 
-	{
-	  printf("%d%d ",i,k);
-	  fprintf(ficlog,"%d%d ",i,k);
-	  fprintf(ficres,"%1d%1d ",i,k);
-	  for(j=1; j <=ncovmodel; j++){
-	    printf("%f ",p[jk]);
-	    fprintf(ficlog,"%f ",p[jk]);
-	    fprintf(ficres,"%f ",p[jk]);
-	    jk++; 
-	  }
-	  printf("\n");
-	  fprintf(ficlog,"\n");
-	  fprintf(ficres,"\n");
+      if (k != i) {
+	printf("%d%d ",i,k);
+	fprintf(ficlog,"%d%d ",i,k);
+	fprintf(ficres,"%1d%1d ",i,k);
+	for(j=1; j <=ncovmodel; j++){
+	  printf("%f ",p[jk]);
+	  fprintf(ficlog,"%f ",p[jk]);
+	  fprintf(ficres,"%f ",p[jk]);
+	  jk++; 
 	}
+	printf("\n");
+	fprintf(ficlog,"\n");
+	fprintf(ficres,"\n");
+      }
     }
   }
   if(mle!=0){
@@ -4653,32 +4683,91 @@ Interval (in months) between two waves:
   }
    
   fprintf(ficres,"# Covariance matrix \n# 121 Var(a12)\n# 122 Cov(b12,a12) Var(b12)\n#   ...\n# 232 Cov(b23,a12)  Cov(b23,b12) ... Var (b23)\n");
-  if(mle==1)
+  if(mle>=1)
     printf("# Covariance matrix \n# 121 Var(a12)\n# 122 Cov(b12,a12) Var(b12)\n#   ...\n# 232 Cov(b23,a12)  Cov(b23,b12) ... Var (b23)\n");
   fprintf(ficlog,"# Covariance matrix \n# 121 Var(a12)\n# 122 Cov(b12,a12) Var(b12)\n#   ...\n# 232 Cov(b23,a12)  Cov(b23,b12) ... Var (b23)\n");
-  for(i=1,k=1;i<=npar;i++){
-    /*  if (k>nlstate) k=1;
-	i1=(i-1)/(ncovmodel*nlstate)+1; 
-	fprintf(ficres,"%s%d%d",alph[k],i1,tab[i]);
-	printf("%s%d%d",alph[k],i1,tab[i]);
-    */
-    fprintf(ficres,"%3d",i);
-    if(mle==1)
-      printf("%3d",i);
-    fprintf(ficlog,"%3d",i);
-    for(j=1; j<=i;j++){
-      fprintf(ficres," %.5e",matcov[i][j]);
-      if(mle==1)
-	printf(" %.5e",matcov[i][j]);
-      fprintf(ficlog," %.5e",matcov[i][j]);
-    }
-    fprintf(ficres,"\n");
-    if(mle==1)
-      printf("\n");
-    fprintf(ficlog,"\n");
-    k++;
-  }
-   
+/* # 121 Var(a12)\n\ */
+/* # 122 Cov(b12,a12) Var(b12)\n\ */
+/* # 131 Cov(a13,a12) Cov(a13,b12, Var(a13)\n\ */
+/* # 132 Cov(b13,a12) Cov(b13,b12, Cov(b13,a13) Var(b13)\n\ */
+/* # 212 Cov(a21,a12) Cov(a21,b12, Cov(a21,a13) Cov(a21,b13) Var(a21)\n\ */
+/* # 212 Cov(b21,a12) Cov(b21,b12, Cov(b21,a13) Cov(b21,b13) Cov(b21,a21) Var(b21)\n\ */
+/* # 232 Cov(a23,a12) Cov(a23,b12, Cov(a23,a13) Cov(a23,b13) Cov(a23,a21) Cov(a23,b21) Var(a23)\n\ */
+/* # 232 Cov(b23,a12) Cov(b23,b12) ... Var (b23)\n" */
+
+
+/* Just to have a covariance matrix which will be more understandable
+   even is we still don't want to manage dictionary of variables
+*/
+  for(itimes=1;itimes<=2;itimes++){
+    jj=0;
+    for(i=1; i <=nlstate; i++){
+      for(j=1; j <=nlstate+ndeath; j++){
+	if(j==i) continue;
+	for(k=1; k<=ncovmodel;k++){
+	  jj++;
+	  ca[0]= k+'a'-1;ca[1]='\0';
+	  if(itimes==1){
+	    if(mle>=1)
+	      printf("#%1d%1d%d",i,j,k);
+	    fprintf(ficlog,"#%1d%1d%d",i,j,k);
+	    fprintf(ficres,"#%1d%1d%d",i,j,k);
+	  }else{
+	    if(mle>=1)
+	      printf("%1d%1d%d",i,j,k);
+	    fprintf(ficlog,"%1d%1d%d",i,j,k);
+	    fprintf(ficres,"%1d%1d%d",i,j,k);
+	  }
+	  ll=0;
+	  for(li=1;li <=nlstate; li++){
+	    for(lj=1;lj <=nlstate+ndeath; lj++){
+	      if(lj==li) continue;
+	      for(lk=1;lk<=ncovmodel;lk++){
+		ll++;
+		if(ll<=jj){
+		  cb[0]= lk +'a'-1;cb[1]='\0';
+		  if(ll<jj){
+		    if(itimes==1){
+		      if(mle>=1)
+			printf(" Cov(%s%1d%1d,%s%1d%1d)",ca,i,j,cb, li,lj);
+		      fprintf(ficlog," Cov(%s%1d%1d,%s%1d%1d)",ca,i,j,cb, li,lj);
+		      fprintf(ficres," Cov(%s%1d%1d,%s%1d%1d)",ca,i,j,cb, li,lj);
+		    }else{
+		      if(mle>=1)
+			printf(" %.5e",matcov[jj][ll]); 
+		      fprintf(ficlog," %.5e",matcov[jj][ll]); 
+		      fprintf(ficres," %.5e",matcov[jj][ll]); 
+		    }
+		  }else{
+		    if(itimes==1){
+		      if(mle>=1)
+			printf(" Var(%s%1d%1d)",ca,i,j);
+		      fprintf(ficlog," Var(%s%1d%1d)",ca,i,j);
+		      fprintf(ficres," Var(%s%1d%1d)",ca,i,j);
+		    }else{
+		      if(mle>=1)
+			printf(" %.5e",matcov[jj][ll]); 
+		      fprintf(ficlog," %.5e",matcov[jj][ll]); 
+		      fprintf(ficres," %.5e",matcov[jj][ll]); 
+		    }
+		  }
+		}
+	      } /* end lk */
+	    } /* end lj */
+	  } /* end li */
+	  if(mle>=1)
+	    printf("\n");
+	  fprintf(ficlog,"\n");
+	  fprintf(ficres,"\n");
+	  numlinepar++;
+	} /* end k*/
+      } /*end j */
+    } /* end i */
+  } /* end itimes */
+
+  fflush(ficlog);
+  fflush(ficres);
+
   while((c=getc(ficpar))=='#' && c!= EOF){
     ungetc(c,ficpar);
     fgets(line, MAXLINE, ficpar);
@@ -4752,9 +4841,9 @@ Interval (in months) between two waves:
   }
   ungetc(c,ficpar);
 
-  fscanf(ficpar,"popforecast=%d popfile=%s popfiledate=%lf/%lf/%lf last-popfiledate=%lf/%lf/%lf\n",&popforecast,popfile,&jpyram,&mpyram,&anpyram,&jpyram1,&mpyram1,&anpyram1);
+  /*  fscanf(ficpar,"popforecast=%d popfile=%s popfiledate=%lf/%lf/%lf last-popfiledate=%lf/%lf/%lf\n",&popforecast,popfile,&jpyram,&mpyram,&anpyram,&jpyram1,&mpyram1,&anpyram1);
   fprintf(ficparo,"popforecast=%d popfile=%s popfiledate=%.lf/%.lf/%.lf last-popfiledate=%.lf/%.lf/%.lf\n",popforecast,popfile,jpyram,mpyram,anpyram,jpyram1,mpyram1,anpyram1);
-  fprintf(ficres,"popforecast=%d popfile=%s popfiledate=%.lf/%.lf/%.lf last-popfiledate=%.lf/%.lf/%.lf\n",popforecast,popfile,jpyram,mpyram,anpyram,jpyram1,mpyram1,anpyram1);
+  fprintf(ficres,"popforecast=%d popfile=%s popfiledate=%.lf/%.lf/%.lf last-popfiledate=%.lf/%.lf/%.lf\n",popforecast,popfile,jpyram,mpyram,anpyram,jpyram1,mpyram1,anpyram1);*/
 
   /*  freqsummary(fileres, agemin, agemax, s, agev, nlstate, imx,Tvaraff,nbcode, ncodemax,mint,anint);*/
   /*,dateprev1,dateprev2,jprev1, mprev1,anprev1,jprev2, mprev2,anprev2);*/
@@ -4893,6 +4982,10 @@ Interval (in months) between two waves:
   fclose(ficrespij);
 
   probs= ma3x(1,AGESUP,1,NCOVMAX, 1,NCOVMAX);
+  for(i=1;i<=AGESUP;i++)
+    for(j=1;j<=NCOVMAX;j++)
+      for(k=1;k<=NCOVMAX;k++)
+	probs[i][j][k]=0.;
 
   /*---------- Forecasting ------------------*/
   /*if((stepm == 1) && (strcmp(model,".")==0)){*/
@@ -5089,7 +5182,7 @@ ageminpar, agemax, s[lastpass][imx], age
   
 
   if((nberr >0) || (nbwarn>0)){
-    printf("End of Imach with %d errors and/or warnings %d\n",nberr,nbwarn);
+    printf("End of Imach with %d errors and/or %d warnings\n",nberr,nbwarn);
     fprintf(ficlog,"End of Imach with %d errors and/or warnings %d\n",nberr,nbwarn);
   }else{
     printf("End of Imach\n");