Commit 2d26e4a3 authored by ronesy's avatar ronesy
Browse files

Made scaling of prior error covariance matrix optional in SETTINGS_config

parent 1278c7be
......@@ -121,5 +121,6 @@ sigma_ocean: 1000.
sigmatime: 7.
# Total error for domain (Tg/y)
# if globerr <= 0 prior error covariance matrix not scaled
globerr: 300.
......@@ -121,5 +121,6 @@ sigma_ocean: 1000.
sigmatime: 7.
# Total error for domain (Tg/y)
# if globerr <= 0 prior error covariance matrix not scaled
globerr: 300.
......@@ -113,5 +113,6 @@ sigma_ocean: 1000.
sigmatime: 30.
# Total error for inversion domain (Tg/y)
# if globerr <= 0 prior error covariance matrix not scaled
globerr: 10.
......@@ -113,5 +113,6 @@ sigma_ocean: 1000.
sigmatime: 30.
# Total error for inversion domain (Tg/y)
# if globerr <= 0 prior error covariance matrix not scaled
globerr: 10.
......@@ -93,11 +93,12 @@ subroutine error_cov(config, files, states, covar, corr, xerr)
! scale to total error of nested domain (unit Tg/y)
toterr = sqrt(covsum)*3600.*24.*365./1.e9/numscale
erscalar = config%globerr/toterr
cov = cov * dble(erscalar)**2
write(logid,*) 'Total error: ',toterr
write(logid,*) 'Total error scaled by: ',erscalar
erscalar = config%globerr/toterr
if ( config%globerr > 0. ) then
cov = cov * dble(erscalar)**2
write(logid,*) 'Total error scaled by: ',erscalar
end if
! calculate eigendecomposition of cov
! -----------------------------------
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment