Update calc_SNR2

Added error if sig<=0
Now SNR is calculated for each obsravtion
Added condition that SNR>=0
Update model1

Added initialization of A_mod and fi_mod