Skip to content

Add obs. and ensemble at initial time in 4D tests #23

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
benjaminmenetrier wants to merge 1 commit into
base: develop
Choose a base branch
from
Open

Add obs. and ensemble at initial time in 4D tests #23

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion test/glb/testinput/4densvar.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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]
Expand Down
2 changes: 1 addition & 1 deletion test/glb/testinput/4densvar_parallel.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
2 changes: 1 addition & 1 deletion test/glb/testinput/4dvar.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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
7 changes: 7 additions & 0 deletions test/glb/testinput/letkf_linear_4d.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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%
Expand Down Expand Up @@ -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
Expand All @@ -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
Expand Down
7 changes: 7 additions & 0 deletions test/glb/testinput/letkf_nonlinear_4d.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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%
Expand Down Expand Up @@ -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
Expand All @@ -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
Expand Down
6 changes: 5 additions & 1 deletion test/glb/testinput/makeobs_4d.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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]
Expand Down
36 changes: 18 additions & 18 deletions test/glb/testref/4densvar.ref
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -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
36 changes: 18 additions & 18 deletions test/glb/testref/4densvar_parallel.ref
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -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
20 changes: 10 additions & 10 deletions test/glb/testref/4dvar.ref
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -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
Loading