Download - The current status
The current status
• Build the xml and ksh scripts for ENKF and GSI(40km)
• 40 members was divided to 8 groups. GSI: 5~10 minutes/groupENKF: <5 minutesWRF(3hours forecast): 5~10 minutes/group
What I have tested---GSI analysis only?• Configuration in WRF using the 3.1 and 3.2 3.1damp_opt = 3,dampcoef = 0.05moist_adv_dfi_opt = 0
Vs 3.2ishallow = 1,cu_rad_feedback = .true.,prec_acc_dt = 60.damp_opt = 1,dampcoef = 0.02moist_adv_dfi_opt = 2
What I have tested---ENKF analysis using the WRF3.2 confiuration?
• Configuration in ENKFCorrelation length corrlengthnh = 1000 .or. 600corrlengthtr = 1000 .or. 600corrlengthsh = 1000 .or. 600
Time windowsobtimelnh = 10.0 .or 6.0obtimeltr = 10.0 .or 6.0obtimelsh = 10.0 .or 6.0
Inflation coefficientanalpertwt = 0.5 .or. 0.2covinflatenh = 0.2 .or. 0.1covinflatetr = 0.2 .or. 0.1covinflatesh = 0.2 .or. 0.1
Observation Qcsprd_tol = 9.9e31 .or. 1.0
Inflation limitcovinflatemax = 1.e30 .or. 1.5
casename Control variable
corrlength
obtimel Sprd_tol(ob qc)
Input Obs qc
analperwt
covinflat covinflatmax
cenkfgsih3 U,v,t,ps,q 1000 10 9.9e31 no 0.5 0.2 1.e30
cenkfgsih3wa5
U,v,t,ps,q 600 6 1.0 no 0.5 0.1 1.e30
cenkfgsih3wa2
U,v,t,ps,q 600 6 1.0 no 0.2 0.1 1.e30
cenkfgsih3wa2qc
U,v,t,ps,q 600 6 1.0 Yes 0.2 0.1 1.e30
cenkfgsih3adjprsul
U,v,t,p,q 600 6 1.0 Yes 0.8 0.1 1.e30
cenkfgsih3adjprsl
U,v,t,p,q 600 6 1.0 Yes 0.8 0.1 1.5
cenkfgsih3adjprsoffl
U,v,t,p,q 600 6 1.0 Yes 0 0 1.5
Surface pressure, obtype=180, SFCSHP
Surface pressure, obtype=181, ADPSFC
Surface pressure, obtype=187, ADPSFC
Surface pressure, obtype=120, ADPUPA
Surface pressure, obtype=183, ADPSFC,SFCSHP
Obtype=153, GPSIPW
Observation check
Qvapor, obtype=120, ADPUPA Qvapor, obtype=133 AIRCAR Qvapor, obtype=134, AIRCFT
Qvapor, obtype=180, SFCSHP Qvapor, obtype=181, ADPSFC
Qvapor, obtype=183, ADPSFC,SFCSHPQvapor, obtype=187, ADPSFC
temperature, obtype=120, ADPUPAtemperature, obtype=120, RASSDA
temperature, obtype=133 AIRCAR
temperature, obtype=134, AIRCFT temperature, obtype=180, SFCSHP
temperature, obtype=130 AIRCFT
temperature, obtype=131 AIRCFT
temperature, obtype=135 AIRCFT
temperature, obtype=181, ADPSFCtemperature, obtype=187, ADPSFC
temperature, obtype=183, ADPSFC,SFCSHP
U, obtype=220, ADPUPA U, obtype=223, PROFLR U, obtype=224, VADWND
U, obtype=243, SATWND U, obtype=245, SATWND U, obtype=246, SATWND
U, obtype=233, AIRCAR U, obtype=234, AIRCFT U, obtype=235, AIRCFT
U, obtype=230, AIRCFT U, obtype=231, AIRCFT
U, obtype=247, SATWND U, obtype=251, SATWND U, obtype=253, SATWND
U, obtype=254, SATWND U, obtype=280, SFCSHP U, obtype=281, ADPSFC
U, obtype=282, SFCSHP U, obtype=284, ADPSFC,SFCSHP U, obtype=287, ADPSFC
Add constraint in the ‘get_convobs_data’
• Temperature: If abs(Hx-O) > 5.0 Set x_qc(nob) = -999
• U and V: If abs(Hx-O) > 20.0 Set x_qc(nob) = -999• Surface pressure: If abs(Hx-O) > 5.0 Set
x_qc(nob) = -999• Qvapor: If abs(Hx-O) > 2.0 Set x_qc(nob) = -999
• Max and min ensemble mean analysis increment (ug,vg,tempg,qg,psg): at 00UTC -4.556772 4.641273 -5.453091 5.320977 -3.193622 2.471778 -2.7511111E-03 1.9388166E-03 -2.564441 2.609161
• Max and min ensemble mean analysis increment (ug,vg,tempg,qg,psg): at 03UTC -10.69518 7.964006 -6.617310 8.373474 -4.142322 4.528957 -2.0709885E-03 3.2286725E-03 -2.207622 2.785254• Max and min ensemble mean analysis increment (ug,vg,tempg,qg,psg): at 06UTC-14.55259 13.19032 -11.47718 11.15390 -5.021232 6.807072 -2.8613484E-03 4.4340347E-03 -4.340828 3.731035• Max and min ensemble mean analysis increment (ug,vg,tempg,qg,psg): at 09UTC -26.25758 20.79196 -19.61047 14.91415 -8.054179 11.22992 -3.4274021E-03 3.6110494E-03 -3.366595 4.557354• Max and min ensemble mean analysis increment (ug,vg,tempg,qg,psg): at 12UTC-27.53129 25.76617 -26.61810 21.47729 -9.255808 15.93376 -3.0404832E-03 4.1383011E-03 -4.060146 9.327016• Max and min ensemble mean analysis increment (ug,vg,tempg,qg,psg): at 15UTC -35.08931 23.87338 -31.59468 19.87436 -11.49796 12.52381 -3.2191891E-03 3.6776178E-03 -4.257873 10.33907• Max and min ensemble mean analysis increment (ug,vg,tempg,qg,psg): at 18UTC -42.62371 34.48653 -31.65698 29.22342 -14.25782 13.88359 -2.9138320E-03 4.0013785E-03 -10.70951 36.57909
inflation' 2 ' 2
' 2 ' 2
analpertwt *( analpertwt , analpertwt 1r1=
analpertwt *( , analpertwt 1
r = abs(logsigma/lnsigcovinfcutoff)R=r1+taper(r)
) / ( ) 1
) / (
*covinflateanalpertwt 0.8;covinflate 0.1; lnsig
)
covi
b a
b a
x x
x x
nfcutoff 6
taper(r)
Inflation check• With inflation: at 00UTC,May 9no covinflate.dat file available for previous analysis min/max inflation = 1.001976 358.3088 min/max ps inflation = 1.100000 302.6647nanal,psg , max/min increment: 40 91.83405 -43.33380 nanal,ug , max/min increment: 40 155.1773 -116.7216 nanal,vg , max/min increment: 40 100.2846 -121.5881 nanal,qg , max/min increment: 40 7.3264516E-03 -7.2183665E-03 nanal,tempg, max/min increment: 40 53.34357 -62.76953
• Turned off the inflation: at 00UTC,May 9no covinflate.dat file available for previous analysis min/max inflation = 1.000000 1.000000 min/max ps inflation = 1.000000 1.000000nanal,psg , max/min increment: 40 60.55426 -17.94000 nanal,ug , max/min increment: 40 76.97275 -82.57754 nanal,vg , max/min increment: 40 89.04576 -62.84710 nanal,qg , max/min increment: 40 9.3078632E-03 -7.4219871E-03 nanal,tempg, max/min increment: 40 55.46576 -57.95239
P inflation coefficient
Q inflation coefficient
T inflation coefficient
U inflation coefficient
V inflation coefficient
P Spread prior (upper) , post (lower)
Q Spread prior (upper) , post (lower)
T Spread prior (upper) , post (lower)
U Spread prior (upper) , post (lower)
V Spread prior (upper) , post (lower)
(1.0 )1 263.15
233.15 233.15 263.15263.15 233.15
0 233.15
cw ct ci ctq weight q q weight q
TTweight T
T
,
RR straitform weight function
, (1.0 )
1 268.150.05*( 248.15) 248.15 268.15
0 248.15
cw ct ci ctq weight q q weight q
Tweight t T
T
ARPS convective weight function
Thompson, G., R. M. Rasmussen, et al. (2004). "Explicit forecasts of winter precipitation using an improved bulk microphysics scheme. Part I: Description and sensitivity analysis." Monthly Weather Review 132(2): 519-542.
10h forecast
1 *exp( / )
29.263 /meth sfc sfc sfcslp p z KT
K m K
2 *exp( / )
29.2898 / ; (2 )/2meth sfc sfc
sfc sfc
slp p z KT
K m K T T z
3 / 9.2meth sfc sfcslp p z
24 *exp(0.119 0.0013 )meth sfc sfc sfcslp p z z
king
elre
guth
okcn
wato
hint