--- imach/src/imach.c	2015/12/16 06:57:54	1.214
+++ imach/src/imach.c	2015/12/16 08:52:24	1.215
@@ -1,6 +1,9 @@
-/* $Id: imach.c,v 1.214 2015/12/16 06:57:54 brouard Exp $
+/* $Id: imach.c,v 1.215 2015/12/16 08:52:24 brouard Exp $
   $State: Exp $
   $Log: imach.c,v $
+  Revision 1.215  2015/12/16 08:52:24  brouard
+  Summary: 0.98r4 working
+
   Revision 1.214  2015/12/16 06:57:54  brouard
   Summary: temporary not working
 
@@ -789,12 +792,12 @@ typedef struct {
 #define ODIRSEPARATOR '\\'
 #endif
 
-/* $Id: imach.c,v 1.214 2015/12/16 06:57:54 brouard Exp $ */
+/* $Id: imach.c,v 1.215 2015/12/16 08:52:24 brouard Exp $ */
 /* $State: Exp $ */
 #include "version.h"
 char version[]=__IMACH_VERSION__;
 char copyright[]="October 2015,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015";
-char fullversion[]="$Revision: 1.214 $ $Date: 2015/12/16 06:57:54 $"; 
+char fullversion[]="$Revision: 1.215 $ $Date: 2015/12/16 08:52:24 $"; 
 char strstart[80];
 char optionfilext[10], optionfilefiname[FILENAMELENGTH];
 int erreur=0, nberr=0, nbwarn=0; /* Error number, number of errors number of warnings  */
@@ -5704,8 +5707,8 @@ void prevforecast(char fileres[], double
     printf("Problem with forecast resultfile: %s\n", fileresf);
     fprintf(ficlog,"Problem with forecast resultfile: %s\n", fileresf);
   }
-  printf("Computing forecasting: result on file '%s' \n", fileresf);
-  fprintf(ficlog,"Computing forecasting: result on file '%s' \n", fileresf);
+  printf("Computing forecasting: result on file '%s', please wait... \n", fileresf);
+  fprintf(ficlog,"Computing forecasting: result on file '%s', please wait... \n", fileresf);
 
   if (cptcoveff==0) ncodemax[cptcoveff]=1;
 
@@ -5800,6 +5803,9 @@ void prevforecast(char fileres[], double
   if (mobilav!=0) free_ma3x(mobaverage,1, AGESUP,1,NCOVMAX, 1,NCOVMAX);
 
   fclose(ficresf);
+  printf("End of Computing forecasting \n");
+  fprintf(ficlog,"End of Computing forecasting\n");
+
 }
 
 /************** Forecasting *****not tested NB*************/
@@ -7685,8 +7691,8 @@ Please run with mle=-1 to get a correct
   free_vector(annais,1,n);
   /* free_matrix(mint,1,maxwav,1,n);
      free_matrix(anint,1,maxwav,1,n);*/
-  free_vector(moisdc,1,n);
-  free_vector(andc,1,n);
+  /* free_vector(moisdc,1,n); */
+  /* free_vector(andc,1,n); */
   /* */
   
   wav=ivector(1,imx);
@@ -7708,6 +7714,9 @@ Please run with mle=-1 to get a correct
   concatwav(wav, dh, bh, mw, s, agedc, agev,  firstpass, lastpass, imx, nlstate, stepm);
   /* */
  
+  free_vector(moisdc,1,n);
+  free_vector(andc,1,n);
+
   /* Routine tricode is to calculate cptcoveff (real number of unique covariates) and to associate covariable number and modality */
 
   nbcode=imatrix(0,NCOVMAX,0,NCOVMAX); 
@@ -8448,10 +8457,10 @@ Please run with mle=-1 to get a correct
    /*------------ free_vector  -------------*/
    /*  chdir(path); */
  
-    free_ivector(wav,1,imx);
-    free_imatrix(dh,1,lastpass-firstpass+2,1,imx);
-    free_imatrix(bh,1,lastpass-firstpass+2,1,imx);
-    free_imatrix(mw,1,lastpass-firstpass+2,1,imx);   
+    /* free_ivector(wav,1,imx); */  /* Moved after last prevalence call */
+    /* free_imatrix(dh,1,lastpass-firstpass+2,1,imx); */
+    /* free_imatrix(bh,1,lastpass-firstpass+2,1,imx); */
+    /* free_imatrix(mw,1,lastpass-firstpass+2,1,imx);    */
     free_lvector(num,1,n);
     free_vector(agedc,1,n);
     /*free_matrix(covar,0,NCOVMAX,1,n);*/
@@ -8511,6 +8520,11 @@ Please run with mle=-1 to get a correct
     /*  printf("ageminpar=%f, agemax=%f, s[lastpass][imx]=%d, agev[lastpass][imx]=%f, nlstate=%d, imx=%d,  mint[lastpass][imx]=%f, anint[lastpass][imx]=%f,dateprev1=%f, dateprev2=%f, firstpass=%d, lastpass=%d\n",\
 	ageminpar, agemax, s[lastpass][imx], agev[lastpass][imx], nlstate, imx, mint[lastpass][imx],anint[lastpass][imx], dateprev1, dateprev2, firstpass, lastpass);
     */
+    free_ivector(wav,1,imx);
+    free_imatrix(dh,1,lastpass-firstpass+2,1,imx);
+    free_imatrix(bh,1,lastpass-firstpass+2,1,imx);
+    free_imatrix(mw,1,lastpass-firstpass+2,1,imx);   
+
 
     if (mobilav!=0) {
       mobaverage= ma3x(1, AGESUP,1,NCOVMAX, 1,NCOVMAX);