Check sigma and sigma2 >0

src/readBayesModels.c

@@ -294,6 +294,11 @@ void    setBayesModel( long int iVal, long int nVal, double **array)
                 if( sblock[i][ipx] != 0 )
                     o_set_source( &source[i], ipx, method(nVal, iParam++, array));
 
+    	    if((sblock[i][SA] !=0) & (source[i].E.a<0))
+		    source[i].E.a=-source[i].E.a;
+    	    if((sblock[i][SA2] !=0) & (source[i].E2.a<0))
+		    source[i].E2.a=-source[i].E2.a;
+
             if((sblock[i][SEPS] != 0 )|| (sblock[i][SA] !=0))
                 source[i].E.b = source[i].E.a * (1. - source[i].eps) / (1. + source[i].eps);
             if((sblock[i][SEPS2] != 0 )|| (sblock[i][SA2] !=0))