bug fix in cost calculation for using satellite observations with the analytic method

source/mod_analytic.f90

@@ -587,10 +587,12 @@ module mod_analytic
     cost_o = 0.
     do i = 1, nobs
       if ( trim(config%spec).eq.'co2' ) then
-        obs%delta(i) = obs%model(i) + obs%nee(i) + obs%fff(i) + obs%ocn(i) - &
+        obs%delta(i) = obs%cpri(i) - obs%cakpri(i) + &
+                        obs%model(i) + obs%nee(i) + obs%fff(i) + obs%ocn(i) - &
                         obs%conc(i) + obs%bkg(i) + obs%cinipos(i)
       else
-        obs%delta(i) = obs%model(i) + obs%ghg(i) - obs%conc(i) + obs%bkg(i) + obs%cinipos(i)
+        obs%delta(i) = obs%cpri(i) - obs%cakpri(i) + &
+                        obs%model(i) + obs%ghg(i) - obs%conc(i) + obs%bkg(i) + obs%cinipos(i)
       endif
       cost_o = cost_o + obs%delta(i)**2/obs%err(i)**2
     end do