diff --git a/test/glb/testinput/4densvar.yaml b/test/glb/testinput/4densvar.yaml index f69aa04..619e18a 100644 --- a/test/glb/testinput/4densvar.yaml +++ b/test/glb/testinput/4densvar.yaml @@ -2,7 +2,7 @@ cost function: cost type: 4D-Ens-Var time window: begin: 2010-01-01T00:00:00Z - length: PT06H + length: PT6H subwindow: PT3H parallel subwindows: false analysis variables: [air_temperature] diff --git a/test/glb/testinput/4densvar_parallel.yaml b/test/glb/testinput/4densvar_parallel.yaml index 27a7b7e..a8eee0d 100644 --- a/test/glb/testinput/4densvar_parallel.yaml +++ b/test/glb/testinput/4densvar_parallel.yaml @@ -2,7 +2,7 @@ cost function: cost type: 4D-Ens-Var time window: begin: 2010-01-01T00:00:00Z - length: PT06H + length: PT6H subwindow: PT3H analysis variables: [air_temperature] geometry: &geom diff --git a/test/glb/testinput/4dvar.yaml b/test/glb/testinput/4dvar.yaml index 68f97e2..6a254dd 100644 --- a/test/glb/testinput/4dvar.yaml +++ b/test/glb/testinput/4dvar.yaml @@ -89,6 +89,6 @@ final: output: filepath: testdata/4dvar_analysis__MPI___DATE_ date pattern: _DATE_ - frequency: PT06H + frequency: PT6H test: reference filename: testref/4dvar.ref diff --git a/test/glb/testinput/letkf_linear_4d.yaml b/test/glb/testinput/letkf_linear_4d.yaml index 1ce1901..2b3259f 100644 --- a/test/glb/testinput/letkf_linear_4d.yaml +++ b/test/glb/testinput/letkf_linear_4d.yaml @@ -16,6 +16,9 @@ background: members from template: template: states: + - date: 2010-01-01T00:00:00Z + state variables: [air_temperature] + filepath: testdata/ens__MPI__PT0S_%mem% - date: 2010-01-01T03:00:00Z state variables: [air_temperature] filepath: testdata/ens__MPI__PT3H_%mem% @@ -69,6 +72,8 @@ local ensemble DA: use linear observer: true output: states: + - date: 2010-01-01T00:00:00Z + filepath: testdata/letkf_linear_4d_state__MPI__PT0S - date: 2010-01-01T03:00:00Z filepath: testdata/letkf_linear_4d_state__MPI__PT3H - date: 2010-01-01T06:00:00Z @@ -79,6 +84,8 @@ output ensemble increments: filepath: testdata/letkf_linear_4d_increment__MPI_ output mean prior: states: + - date: 2010-01-01T00:00:00Z + filepath: testdata/letkf_linear_4d_mean_prior__MPI__PT0S - date: 2010-01-01T03:00:00Z filepath: testdata/letkf_linear_4d_mean_prior__MPI__PT3H - date: 2010-01-01T06:00:00Z diff --git a/test/glb/testinput/letkf_nonlinear_4d.yaml b/test/glb/testinput/letkf_nonlinear_4d.yaml index 88308c0..26a9e88 100644 --- a/test/glb/testinput/letkf_nonlinear_4d.yaml +++ b/test/glb/testinput/letkf_nonlinear_4d.yaml @@ -16,6 +16,9 @@ background: members from template: template: states: + - date: 2010-01-01T00:00:00Z + state variables: [air_temperature] + filepath: testdata/ens__MPI__PT0S_%mem% - date: 2010-01-01T03:00:00Z state variables: [air_temperature] filepath: testdata/ens__MPI__PT3H_%mem% @@ -68,6 +71,8 @@ local ensemble DA: mult: 1.1 output: states: + - date: 2010-01-01T00:00:00Z + filepath: testdata/letkf_nonlinear_4d_state__MPI__PT0S - date: 2010-01-01T03:00:00Z filepath: testdata/letkf_nonlinear_4d_state__MPI__PT3H - date: 2010-01-01T06:00:00Z @@ -78,6 +83,8 @@ output ensemble increments: filepath: testdata/letkf_nonlinear_4d_increment__MPI_ output mean prior: states: + - date: 2010-01-01T00:00:00Z + filepath: testdata/letkf_nonlinear_4d_mean_prior__MPI__PT0S - date: 2010-01-01T03:00:00Z filepath: testdata/letkf_nonlinear_4d_mean_prior__MPI__PT3H - date: 2010-01-01T06:00:00Z diff --git a/test/glb/testinput/makeobs_4d.yaml b/test/glb/testinput/makeobs_4d.yaml index 40e9965..79e9a0f 100644 --- a/test/glb/testinput/makeobs_4d.yaml +++ b/test/glb/testinput/makeobs_4d.yaml @@ -31,13 +31,17 @@ observations: engine: type: GenList lats: [-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0, + -60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0, -60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0,-60.0,-30.0,0.0,30.0,60.0] lons: [0.0,0.0,0.0,0.0,0.0,60.0,60.0,60.0,60.0,60.0,120.0,120.0,120.0,120.0,120.0,180.0,180.0,180.0,180.0,180.0,240.0,240.0,240.0,240.0,240.0,300.0,300.0,300.0,300.0,300.0, + 0.0,0.0,0.0,0.0,0.0,60.0,60.0,60.0,60.0,60.0,120.0,120.0,120.0,120.0,120.0,180.0,180.0,180.0,180.0,180.0,240.0,240.0,240.0,240.0,240.0,300.0,300.0,300.0,300.0,300.0, 0.0,0.0,0.0,0.0,0.0,60.0,60.0,60.0,60.0,60.0,120.0,120.0,120.0,120.0,120.0,180.0,180.0,180.0,180.0,180.0,240.0,240.0,240.0,240.0,240.0,300.0,300.0,300.0,300.0,300.0] - dateTimes: [10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800, + dateTimes: [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, + 10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800, 21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600] vert coord type: height vert coords: [1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3, + 1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3, 1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3] epoch: seconds since 2010-01-01T00:00:00Z obs errors: [0.1] diff --git a/test/glb/testref/4densvar.ref b/test/glb/testref/4densvar.ref index a6028e9..113d75f 100644 --- a/test/glb/testref/4densvar.ref +++ b/test/glb/testref/4densvar.ref @@ -1,32 +1,32 @@ -CostJo : Nonlinear Jo(air_temperature) = 4.4054587285983096e+01, nobs = 60, Jo/n = 7.3424312143305159e-01, err = 1.0000000149011609e-01 +CostJo : Nonlinear Jo(air_temperature) = 7.1581287932829468e+01, nobs = 60, Jo/n = 1.1930214655471578e+00, err = 1.0000000149011609e-01 CostJb : Nonlinear Jb = 0.0000000000000000e+00 -CostFunction: Nonlinear J = 4.4054587285983096e+01 -PCGMinimizer: reduction in residual norm = 1.6618873406073143e-02 +CostFunction: Nonlinear J = 7.1581287932829468e+01 +PCGMinimizer: reduction in residual norm = 1.6111967195909689e-02 CostFunction::addIncrement: Analysis: - Valid time: 2010-01-01T00:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -2.9638508919518216e+00 - + max = 1.0452529273299913e+01 - + mean = 2.9096803545872811e+00 - + stddev = 2.4292445733327712e+00 + + min = -2.9644582222915727e+00 + + max = 1.0452190000356063e+01 + + mean = 2.9109356463691181e+00 + + stddev = 2.4293809396571047e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7350390690600599e+00 - + max = 1.2026914686261097e+01 - + mean = 2.9241588464668880e+00 - + stddev = 2.4742857781092926e+00 + + min = -3.7369111821520757e+00 + + max = 1.2025150523460592e+01 + + mean = 2.9254141382487311e+00 + + stddev = 2.4754122645352763e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.1923878727186010e+00 - + max = 1.1772997316842755e+01 - + mean = 2.8998842447663544e+00 - + stddev = 2.5147091409180176e+00 -CostJo : Nonlinear Jo(air_temperature) = 2.4054842703565615e+01, nobs = 60, Jo/n = 4.0091404505942690e-01, err = 1.0000000149011609e-01 -CostJb : Nonlinear Jb = 6.1253934678079780e+00 -CostFunction: Nonlinear J = 3.0180236171373593e+01 + + min = -3.1911907597327307e+00 + + max = 1.1771310192660550e+01 + + mean = 2.9011395365481851e+00 + + stddev = 2.5142939478344926e+00 +CostJo : Nonlinear Jo(air_temperature) = 2.9177965344398764e+01, nobs = 60, Jo/n = 4.8629942240664609e-01, err = 1.0000000149011609e-01 +CostJb : Nonlinear Jb = 1.2912994390392910e+01 +CostFunction: Nonlinear J = 4.2090959734791674e+01 diff --git a/test/glb/testref/4densvar_parallel.ref b/test/glb/testref/4densvar_parallel.ref index ea85d5f..516949d 100644 --- a/test/glb/testref/4densvar_parallel.ref +++ b/test/glb/testref/4densvar_parallel.ref @@ -1,32 +1,32 @@ -CostJo : Nonlinear Jo(air_temperature) = 4.4054587285983104e+01, nobs = 60, Jo/n = 7.3424312143305170e-01, err = 1.0000000149011612e-01 +CostJo : Nonlinear Jo(air_temperature) = 7.1581287932829468e+01, nobs = 60, Jo/n = 1.1930214655471578e+00, err = 1.0000000149011612e-01 CostJb : Nonlinear Jb = 0.0000000000000000e+00 -CostFunction: Nonlinear J = 4.4054587285983104e+01 -PCGMinimizer: reduction in residual norm = 1.6618873406073143e-02 +CostFunction: Nonlinear J = 7.1581287932829468e+01 +PCGMinimizer: reduction in residual norm = 1.6111967195909689e-02 CostFunction::addIncrement: Analysis: - Valid time: 2010-01-01T00:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -2.9638508919518216e+00 - + max = 1.0452529273299913e+01 - + mean = 2.9096803545872811e+00 - + stddev = 2.4292445733327712e+00 + + min = -2.9644582222915727e+00 + + max = 1.0452190000356063e+01 + + mean = 2.9109356463691181e+00 + + stddev = 2.4293809396571047e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7350390690600599e+00 - + max = 1.2026914686261097e+01 - + mean = 2.9241588464668880e+00 - + stddev = 2.4742857781092926e+00 + + min = -3.7369111821520757e+00 + + max = 1.2025150523460592e+01 + + mean = 2.9254141382487311e+00 + + stddev = 2.4754122645352763e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.1923878727186010e+00 - + max = 1.1772997316842755e+01 - + mean = 2.8998842447663544e+00 - + stddev = 2.5147091409180176e+00 -CostJo : Nonlinear Jo(air_temperature) = 2.4054842703565615e+01, nobs = 60, Jo/n = 4.0091404505942690e-01, err = 1.0000000149011612e-01 -CostJb : Nonlinear Jb = 6.1253934678079780e+00 -CostFunction: Nonlinear J = 3.0180236171373593e+01 + + min = -3.1911907597327307e+00 + + max = 1.1771310192660550e+01 + + mean = 2.9011395365481851e+00 + + stddev = 2.5142939478344926e+00 +CostJo : Nonlinear Jo(air_temperature) = 2.9177965344398761e+01, nobs = 60, Jo/n = 4.8629942240664603e-01, err = 1.0000000149011612e-01 +CostJb : Nonlinear Jb = 1.2912994390392910e+01 +CostFunction: Nonlinear J = 4.2090959734791667e+01 diff --git a/test/glb/testref/4dvar.ref b/test/glb/testref/4dvar.ref index 81cd9af..94ca7e2 100644 --- a/test/glb/testref/4dvar.ref +++ b/test/glb/testref/4dvar.ref @@ -1,17 +1,17 @@ -CostJo : Nonlinear Jo(air_temperature) = 5.2823922975696632e+01, nobs = 60, Jo/n = 8.8039871626161059e-01, err = 1.0000000149011609e-01 +CostJo : Nonlinear Jo(air_temperature) = 5.7726742375330623e+01, nobs = 60, Jo/n = 9.6211237292217711e-01, err = 1.0000000149011609e-01 Norm of input parameter StdDev: 3.9038442591885909e+00 CostJb : Nonlinear Jb = 0.0000000000000000e+00 -CostFunction: Nonlinear J = 5.2823922975696632e+01 -PCGMinimizer: reduction in residual norm = 1.0461621924776851e-02 +CostFunction: Nonlinear J = 5.7726742375330623e+01 +PCGMinimizer: reduction in residual norm = 1.1374486331151328e-02 CostFunction::addIncrement: Analysis: - Valid time: 2010-01-01T00:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -2.9464302173552590e+00 - + max = 1.0512453122031689e+01 - + mean = 2.9117600857746280e+00 - + stddev = 2.4307363499441710e+00 -CostJo : Nonlinear Jo(air_temperature) = 2.3645017642682241e+01, nobs = 60, Jo/n = 3.9408362737803737e-01, err = 1.0000000149011609e-01 -CostJb : Nonlinear Jb = 9.8486776302394912e+00 -CostFunction: Nonlinear J = 3.3493695272921734e+01 + + min = -2.9410216859986762e+00 + + max = 1.0470338268671487e+01 + + mean = 2.9005605440000939e+00 + + stddev = 2.4293505392687185e+00 +CostJo : Nonlinear Jo(air_temperature) = 2.7980203587538721e+01, nobs = 60, Jo/n = 4.6633672645897867e-01, err = 1.0000000149011609e-01 +CostJb : Nonlinear Jb = 9.7315176683962754e+00 +CostFunction: Nonlinear J = 3.7711721255934997e+01 diff --git a/test/glb/testref/letkf_linear_4d.ref b/test/glb/testref/letkf_linear_4d.ref index 04082af..839520d 100644 --- a/test/glb/testref/letkf_linear_4d.ref +++ b/test/glb/testref/letkf_linear_4d.ref @@ -1,249 +1,363 @@ Initial state for member 1: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9070402095754342e+00 + + max = 1.0443230228491029e+01 + + mean = 2.8890224246311544e+00 + + stddev = 2.4337264386471138e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7984786781412843e+00 - + max = 1.2105572354693958e+01 + + min = -3.7984786781412785e+00 + + max = 1.2105572354693960e+01 + mean = 2.9035009165107679e+00 - + stddev = 2.4802542460785419e+00 + + stddev = 2.4802542460785428e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.2324168099554287e+00 + + min = -3.2324168099554305e+00 + max = 1.1684028286658434e+01 + mean = 2.8792263148102166e+00 + stddev = 2.5189832758053812e+00 Initial state for member 2: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9169112429360231e+00 + + max = 1.0640558928014917e+01 + + mean = 2.9050687446303867e+00 + + stddev = 2.4304624544932087e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.4686154982654220e+00 - + max = 1.1912098410366392e+01 + + min = -3.4686154982654163e+00 + + max = 1.1912098410366394e+01 + mean = 2.9195472365100055e+00 + stddev = 2.4687036163275997e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.2664917313325974e+00 + + min = -3.2664917313325992e+00 + max = 1.1761288355765901e+01 + mean = 2.8952726348094577e+00 + stddev = 2.5149762090011363e+00 Initial state for member 3: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.8439500143940992e+00 + + max = 1.0440496906472719e+01 + + mean = 2.9144223410218180e+00 + + stddev = 2.4306939104918670e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.8669960664368079e+00 - + max = 1.2115786341942979e+01 - + mean = 2.9289008329014217e+00 - + stddev = 2.4776428547666072e+00 + + min = -3.8669960664368022e+00 + + max = 1.2115786341942981e+01 + + mean = 2.9289008329014226e+00 + + stddev = 2.4776428547666063e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.2850289535587947e+00 + + min = -3.2850289535587964e+00 + max = 1.1990238825003495e+01 - + mean = 2.9046262312008806e+00 - + stddev = 2.5193094400446165e+00 + + mean = 2.9046262312008815e+00 + + stddev = 2.5193094400446161e+00 Initial state for member 4: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9719403295366162e+00 + + max = 1.0515799598376681e+01 + + mean = 2.9227258152368534e+00 + + stddev = 2.4350683503339825e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7182596232436875e+00 - + max = 1.2195305011408491e+01 - + mean = 2.9372043071164575e+00 - + stddev = 2.4773099230201883e+00 + + min = -3.7182596232436844e+00 + + max = 1.2195305011408493e+01 + + mean = 2.9372043071164584e+00 + + stddev = 2.4773099230201887e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.1121660976460981e+00 + + min = -3.1121660976460999e+00 + max = 1.1702623812413986e+01 - + mean = 2.9129297054159191e+00 + + mean = 2.9129297054159196e+00 + stddev = 2.5169431213654003e+00 Initial state for member 5: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9612817175857606e+00 + + max = 1.0551713700125084e+01 + + mean = 2.9126327750528640e+00 + + stddev = 2.4337589378004414e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.5956262740496951e+00 - + max = 1.2171329064345073e+01 - + mean = 2.9271112669324744e+00 + + min = -3.5956262740496894e+00 + + max = 1.2171329064345075e+01 + + mean = 2.9271112669324730e+00 + stddev = 2.4804768016832432e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -2.9774585234633828e+00 + + min = -2.9774585234633846e+00 + max = 1.1723207432353661e+01 - + mean = 2.9028366652319368e+00 - + stddev = 2.5131412139072196e+00 + + mean = 2.9028366652319351e+00 + + stddev = 2.5131412139072200e+00 Initial state for member 6: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9269592310066841e+00 + + max = 1.0369847700790473e+01 + + mean = 2.9000458775643962e+00 + + stddev = 2.4290801634301578e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7978844217154695e+00 - + max = 1.1989175148280928e+01 + + min = -3.7978844217154637e+00 + + max = 1.1989175148280930e+01 + mean = 2.9145243694440039e+00 + stddev = 2.4805039059431024e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.3171213032494618e+00 + + min = -3.3171213032494635e+00 + max = 1.1573971625704308e+01 - + mean = 2.8902497677434691e+00 + + mean = 2.8902497677434695e+00 + stddev = 2.5187596854408216e+00 Initial state for member 7: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9082631390290494e+00 + + max = 1.0454101188243557e+01 + + mean = 2.8798595233266093e+00 + + stddev = 2.4299793068145683e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7975749390902465e+00 - + max = 1.2251850212657780e+01 + + min = -3.7975749390902407e+00 + + max = 1.2251850212657782e+01 + mean = 2.8943380152062188e+00 + stddev = 2.4812872150538610e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.0378682439934033e+00 + + min = -3.0378682439934050e+00 + max = 1.1826292199819925e+01 - + mean = 2.8700634135056800e+00 - + stddev = 2.5180079590924711e+00 + + mean = 2.8700634135056804e+00 + + stddev = 2.5180079590924707e+00 Initial state for member 8: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -3.0108771430918098e+00 + + max = 1.0392209617634274e+01 + + mean = 2.9186225207429057e+00 + + stddev = 2.4319235862110085e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7743073814984451e+00 - + max = 1.2143178226825979e+01 + + min = -3.7743073814984394e+00 + + max = 1.2143178226825981e+01 + mean = 2.9331010126225117e+00 - + stddev = 2.4816598891775685e+00 + + stddev = 2.4816598891775690e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.2018401754241959e+00 + + min = -3.2018401754241976e+00 + max = 1.1716628065200444e+01 - + mean = 2.9088264109219706e+00 - + stddev = 2.5152523912857929e+00 + + mean = 2.9088264109219693e+00 + + stddev = 2.5152523912857925e+00 Initial state for member 9: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -3.0167976800036391e+00 + + max = 1.0342283286192707e+01 + + mean = 2.9140939531907932e+00 + + stddev = 2.4275906316159497e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.5964798492363159e+00 - + max = 1.2102612227690715e+01 - + mean = 2.9285724450704063e+00 - + stddev = 2.4741424832790200e+00 + + min = -3.5964798492363101e+00 + + max = 1.2102612227690717e+01 + + mean = 2.9285724450704054e+00 + + stddev = 2.4741424832790195e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.0891739805595986e+00 + + min = -3.0891739805596004e+00 + max = 1.1882456957369239e+01 - + mean = 2.9042978433698612e+00 - + stddev = 2.5188529136466817e+00 + + mean = 2.9042978433698621e+00 + + stddev = 2.5188529136466813e+00 Initial state for member 10: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.8369190735884664e+00 + + max = 1.0516903689881381e+01 + + mean = 2.9112468504917692e+00 + + stddev = 2.4261568477618400e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.8996379141304436e+00 - + max = 1.2062967451595041e+01 - + mean = 2.9257253423713774e+00 + + min = -3.8996379141304405e+00 + + max = 1.2062967451595043e+01 + + mean = 2.9257253423713792e+00 + stddev = 2.4751851755781749e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.1391896996005784e+00 + + min = -3.1391896996005801e+00 + max = 1.1657705040131665e+01 - + mean = 2.9014507406708465e+00 - + stddev = 2.5208945216188483e+00 + + mean = 2.9014507406708456e+00 + + stddev = 2.5208945216188479e+00 H(x) for member 1: -air_temperature nobs= 60 Min=-7.2610476813309921e-01, Max=6.9439433725053235e+00, RMS=2.9491597734195136e+00 +air_temperature nobs= 60 Min=-7.2610476813309721e-01, Max=6.9439433725053235e+00, RMS=2.9491597734195141e+00 H(x) for member 2: -air_temperature nobs= 60 Min=-6.8682889002406133e-01, Max=7.0707356969794048e+00, RMS=2.9863156060296117e+00 +air_temperature nobs= 60 Min=-6.8682889002405934e-01, Max=7.0707356969794048e+00, RMS=2.9863156060296117e+00 H(x) for member 3: -air_temperature nobs= 60 Min=-7.6614106926166547e-01, Max=7.1354212151342793e+00, RMS=2.9847056166657615e+00 +air_temperature nobs= 60 Min=-7.6614106926166348e-01, Max=7.1354212151342793e+00, RMS=2.9847056166657615e+00 H(x) for member 4: -air_temperature nobs= 60 Min=-5.9029356124841703e-01, Max=7.0909195100391074e+00, RMS=2.9961751265341521e+00 +air_temperature nobs= 60 Min=-5.9029356124841503e-01, Max=7.0909195100391074e+00, RMS=2.9961751265341516e+00 H(x) for member 5: -air_temperature nobs= 60 Min=-6.1535834255659116e-01, Max=7.1772045032555027e+00, RMS=3.0153926854057795e+00 +air_temperature nobs= 60 Min=-6.1535834255658917e-01, Max=7.1772045032555027e+00, RMS=3.0153926854057795e+00 H(x) for member 6: -air_temperature nobs= 60 Min=-7.7678864183866514e-01, Max=7.1493712650114460e+00, RMS=2.9928362819356846e+00 +air_temperature nobs= 60 Min=-7.7678864183866314e-01, Max=7.1493712650114460e+00, RMS=2.9928362819356860e+00 H(x) for member 7: -air_temperature nobs= 60 Min=-8.6108647766553892e-01, Max=6.9505081982785626e+00, RMS=2.9549896130109961e+00 +air_temperature nobs= 60 Min=-8.6108647766553692e-01, Max=6.9505081982785626e+00, RMS=2.9549896130109961e+00 H(x) for member 8: -air_temperature nobs= 60 Min=-7.1565823247687055e-01, Max=7.2152245999747677e+00, RMS=2.9903270766404169e+00 +air_temperature nobs= 60 Min=-7.1565823247686855e-01, Max=7.2152245999747677e+00, RMS=2.9903270766404169e+00 H(x) for member 9: -air_temperature nobs= 60 Min=-7.3895062613331808e-01, Max=6.8982934416705532e+00, RMS=2.9806028563771765e+00 +air_temperature nobs= 60 Min=-7.3895062613331608e-01, Max=6.8982934416705532e+00, RMS=2.9806028563771765e+00 H(x) for member 10: -air_temperature nobs= 60 Min=-7.0307266539924873e-01, Max=7.0775172184655233e+00, RMS=2.9671364758528869e+00 +air_temperature nobs= 60 Min=-7.0307266539924673e-01, Max=7.0775172184655233e+00, RMS=2.9671364758528869e+00 H(x) ensemble background mean: -air_temperature nobs= 60 Min=-7.1802832747374767e-01, Max=7.0709139021314478e+00, RMS=2.9806951511382449e+00 +air_temperature nobs= 60 Min=-7.1802832747374545e-01, Max=7.0709139021314478e+00, RMS=2.9806951511382449e+00 background y - H(x): -air_temperature nobs= 60 Min=-2.7821336125458485e-01, Max=2.9035468632568628e-01, RMS=1.2909270475559509e-01 +air_temperature nobs= 60 Min=-3.0100697971122448e-01, Max=4.3984135999299956e-01, RMS=1.6129393728868352e-01 Background mean : +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9073846086571313e+00 + + max = 1.0466714484422283e+01 + + mean = 2.9067740825889588e+00 + + stddev = 2.4291217731685397e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7053062732152808e+00 + + min = -3.7053062732152755e+00 + max = 1.2104987444980736e+01 + mean = 2.9212525744685660e+00 - + stddev = 2.4760285070998300e+00 + + stddev = 2.4760285070998291e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.1658755518783543e+00 + + min = -3.1658755518783561e+00 + max = 1.1751844060042107e+01 - + mean = 2.8969779727680245e+00 - + stddev = 2.5158486977911565e+00 + + mean = 2.8969779727680240e+00 + + stddev = 2.5158486977911561e+00 Analysis mean : +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9073846086571313e+00 + + max = 1.0466714484422283e+01 + + mean = 2.9076689992076608e+00 + + stddev = 2.4290761916962622e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7053062732152808e+00 + + min = -3.7053062732152755e+00 + max = 1.2104987444980736e+01 - + mean = 2.9221231101451690e+00 - + stddev = 2.4761618318061847e+00 + + mean = 2.9221474910872711e+00 + + stddev = 2.4764948117072509e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.1685441889181298e+00 + + min = -3.1686694859290410e+00 + max = 1.1751844060042107e+01 - + mean = 2.8978485084446270e+00 - + stddev = 2.5155427507588950e+00 + + mean = 2.8978728893867300e+00 + + stddev = 2.5153541846932703e+00 +Analysis mean increment : +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -8.4930894365067733e-02 + + max = 1.4829482202157740e-01 + + mean = 8.9491661870372025e-04 + + stddev = 1.2963905605867066e-02 Analysis mean increment : - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.9504226916897700e-02 - + max = 1.0704580639539252e-01 - + mean = 8.7053567660182846e-04 - + stddev = 8.8192004178665497e-03 + + min = -8.4930894365069065e-02 + + max = 1.4829482202157518e-01 + + mean = 8.9491661870371938e-04 + + stddev = 1.2963905605867090e-02 Analysis mean increment : - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.9504226916898588e-02 - + max = 1.0704580639539207e-01 - + mean = 8.7053567660183692e-04 - + stddev = 8.8192004178665340e-03 + + min = -8.4930894365068510e-02 + + max = 1.4829482202157740e-01 + + mean = 8.9491661870372491e-04 + + stddev = 1.2963905605867103e-02 +Forecast variance : +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = 1.2171120495856585e-03 + + max = 2.8858115544075500e-02 + + mean = 9.4854503584923446e-03 + + stddev = 4.6904630703599563e-03 Forecast variance : - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] @@ -251,7 +365,7 @@ Forecast variance : - air_temperature (2 levels): + min = 1.2171120495856600e-03 + max = 2.8858115544075549e-02 - + mean = 9.4854503584923446e-03 + + mean = 9.4854503584923464e-03 + stddev = 4.6904630703599546e-03 Forecast variance : - Valid time: 2010-01-01T06:00:00Z @@ -259,50 +373,59 @@ Forecast variance : Fields: - air_temperature (2 levels): + min = 1.2171120495856624e-03 - + max = 2.8858115544075472e-02 + + max = 2.8858115544075500e-02 + mean = 9.4854503584923446e-03 + stddev = 4.6904630703599563e-03 Analysis variance : +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = 1.2171120495856585e-03 + + max = 2.8858115544075500e-02 + + mean = 9.5013453418191837e-03 + + stddev = 4.7042078663046671e-03 +Analysis variance : - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): + min = 1.2171120495856600e-03 + max = 2.8858115544075549e-02 - + mean = 9.5013453249963038e-03 - + stddev = 4.7042078149935641e-03 + + mean = 9.5013453418191837e-03 + + stddev = 4.7042078663046688e-03 Analysis variance : - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): + min = 1.2171120495856624e-03 - + max = 2.8858115544075472e-02 - + mean = 9.5013453249963055e-03 - + stddev = 4.7042078149935650e-03 + + max = 2.8858115544075500e-02 + + mean = 9.5013453418191837e-03 + + stddev = 4.7042078663046662e-03 H(x) for member 1: -air_temperature nobs= 60 Min=-7.5723120345877337e-01, Max=6.9632137008589163e+00, RMS=2.9569398698299780e+00 +air_temperature nobs= 60 Min=-7.5915611667509786e-01, Max=6.9597031075864990e+00, RMS=2.9576089300074218e+00 H(x) for member 2: -air_temperature nobs= 60 Min=-7.2255607470081973e-01, Max=7.0820201747988287e+00, RMS=2.9938030815438932e+00 +air_temperature nobs= 60 Min=-7.2448098465751520e-01, Max=7.0785096005021089e+00, RMS=2.9946241505817182e+00 H(x) for member 3: -air_temperature nobs= 60 Min=-7.9188344589577842e-01, Max=7.1428012103420153e+00, RMS=2.9922844122795844e+00 +air_temperature nobs= 60 Min=-7.9380836330305449e-01, Max=7.1392906356960495e+00, RMS=2.9930685402167252e+00 H(x) for member 4: -air_temperature nobs= 60 Min=-6.3514482191907096e-01, Max=7.0990774073284522e+00, RMS=3.0029534613315598e+00 +air_temperature nobs= 60 Min=-6.3706973560105673e-01, Max=7.0955668103307445e+00, RMS=3.0040227709466265e+00 H(x) for member 5: -air_temperature nobs= 60 Min=-6.5777879170285392e-01, Max=7.1832733192127600e+00, RMS=3.0209730891180322e+00 +air_temperature nobs= 60 Min=-6.5970370538483969e-01, Max=7.1797627333909233e+00, RMS=3.0218188088708509e+00 H(x) for member 6: -air_temperature nobs= 60 Min=-8.0290583483086753e-01, Max=7.1566844843905821e+00, RMS=2.9998461395174276e+00 +air_temperature nobs= 60 Min=-8.0483074851285330e-01, Max=7.1531738985687454e+00, RMS=3.0006193457062751e+00 H(x) for member 7: -air_temperature nobs= 60 Min=-8.7859109751092124e-01, Max=6.9701160081308737e+00, RMS=2.9639572451847611e+00 +air_temperature nobs= 60 Min=-8.8051602609406820e-01, Max=6.9666054223090370e+00, RMS=2.9648193636886280e+00 H(x) for member 8: -air_temperature nobs= 60 Min=-7.4735064254285621e-01, Max=7.2167436890881911e+00, RMS=2.9971953239377438e+00 +air_temperature nobs= 60 Min=-7.4927555645767263e-01, Max=7.2132331032663544e+00, RMS=2.9980787828684474e+00 H(x) for member 9: -air_temperature nobs= 60 Min=-7.6866369656430411e-01, Max=6.9206616692822829e+00, RMS=2.9880902709663997e+00 +air_temperature nobs= 60 Min=-7.7058861024628988e-01, Max=6.9171510685592850e+00, RMS=2.9890677884855115e+00 H(x) for member 10: -air_temperature nobs= 60 Min=-7.3686359319047856e-01, Max=7.0884226868506923e+00, RMS=2.9754692807636975e+00 +air_temperature nobs= 60 Min=-7.3878850780378691e-01, Max=7.0849121010288556e+00, RMS=2.9763713219297534e+00 H(x) ensemble analysis mean: -air_temperature nobs= 60 Min=-7.4989692023167243e-01, Max=7.0823014350283602e+00, RMS=2.9881725838975504e+00 +air_temperature nobs= 60 Min=-7.5182183547362358e-01, Max=7.0787908481238606e+00, RMS=2.9890318530437305e+00 analysis y - H(x): -air_temperature nobs= 60 Min=-2.4634476849666009e-01, Max=2.4145483611741847e-01, RMS=1.1138454603568704e-01 -ombg RMS: 1.2909270475559509e-01 -oman RMS: 1.1138454603568704e-01 +air_temperature nobs= 60 Min=-2.5219727273170456e-01, Max=3.7279251710505634e-01, RMS=1.3189979301567081e-01 +ombg RMS: 1.6129393728868352e-01 +oman RMS: 1.3189979301567081e-01 diff --git a/test/glb/testref/letkf_nonlinear_4d.ref b/test/glb/testref/letkf_nonlinear_4d.ref index 41cbd17..22c1d5d 100644 --- a/test/glb/testref/letkf_nonlinear_4d.ref +++ b/test/glb/testref/letkf_nonlinear_4d.ref @@ -1,249 +1,363 @@ Initial state for member 1: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9070402095754342e+00 + + max = 1.0443230228491029e+01 + + mean = 2.8890224246311544e+00 + + stddev = 2.4337264386471138e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7984786781412843e+00 - + max = 1.2105572354693958e+01 + + min = -3.7984786781412785e+00 + + max = 1.2105572354693960e+01 + mean = 2.9035009165107679e+00 - + stddev = 2.4802542460785419e+00 + + stddev = 2.4802542460785428e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.2324168099554287e+00 + + min = -3.2324168099554305e+00 + max = 1.1684028286658434e+01 + mean = 2.8792263148102166e+00 + stddev = 2.5189832758053812e+00 Initial state for member 2: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9169112429360231e+00 + + max = 1.0640558928014917e+01 + + mean = 2.9050687446303867e+00 + + stddev = 2.4304624544932087e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.4686154982654220e+00 - + max = 1.1912098410366392e+01 + + min = -3.4686154982654163e+00 + + max = 1.1912098410366394e+01 + mean = 2.9195472365100055e+00 + stddev = 2.4687036163275997e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.2664917313325974e+00 + + min = -3.2664917313325992e+00 + max = 1.1761288355765901e+01 + mean = 2.8952726348094577e+00 + stddev = 2.5149762090011363e+00 Initial state for member 3: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.8439500143940992e+00 + + max = 1.0440496906472719e+01 + + mean = 2.9144223410218180e+00 + + stddev = 2.4306939104918670e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.8669960664368079e+00 - + max = 1.2115786341942979e+01 - + mean = 2.9289008329014217e+00 - + stddev = 2.4776428547666072e+00 + + min = -3.8669960664368022e+00 + + max = 1.2115786341942981e+01 + + mean = 2.9289008329014226e+00 + + stddev = 2.4776428547666063e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.2850289535587947e+00 + + min = -3.2850289535587964e+00 + max = 1.1990238825003495e+01 - + mean = 2.9046262312008806e+00 - + stddev = 2.5193094400446165e+00 + + mean = 2.9046262312008815e+00 + + stddev = 2.5193094400446161e+00 Initial state for member 4: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9719403295366162e+00 + + max = 1.0515799598376681e+01 + + mean = 2.9227258152368534e+00 + + stddev = 2.4350683503339825e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7182596232436875e+00 - + max = 1.2195305011408491e+01 - + mean = 2.9372043071164575e+00 - + stddev = 2.4773099230201883e+00 + + min = -3.7182596232436844e+00 + + max = 1.2195305011408493e+01 + + mean = 2.9372043071164584e+00 + + stddev = 2.4773099230201887e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.1121660976460981e+00 + + min = -3.1121660976460999e+00 + max = 1.1702623812413986e+01 - + mean = 2.9129297054159191e+00 + + mean = 2.9129297054159196e+00 + stddev = 2.5169431213654003e+00 Initial state for member 5: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9612817175857606e+00 + + max = 1.0551713700125084e+01 + + mean = 2.9126327750528640e+00 + + stddev = 2.4337589378004414e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.5956262740496951e+00 - + max = 1.2171329064345073e+01 - + mean = 2.9271112669324744e+00 + + min = -3.5956262740496894e+00 + + max = 1.2171329064345075e+01 + + mean = 2.9271112669324730e+00 + stddev = 2.4804768016832432e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -2.9774585234633828e+00 + + min = -2.9774585234633846e+00 + max = 1.1723207432353661e+01 - + mean = 2.9028366652319368e+00 - + stddev = 2.5131412139072196e+00 + + mean = 2.9028366652319351e+00 + + stddev = 2.5131412139072200e+00 Initial state for member 6: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9269592310066841e+00 + + max = 1.0369847700790473e+01 + + mean = 2.9000458775643962e+00 + + stddev = 2.4290801634301578e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7978844217154695e+00 - + max = 1.1989175148280928e+01 + + min = -3.7978844217154637e+00 + + max = 1.1989175148280930e+01 + mean = 2.9145243694440039e+00 + stddev = 2.4805039059431024e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.3171213032494618e+00 + + min = -3.3171213032494635e+00 + max = 1.1573971625704308e+01 - + mean = 2.8902497677434691e+00 + + mean = 2.8902497677434695e+00 + stddev = 2.5187596854408216e+00 Initial state for member 7: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9082631390290494e+00 + + max = 1.0454101188243557e+01 + + mean = 2.8798595233266093e+00 + + stddev = 2.4299793068145683e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7975749390902465e+00 - + max = 1.2251850212657780e+01 + + min = -3.7975749390902407e+00 + + max = 1.2251850212657782e+01 + mean = 2.8943380152062188e+00 + stddev = 2.4812872150538610e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.0378682439934033e+00 + + min = -3.0378682439934050e+00 + max = 1.1826292199819925e+01 - + mean = 2.8700634135056800e+00 - + stddev = 2.5180079590924711e+00 + + mean = 2.8700634135056804e+00 + + stddev = 2.5180079590924707e+00 Initial state for member 8: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -3.0108771430918098e+00 + + max = 1.0392209617634274e+01 + + mean = 2.9186225207429057e+00 + + stddev = 2.4319235862110085e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7743073814984451e+00 - + max = 1.2143178226825979e+01 + + min = -3.7743073814984394e+00 + + max = 1.2143178226825981e+01 + mean = 2.9331010126225117e+00 - + stddev = 2.4816598891775685e+00 + + stddev = 2.4816598891775690e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.2018401754241959e+00 + + min = -3.2018401754241976e+00 + max = 1.1716628065200444e+01 - + mean = 2.9088264109219706e+00 - + stddev = 2.5152523912857929e+00 + + mean = 2.9088264109219693e+00 + + stddev = 2.5152523912857925e+00 Initial state for member 9: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -3.0167976800036391e+00 + + max = 1.0342283286192707e+01 + + mean = 2.9140939531907932e+00 + + stddev = 2.4275906316159497e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.5964798492363159e+00 - + max = 1.2102612227690715e+01 - + mean = 2.9285724450704063e+00 - + stddev = 2.4741424832790200e+00 + + min = -3.5964798492363101e+00 + + max = 1.2102612227690717e+01 + + mean = 2.9285724450704054e+00 + + stddev = 2.4741424832790195e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.0891739805595986e+00 + + min = -3.0891739805596004e+00 + max = 1.1882456957369239e+01 - + mean = 2.9042978433698612e+00 - + stddev = 2.5188529136466817e+00 + + mean = 2.9042978433698621e+00 + + stddev = 2.5188529136466813e+00 Initial state for member 10: +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.8369190735884664e+00 + + max = 1.0516903689881381e+01 + + mean = 2.9112468504917692e+00 + + stddev = 2.4261568477618400e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.8996379141304436e+00 - + max = 1.2062967451595041e+01 - + mean = 2.9257253423713774e+00 + + min = -3.8996379141304405e+00 + + max = 1.2062967451595043e+01 + + mean = 2.9257253423713792e+00 + stddev = 2.4751851755781749e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.1391896996005784e+00 + + min = -3.1391896996005801e+00 + max = 1.1657705040131665e+01 - + mean = 2.9014507406708465e+00 - + stddev = 2.5208945216188483e+00 + + mean = 2.9014507406708456e+00 + + stddev = 2.5208945216188479e+00 H(x) for member 1: -air_temperature nobs= 60 Min=-7.2610476795019063e-01, Max=6.9439433722032051e+00, RMS=2.9491597730964592e+00 +air_temperature nobs= 60 Min=-7.2610476795018863e-01, Max=6.9439433722032042e+00, RMS=2.9491597730964600e+00 H(x) for member 2: -air_temperature nobs= 60 Min=-6.8682888958949873e-01, Max=7.0707356969791579e+00, RMS=2.9863156062997747e+00 +air_temperature nobs= 60 Min=-6.8682888958949673e-01, Max=7.0707356969791579e+00, RMS=2.9863156062997747e+00 H(x) for member 3: -air_temperature nobs= 60 Min=-7.6614106824035944e-01, Max=7.1354212132202015e+00, RMS=2.9847056168956945e+00 +air_temperature nobs= 60 Min=-7.6614106824035755e-01, Max=7.1354212132202015e+00, RMS=2.9847056168956945e+00 H(x) for member 4: -air_temperature nobs= 60 Min=-5.9029356345726924e-01, Max=7.0909195092985735e+00, RMS=2.9961751264560275e+00 +air_temperature nobs= 60 Min=-5.9029356345726725e-01, Max=7.0909195092985726e+00, RMS=2.9961751264560275e+00 H(x) for member 5: -air_temperature nobs= 60 Min=-6.1535833923841787e-01, Max=7.1772045067809938e+00, RMS=3.0153926850983885e+00 +air_temperature nobs= 60 Min=-6.1535833923841599e-01, Max=7.1772045067809938e+00, RMS=3.0153926850983885e+00 H(x) for member 6: -air_temperature nobs= 60 Min=-7.7678864189652597e-01, Max=7.1493712632461115e+00, RMS=2.9928362822642169e+00 +air_temperature nobs= 60 Min=-7.7678864189652397e-01, Max=7.1493712632461115e+00, RMS=2.9928362822642174e+00 H(x) for member 7: -air_temperature nobs= 60 Min=-8.6108647142780237e-01, Max=6.9505081985040729e+00, RMS=2.9549896128035780e+00 +air_temperature nobs= 60 Min=-8.6108647142780059e-01, Max=6.9505081985040720e+00, RMS=2.9549896128035780e+00 H(x) for member 8: -air_temperature nobs= 60 Min=-7.1565823258985661e-01, Max=7.2152246002289973e+00, RMS=2.9903270767258325e+00 +air_temperature nobs= 60 Min=-7.1565823258985484e-01, Max=7.2152246002289973e+00, RMS=2.9903270767258325e+00 H(x) for member 9: -air_temperature nobs= 60 Min=-7.3895062696432534e-01, Max=6.8982934418590860e+00, RMS=2.9806028562058509e+00 +air_temperature nobs= 60 Min=-7.3895062696432334e-01, Max=6.8982934418590860e+00, RMS=2.9806028562058513e+00 H(x) for member 10: -air_temperature nobs= 60 Min=-7.0307266523415513e-01, Max=7.0775172182955686e+00, RMS=2.9671364757248346e+00 +air_temperature nobs= 60 Min=-7.0307266523415324e-01, Max=7.0775172182955686e+00, RMS=2.9671364757248346e+00 H(x) ensemble background mean: -air_temperature nobs= 60 Min=-7.1802832665884020e-01, Max=7.0709139020615970e+00, RMS=2.9806951510999538e+00 +air_temperature nobs= 60 Min=-7.1802832665883809e-01, Max=7.0709139020615970e+00, RMS=2.9806951510999542e+00 background y - H(x): -air_temperature nobs= 60 Min=-2.7821336206949232e-01, Max=2.9035468786236862e-01, RMS=1.2909270488859084e-01 +air_temperature nobs= 60 Min=-3.0100697969376222e-01, Max=4.3984136152968212e-01, RMS=1.6129393736688719e-01 Background mean : +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9073846086571313e+00 + + max = 1.0466714484422283e+01 + + mean = 2.9067740825889588e+00 + + stddev = 2.4291217731685397e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7053062732152808e+00 + + min = -3.7053062732152755e+00 + max = 1.2104987444980736e+01 + mean = 2.9212525744685660e+00 - + stddev = 2.4760285070998300e+00 + + stddev = 2.4760285070998291e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.1658755518783543e+00 + + min = -3.1658755518783561e+00 + max = 1.1751844060042107e+01 - + mean = 2.8969779727680245e+00 - + stddev = 2.5158486977911565e+00 + + mean = 2.8969779727680240e+00 + + stddev = 2.5158486977911561e+00 Analysis mean : +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -2.9073846086571313e+00 + + max = 1.0466714484422283e+01 + + mean = 2.9076689992449807e+00 + + stddev = 2.4290761917065598e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.7053062732152808e+00 + + min = -3.7053062732152755e+00 + max = 1.2104987444980736e+01 - + mean = 2.9221231101474494e+00 - + stddev = 2.4761618318054239e+00 + + mean = 2.9221474911245902e+00 + + stddev = 2.4764948117212917e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.1685441889315054e+00 + + min = -3.1686694859290410e+00 + max = 1.1751844060042107e+01 - + mean = 2.8978485084469061e+00 - + stddev = 2.5155427507595931e+00 + + mean = 2.8978728894240491e+00 + + stddev = 2.5153541846877907e+00 +Analysis mean increment : +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = -8.4930894365067733e-02 + + max = 1.4829482202157740e-01 + + mean = 8.9491665602318450e-04 + + stddev = 1.2963905717704022e-02 Analysis mean increment : - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.9504226991688540e-02 - + max = 1.0704580721221202e-01 - + mean = 8.7053567888188167e-04 - + stddev = 8.8192004424942263e-03 + + min = -8.4930894365069065e-02 + + max = 1.4829482202157518e-01 + + mean = 8.9491665602318504e-04 + + stddev = 1.2963905717704054e-02 Analysis mean increment : - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.9504226991689650e-02 - + max = 1.0704580721221135e-01 - + mean = 8.7053567888189447e-04 - + stddev = 8.8192004424942367e-03 + + min = -8.4930894365068510e-02 + + max = 1.4829482202157740e-01 + + mean = 8.9491665602318537e-04 + + stddev = 1.2963905717704052e-02 +Forecast variance : +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = 1.2171120495856585e-03 + + max = 2.8858115544075500e-02 + + mean = 9.4854503584923446e-03 + + stddev = 4.6904630703599563e-03 Forecast variance : - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] @@ -251,7 +365,7 @@ Forecast variance : - air_temperature (2 levels): + min = 1.2171120495856600e-03 + max = 2.8858115544075549e-02 - + mean = 9.4854503584923446e-03 + + mean = 9.4854503584923464e-03 + stddev = 4.6904630703599546e-03 Forecast variance : - Valid time: 2010-01-01T06:00:00Z @@ -259,50 +373,59 @@ Forecast variance : Fields: - air_temperature (2 levels): + min = 1.2171120495856624e-03 - + max = 2.8858115544075472e-02 + + max = 2.8858115544075500e-02 + mean = 9.4854503584923446e-03 + stddev = 4.6904630703599563e-03 Analysis variance : +- Valid time: 2010-01-01T00:00:00Z + Geometry: L40x21 [840] + Fields: + - air_temperature (2 levels): + + min = 1.2171120495856585e-03 + + max = 2.8858115544075500e-02 + + mean = 9.5013453418191837e-03 + + stddev = 4.7042078663046671e-03 +Analysis variance : - Valid time: 2010-01-01T03:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): + min = 1.2171120495856600e-03 + max = 2.8858115544075549e-02 - + mean = 9.5013453249963038e-03 - + stddev = 4.7042078149935641e-03 + + mean = 9.5013453418191837e-03 + + stddev = 4.7042078663046688e-03 Analysis variance : - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): + min = 1.2171120495856624e-03 - + max = 2.8858115544075472e-02 - + mean = 9.5013453249963055e-03 - + stddev = 4.7042078149935659e-03 + + max = 2.8858115544075500e-02 + + mean = 9.5013453418191837e-03 + + stddev = 4.7042078663046662e-03 H(x) for member 1: -air_temperature nobs= 60 Min=-7.5723120387519338e-01, Max=6.9632137026234044e+00, RMS=2.9569398692750815e+00 +air_temperature nobs= 60 Min=-7.5915611668703398e-01, Max=6.9597031093445718e+00, RMS=2.9576089305744047e+00 H(x) for member 2: -air_temperature nobs= 60 Min=-7.2255607465702409e-01, Max=7.0820201748351934e+00, RMS=2.9938030816635370e+00 +air_temperature nobs= 60 Min=-7.2448098545098305e-01, Max=7.0785096005021728e+00, RMS=2.9946241505496163e+00 H(x) for member 3: -air_temperature nobs= 60 Min=-7.9188344678414802e-01, Max=7.1428012115166366e+00, RMS=2.9922844123608421e+00 +air_temperature nobs= 60 Min=-7.9380836339824401e-01, Max=7.1392906354828689e+00, RMS=2.9930685408711861e+00 H(x) for member 4: -air_temperature nobs= 60 Min=-6.3514482432621477e-01, Max=7.0990774070306504e+00, RMS=3.0029534613122757e+00 +air_temperature nobs= 60 Min=-6.3706973487661622e-01, Max=7.0955668108904710e+00, RMS=3.0040227717245376e+00 H(x) for member 5: -air_temperature nobs= 60 Min=-6.5777879129903249e-01, Max=7.1832733171709666e+00, RMS=3.0209730894996256e+00 +air_temperature nobs= 60 Min=-6.5970370215326701e-01, Max=7.1797627326876352e+00, RMS=3.0218188089032609e+00 H(x) for member 6: -air_temperature nobs= 60 Min=-8.0290583439235552e-01, Max=7.1566844852751501e+00, RMS=2.9998461399966514e+00 +air_temperature nobs= 60 Min=-8.0483074813441990e-01, Max=7.1531739001390475e+00, RMS=3.0006193459644894e+00 H(x) for member 7: -air_temperature nobs= 60 Min=-8.7859110405364904e-01, Max=6.9701160045639821e+00, RMS=2.9639572451950476e+00 +air_temperature nobs= 60 Min=-8.8051602060308520e-01, Max=6.9666054195817804e+00, RMS=2.9648193643223668e+00 H(x) for member 8: -air_temperature nobs= 60 Min=-7.4735064290465070e-01, Max=7.2167436896909889e+00, RMS=2.9971953241326506e+00 +air_temperature nobs= 60 Min=-7.4927555640853982e-01, Max=7.2132331001114469e+00, RMS=2.9980787832891505e+00 H(x) for member 9: -air_temperature nobs= 60 Min=-7.6866369751858588e-01, Max=6.9206616626202262e+00, RMS=2.9880902707186485e+00 +air_temperature nobs= 60 Min=-7.7058861007825308e-01, Max=6.9171510739617394e+00, RMS=2.9890677894506021e+00 H(x) for member 10: -air_temperature nobs= 60 Min=-7.3686359322522954e-01, Max=7.0884226868516098e+00, RMS=2.9754692807702341e+00 +air_temperature nobs= 60 Min=-7.3878850739973601e-01, Max=7.0849121009524882e+00, RMS=2.9763713225125139e+00 H(x) ensemble analysis mean: -air_temperature nobs= 60 Min=-7.4989692130360841e-01, Max=7.0823014342178814e+00, RMS=2.9881725839404680e+00 +air_temperature nobs= 60 Min=-7.5182183451901796e-01, Max=7.0787908483654221e+00, RMS=2.9890318535353999e+00 analysis y - H(x): -air_temperature nobs= 60 Min=-2.4634476742472411e-01, Max=2.4145483602368234e-01, RMS=1.1138454604690103e-01 -ombg RMS: 1.2909270488859084e-01 -oman RMS: 1.1138454604690103e-01 +air_temperature nobs= 60 Min=-2.5219727235917566e-01, Max=3.7279251480228481e-01, RMS=1.3189979291417500e-01 +ombg RMS: 1.6129393736688719e-01 +oman RMS: 1.3189979291417500e-01 diff --git a/test/glb/testref/makeobs_4d.ref b/test/glb/testref/makeobs_4d.ref index 1981c41..5693c7d 100644 --- a/test/glb/testref/makeobs_4d.ref +++ b/test/glb/testref/makeobs_4d.ref @@ -3,22 +3,22 @@ Initial state: Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.0181907010909312e+00 - + max = 1.0487693671589927e+01 - + mean = 2.9143919578582302e+00 - + stddev = 2.4294364876410390e+00 + + min = -3.0181907010909321e+00 + + max = 1.0487693671589923e+01 + + mean = 2.9143919578582307e+00 + + stddev = 2.4294364876410386e+00 Final state: - Valid time: 2010-01-01T06:00:00Z Geometry: L40x21 [840] Fields: - air_temperature (2 levels): - + min = -3.1224054601450542e+00 + + min = -3.1224054601450559e+00 + max = 1.1864022317124880e+01 - + mean = 2.9045958480372880e+00 + + mean = 2.9045958480372889e+00 + stddev = 2.5173193882399643e+00 H(x): -air_temperature nobs= 60 Min=-6.7905183539386926e-01, Max=7.1542350417439771e+00, RMS=3.0132229461012052e+00 +air_temperature nobs= 60 Min=-6.7905183539386749e-01, Max=7.1542350417439762e+00, RMS=3.0132229461012052e+00 End H(x) Perturbed H(x): -air_temperature nobs= 60 Min=-9.9624168599816576e-01, Max=7.1593288511702990e+00, RMS=3.0166879919848535e+00 +air_temperature nobs= 60 Min=-5.5786034384246941e-01, Max=7.0344469027142074e+00, RMS=3.0012214106026036e+00 End Perturbed H(x) diff --git a/test/reg/testinput/4densvar.yaml b/test/reg/testinput/4densvar.yaml index cbd7acd..f6efed8 100644 --- a/test/reg/testinput/4densvar.yaml +++ b/test/reg/testinput/4densvar.yaml @@ -2,7 +2,7 @@ cost function: cost type: 4D-Ens-Var time window: begin: 2010-01-01T00:00:00Z - length: PT06H + length: PT6H subwindow: PT3H parallel subwindows: false analysis variables: [air_temperature] diff --git a/test/reg/testinput/4densvar_parallel.yaml b/test/reg/testinput/4densvar_parallel.yaml index e00441e..0490dac 100644 --- a/test/reg/testinput/4densvar_parallel.yaml +++ b/test/reg/testinput/4densvar_parallel.yaml @@ -2,7 +2,7 @@ cost function: cost type: 4D-Ens-Var time window: begin: 2010-01-01T00:00:00Z - length: PT06H + length: PT6H subwindow: PT3H analysis variables: [air_temperature] geometry: &geom diff --git a/test/reg/testinput/4dvar.yaml b/test/reg/testinput/4dvar.yaml index 324d9cf..341e863 100644 --- a/test/reg/testinput/4dvar.yaml +++ b/test/reg/testinput/4dvar.yaml @@ -99,6 +99,6 @@ final: output: filepath: testdata/4dvar_analysis__MPI___DATE_ date pattern: _DATE_ - frequency: PT06H + frequency: PT6H test: reference filename: testref/4dvar.ref diff --git a/test/reg/testinput/letkf_linear_4d.yaml b/test/reg/testinput/letkf_linear_4d.yaml index e4806f5..1502733 100644 --- a/test/reg/testinput/letkf_linear_4d.yaml +++ b/test/reg/testinput/letkf_linear_4d.yaml @@ -26,6 +26,9 @@ background: members from template: template: states: + - date: 2010-01-01T00:00:00Z + state variables: [air_temperature] + filepath: testdata/ens__MPI__PT0S_%mem% - date: 2010-01-01T03:00:00Z state variables: [air_temperature] filepath: testdata/ens__MPI__PT3H_%mem% @@ -79,6 +82,8 @@ local ensemble DA: use linear observer: true output: states: + - date: 2010-01-01T00:00:00Z + filepath: testdata/letkf_linear_4d_state__MPI__PT0S - date: 2010-01-01T03:00:00Z filepath: testdata/letkf_linear_4d_state__MPI__PT3H - date: 2010-01-01T06:00:00Z @@ -89,6 +94,8 @@ output ensemble increments: filepath: testdata/letkf_linear_4d_increment__MPI_ output mean prior: states: + - date: 2010-01-01T00:00:00Z + filepath: testdata/letkf_linear_4d_mean_prior__MPI__PT0S - date: 2010-01-01T03:00:00Z filepath: testdata/letkf_linear_4d_mean_prior__MPI__PT3H - date: 2010-01-01T06:00:00Z diff --git a/test/reg/testinput/letkf_nonlinear_4d.yaml b/test/reg/testinput/letkf_nonlinear_4d.yaml index 19a03c4..b17d062 100644 --- a/test/reg/testinput/letkf_nonlinear_4d.yaml +++ b/test/reg/testinput/letkf_nonlinear_4d.yaml @@ -26,6 +26,9 @@ background: members from template: template: states: + - date: 2010-01-01T00:00:00Z + state variables: [air_temperature] + filepath: testdata/ens__MPI__PT0S_%mem% - date: 2010-01-01T03:00:00Z state variables: [air_temperature] filepath: testdata/ens__MPI__PT3H_%mem% @@ -78,6 +81,8 @@ local ensemble DA: mult: 1.1 output: states: + - date: 2010-01-01T00:00:00Z + filepath: testdata/letkf_nonlinear_4d_state__MPI__PT0S - date: 2010-01-01T03:00:00Z filepath: testdata/letkf_nonlinear_4d_state__MPI__PT3H - date: 2010-01-01T06:00:00Z @@ -88,6 +93,8 @@ output ensemble increments: filepath: testdata/letkf_nonlinear_4d_increment__MPI_ output mean prior: states: + - date: 2010-01-01T00:00:00Z + filepath: testdata/letkf_nonlinear_4d_mean_prior__MPI__PT0S - date: 2010-01-01T03:00:00Z filepath: testdata/letkf_nonlinear_4d_mean_prior__MPI__PT3H - date: 2010-01-01T06:00:00Z diff --git a/test/reg/testinput/makeobs_4d.yaml b/test/reg/testinput/makeobs_4d.yaml index 4346d63..909544e 100644 --- a/test/reg/testinput/makeobs_4d.yaml +++ b/test/reg/testinput/makeobs_4d.yaml @@ -41,13 +41,17 @@ observations: engine: type: GenList lats: [55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3, + 55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3, 55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3,55.3,55.8,56.3,56.8,57.3] lons: [8.9,8.9,8.9,8.9,8.9,9.4,9.4,9.4,9.4,9.4,9.9,9.9,9.9,9.9,9.9,10.4,10.4,10.4,10.4,10.4,10.9,10.9,10.9,10.9,10.9, + 8.9,8.9,8.9,8.9,8.9,9.4,9.4,9.4,9.4,9.4,9.9,9.9,9.9,9.9,9.9,10.4,10.4,10.4,10.4,10.4,10.9,10.9,10.9,10.9,10.9, 8.9,8.9,8.9,8.9,8.9,9.4,9.4,9.4,9.4,9.4,9.9,9.9,9.9,9.9,9.9,10.4,10.4,10.4,10.4,10.4,10.9,10.9,10.9,10.9,10.9] - dateTimes: [10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800, + dateTimes: [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, + 10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800,10800, 21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600,21600] vert coord type: height vert coords: [1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3, + 1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3, 1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3] epoch: seconds since 2010-01-01T00:00:00Z obs errors: [0.1] diff --git a/test/reg/testref/4densvar.ref b/test/reg/testref/4densvar.ref index 914a3b6..5f7b828 100644 --- a/test/reg/testref/4densvar.ref +++ b/test/reg/testref/4densvar.ref @@ -1,32 +1,32 @@ -CostJo : Nonlinear Jo(air_temperature) = 3.5705638983357709e+01, nobs = 50, Jo/n = 7.1411277966715414e-01, err = 1.0000000149011612e-01 +CostJo : Nonlinear Jo(air_temperature) = 4.7665289478449964e+01, nobs = 50, Jo/n = 9.5330578956899925e-01, err = 1.0000000149011612e-01 CostJb : Nonlinear Jb = 0.0000000000000000e+00 -CostFunction: Nonlinear J = 3.5705638983357709e+01 -PCGMinimizer: reduction in residual norm = 2.3497013115036190e-02 +CostFunction: Nonlinear J = 4.7665289478449964e+01 +PCGMinimizer: reduction in residual norm = 2.2835528481999670e-02 CostFunction::addIncrement: Analysis: - Valid time: 2010-01-01T00:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0280247877836057e+01 - + max = 1.9629953449148527e+01 - + mean = 2.4609351735314511e+00 - + stddev = 3.1926442078850124e+00 + + min = -1.0280343575761970e+01 + + max = 1.9630068295996622e+01 + + mean = 2.4616193984686672e+00 + + stddev = 3.1920354877830919e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0280247877836057e+01 - + max = 1.8692589583074255e+01 - + mean = 2.1962092817309111e+00 - + stddev = 2.8452107671595677e+00 + + min = -1.0280343575761970e+01 + + max = 1.8692704429922344e+01 + + mean = 2.1968935066681241e+00 + + stddev = 2.8453126294711706e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0280247877836057e+01 - + max = 1.7495902817608833e+01 - + mean = 2.1253872227991235e+00 - + stddev = 2.5150768434627597e+00 -CostJo : Nonlinear Jo(air_temperature) = 1.8098463188595737e+01, nobs = 50, Jo/n = 3.6196926377191474e-01, err = 1.0000000149011612e-01 -CostJb : Nonlinear Jb = 5.3887050259078837e+00 -CostFunction: Nonlinear J = 2.3487168214503619e+01 + + min = -1.0280343575761970e+01 + + max = 1.7496017664456922e+01 + + mean = 2.1260714477363338e+00 + + stddev = 2.5152244495880551e+00 +CostJo : Nonlinear Jo(air_temperature) = 1.9113410491362416e+01, nobs = 50, Jo/n = 3.8226820982724830e-01, err = 1.0000000149011612e-01 +CostJb : Nonlinear Jb = 8.6280766192006126e+00 +CostFunction: Nonlinear J = 2.7741487110563028e+01 diff --git a/test/reg/testref/4densvar_parallel.ref b/test/reg/testref/4densvar_parallel.ref index 7b3af98..5f7b828 100644 --- a/test/reg/testref/4densvar_parallel.ref +++ b/test/reg/testref/4densvar_parallel.ref @@ -1,32 +1,32 @@ -CostJo : Nonlinear Jo(air_temperature) = 3.5705638983357701e+01, nobs = 50, Jo/n = 7.1411277966715403e-01, err = 1.0000000149011612e-01 +CostJo : Nonlinear Jo(air_temperature) = 4.7665289478449964e+01, nobs = 50, Jo/n = 9.5330578956899925e-01, err = 1.0000000149011612e-01 CostJb : Nonlinear Jb = 0.0000000000000000e+00 -CostFunction: Nonlinear J = 3.5705638983357701e+01 -PCGMinimizer: reduction in residual norm = 2.3497013115036190e-02 +CostFunction: Nonlinear J = 4.7665289478449964e+01 +PCGMinimizer: reduction in residual norm = 2.2835528481999670e-02 CostFunction::addIncrement: Analysis: - Valid time: 2010-01-01T00:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0280247877836057e+01 - + max = 1.9629953449148527e+01 - + mean = 2.4609351735314511e+00 - + stddev = 3.1926442078850124e+00 + + min = -1.0280343575761970e+01 + + max = 1.9630068295996622e+01 + + mean = 2.4616193984686672e+00 + + stddev = 3.1920354877830919e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0280247877836057e+01 - + max = 1.8692589583074255e+01 - + mean = 2.1962092817309111e+00 - + stddev = 2.8452107671595677e+00 + + min = -1.0280343575761970e+01 + + max = 1.8692704429922344e+01 + + mean = 2.1968935066681241e+00 + + stddev = 2.8453126294711706e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0280247877836057e+01 - + max = 1.7495902817608833e+01 - + mean = 2.1253872227991235e+00 - + stddev = 2.5150768434627597e+00 -CostJo : Nonlinear Jo(air_temperature) = 1.8098463188595737e+01, nobs = 50, Jo/n = 3.6196926377191474e-01, err = 1.0000000149011612e-01 -CostJb : Nonlinear Jb = 5.3887050259078837e+00 -CostFunction: Nonlinear J = 2.3487168214503619e+01 + + min = -1.0280343575761970e+01 + + max = 1.7496017664456922e+01 + + mean = 2.1260714477363338e+00 + + stddev = 2.5152244495880551e+00 +CostJo : Nonlinear Jo(air_temperature) = 1.9113410491362416e+01, nobs = 50, Jo/n = 3.8226820982724830e-01, err = 1.0000000149011612e-01 +CostJb : Nonlinear Jb = 8.6280766192006126e+00 +CostFunction: Nonlinear J = 2.7741487110563028e+01 diff --git a/test/reg/testref/4dvar.ref b/test/reg/testref/4dvar.ref index cf13f2a..9bef459 100644 --- a/test/reg/testref/4dvar.ref +++ b/test/reg/testref/4dvar.ref @@ -1,17 +1,17 @@ -CostJo : Nonlinear Jo(air_temperature) = 4.3500488593420243e+01, nobs = 50, Jo/n = 8.7000977186840489e-01, err = 1.0000000149011612e-01 +CostJo : Nonlinear Jo(air_temperature) = 5.1692441986392836e+01, nobs = 50, Jo/n = 1.0338488397278567e+00, err = 1.0000000149011612e-01 Norm of input parameter StdDev: 6.0909769331365213e+00 CostJb : Nonlinear Jb = 0.0000000000000000e+00 -CostFunction: Nonlinear J = 4.3500488593420243e+01 -PCGMinimizer: reduction in residual norm = 1.1431739638690886e-02 +CostFunction: Nonlinear J = 5.1692441986392836e+01 +PCGMinimizer: reduction in residual norm = 1.0521704459571303e-02 CostFunction::addIncrement: Analysis: - Valid time: 2010-01-01T00:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0282252083547204e+01 - + max = 1.9585347212804773e+01 - + mean = 2.4604025519081754e+00 - + stddev = 3.1925576048880280e+00 -CostJo : Nonlinear Jo(air_temperature) = 2.0830904522787939e+01, nobs = 50, Jo/n = 4.1661809045575882e-01, err = 1.0000000149011612e-01 -CostJb : Nonlinear Jb = 7.1849842005136626e+00 -CostFunction: Nonlinear J = 2.8015888723301600e+01 + + min = -1.0280933961157031e+01 + + max = 1.9605094074831282e+01 + + mean = 2.4642846081132270e+00 + + stddev = 3.1926502639057088e+00 +CostJo : Nonlinear Jo(air_temperature) = 2.2358465087978473e+01, nobs = 50, Jo/n = 4.4716930175956948e-01, err = 1.0000000149011612e-01 +CostJb : Nonlinear Jb = 9.2030594168369557e+00 +CostFunction: Nonlinear J = 3.1561524504815431e+01 diff --git a/test/reg/testref/getkf_nonlinear.ref b/test/reg/testref/getkf_nonlinear.ref index 312a665..6437efa 100644 --- a/test/reg/testref/getkf_nonlinear.ref +++ b/test/reg/testref/getkf_nonlinear.ref @@ -4,141 +4,141 @@ Initial state for member 1: Fields: - air_temperature (2 levels): + min = -1.0272916804277608e+01 - + max = 1.8740481760074051e+01 - + mean = 2.1961001758681213e+00 - + stddev = 2.8514424414275474e+00 + + max = 1.8740481760074047e+01 + + mean = 2.1961001758681178e+00 + + stddev = 2.8514424414275528e+00 Initial state for member 2: - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0345315914106127e+01 - + max = 1.8571805261116310e+01 - + mean = 2.2192193482732367e+00 - + stddev = 2.8478755388573083e+00 + + max = 1.8571805261116307e+01 + + mean = 2.2192193482732310e+00 + + stddev = 2.8478755388573123e+00 Initial state for member 3: - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0279656748962816e+01 - + max = 1.8644069152415483e+01 + + max = 1.8644069152415479e+01 + mean = 2.1877262258572960e+00 - + stddev = 2.8441569245650746e+00 + + stddev = 2.8441569245650706e+00 Initial state for member 4: - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0240502338137325e+01 - + max = 1.8675853918837795e+01 - + mean = 2.2036774297260631e+00 - + stddev = 2.8472497796821701e+00 + + max = 1.8675853918837792e+01 + + mean = 2.2036774297260626e+00 + + stddev = 2.8472497796821692e+00 Initial state for member 5: - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0283807777570551e+01 - + max = 1.8644989920588163e+01 - + mean = 2.1932883221742983e+00 - + stddev = 2.8473063130026253e+00 + + max = 1.8644989920588159e+01 + + mean = 2.1932883221742965e+00 + + stddev = 2.8473063130026275e+00 Initial state for member 6: - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0406692783814259e+01 - + max = 1.8819427931232656e+01 - + mean = 2.2029788605963052e+00 - + stddev = 2.8445180090287181e+00 + + max = 1.8819427931232653e+01 + + mean = 2.2029788605963097e+00 + + stddev = 2.8445180090287163e+00 Initial state for member 7: - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0220623189263415e+01 - + max = 1.8876841842386153e+01 - + mean = 2.1960289287861086e+00 - + stddev = 2.8482208026553208e+00 + + max = 1.8876841842386149e+01 + + mean = 2.1960289287861112e+00 + + stddev = 2.8482208026553186e+00 Initial state for member 8: - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0272136420138436e+01 - + max = 1.8664322482275654e+01 - + mean = 2.1887556044999519e+00 - + stddev = 2.8412451259362306e+00 + + max = 1.8664322482275651e+01 + + mean = 2.1887556044999452e+00 + + stddev = 2.8412451259362301e+00 Initial state for member 9: - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0137886041649526e+01 - + max = 1.8497301888272666e+01 - + mean = 2.1985712328394129e+00 - + stddev = 2.8474739075373883e+00 + + max = 1.8497301888272663e+01 + + mean = 2.1985712328394080e+00 + + stddev = 2.8474739075373869e+00 Initial state for member 10: - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0323479917719220e+01 - + max = 1.8652489165589472e+01 - + mean = 2.1986593782124051e+00 - + stddev = 2.8505535212368502e+00 + + max = 1.8652489165589468e+01 + + mean = 2.1986593782124038e+00 + + stddev = 2.8505535212368454e+00 H(x) for member 1: -air_temperature nobs= 25 Min=-8.0033176611693579e-01, Max=4.8198743621834392e+00, RMS=2.4503244194634628e+00 +air_temperature nobs= 25 Min=-8.0033176611693557e-01, Max=4.8198743621834312e+00, RMS=2.4503244194634619e+00 H(x) for member 2: -air_temperature nobs= 25 Min=-7.7341243668325133e-01, Max=4.7427526471731349e+00, RMS=2.4493810775645892e+00 +air_temperature nobs= 25 Min=-7.7341243668325110e-01, Max=4.7427526471731278e+00, RMS=2.4493810775645883e+00 H(x) for member 3: -air_temperature nobs= 25 Min=-7.3237270521936493e-01, Max=4.8281233727289425e+00, RMS=2.4161000289397401e+00 +air_temperature nobs= 25 Min=-7.3237270521936482e-01, Max=4.8281233727289354e+00, RMS=2.4161000289397396e+00 H(x) for member 4: -air_temperature nobs= 25 Min=-8.3278958305467710e-01, Max=4.6785267644627782e+00, RMS=2.4420403669141129e+00 +air_temperature nobs= 25 Min=-8.3278958305467699e-01, Max=4.6785267644627710e+00, RMS=2.4420403669141120e+00 H(x) for member 5: -air_temperature nobs= 25 Min=-8.3424927599662113e-01, Max=4.5588086888390489e+00, RMS=2.4341937566450422e+00 +air_temperature nobs= 25 Min=-8.3424927599662102e-01, Max=4.5588086888390409e+00, RMS=2.4341937566450413e+00 H(x) for member 6: -air_temperature nobs= 25 Min=-7.1947793240223246e-01, Max=4.6839379579677320e+00, RMS=2.4422220408939115e+00 +air_temperature nobs= 25 Min=-7.1947793240223246e-01, Max=4.6839379579677249e+00, RMS=2.4422220408939097e+00 H(x) for member 7: -air_temperature nobs= 25 Min=-8.9702681258623451e-01, Max=4.7376379880923762e+00, RMS=2.4503699734336322e+00 +air_temperature nobs= 25 Min=-8.9702681258623462e-01, Max=4.7376379880923682e+00, RMS=2.4503699734336313e+00 H(x) for member 8: -air_temperature nobs= 25 Min=-8.5064128616631240e-01, Max=4.8176875859650128e+00, RMS=2.4450249099270289e+00 +air_temperature nobs= 25 Min=-8.5064128616631252e-01, Max=4.8176875859650057e+00, RMS=2.4450249099270280e+00 H(x) for member 9: -air_temperature nobs= 25 Min=-8.4394015351321827e-01, Max=4.6950845406414174e+00, RMS=2.4337539915153563e+00 +air_temperature nobs= 25 Min=-8.4394015351321827e-01, Max=4.6950845406414103e+00, RMS=2.4337539915153554e+00 H(x) for member 10: -air_temperature nobs= 25 Min=-7.9330489397198778e-01, Max=4.7083014262153355e+00, RMS=2.4339895614870266e+00 +air_temperature nobs= 25 Min=-7.9330489397198756e-01, Max=4.7083014262153284e+00, RMS=2.4339895614870257e+00 H(x) ensemble background mean: -air_temperature nobs= 25 Min=-8.0775468457108346e-01, Max=4.7270735334269220e+00, RMS=2.4384935788933442e+00 +air_temperature nobs= 25 Min=-8.0775468457108346e-01, Max=4.7270735334269141e+00, RMS=2.4384935788933428e+00 background y - H(x): -air_temperature nobs= 25 Min=-2.3507591848416887e-01, Max=2.6369343316652438e-01, RMS=1.4649775063202392e-01 +air_temperature nobs= 25 Min=-2.3507591848417064e-01, Max=2.6369343316652433e-01, RMS=1.4649775063202405e-01 Background mean : - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0278301793563930e+01 - + max = 1.8678758332278843e+01 - + mean = 2.1985005506833204e+00 - + stddev = 2.8454484711447656e+00 + + max = 1.8678758332278839e+01 + + mean = 2.1985005506833186e+00 + + stddev = 2.8454484711447687e+00 Analysis mean : - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0267858117053725e+01 - + max = 1.8591510736129553e+01 - + mean = 2.2020219314428973e+00 - + stddev = 2.8440813362997277e+00 + + min = -1.0267894625343049e+01 + + max = 1.8591514416331925e+01 + + mean = 2.2020223360665776e+00 + + stddev = 2.8440807683617217e+00 Analysis mean increment : - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -3.0869792050202260e-01 - + max = 3.5868891971667527e-01 - + mean = 3.5213807595765238e-03 - + stddev = 8.1150620789074049e-02 + + min = -3.0869790492080984e-01 + + max = 3.5870935083111877e-01 + + mean = 3.5217853832612955e-03 + + stddev = 8.1148641857371015e-02 Forecast variance : - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] @@ -146,40 +146,40 @@ Forecast variance : - air_temperature (2 levels): + min = 4.6034591952436623e-04 + max = 3.6294641722342598e-02 - + mean = 9.9273235952287049e-03 + + mean = 9.9273235952286910e-03 + stddev = 4.6717843128316842e-03 Analysis variance : - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = 4.0482264433241598e-04 - + max = 3.0755248308003937e-02 - + mean = 8.5478016138302200e-03 - + stddev = 4.0085150278443959e-03 + + min = 4.0556394300518707e-04 + + max = 3.0755254011557241e-02 + + mean = 8.5478809181325863e-03 + + stddev = 4.0083051286260146e-03 H(x) for member 1: -air_temperature nobs= 25 Min=-7.6382821093977382e-01, Max=4.7584468335224930e+00, RMS=2.4371538169704690e+00 +air_temperature nobs= 25 Min=-7.6381830724515254e-01, Max=4.7588040630924775e+00, RMS=2.4372715772822771e+00 H(x) for member 2: -air_temperature nobs= 25 Min=-7.4169735947068582e-01, Max=4.6895000369642013e+00, RMS=2.4395228941307767e+00 +air_temperature nobs= 25 Min=-7.4165676996923380e-01, Max=4.6897318211180359e+00, RMS=2.4394798820975248e+00 H(x) for member 3: -air_temperature nobs= 25 Min=-7.1058363480398890e-01, Max=4.7650046486989757e+00, RMS=2.4088898381551140e+00 +air_temperature nobs= 25 Min=-7.1058466893163275e-01, Max=4.7647778528161382e+00, RMS=2.4089279634785159e+00 H(x) for member 4: -air_temperature nobs= 25 Min=-7.9726958150301974e-01, Max=4.6316951483191549e+00, RMS=2.4320503612980264e+00 +air_temperature nobs= 25 Min=-7.9725068554990919e-01, Max=4.6315817818069753e+00, RMS=2.4319935139087465e+00 H(x) for member 5: -air_temperature nobs= 25 Min=-7.9965882187205384e-01, Max=4.5234482202146333e+00, RMS=2.4240410574223472e+00 +air_temperature nobs= 25 Min=-7.9965827072765450e-01, Max=4.5233251454884336e+00, RMS=2.4239947043525762e+00 H(x) for member 6: -air_temperature nobs= 25 Min=-6.9677976311562229e-01, Max=4.6304895704087077e+00, RMS=2.4329311773746953e+00 +air_temperature nobs= 25 Min=-6.9679561041505078e-01, Max=4.6305445759721469e+00, RMS=2.4328376414149369e+00 H(x) for member 7: -air_temperature nobs= 25 Min=-8.5699789783457425e-01, Max=4.6796877508953294e+00, RMS=2.4382782808050636e+00 +air_temperature nobs= 25 Min=-8.5695521796464491e-01, Max=4.6798966092286802e+00, RMS=2.4384236875628247e+00 H(x) for member 8: -air_temperature nobs= 25 Min=-8.1799931583250385e-01, Max=4.7597893508497524e+00, RMS=2.4347344820705130e+00 +air_temperature nobs= 25 Min=-8.1805232332378375e-01, Max=4.7599989526028423e+00, RMS=2.4347344764270415e+00 H(x) for member 9: -air_temperature nobs= 25 Min=-8.0427379608589944e-01, Max=4.6417688370257277e+00, RMS=2.4238088624570837e+00 +air_temperature nobs= 25 Min=-8.0428843488305346e-01, Max=4.6419443145244212e+00, RMS=2.4237584579131299e+00 H(x) for member 10: -air_temperature nobs= 25 Min=-7.6231738256458748e-01, Max=4.6594217363845045e+00, RMS=2.4245658557373448e+00 +air_temperature nobs= 25 Min=-7.6227761651844506e-01, Max=4.6587289698892125e+00, RMS=2.4245629681804752e+00 H(x) ensemble analysis mean: -air_temperature nobs= 25 Min=-7.7514057640227108e-01, Max=4.6739252133283475e+00, RMS=2.4285880851102810e+00 +air_temperature nobs= 25 Min=-7.7513379055285614e-01, Max=4.6739334086539364e+00, RMS=2.4285894901670666e+00 analysis y - H(x): -air_temperature nobs= 25 Min=-1.1650182525202135e-01, Max=1.7882287419197751e-01, RMS=9.2710257122365330e-02 -ombg RMS: 1.4649775063202392e-01 -oman RMS: 9.2710257122365330e-02 +air_temperature nobs= 25 Min=-1.1649017756368174e-01, Max=1.7882796032666315e-01, RMS=9.2697725826235158e-02 +ombg RMS: 1.4649775063202405e-01 +oman RMS: 9.2697725826235158e-02 diff --git a/test/reg/testref/letkf_linear_4d.ref b/test/reg/testref/letkf_linear_4d.ref index 6d258bb..7c21e62 100644 --- a/test/reg/testref/letkf_linear_4d.ref +++ b/test/reg/testref/letkf_linear_4d.ref @@ -1,44 +1,68 @@ Initial state for member 1: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0272916804277608e+01 + + max = 1.9677845626148326e+01 + + mean = 2.4608260676686631e+00 + + stddev = 3.1983604806334918e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0272916804277608e+01 - + max = 1.8740481760074051e+01 - + mean = 2.1961001758681182e+00 - + stddev = 2.8514424414275532e+00 + + max = 1.8740481760074047e+01 + + mean = 2.1961001758681178e+00 + + stddev = 2.8514424414275528e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0272916804277608e+01 - + max = 1.7543794994608628e+01 - + mean = 2.1252781169363306e+00 - + stddev = 2.5242824652600757e+00 + + max = 1.7543794994608625e+01 + + mean = 2.1252781169363302e+00 + + stddev = 2.5242824652600762e+00 Initial state for member 2: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0345315914106127e+01 + + max = 1.9509169127190585e+01 + + mean = 2.4839452400737745e+00 + + stddev = 3.1910382432457309e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0345315914106127e+01 - + max = 1.8571805261116310e+01 - + mean = 2.2192193482732296e+00 - + stddev = 2.8478755388573127e+00 + + max = 1.8571805261116307e+01 + + mean = 2.2192193482732310e+00 + + stddev = 2.8478755388573123e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0345315914106127e+01 - + max = 1.7375118495650888e+01 - + mean = 2.1483972893414474e+00 + + max = 1.7375118495650884e+01 + + mean = 2.1483972893414469e+00 + stddev = 2.5201721677971483e+00 Initial state for member 3: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0279656748962816e+01 + + max = 1.9581433018489758e+01 + + mean = 2.4524521176578351e+00 + + stddev = 3.1928803598372117e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0279656748962816e+01 - + max = 1.8644069152415483e+01 + + max = 1.8644069152415479e+01 + mean = 2.1877262258572960e+00 + stddev = 2.8441569245650706e+00 - Valid time: 2010-01-01T06:00:00Z @@ -46,204 +70,294 @@ Initial state for member 3: Fields: - air_temperature (2 levels): + min = -1.0279656748962816e+01 - + max = 1.7447382386950061e+01 - + mean = 2.1169041669255124e+00 + + max = 1.7447382386950057e+01 + + mean = 2.1169041669255129e+00 + stddev = 2.5177118723945289e+00 Initial state for member 4: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0240502338137325e+01 + + max = 1.9613217784912070e+01 + + mean = 2.4684033215266141e+00 + + stddev = 3.1901143631877100e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0240502338137325e+01 - + max = 1.8675853918837795e+01 - + mean = 2.2036774297260622e+00 - + stddev = 2.8472497796821687e+00 + + max = 1.8675853918837792e+01 + + mean = 2.2036774297260626e+00 + + stddev = 2.8472497796821692e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0240502338137325e+01 - + max = 1.7479167153372373e+01 - + mean = 2.1328553707942635e+00 - + stddev = 2.5186251374861062e+00 + + max = 1.7479167153372369e+01 + + mean = 2.1328553707942643e+00 + + stddev = 2.5186251374861057e+00 Initial state for member 5: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0283807777570551e+01 + + max = 1.9582353786662438e+01 + + mean = 2.4580142139748378e+00 + + stddev = 3.1966066810668345e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0283807777570551e+01 - + max = 1.8644989920588163e+01 - + mean = 2.1932883221742960e+00 + + max = 1.8644989920588159e+01 + + mean = 2.1932883221742965e+00 + stddev = 2.8473063130026275e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0283807777570551e+01 - + max = 1.7448303155122741e+01 + + max = 1.7448303155122737e+01 + mean = 2.1224662632425075e+00 + stddev = 2.5186490031501503e+00 Initial state for member 6: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0406692783814259e+01 + + max = 1.9756791797306931e+01 + + mean = 2.4677047523968532e+00 + + stddev = 3.1943041888207104e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0406692783814259e+01 - + max = 1.8819427931232656e+01 - + mean = 2.2029788605963101e+00 - + stddev = 2.8445180090287168e+00 + + max = 1.8819427931232653e+01 + + mean = 2.2029788605963097e+00 + + stddev = 2.8445180090287163e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0406692783814259e+01 - + max = 1.7622741165767234e+01 + + max = 1.7622741165767231e+01 + mean = 2.1321568016645167e+00 + stddev = 2.5168901199567433e+00 Initial state for member 7: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0220623189263415e+01 + + max = 1.9814205708460428e+01 + + mean = 2.4607548205866490e+00 + + stddev = 3.1952640926878892e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0220623189263415e+01 - + max = 1.8876841842386153e+01 + + max = 1.8876841842386149e+01 + mean = 2.1960289287861112e+00 - + stddev = 2.8482208026553191e+00 + + stddev = 2.8482208026553186e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0220623189263415e+01 - + max = 1.7680155076920730e+01 - + mean = 2.1252068698543223e+00 - + stddev = 2.5217580620433693e+00 + + max = 1.7680155076920727e+01 + + mean = 2.1252068698543214e+00 + + stddev = 2.5217580620433697e+00 Initial state for member 8: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0272136420138436e+01 + + max = 1.9601686348349929e+01 + + mean = 2.4534814963004927e+00 + + stddev = 3.1936019692500373e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0272136420138436e+01 - + max = 1.8664322482275654e+01 - + mean = 2.1887556044999461e+00 - + stddev = 2.8412451259362306e+00 + + max = 1.8664322482275651e+01 + + mean = 2.1887556044999452e+00 + + stddev = 2.8412451259362301e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0272136420138436e+01 - + max = 1.7467635716810232e+01 - + mean = 2.1179335455681629e+00 + + max = 1.7467635716810229e+01 + + mean = 2.1179335455681625e+00 + stddev = 2.5179398372939685e+00 Initial state for member 9: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0137886041649526e+01 + + max = 1.9434665754346941e+01 + + mean = 2.4632971246399511e+00 + + stddev = 3.1949428489924161e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0137886041649526e+01 - + max = 1.8497301888272666e+01 - + mean = 2.1985712328394085e+00 - + stddev = 2.8474739075373878e+00 + + max = 1.8497301888272663e+01 + + mean = 2.1985712328394080e+00 + + stddev = 2.8474739075373869e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0137886041649526e+01 - + max = 1.7300615122807244e+01 - + mean = 2.1277491739076266e+00 - + stddev = 2.5172226982177821e+00 + + max = 1.7300615122807240e+01 + + mean = 2.1277491739076257e+00 + + stddev = 2.5172226982177826e+00 Initial state for member 10: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0323479917719220e+01 + + max = 1.9589853031663747e+01 + + mean = 2.4633852700129450e+00 + + stddev = 3.1983965297893397e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0323479917719220e+01 - + max = 1.8652489165589472e+01 - + mean = 2.1986593782124042e+00 + + max = 1.8652489165589468e+01 + + mean = 2.1986593782124038e+00 + stddev = 2.8505535212368454e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0323479917719220e+01 - + max = 1.7455802400124050e+01 - + mean = 2.1278373192806099e+00 - + stddev = 2.5151283324109097e+00 + + max = 1.7455802400124046e+01 + + mean = 2.1278373192806095e+00 + + stddev = 2.5151283324109093e+00 H(x) for member 1: -air_temperature nobs= 50 Min=-1.2695618997552252e+00, Max=5.0514508890641263e+00, RMS=2.4041594135950661e+00 +air_temperature nobs= 50 Min=-1.2695618997552245e+00, Max=5.0514508890641245e+00, RMS=2.4041594135950661e+00 H(x) for member 2: -air_temperature nobs= 50 Min=-1.3208795226403331e+00, Max=5.1054957476166418e+00, RMS=2.4164010784150358e+00 +air_temperature nobs= 50 Min=-1.3208795226403325e+00, Max=5.1054957476166400e+00, RMS=2.4164010784150358e+00 H(x) for member 3: -air_temperature nobs= 50 Min=-1.3365986141600465e+00, Max=5.1050667757821611e+00, RMS=2.3922068649386721e+00 +air_temperature nobs= 50 Min=-1.3365986141600459e+00, Max=5.1050667757821593e+00, RMS=2.3922068649386721e+00 H(x) for member 4: -air_temperature nobs= 50 Min=-1.1680394913294172e+00, Max=5.2304605650437406e+00, RMS=2.3951684520307550e+00 +air_temperature nobs= 50 Min=-1.1680394913294165e+00, Max=5.2304605650437388e+00, RMS=2.3951684520307550e+00 H(x) for member 5: -air_temperature nobs= 50 Min=-1.3527164104618883e+00, Max=5.1972187640798619e+00, RMS=2.4058700195403713e+00 +air_temperature nobs= 50 Min=-1.3527164104618876e+00, Max=5.1972187640798602e+00, RMS=2.4058700195403713e+00 H(x) for member 6: -air_temperature nobs= 50 Min=-1.2222796957232809e+00, Max=5.1090079401341608e+00, RMS=2.3988762858953909e+00 +air_temperature nobs= 50 Min=-1.2222796957232802e+00, Max=5.1090079401341590e+00, RMS=2.3988762858953909e+00 H(x) for member 7: -air_temperature nobs= 50 Min=-1.2717312946208572e+00, Max=5.0776092740428975e+00, RMS=2.4115399844833574e+00 +air_temperature nobs= 50 Min=-1.2717312946208565e+00, Max=5.0776092740428957e+00, RMS=2.4115399844833574e+00 H(x) for member 8: -air_temperature nobs= 50 Min=-1.2938306498565531e+00, Max=5.0743719744217923e+00, RMS=2.3991798013467434e+00 +air_temperature nobs= 50 Min=-1.2938306498565524e+00, Max=5.0743719744217906e+00, RMS=2.3991798013467438e+00 H(x) for member 9: -air_temperature nobs= 50 Min=-1.5175354908444738e+00, Max=5.1994875105751088e+00, RMS=2.4039977272306801e+00 +air_temperature nobs= 50 Min=-1.5175354908444731e+00, Max=5.1994875105751071e+00, RMS=2.4039977272306796e+00 H(x) for member 10: -air_temperature nobs= 50 Min=-1.1628020684459066e+00, Max=5.0503043415678075e+00, RMS=2.3910141656603314e+00 +air_temperature nobs= 50 Min=-1.1628020684459059e+00, Max=5.0503043415678057e+00, RMS=2.3910141656603314e+00 H(x) ensemble background mean: -air_temperature nobs= 50 Min=-1.2915975137837981e+00, Max=5.1200473782328313e+00, RMS=2.4005666633340832e+00 +air_temperature nobs= 50 Min=-1.2915975137837976e+00, Max=5.1200473782328295e+00, RMS=2.4005666633340836e+00 background y - H(x): -air_temperature nobs= 50 Min=-3.0904601954408850e-01, Max=2.4788732602700980e-01, RMS=1.2718289261558888e-01 +air_temperature nobs= 50 Min=-3.0904607403139422e-01, Max=4.2503948242390877e-01, RMS=1.4489786187411210e-01 Background mean : +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0278301793563930e+01 + + max = 1.9616122198353118e+01 + + mean = 2.4632264424838528e+00 + + stddev = 3.1931643516640591e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0278301793563930e+01 - + max = 1.8678758332278843e+01 - + mean = 2.1985005506833173e+00 - + stddev = 2.8454484711447692e+00 + + max = 1.8678758332278839e+01 + + mean = 2.1985005506833186e+00 + + stddev = 2.8454484711447687e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0278301793563930e+01 - + max = 1.7482071566813421e+01 + + max = 1.7482071566813417e+01 + mean = 2.1276784917515221e+00 + stddev = 2.5170789792474042e+00 Analysis mean : +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0344294798986292e+01 + + max = 1.9576028887721325e+01 + + mean = 2.4640237350325238e+00 + + stddev = 3.1929130257880263e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0349656866548454e+01 - + max = 1.8630768031596617e+01 - + mean = 2.2059727047276412e+00 - + stddev = 2.8435166364882152e+00 + + min = -1.0344294798986292e+01 + + max = 1.8638665021647046e+01 + + mean = 2.1992978432319803e+00 + + stddev = 2.8427786834913840e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0349656866548454e+01 - + max = 1.7434081266131194e+01 - + mean = 2.1351506457958465e+00 - + stddev = 2.5152958790498405e+00 + + min = -1.0344294798986292e+01 + + max = 1.7441978256181624e+01 + + mean = 2.1284757843001882e+00 + + stddev = 2.5136874391858512e+00 +Analysis mean increment : +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -2.8773179864033316e-01 + + max = 2.7945463804838555e-01 + + mean = 7.9729254865935884e-04 + + stddev = 8.8521587291884321e-02 Analysis mean increment : - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -2.3180136230557624e-01 - + max = 2.4052636708797248e-01 - + mean = 7.4721540443193736e-03 - + stddev = 7.9862988081685821e-02 + + min = -2.8773179864033316e-01 + + max = 2.7945463804838522e-01 + + mean = 7.9729254865939137e-04 + + stddev = 8.8521587291884279e-02 Analysis mean increment : - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -2.3180136230557619e-01 - + max = 2.4052636708796982e-01 - + mean = 7.4721540443193832e-03 - + stddev = 7.9862988081685807e-02 + + min = -2.8773179864033316e-01 + + max = 2.7945463804838511e-01 + + mean = 7.9729254865939039e-04 + + stddev = 8.8521587291884279e-02 +Forecast variance : +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = 4.6034591952436564e-04 + + max = 3.6294641722342640e-02 + + mean = 9.9273235952286910e-03 + + stddev = 4.6717843128316850e-03 Forecast variance : - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] @@ -252,7 +366,7 @@ Forecast variance : + min = 4.6034591952436623e-04 + max = 3.6294641722342598e-02 + mean = 9.9273235952286910e-03 - + stddev = 4.6717843128316850e-03 + + stddev = 4.6717843128316842e-03 Forecast variance : - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] @@ -263,46 +377,55 @@ Forecast variance : + mean = 9.9273235952286910e-03 + stddev = 4.6717843128316859e-03 Analysis variance : +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = 3.5944749374122812e-04 + + max = 2.8343292475071127e-02 + + mean = 7.6767697298578703e-03 + + stddev = 3.6123653670426126e-03 +Analysis variance : - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = 3.5944746468038726e-04 - + max = 2.8343289605209567e-02 - + mean = 7.6767696481328503e-03 - + stddev = 3.6123653099917569e-03 + + min = 3.5944749374122850e-04 + + max = 2.8343292475071100e-02 + + mean = 7.6767697298578711e-03 + + stddev = 3.6123653670426122e-03 Analysis variance : - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = 3.5944746468038682e-04 - + max = 2.8343289605209567e-02 - + mean = 7.6767696481328520e-03 - + stddev = 3.6123653099917560e-03 + + min = 3.5944749374122850e-04 + + max = 2.8343292475071100e-02 + + mean = 7.6767697298578729e-03 + + stddev = 3.6123653670426109e-03 H(x) for member 1: -air_temperature nobs= 50 Min=-1.2602201721165966e+00, Max=5.1330192634701417e+00, RMS=2.4059920585876018e+00 +air_temperature nobs= 50 Min=-1.2676489427656799e+00, Max=5.1403154034170990e+00, RMS=2.4016023881337074e+00 H(x) for member 2: -air_temperature nobs= 50 Min=-1.2923765501473736e+00, Max=5.1696522315829680e+00, RMS=2.4198597548861955e+00 +air_temperature nobs= 50 Min=-1.2998053347662955e+00, Max=5.1769483752552157e+00, RMS=2.4152747233980052e+00 H(x) for member 3: -air_temperature nobs= 50 Min=-1.3156692332123112e+00, Max=5.1812981344625042e+00, RMS=2.3971151033389173e+00 +air_temperature nobs= 50 Min=-1.3230980187625556e+00, Max=5.1885942725468164e+00, RMS=2.3932075704193410e+00 H(x) for member 4: -air_temperature nobs= 50 Min=-1.1658482125495266e+00, Max=5.2888126248359368e+00, RMS=2.4000087731710060e+00 +air_temperature nobs= 50 Min=-1.1732769757480292e+00, Max=5.2961087647828942e+00, RMS=2.3955793525083600e+00 H(x) for member 5: -air_temperature nobs= 50 Min=-1.3282343468163800e+00, Max=5.2576914304732965e+00, RMS=2.4097572692891300e+00 +air_temperature nobs= 50 Min=-1.3356631286413341e+00, Max=5.2649875629696732e+00, RMS=2.4056895272156376e+00 H(x) for member 6: -air_temperature nobs= 50 Min=-1.2278214527760338e+00, Max=5.1828877335279868e+00, RMS=2.4046086954634092e+00 +air_temperature nobs= 50 Min=-1.2352502308756976e+00, Max=5.1901838660243635e+00, RMS=2.4006474862168012e+00 H(x) for member 7: -air_temperature nobs= 50 Min=-1.2616896557156634e+00, Max=5.1563734473526166e+00, RMS=2.4141614303676544e+00 +air_temperature nobs= 50 Min=-1.2691184338153272e+00, Max=5.1636695872995739e+00, RMS=2.4098034739888710e+00 H(x) for member 8: -air_temperature nobs= 50 Min=-1.2723836621735882e+00, Max=5.1497118124603913e+00, RMS=2.4033651427839540e+00 +air_temperature nobs= 50 Min=-1.2798124495864778e+00, Max=5.1570079561326390e+00, RMS=2.3992823652166573e+00 H(x) for member 9: -air_temperature nobs= 50 Min=-1.4768283909653019e+00, Max=5.2600023249029801e+00, RMS=2.4074343523044610e+00 +air_temperature nobs= 50 Min=-1.4842571541638045e+00, Max=5.2672984723005181e+00, RMS=2.4032589466082093e+00 H(x) for member 10: -air_temperature nobs= 50 Min=-1.1730448177550625e+00, Max=5.1348378570615933e+00, RMS=2.3968069449151539e+00 +air_temperature nobs= 50 Min=-1.1804736033053069e+00, Max=5.1421339970085507e+00, RMS=2.3933314993484078e+00 H(x) ensemble analysis mean: -air_temperature nobs= 50 Min=-1.2774116494227838e+00, Max=5.1914286860130412e+00, RMS=2.4049785937466508e+00 +air_temperature nobs= 50 Min=-1.2848404272430509e+00, Max=5.1987248257737342e+00, RMS=2.4008331781248939e+00 analysis y - H(x): -air_temperature nobs= 50 Min=-3.0186856422944142e-01, Max=1.7205103741094346e-01, RMS=9.0847493539626192e-02 -ombg RMS: 1.2718289261558888e-01 -oman RMS: 9.0847493539626192e-02 +air_temperature nobs= 50 Min=-2.4469863608637291e-01, Max=3.4636203488300410e-01, RMS=1.0628607375331114e-01 +ombg RMS: 1.4489786187411210e-01 +oman RMS: 1.0628607375331114e-01 diff --git a/test/reg/testref/letkf_nonlinear_4d.ref b/test/reg/testref/letkf_nonlinear_4d.ref index d6abd46..a571ef5 100644 --- a/test/reg/testref/letkf_nonlinear_4d.ref +++ b/test/reg/testref/letkf_nonlinear_4d.ref @@ -1,44 +1,68 @@ Initial state for member 1: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0272916804277608e+01 + + max = 1.9677845626148326e+01 + + mean = 2.4608260676686631e+00 + + stddev = 3.1983604806334918e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0272916804277608e+01 - + max = 1.8740481760074051e+01 - + mean = 2.1961001758681182e+00 - + stddev = 2.8514424414275532e+00 + + max = 1.8740481760074047e+01 + + mean = 2.1961001758681178e+00 + + stddev = 2.8514424414275528e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0272916804277608e+01 - + max = 1.7543794994608628e+01 - + mean = 2.1252781169363306e+00 - + stddev = 2.5242824652600757e+00 + + max = 1.7543794994608625e+01 + + mean = 2.1252781169363302e+00 + + stddev = 2.5242824652600762e+00 Initial state for member 2: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0345315914106127e+01 + + max = 1.9509169127190585e+01 + + mean = 2.4839452400737745e+00 + + stddev = 3.1910382432457309e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0345315914106127e+01 - + max = 1.8571805261116310e+01 - + mean = 2.2192193482732296e+00 - + stddev = 2.8478755388573127e+00 + + max = 1.8571805261116307e+01 + + mean = 2.2192193482732310e+00 + + stddev = 2.8478755388573123e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0345315914106127e+01 - + max = 1.7375118495650888e+01 - + mean = 2.1483972893414474e+00 + + max = 1.7375118495650884e+01 + + mean = 2.1483972893414469e+00 + stddev = 2.5201721677971483e+00 Initial state for member 3: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0279656748962816e+01 + + max = 1.9581433018489758e+01 + + mean = 2.4524521176578351e+00 + + stddev = 3.1928803598372117e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0279656748962816e+01 - + max = 1.8644069152415483e+01 + + max = 1.8644069152415479e+01 + mean = 2.1877262258572960e+00 + stddev = 2.8441569245650706e+00 - Valid time: 2010-01-01T06:00:00Z @@ -46,204 +70,294 @@ Initial state for member 3: Fields: - air_temperature (2 levels): + min = -1.0279656748962816e+01 - + max = 1.7447382386950061e+01 - + mean = 2.1169041669255124e+00 + + max = 1.7447382386950057e+01 + + mean = 2.1169041669255129e+00 + stddev = 2.5177118723945289e+00 Initial state for member 4: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0240502338137325e+01 + + max = 1.9613217784912070e+01 + + mean = 2.4684033215266141e+00 + + stddev = 3.1901143631877100e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0240502338137325e+01 - + max = 1.8675853918837795e+01 - + mean = 2.2036774297260622e+00 - + stddev = 2.8472497796821687e+00 + + max = 1.8675853918837792e+01 + + mean = 2.2036774297260626e+00 + + stddev = 2.8472497796821692e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0240502338137325e+01 - + max = 1.7479167153372373e+01 - + mean = 2.1328553707942635e+00 - + stddev = 2.5186251374861062e+00 + + max = 1.7479167153372369e+01 + + mean = 2.1328553707942643e+00 + + stddev = 2.5186251374861057e+00 Initial state for member 5: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0283807777570551e+01 + + max = 1.9582353786662438e+01 + + mean = 2.4580142139748378e+00 + + stddev = 3.1966066810668345e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0283807777570551e+01 - + max = 1.8644989920588163e+01 - + mean = 2.1932883221742960e+00 + + max = 1.8644989920588159e+01 + + mean = 2.1932883221742965e+00 + stddev = 2.8473063130026275e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0283807777570551e+01 - + max = 1.7448303155122741e+01 + + max = 1.7448303155122737e+01 + mean = 2.1224662632425075e+00 + stddev = 2.5186490031501503e+00 Initial state for member 6: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0406692783814259e+01 + + max = 1.9756791797306931e+01 + + mean = 2.4677047523968532e+00 + + stddev = 3.1943041888207104e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0406692783814259e+01 - + max = 1.8819427931232656e+01 - + mean = 2.2029788605963101e+00 - + stddev = 2.8445180090287168e+00 + + max = 1.8819427931232653e+01 + + mean = 2.2029788605963097e+00 + + stddev = 2.8445180090287163e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0406692783814259e+01 - + max = 1.7622741165767234e+01 + + max = 1.7622741165767231e+01 + mean = 2.1321568016645167e+00 + stddev = 2.5168901199567433e+00 Initial state for member 7: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0220623189263415e+01 + + max = 1.9814205708460428e+01 + + mean = 2.4607548205866490e+00 + + stddev = 3.1952640926878892e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0220623189263415e+01 - + max = 1.8876841842386153e+01 + + max = 1.8876841842386149e+01 + mean = 2.1960289287861112e+00 - + stddev = 2.8482208026553191e+00 + + stddev = 2.8482208026553186e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0220623189263415e+01 - + max = 1.7680155076920730e+01 - + mean = 2.1252068698543223e+00 - + stddev = 2.5217580620433693e+00 + + max = 1.7680155076920727e+01 + + mean = 2.1252068698543214e+00 + + stddev = 2.5217580620433697e+00 Initial state for member 8: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0272136420138436e+01 + + max = 1.9601686348349929e+01 + + mean = 2.4534814963004927e+00 + + stddev = 3.1936019692500373e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0272136420138436e+01 - + max = 1.8664322482275654e+01 - + mean = 2.1887556044999461e+00 - + stddev = 2.8412451259362306e+00 + + max = 1.8664322482275651e+01 + + mean = 2.1887556044999452e+00 + + stddev = 2.8412451259362301e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0272136420138436e+01 - + max = 1.7467635716810232e+01 - + mean = 2.1179335455681629e+00 + + max = 1.7467635716810229e+01 + + mean = 2.1179335455681625e+00 + stddev = 2.5179398372939685e+00 Initial state for member 9: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0137886041649526e+01 + + max = 1.9434665754346941e+01 + + mean = 2.4632971246399511e+00 + + stddev = 3.1949428489924161e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0137886041649526e+01 - + max = 1.8497301888272666e+01 - + mean = 2.1985712328394085e+00 - + stddev = 2.8474739075373878e+00 + + max = 1.8497301888272663e+01 + + mean = 2.1985712328394080e+00 + + stddev = 2.8474739075373869e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0137886041649526e+01 - + max = 1.7300615122807244e+01 - + mean = 2.1277491739076266e+00 - + stddev = 2.5172226982177821e+00 + + max = 1.7300615122807240e+01 + + mean = 2.1277491739076257e+00 + + stddev = 2.5172226982177826e+00 Initial state for member 10: +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0323479917719220e+01 + + max = 1.9589853031663747e+01 + + mean = 2.4633852700129450e+00 + + stddev = 3.1983965297893397e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0323479917719220e+01 - + max = 1.8652489165589472e+01 - + mean = 2.1986593782124042e+00 + + max = 1.8652489165589468e+01 + + mean = 2.1986593782124038e+00 + stddev = 2.8505535212368454e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0323479917719220e+01 - + max = 1.7455802400124050e+01 - + mean = 2.1278373192806099e+00 - + stddev = 2.5151283324109097e+00 + + max = 1.7455802400124046e+01 + + mean = 2.1278373192806095e+00 + + stddev = 2.5151283324109093e+00 H(x) for member 1: -air_temperature nobs= 50 Min=-1.2695618989387294e+00, Max=5.0514508890327390e+00, RMS=2.4041594136466431e+00 +air_temperature nobs= 50 Min=-1.2695618989387289e+00, Max=5.0514508890327381e+00, RMS=2.4041594136466435e+00 H(x) for member 2: -air_temperature nobs= 50 Min=-1.3208795227495509e+00, Max=5.1054957471914317e+00, RMS=2.4164010783353675e+00 +air_temperature nobs= 50 Min=-1.3208795227495504e+00, Max=5.1054957471914317e+00, RMS=2.4164010783353675e+00 H(x) for member 3: -air_temperature nobs= 50 Min=-1.3365986129438956e+00, Max=5.1050667755833823e+00, RMS=2.3922068645852597e+00 +air_temperature nobs= 50 Min=-1.3365986129438949e+00, Max=5.1050667755833814e+00, RMS=2.3922068645852592e+00 H(x) for member 4: -air_temperature nobs= 50 Min=-1.1680394925586961e+00, Max=5.2304605643755426e+00, RMS=2.3951684520472365e+00 +air_temperature nobs= 50 Min=-1.1680394925586957e+00, Max=5.2304605643755426e+00, RMS=2.3951684520472365e+00 H(x) for member 5: -air_temperature nobs= 50 Min=-1.3527164115068506e+00, Max=5.1972187664544709e+00, RMS=2.4058700202348420e+00 +air_temperature nobs= 50 Min=-1.3527164115068500e+00, Max=5.1972187664544700e+00, RMS=2.4058700202348415e+00 H(x) for member 6: -air_temperature nobs= 50 Min=-1.2222796970789762e+00, Max=5.1090079401287181e+00, RMS=2.3988762856803163e+00 +air_temperature nobs= 50 Min=-1.2222796970789760e+00, Max=5.1090079401287172e+00, RMS=2.3988762856803159e+00 H(x) for member 7: -air_temperature nobs= 50 Min=-1.2717312950315394e+00, Max=5.0776092731857974e+00, RMS=2.4115399847673236e+00 +air_temperature nobs= 50 Min=-1.2717312950315391e+00, Max=5.0776092731857956e+00, RMS=2.4115399847673231e+00 H(x) for member 8: -air_temperature nobs= 50 Min=-1.2938306497937677e+00, Max=5.0743719762514115e+00, RMS=2.3991798019285939e+00 +air_temperature nobs= 50 Min=-1.2938306497937671e+00, Max=5.0743719762514115e+00, RMS=2.3991798019285939e+00 H(x) for member 9: -air_temperature nobs= 50 Min=-1.5175354879218570e+00, Max=5.1994875108469349e+00, RMS=2.4039977272901951e+00 +air_temperature nobs= 50 Min=-1.5175354879218563e+00, Max=5.1994875108469341e+00, RMS=2.4039977272901947e+00 H(x) for member 10: -air_temperature nobs= 50 Min=-1.1628020655888280e+00, Max=5.0503043448657952e+00, RMS=2.3910141660330040e+00 +air_temperature nobs= 50 Min=-1.1628020655888276e+00, Max=5.0503043448657952e+00, RMS=2.3910141660330044e+00 H(x) ensemble background mean: -air_temperature nobs= 50 Min=-1.2915975134112689e+00, Max=5.1200473787916225e+00, RMS=2.4005666634776146e+00 +air_temperature nobs= 50 Min=-1.2915975134112687e+00, Max=5.1200473787916216e+00, RMS=2.4005666634776146e+00 background y - H(x): -air_temperature nobs= 50 Min=-3.0904601861276593e-01, Max=2.4788732616670828e-01, RMS=1.2718289260451002e-01 +air_temperature nobs= 50 Min=-3.0904607310007171e-01, Max=4.2503948186511664e-01, RMS=1.4489786175292821e-01 Background mean : +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0278301793563930e+01 + + max = 1.9616122198353118e+01 + + mean = 2.4632264424838528e+00 + + stddev = 3.1931643516640591e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0278301793563930e+01 - + max = 1.8678758332278843e+01 - + mean = 2.1985005506833173e+00 - + stddev = 2.8454484711447692e+00 + + max = 1.8678758332278839e+01 + + mean = 2.1985005506833186e+00 + + stddev = 2.8454484711447687e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): + min = -1.0278301793563930e+01 - + max = 1.7482071566813421e+01 + + max = 1.7482071566813417e+01 + mean = 2.1276784917515221e+00 + stddev = 2.5170789792474042e+00 Analysis mean : +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -1.0344294796037092e+01 + + max = 1.9576028887577642e+01 + + mean = 2.4640237350847283e+00 + + stddev = 3.1929130256293661e+00 - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0349656866209049e+01 - + max = 1.8630768031331318e+01 - + mean = 2.2059727047269719e+00 - + stddev = 2.8435166364902660e+00 + + min = -1.0344294796037092e+01 + + max = 1.8638665021503368e+01 + + mean = 2.1992978432841888e+00 + + stddev = 2.8427786834232593e+00 - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -1.0349656866209049e+01 - + max = 1.7434081265865895e+01 - + mean = 2.1351506457951785e+00 - + stddev = 2.5152958790462518e+00 + + min = -1.0344294796037092e+01 + + max = 1.7441978256037945e+01 + + mean = 2.1284757843523976e+00 + + stddev = 2.5136874390805120e+00 +Analysis mean increment : +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = -2.8773181432599149e-01 + + max = 2.7945463804838555e-01 + + mean = 7.9729260086807818e-04 + + stddev = 8.8521587414175915e-02 Analysis mean increment : - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -2.3180136250854516e-01 - + max = 2.4052636772270386e-01 - + mean = 7.4721540436532138e-03 - + stddev = 7.9862988144120822e-02 + + min = -2.8773181432599149e-01 + + max = 2.7945463804838522e-01 + + mean = 7.9729260086810745e-04 + + stddev = 8.8521587414175873e-02 Analysis mean increment : - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = -2.3180136250854513e-01 - + max = 2.4052636772270475e-01 - + mean = 7.4721540436532025e-03 - + stddev = 7.9862988144120781e-02 + + min = -2.8773181432599149e-01 + + max = 2.7945463804838511e-01 + + mean = 7.9729260086811471e-04 + + stddev = 8.8521587414175928e-02 +Forecast variance : +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = 4.6034591952436564e-04 + + max = 3.6294641722342640e-02 + + mean = 9.9273235952286910e-03 + + stddev = 4.6717843128316850e-03 Forecast variance : - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] @@ -252,7 +366,7 @@ Forecast variance : + min = 4.6034591952436623e-04 + max = 3.6294641722342598e-02 + mean = 9.9273235952286910e-03 - + stddev = 4.6717843128316850e-03 + + stddev = 4.6717843128316842e-03 Forecast variance : - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] @@ -263,46 +377,55 @@ Forecast variance : + mean = 9.9273235952286910e-03 + stddev = 4.6717843128316859e-03 Analysis variance : +- Valid time: 2010-01-01T00:00:00Z + Geometry: structured [1855] + Fields: + - air_temperature (2 levels): + + min = 3.5944749374122823e-04 + + max = 2.8343292475071127e-02 + + mean = 7.6767697298578720e-03 + + stddev = 3.6123653670426126e-03 +Analysis variance : - Valid time: 2010-01-01T03:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = 3.5944746468038726e-04 - + max = 2.8343289605209567e-02 - + mean = 7.6767696481328503e-03 - + stddev = 3.6123653099917569e-03 + + min = 3.5944749374122866e-04 + + max = 2.8343292475071100e-02 + + mean = 7.6767697298578711e-03 + + stddev = 3.6123653670426122e-03 Analysis variance : - Valid time: 2010-01-01T06:00:00Z Geometry: structured [1855] Fields: - air_temperature (2 levels): - + min = 3.5944746468038699e-04 - + max = 2.8343289605209574e-02 - + mean = 7.6767696481328538e-03 - + stddev = 3.6123653099917564e-03 + + min = 3.5944749374122850e-04 + + max = 2.8343292475071100e-02 + + mean = 7.6767697298578729e-03 + + stddev = 3.6123653670426109e-03 H(x) for member 1: -air_temperature nobs= 50 Min=-1.2602201729385072e+00, Max=5.1330192642603407e+00, RMS=2.4059920585553405e+00 +air_temperature nobs= 50 Min=-1.2676489436576936e+00, Max=5.1403154011636119e+00, RMS=2.4016023872185390e+00 H(x) for member 2: -air_temperature nobs= 50 Min=-1.2923765506351121e+00, Max=5.1696522313741102e+00, RMS=2.4198597547780278e+00 +air_temperature nobs= 50 Min=-1.2998053359168009e+00, Max=5.1769483742064946e+00, RMS=2.4152747218592490e+00 H(x) for member 3: -air_temperature nobs= 50 Min=-1.3156692336237916e+00, Max=5.1812981343047744e+00, RMS=2.3971151038107212e+00 +air_temperature nobs= 50 Min=-1.3230980201509781e+00, Max=5.1885942719874372e+00, RMS=2.3932075694059058e+00 H(x) for member 4: -air_temperature nobs= 50 Min=-1.1658482096884737e+00, Max=5.2888126217124434e+00, RMS=2.4000087732214772e+00 +air_temperature nobs= 50 Min=-1.1732769775567262e+00, Max=5.2961087626661243e+00, RMS=2.3955793518392379e+00 H(x) for member 5: -air_temperature nobs= 50 Min=-1.3282343466485462e+00, Max=5.2576914294883945e+00, RMS=2.4097572687709898e+00 +air_temperature nobs= 50 Min=-1.3356631285361307e+00, Max=5.2649875622504041e+00, RMS=2.4056895265487062e+00 H(x) for member 6: -air_temperature nobs= 50 Min=-1.2278214532690122e+00, Max=5.1828877336762202e+00, RMS=2.4046086959054307e+00 +air_temperature nobs= 50 Min=-1.2352502330247848e+00, Max=5.1901838650364462e+00, RMS=2.4006474851364619e+00 H(x) for member 7: -air_temperature nobs= 50 Min=-1.2616896557694877e+00, Max=5.1563734484991315e+00, RMS=2.4141614298168408e+00 +air_temperature nobs= 50 Min=-1.2691184345338069e+00, Max=5.1636695868376385e+00, RMS=2.4098034732776470e+00 H(x) for member 8: -air_temperature nobs= 50 Min=-1.2723836626175438e+00, Max=5.1497118137537043e+00, RMS=2.4033651431426590e+00 +air_temperature nobs= 50 Min=-1.2798124511887341e+00, Max=5.1570079563053808e+00, RMS=2.3992823650203872e+00 H(x) for member 9: -air_temperature nobs= 50 Min=-1.4768283899871020e+00, Max=5.2600023269201088e+00, RMS=2.4074343522937194e+00 +air_temperature nobs= 50 Min=-1.4842571615768250e+00, Max=5.2672984732544066e+00, RMS=2.4032589457605251e+00 H(x) for member 10: -air_temperature nobs= 50 Min=-1.1730448204906045e+00, Max=5.1348378588795995e+00, RMS=2.3968069443557942e+00 +air_temperature nobs= 50 Min=-1.1804736014879489e+00, Max=5.1421339968118609e+00, RMS=2.3933314986009901e+00 H(x) ensemble analysis mean: -air_temperature nobs= 50 Min=-1.2774116495668182e+00, Max=5.1914286862868835e+00, RMS=2.4049785937079569e+00 +air_temperature nobs= 50 Min=-1.2848404287630428e+00, Max=5.1987248250519800e+00, RMS=2.4008331772890501e+00 analysis y - H(x): -air_temperature nobs= 50 Min=-3.0186856435971698e-01, Max=1.7205103728066801e-01, RMS=9.0847493643301511e-02 -ombg RMS: 1.2718289260451002e-01 -oman RMS: 9.0847493643301511e-02 +air_temperature nobs= 50 Min=-2.4469863644279261e-01, Max=3.4636203560475831e-01, RMS=1.0628607412638601e-01 +ombg RMS: 1.4489786175292821e-01 +oman RMS: 1.0628607412638601e-01 diff --git a/test/reg/testref/makeobs_4d.ref b/test/reg/testref/makeobs_4d.ref index 4bf16c7..97b23c5 100644 --- a/test/reg/testref/makeobs_4d.ref +++ b/test/reg/testref/makeobs_4d.ref @@ -4,8 +4,8 @@ Initial state: Fields: - air_temperature (2 levels): + min = -1.0287764472676125e+01 - + max = 1.9582387793721100e+01 - + mean = 2.4612641179440802e+00 + + max = 1.9582387793721104e+01 + + mean = 2.4612641179440797e+00 + stddev = 3.1910365906225642e+00 Final state: - Valid time: 2010-01-01T06:00:00Z @@ -13,12 +13,12 @@ Final state: Fields: - air_temperature (2 levels): + min = -1.0287764472676125e+01 - + max = 1.7448337162181407e+01 - + mean = 2.1257161672117437e+00 - + stddev = 2.5109522301847389e+00 + + max = 1.7448337162181403e+01 + + mean = 2.1257161672117442e+00 + + stddev = 2.5109522301847385e+00 H(x): -air_temperature nobs= 50 Min=-1.3311231528550436e+00, Max=5.2747832452627996e+00, RMS=2.4044863832455081e+00 +air_temperature nobs= 50 Min=-1.3311231528550431e+00, Max=5.2747832452627987e+00, RMS=2.4044863832455081e+00 End H(x) Perturbed H(x): -air_temperature nobs= 50 Min=-1.3028111157343716e+00, Max=5.3402882786570558e+00, RMS=2.4030954753449647e+00 +air_temperature nobs= 50 Min=-1.2467055664814741e+00, Max=5.5450869038783122e+00, RMS=2.4209308316571185e+00 End Perturbed H(x)