Skip to content

Commit d9efa12

Browse files
gmao-rreichlegmao-rreichle
committed
edited LDASsa_SPECIAL_inputs_ensupd.nml for coupled land-atm DAS:
- added comment about purpose of file - made consistent with latest version of LDASsa_DEFAULT_inputs_ensupd.nml - templated resolution-dependency of %scalepath and %scalename - replaced Discover physical path (/gpfsm/dnb..) with logical path
1 parent 95e312c commit d9efa12

File tree

1 file changed

+70
-38
lines changed

1 file changed

+70
-38
lines changed

‎GEOSldas_App/SPECnml/LDASsa_SPECIAL_inputs_ensupd.nml‎

Lines changed: 70 additions & 38 deletions
Original file line numberDiff line numberDiff line change
@@ -1,39 +1,44 @@
11
!
22
! namelist of EnKF inputs for land EnKF update
33
!
4-
! reichle, 28 Jan 2005
5-
! reichle, 13 Jun 2011 - updated for SMOS angles and downscaling ("FOV")
6-
! reichle, 8 Jun 2017 - added "%flistpath" and "%flistname"; updated comments
4+
! CUSTOM version for coupled land-atm DAS with soil moisture analysis based on SMAP Tb obs
75
!
8-
! --------------------------------------------------------------------
6+
! ----------------------------------------------------------------------
97

108
&ens_upd_inputs
119

12-
1310
! ----------------------------------------------------------------------
1411
!
1512
! update type - for details see subroutine cat_enkf_update()
1613
! (note: all 3d updates use compact support)
1714
!
18-
! local = "1d", regional = "3d"
19-
!
20-
! update_type = 0: NO assimilation, NO bias correction
21-
! update_type = 1: 1d soil moisture analysis; sfmc obs
22-
! update_type = 2: 3d soil moisture analysis; sfmc obs
23-
! update_type = 3: 1d Tskin (assim incr NOT applied, use w/ bias corr) analysis; Tskin obs
24-
! update_type = 4: 1d Tskin/ght1 (assim incr applied, use w/ or w/o bias corr) analysis; Tskin obs
25-
! update_type = 5: 1d Tskin/ght1 (assim incr NOT applied, use w/ bias corr) analysis; Tskin obs
26-
! update_type = 6: 1d soil moisture/Tskin/ght(1); TB obs
27-
! update_type = 7: 3d Tskin/ght1 update; Tskin obs
28-
! update_type = 8: 3d soil moisture/Tskin/ght(1); TB obs
29-
! update_type = 9: 1d Tskin/ght1 update; FT obs
30-
! update_type = 10: 3d soil moisture excl catdef/Tskin/ght(1); TB obs
15+
! local = "1d", regional = "3d"
16+
!
17+
! # = no longer supported
18+
!
19+
! update_type | analysis state vector | assimilated obs
20+
! -------------------------------------------------------------------------------------------------
21+
! 0 | NO assimilation, NO bias correction | n/a
22+
! # 1 | 1d soil moisture | sfmc
23+
! # 2 | 3d soil moisture | sfmc
24+
! 3 | 1d Tskin (assim incr NOT applied, use w/ bias corr) | Tskin
25+
! 4 | 1d Tskin/ght1 (assim incr applied, use w/ or w/o bias corr) | Tskin
26+
! 5 | 1d Tskin/ght1 (assim incr NOT applied, use w/ bias corr) | Tskin
27+
! 6 | 1d soil moisture/Tskin/ght(1) | Tb
28+
! 7 | 3d Tskin/ght1 update | Tskin
29+
! 8 | 3d soil moisture/Tskin/ght(1) | Tb
30+
! 9 | 1d Tskin/ght1 update | FT
31+
! 10 | 3d soil moisture/Tskin/ght(1) excl. catdef unless PEATCLSM tile | Tb
32+
! 11 | 1d snow analysis (Toure et al. 2018 empirical gain) | SCF
33+
! 12 | 3d soil moisture/Tskin/ght(1) excl. catdef unless PEATCLSM tile | sfmc/sfds, Tb, SCF
34+
! | & 1d snow analysis (Toure et al. 2018 empirical gain) |
35+
! 13 | 3d soil moisture/Tskin/ght(1) excl. catdef unless PEATCLSM tile | sfmc/sfds, Tb
3136

