From b28c96f4923dc26e249d61c9fe036cd49843e5e4 Mon Sep 17 00:00:00 2001
From: ronesy <rlt@nilu.no>
Date: Wed, 8 Dec 2021 15:57:25 +0100
Subject: [PATCH] Correction to posterior error calculation for scalars of
initial concentration in congrad
---
source/congrad.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/source/congrad.f90 b/source/congrad.f90
index ca656a1..e88fd53 100644
--- a/source/congrad.f90
+++ b/source/congrad.f90
@@ -403,7 +403,7 @@ subroutine congrad(iter, grad, files, config, fluxes, obs, states, covar)
end do
end do
if ( config%opt_cini ) then
- tmp(npvar+1:nvar) = tmp(npvar+1:nvar) + config%cinierr * pevecs(jk,npvar+1:nvar)
+ tmp(npvar+1:nvar) = tmp(npvar+1:nvar) + (config%cinierr * pevecs(jk,npvar+1:nvar))**2
endif
end do
! calculate last part
--
GitLab