From ece87145ed016ab0a28c1f0b5df7c8cfe7c9441b Mon Sep 17 00:00:00 2001
From: ronesy <rlt@nilu.no>
Date: Wed, 30 Oct 2024 14:49:19 +0100
Subject: [PATCH] change to observation space uncertainty for when background
is optimized
---
source/calc_conc.f90 | 12 ++++++++++--
1 file changed, 10 insertions(+), 2 deletions(-)
diff --git a/source/calc_conc.f90 b/source/calc_conc.f90
index b70c602..98f266c 100644
--- a/source/calc_conc.f90
+++ b/source/calc_conc.f90
@@ -150,10 +150,18 @@ subroutine calc_conc(config, fluxes, obs, ngrid, gtime, hnest, hbkg, iobs, ix1,
if ( trim(config%spec).eq.'co2' ) then
! CO2
- obs%err(iobs) = sqrt(obs%measerr(iobs)**2 + (sum(obs%cini(iobs,:))*config%cinierr)**2 + bkgerr + ffferr)
+ if ( config%opt_cini ) then
+ obs%err(iobs) = sqrt(obs%measerr(iobs)**2 + bkgerr + ffferr)
+ else
+ obs%err(iobs) = sqrt(obs%measerr(iobs)**2 + (sum(obs%cini(iobs,:))*config%cinierr)**2 + bkgerr + ffferr)
+ endif
else if ( trim(config%spec).eq.'ghg' ) then
! GHG
- obs%err(iobs) = sqrt(obs%measerr(iobs)**2 + (sum(obs%cini(iobs,:))*config%cinierr)**2 + bkgerr)
+ if ( config%opt_cini ) then
+ obs%err(iobs) = sqrt(obs%measerr(iobs)**2 + bkgerr)
+ else
+ obs%err(iobs) = sqrt(obs%measerr(iobs)**2 + (sum(obs%cini(iobs,:))*config%cinierr)**2 + bkgerr)
+ endif
else if ( trim(config%spec).eq.'aero' ) then
! aerosols
obs%err(iobs) = sqrt(obs%measerr(iobs)**2 + bkgerr)
--
GitLab