freekiller
New member
Dear All,
I got strange errors form gen_be_etkf.exe in WRFV4.0.3 and WRF4.1. However, it works fine in WRF3.9.1.1.
The variable of "proj_mn" gets a value of "NaN", indicating that something is wrong. Besides, a error "NCVPT: : NetCDF: Start+count exceeds dimension bound" occurred at the end
Full information is given as follows. T
Thanks in advance.
Jerf
====================================
[1] Initialize information.
Number of ensemble members = 24
Number of prognostic variables = 13
List of prognostic variables = U V W PH T MU QCLOUD QRAIN QVAPOR U10 V10 T2 Q2
naccumt1 = 20
naccumt2 = 20
nstartaccum1 = 1
nstartaccum2 = 1
nout = 1
tainflatinput = 1.00000
rhoinput = 1.00000
infl_fac_file = ./inflation_factor.dat
infl_let_file = ./inflation_letkf.dat
eigen_val_file = ./eigen_value.dat
inno2_val_file = ./innovation_value.dat
proj2_val_file = ./projection_value.dat
infl_fac_TRNK = F
infl_fac_WG03 = F
infl_fac_WG07 = T
infl_fac_BOWL = F
letkf_flag = F
rnd_seed = 2006102800
rnd_nobs = 5000
etkf_erro_max = 20.00000
etkf_erro_min = 0.00001
etkf_inno_max = 20.00000
etkf_inno_min = 0.00001
etkf_erro_flg = T
etkf_inno_flg = T
etkf_wrfda = T
[2] Read observation information.
Number of unfiltered observations = 1193 1
Number of unfiltered observations = 1193 2
Number of unfiltered observations = 1193 3
Number of unfiltered observations = 1193 4
Number of unfiltered observations = 1193 5
Number of unfiltered observations = 1193 6
Number of unfiltered observations = 1193 7
Number of unfiltered observations = 1193 8
Number of unfiltered observations = 1193 9
Number of unfiltered observations = 1193 10
Number of unfiltered observations = 1193 11
Number of unfiltered observations = 1193 12
Number of unfiltered observations = 1193 13
Number of unfiltered observations = 1193 14
Number of unfiltered observations = 1193 15
Number of unfiltered observations = 1193 16
Number of unfiltered observations = 1193 17
Number of unfiltered observations = 1193 18
Number of unfiltered observations = 1193 19
Number of unfiltered observations = 1193 20
Number of unfiltered observations = 1193 21
Number of unfiltered observations = 1193 22
Number of unfiltered observations = 1193 23
Number of unfiltered observations = 1193 24
etkf_wrfda is true: NO ETKF OBS FILTERING
Number of filtered observations = 1193
[3] Set up arrays using input ensemble mean forecast
[4] Extract necessary fields from WRF ensemble forecasts.
[5] Call ETKF:
Enter ETKF with nout = 1
infl_fac_TRNK F
infl_fac_WG03 F
infl_fac_WG07 T
infl_fac_BOWL F
[1] Compute mean(H(xf)) and H(xf)'
[2] Calculate Transform(HZ)(HZ) in Bishop et al. 2001 Eq.(11)
[3] Calculate C and Gamma in Bishop et al. 2001 Eq.(11)
Eigenvalue 1 46.3606521891091
Eigenvalue 2 6.909393081664143E-015
Eigenvalue 3 1.311552345351038E-016
Eigenvalue 4 8.901549475247862E-032
Eigenvalue 5 6.329407068480390E-032
Eigenvalue 6 8.357691306933130E-033
Eigenvalue 7 5.904359660631437E-033
Eigenvalue 8 4.691710728561051E-033
Eigenvalue 9 3.667837110075054E-034
Eigenvalue 10 2.490865889923581E-036
Eigenvalue 11 4.320788832644501E-048
Eigenvalue 12 2.795801049222292E-051
Eigenvalue 13 9.495567745759687E-066
Eigenvalue 14 -9.882498344585418E-049
Eigenvalue 15 -6.023343925557962E-048
Eigenvalue 16 -7.424171949580354E-038
Eigenvalue 17 -8.795444086665627E-034
Eigenvalue 18 -7.568981456494932E-033
Eigenvalue 19 -1.104875173701199E-032
Eigenvalue 20 -1.275584044804543E-032
Eigenvalue 21 -2.356343780986948E-032
Eigenvalue 22 -5.330713169259884E-032
Eigenvalue 23 -4.494450149572382E-019
Eigenvalue 24 -6.563164795896388E-016
Eigenvalue-vector check 1 1.620925615952729E-014
Eigenvalue-vector check 2 8.392592176775793E-015
Eigenvalue-vector check 3 1.017241846312800E-014
Eigenvalue-vector check 4 1.033895191682177E-014
Eigenvalue-vector check 5 9.998078753792328E-015
Eigenvalue-vector check 6 1.027650187168661E-014
Eigenvalue-vector check 7 1.025134838128494E-014
Eigenvalue-vector check 8 1.010649897104088E-014
Eigenvalue-vector check 9 1.051762843484738E-014
Eigenvalue-vector check 10 1.071137536307054E-014
Eigenvalue-vector check 11 1.043316098109658E-014
Eigenvalue-vector check 12 1.059096399323806E-014
Eigenvalue-vector check 13 1.049854647661164E-014
Eigenvalue-vector check 14 1.063680028589553E-014
Eigenvalue-vector check 15 1.057114343463487E-014
Eigenvalue-vector check 16 1.067166645850326E-014
Eigenvalue-vector check 17 1.066161048335346E-014
Eigenvalue-vector check 18 1.065120214249760E-014
Eigenvalue-vector check 19 1.067028410073334E-014
Eigenvalue-vector check 20 1.069283550592104E-014
Eigenvalue-vector check 21 1.065120214249760E-014
Eigenvalue-vector check 22 1.009609063018502E-014
Eigenvalue-vector check 23 9.825473767932635E-015
Eigenvalue-vector check 24 1.053324094613117E-014
[4] Calculate inflation factor
USING WG07 INFLATION FACTOR ALGORITHM
squareinno_now 3626.18604426440
squareinno_mn 3626.18604426440
nobs 1193.00000000000
trace 46.3606521891091
ainflat 52.4838614077133
ainflat_old 1.00000000000000
ainflat_tmp 52.4838614077133
ainflat_mean 52.4838614077133
squareinno_mn 3626.18604426440
proj_mn NaN
nanals-1 23
squareinno_mn 3626.18604426440
nobs 1193
rho NaN
[5] Calculate the grand transformation matrix
ENS
MAX ENS 5.149000131689819E+022
MIN ENS -2128.16205501556
EIGEN/(EIGEN*RHO+1)
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
SUM_EIGNV= 46.3606521891091
SUM_AFTEREIGNV= NaN
AVG_ENS= 1.971353386103933E+015
STD_ENS= 1.015049884464819E+038
[6] Calculate the rescaled ETKF perturbations
ENS
MAX ENS -Infinity
MIN ENS Infinity
AFTER AVG_ENS= NaN
AFTER STD_ENS= NaN
v Variable Ensemble Mean Prior Pert StDv Post. Pert. StDv
1 U 4.1303950***************** NaN
2 V -2.6430138***************** NaN
3 W 0.0001397 3097.3256580 NaN
4 PH 2096.4369864 110570.0798876 NaN
5 T 55.2321618 0.7939410 NaN
6 MU 1847.1413360 1862.0106721 NaN
7 QCLOUD 0.0000000 0.0016617 NaN
8 QRAIN 0.0000000 0.0000005 NaN
9 QVAPOR 0.0008711 0.0016616 NaN
10 U10 0.2617637 90039.0762775 NaN
11 V10 0.1749459 2.6641576 NaN
12 T2 280.9342441 280.9931456 NaN
13 Q2 0.0039820 0.0000472 NaN
[6] Output ETKF analysis ensemble:
NCVPT: : NetCDF: Start+count exceeds dimension bound
I got strange errors form gen_be_etkf.exe in WRFV4.0.3 and WRF4.1. However, it works fine in WRF3.9.1.1.
The variable of "proj_mn" gets a value of "NaN", indicating that something is wrong. Besides, a error "NCVPT: : NetCDF: Start+count exceeds dimension bound" occurred at the end
Full information is given as follows. T
Thanks in advance.
Jerf
====================================
[1] Initialize information.
Number of ensemble members = 24
Number of prognostic variables = 13
List of prognostic variables = U V W PH T MU QCLOUD QRAIN QVAPOR U10 V10 T2 Q2
naccumt1 = 20
naccumt2 = 20
nstartaccum1 = 1
nstartaccum2 = 1
nout = 1
tainflatinput = 1.00000
rhoinput = 1.00000
infl_fac_file = ./inflation_factor.dat
infl_let_file = ./inflation_letkf.dat
eigen_val_file = ./eigen_value.dat
inno2_val_file = ./innovation_value.dat
proj2_val_file = ./projection_value.dat
infl_fac_TRNK = F
infl_fac_WG03 = F
infl_fac_WG07 = T
infl_fac_BOWL = F
letkf_flag = F
rnd_seed = 2006102800
rnd_nobs = 5000
etkf_erro_max = 20.00000
etkf_erro_min = 0.00001
etkf_inno_max = 20.00000
etkf_inno_min = 0.00001
etkf_erro_flg = T
etkf_inno_flg = T
etkf_wrfda = T
[2] Read observation information.
Number of unfiltered observations = 1193 1
Number of unfiltered observations = 1193 2
Number of unfiltered observations = 1193 3
Number of unfiltered observations = 1193 4
Number of unfiltered observations = 1193 5
Number of unfiltered observations = 1193 6
Number of unfiltered observations = 1193 7
Number of unfiltered observations = 1193 8
Number of unfiltered observations = 1193 9
Number of unfiltered observations = 1193 10
Number of unfiltered observations = 1193 11
Number of unfiltered observations = 1193 12
Number of unfiltered observations = 1193 13
Number of unfiltered observations = 1193 14
Number of unfiltered observations = 1193 15
Number of unfiltered observations = 1193 16
Number of unfiltered observations = 1193 17
Number of unfiltered observations = 1193 18
Number of unfiltered observations = 1193 19
Number of unfiltered observations = 1193 20
Number of unfiltered observations = 1193 21
Number of unfiltered observations = 1193 22
Number of unfiltered observations = 1193 23
Number of unfiltered observations = 1193 24
etkf_wrfda is true: NO ETKF OBS FILTERING
Number of filtered observations = 1193
[3] Set up arrays using input ensemble mean forecast
[4] Extract necessary fields from WRF ensemble forecasts.
[5] Call ETKF:
Enter ETKF with nout = 1
infl_fac_TRNK F
infl_fac_WG03 F
infl_fac_WG07 T
infl_fac_BOWL F
[1] Compute mean(H(xf)) and H(xf)'
[2] Calculate Transform(HZ)(HZ) in Bishop et al. 2001 Eq.(11)
[3] Calculate C and Gamma in Bishop et al. 2001 Eq.(11)
Eigenvalue 1 46.3606521891091
Eigenvalue 2 6.909393081664143E-015
Eigenvalue 3 1.311552345351038E-016
Eigenvalue 4 8.901549475247862E-032
Eigenvalue 5 6.329407068480390E-032
Eigenvalue 6 8.357691306933130E-033
Eigenvalue 7 5.904359660631437E-033
Eigenvalue 8 4.691710728561051E-033
Eigenvalue 9 3.667837110075054E-034
Eigenvalue 10 2.490865889923581E-036
Eigenvalue 11 4.320788832644501E-048
Eigenvalue 12 2.795801049222292E-051
Eigenvalue 13 9.495567745759687E-066
Eigenvalue 14 -9.882498344585418E-049
Eigenvalue 15 -6.023343925557962E-048
Eigenvalue 16 -7.424171949580354E-038
Eigenvalue 17 -8.795444086665627E-034
Eigenvalue 18 -7.568981456494932E-033
Eigenvalue 19 -1.104875173701199E-032
Eigenvalue 20 -1.275584044804543E-032
Eigenvalue 21 -2.356343780986948E-032
Eigenvalue 22 -5.330713169259884E-032
Eigenvalue 23 -4.494450149572382E-019
Eigenvalue 24 -6.563164795896388E-016
Eigenvalue-vector check 1 1.620925615952729E-014
Eigenvalue-vector check 2 8.392592176775793E-015
Eigenvalue-vector check 3 1.017241846312800E-014
Eigenvalue-vector check 4 1.033895191682177E-014
Eigenvalue-vector check 5 9.998078753792328E-015
Eigenvalue-vector check 6 1.027650187168661E-014
Eigenvalue-vector check 7 1.025134838128494E-014
Eigenvalue-vector check 8 1.010649897104088E-014
Eigenvalue-vector check 9 1.051762843484738E-014
Eigenvalue-vector check 10 1.071137536307054E-014
Eigenvalue-vector check 11 1.043316098109658E-014
Eigenvalue-vector check 12 1.059096399323806E-014
Eigenvalue-vector check 13 1.049854647661164E-014
Eigenvalue-vector check 14 1.063680028589553E-014
Eigenvalue-vector check 15 1.057114343463487E-014
Eigenvalue-vector check 16 1.067166645850326E-014
Eigenvalue-vector check 17 1.066161048335346E-014
Eigenvalue-vector check 18 1.065120214249760E-014
Eigenvalue-vector check 19 1.067028410073334E-014
Eigenvalue-vector check 20 1.069283550592104E-014
Eigenvalue-vector check 21 1.065120214249760E-014
Eigenvalue-vector check 22 1.009609063018502E-014
Eigenvalue-vector check 23 9.825473767932635E-015
Eigenvalue-vector check 24 1.053324094613117E-014
[4] Calculate inflation factor
USING WG07 INFLATION FACTOR ALGORITHM
squareinno_now 3626.18604426440
squareinno_mn 3626.18604426440
nobs 1193.00000000000
trace 46.3606521891091
ainflat 52.4838614077133
ainflat_old 1.00000000000000
ainflat_tmp 52.4838614077133
ainflat_mean 52.4838614077133
squareinno_mn 3626.18604426440
proj_mn NaN
nanals-1 23
squareinno_mn 3626.18604426440
nobs 1193
rho NaN
[5] Calculate the grand transformation matrix
ENS
MAX ENS 5.149000131689819E+022
MIN ENS -2128.16205501556
EIGEN/(EIGEN*RHO+1)
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
NaN NaN NaN
SUM_EIGNV= 46.3606521891091
SUM_AFTEREIGNV= NaN
AVG_ENS= 1.971353386103933E+015
STD_ENS= 1.015049884464819E+038
[6] Calculate the rescaled ETKF perturbations
ENS
MAX ENS -Infinity
MIN ENS Infinity
AFTER AVG_ENS= NaN
AFTER STD_ENS= NaN
v Variable Ensemble Mean Prior Pert StDv Post. Pert. StDv
1 U 4.1303950***************** NaN
2 V -2.6430138***************** NaN
3 W 0.0001397 3097.3256580 NaN
4 PH 2096.4369864 110570.0798876 NaN
5 T 55.2321618 0.7939410 NaN
6 MU 1847.1413360 1862.0106721 NaN
7 QCLOUD 0.0000000 0.0016617 NaN
8 QRAIN 0.0000000 0.0000005 NaN
9 QVAPOR 0.0008711 0.0016616 NaN
10 U10 0.2617637 90039.0762775 NaN
11 V10 0.1749459 2.6641576 NaN
12 T2 280.9342441 280.9931456 NaN
13 Q2 0.0039820 0.0000472 NaN
[6] Output ETKF analysis ensemble:
NCVPT: : NetCDF: Start+count exceeds dimension bound