From 5eed8c43c387d9b00dd2c1712d06aff8e9e1ea0f Mon Sep 17 00:00:00 2001
From: ronesy <rlt@nilu.no>
Date: Wed, 30 Oct 2024 14:39:01 +0100
Subject: [PATCH] bug fix in cost calculation for using satellite observations
with the analytic method
---
source/mod_analytic.f90 | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/source/mod_analytic.f90 b/source/mod_analytic.f90
index b625ba1..e79d9f7 100644
--- a/source/mod_analytic.f90
+++ b/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
--
GitLab