3237
update_type = 10
3338

39+
out_obslog = .true.
3440
out_ObsFcstAna = .true.
3541
out_smapL4SMaup = .false.
36-
out_obslog = .true.
3742

3843
! ---------------------------------------------------------------------
3944
!
@@ -45,6 +50,17 @@ out_obslog = .true.
4550
xcompact = 1.25 ! [deg] longitude
4651
ycompact = 1.25 ! [deg] latitude
4752

53+
! ---------------------------------------------------------------------
54+
!
55+
! forecast error covariance inflaction factor
56+
!
57+
! - assigns more weight to observations in analysis by inflating forecast error covariance
58+
! - works on std-dev, i.e., var_inflated = var * inflation_fac**2
59+
! - typical values: 1 <= inflation_fac <= 1.5
60+
! - to turn off, set to any negative real number
61+
62+
fcsterr_inflation_fac = -9999.
63+
4864
! ---------------------------------------------------------------------
4965
!
5066
! Definition of measurement species and parameters
@@ -68,15 +84,15 @@ ycompact = 1.25 ! [deg] latitude
6884
! 0 = n/a [eg., in situ obs]
6985
! 1 = ascending
7086
! 2 = descending
71-
! 3 = ascending and descending
87+
! 3 = ascending or descending
7288
! 4 = geostationary
7389
! %pol = polarization
7490
! 0 = n/a [eg., multi-pol. retrieval]
7591
! 1 = horizontal
7692
! 2 = vertical
7793
! 3 = ...
78-
! %N_ang = # angles in species
79-
! %ang = vector of angles
94+
! %N_ang = # satellite viewing angles in species (radiance obs only)
95+
! %ang = vector of satellite viewing angles
8096
! %freq = frequency [Hz]
8197
! %FOV = field-of-view *radius*, see NOTES below
8298
! (if FOV==0. equate obs footprint w/ tile)
@@ -88,9 +104,10 @@ ycompact = 1.25 ! [deg] latitude
88104
! %RTM_ID = ID of radiative transfer model to use for Tb forward modeling
89105
! (subroutine get_obs_pred())
90106
! 0 = none
91-
! 1 = tau-omega model as in De Lannoy et al. 2013 (doi:10.1175/JHM-D-12-092.1)
92-
! 2 = same as 1 but without Pellarin atmospheric corrections
93-
! 3 = ...
107+
! 1 = L-band tau-omega model as in De Lannoy et al. 2013 (doi:10.1175/JHM-D-12-092.1) (SMOS)
108+
! 2 = same as 1 but without Pellarin atm corr (SMAP)
109+
! 3 = same as 1 but with Mironov and SMAP L2_SM pol mixing (SMOS)
110+
! 4 = same as 3 but without Pellarin atm corr (targeted for SMAP L4_SM Version 8)
94111
! %bias_Npar = number of obs bias states tracked per day (integer)
95112
! %bias_trel = e-folding time scale of obs bias memory [s]
96113
! %bias_tcut = cutoff time for confident obs bias estimate [s]
@@ -99,10 +116,12 @@ ycompact = 1.25 ! [deg] latitude
99116
! %units = units (eg., 'K' or 'm3/m3')
100117
! %path = path to measurement files
101118
! %name = name identifier for file containing measurements
119+
! %maskpath = path to obs mask file
120+
! %maskname = filename for obs mask
102121
! %scalepath = path to file(s) with scaling parameters
103122
! %scalename = filename for scaling parameters
104123
! %flistpath = path to file with list of obs file names
105-
! %flistname = name of file with lis of obs file names
124+
! %flistname = name of file with list of obs file names
106125
! %errstd = default obs error std
107126
! %std_normal_max = maximum allowed perturbation (relative to N(0,1))
108127
! %zeromean = enforce zero mean across ensemble
@@ -211,6 +230,9 @@ ycompact = 1.25 ! [deg] latitude
211230
! (see subroutine read_ens_upd_inputs())
212231
!
213232
!
233+
! ------------------------------------------------------------------------
234+
235+
214236
! --------------------------------------------------------------------
215237
!
216238
! SMAP L1C_TB brightness temperature (36 km EASE grid)
@@ -244,10 +266,12 @@ obs_param_nml(31)%nodata = -9999.
244266
obs_param_nml(31)%varname = 'Tb'
245267
obs_param_nml(31)%units = 'K'
246268
obs_param_nml(31)%path = '/discover/nobackup/projects/gmao/smap/SMAP_L4/SMAP/OPS/L1C_TB/'
247-
obs_param_nml(31)%flistpath = '/gpfsm/dnb05/projects/p51/Operations/SPL4SM/SPL4SM_OL7000/import/obs/SMAP/L1C_TB/'
248269
obs_param_nml(31)%name = ''
249-
obs_param_nml(31)%scalepath = '/discover/nobackup/qzhang/LDAS_INPUTS/ZSCOREstats/CF0180x6C_GLOBAL/'
250-
obs_param_nml(31)%scalename = 'ScMO_SMAPc180e24_zscore_stats_2015_p19_2020_p18_hscale_0.00_W_9p_Nmin_20'
270+
obs_param_nml(31)%maskpath = ''
271+
obs_param_nml(31)%maskname = ''
272+
obs_param_nml(31)%scalepath = '/discover/nobackup/qzhang/LDAS_INPUTS/ZSCOREstats/<CFnnnn>x6C_GLOBAL/'
273+
obs_param_nml(31)%scalename = 'ScMO_SMAP<CFnnnn>e24_zscore_stats_2015_p19_2020_p18_hscale_0.00_W_9p_Nmin_20'
274+
obs_param_nml(31)%flistpath = '/discover/nobackup/projects/gmao/smap/Operations/SPL4SM/SPL4SM_OL7000/import/obs/SMAP/L1C_TB/'
251275
obs_param_nml(31)%flistname = 'SMAP_L1C_TB_A_list.txt'
252276
obs_param_nml(31)%errstd = 4.
253277
obs_param_nml(31)%std_normal_max = 2.5
@@ -280,10 +304,12 @@ obs_param_nml(32)%nodata = -9999.
280304
obs_param_nml(32)%varname = 'Tb'
281305
obs_param_nml(32)%units = 'K'
282306
obs_param_nml(32)%path = '/discover/nobackup/projects/gmao/smap/SMAP_L4/SMAP/OPS/L1C_TB/'
283-
obs_param_nml(32)%flistpath = '/gpfsm/dnb05/projects/p51/Operations/SPL4SM/SPL4SM_OL7000/import/obs/SMAP/L1C_TB/'
284307
obs_param_nml(32)%name = ''
285-
obs_param_nml(32)%scalepath = '/discover/nobackup/qzhang/LDAS_INPUTS/ZSCOREstats/CF0180x6C_GLOBAL/'
286-
obs_param_nml(32)%scalename = 'ScMO_SMAPc180e24_zscore_stats_2015_p19_2020_p18_hscale_0.00_W_9p_Nmin_20'
308+
obs_param_nml(32)%maskpath = ''
309+
obs_param_nml(32)%maskname = ''
310+
obs_param_nml(32)%scalepath = '/discover/nobackup/qzhang/LDAS_INPUTS/ZSCOREstats/<CFnnnn>x6C_GLOBAL/'
311+
obs_param_nml(32)%scalename = 'ScMO_SMAP<CFnnnn>e24_zscore_stats_2015_p19_2020_p18_hscale_0.00_W_9p_Nmin_20'
312+
obs_param_nml(32)%flistpath = '/discover/nobackup/projects/gmao/smap/Operations/SPL4SM/SPL4SM_OL7000/import/obs/SMAP/L1C_TB/'
287313
obs_param_nml(32)%flistname = 'SMAP_L1C_TB_D_list.txt'
288314
obs_param_nml(32)%errstd = 4.
289315
obs_param_nml(32)%std_normal_max = 2.5
@@ -316,10 +342,12 @@ obs_param_nml(33)%nodata = -9999.
316342
obs_param_nml(33)%varname = 'Tb'
317343
obs_param_nml(33)%units = 'K'
318344
obs_param_nml(33)%path = '/discover/nobackup/projects/gmao/smap/SMAP_L4/SMAP/OPS/L1C_TB/'
319-
obs_param_nml(33)%flistpath = '/gpfsm/dnb05/projects/p51/Operations/SPL4SM/SPL4SM_OL7000/import/obs/SMAP/L1C_TB/'
320345
obs_param_nml(33)%name = ''
321-
obs_param_nml(33)%scalepath = '/discover/nobackup/qzhang/LDAS_INPUTS/ZSCOREstats/CF0180x6C_GLOBAL/'
322-
obs_param_nml(33)%scalename = 'ScMO_SMAPc180e24_zscore_stats_2015_p19_2020_p18_hscale_0.00_W_9p_Nmin_20'
346+
obs_param_nml(33)%maskpath = ''
347+
obs_param_nml(33)%maskname = ''
348+
obs_param_nml(33)%scalepath = '/discover/nobackup/qzhang/LDAS_INPUTS/ZSCOREstats/<CFnnnn>x6C_GLOBAL/'
349+
obs_param_nml(33)%scalename = 'ScMO_SMAP<CFnnnn>e24_zscore_stats_2015_p19_2020_p18_hscale_0.00_W_9p_Nmin_20'
350+
obs_param_nml(33)%flistpath = '/discover/nobackup/projects/gmao/smap/Operations/SPL4SM/SPL4SM_OL7000/import/obs/SMAP/L1C_TB/'
323351
obs_param_nml(33)%flistname = 'SMAP_L1C_TB_A_list.txt'
324352
obs_param_nml(33)%errstd = 4.
325353
obs_param_nml(33)%std_normal_max = 2.5
@@ -352,10 +380,12 @@ obs_param_nml(34)%nodata = -9999.
352380
obs_param_nml(34)%varname = 'Tb'
353381
obs_param_nml(34)%units = 'K'
354382
obs_param_nml(34)%path = '/discover/nobackup/projects/gmao/smap/SMAP_L4/SMAP/OPS/L1C_TB/'
355-
obs_param_nml(34)%flistpath = '/gpfsm/dnb05/projects/p51/Operations/SPL4SM/SPL4SM_OL7000/import/obs/SMAP/L1C_TB/'
356383
obs_param_nml(34)%name = ''
357-
obs_param_nml(34)%scalepath = '/discover/nobackup/qzhang/LDAS_INPUTS/ZSCOREstats/CF0180x6C_GLOBAL/'
358-
obs_param_nml(34)%scalename = 'SMAPc180e24_zscore_stats_2015_p19_2020_p18_hscale_0.00_W_9p_Nmin_20'
384+
obs_param_nml(34)%maskpath = ''
385+
obs_param_nml(34)%maskname = ''
386+
obs_param_nml(34)%scalepath = '/discover/nobackup/qzhang/LDAS_INPUTS/ZSCOREstats/<CFnnnn>x6C_GLOBAL/'
387+
obs_param_nml(34)%scalename = 'ScMO_SMAP<CFnnnn>e24_zscore_stats_2015_p19_2020_p18_hscale_0.00_W_9p_Nmin_20'
388+
obs_param_nml(34)%flistpath = '/discover/nobackup/projects/gmao/smap/Operations/SPL4SM/SPL4SM_OL7000/import/obs/SMAP/L1C_TB/'
359389
obs_param_nml(34)%flistname = 'SMAP_L1C_TB_D_list.txt'
360390
obs_param_nml(34)%errstd = 4.
361391
obs_param_nml(34)%std_normal_max = 2.5
@@ -365,8 +395,10 @@ obs_param_nml(34)%xcorr = 0.25
365395
obs_param_nml(34)%ycorr = 0.25
366396
obs_param_nml(34)%adapt = 0
367397

398+
368399
! --------------------------------------------------------------------
369400

401+
370402
/
371403

372404
! =========================== EOF =======================================

0 commit comments

Comments
 (0)