Skip to content
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

add ensemble sdv to diag ( selected variables ) #635

Merged
merged 12 commits into from
Mar 23, 2023
Merged

add ensemble sdv to diag ( selected variables ) #635

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
12 commits
Select commit Hold shift + click to select a range
1f7004b
add ensemble sdv to diag ( selected variables )
saraqzhang Mar 8, 2023
3c360e8
(1) make "std" in variable names for ens. standard deviatio…
saraqzhang Mar 15, 2023
a8df7fc
Merge branch 'develop' into feature/saraqzhang/enssdvfordiag
gmao-rreichle Mar 15, 2023
5d08507
additional changes in variables names for ENSSTD output
gmao-rreichle Mar 16, 2023
4ef7473
added new operators for cat_diagS type; protected against division by…
gmao-rreichle Mar 16, 2023
3a16724
additional protections against division by zero for calc of ens std
gmao-rreichle Mar 16, 2023
435066f
fixed new multiplication operators for cat_diagS (order of args)
gmao-rreichle Mar 16, 2023
37dbf00
ensure non-negative ens std output
gmao-rreichle Mar 17, 2023
42f5e71
facilitate HISTORY output of ensstd even when ensavg is not written o…
gmao-rreichle Mar 20, 2023
f6a1b71
undoing previous commit which failed to facilitate ensstd HISTORY out…
gmao-rreichle Mar 22, 2023
d40091f
Merge branch 'develop' into feature/saraqzhang/enssdvfordiag
gmao-rreichle Mar 22, 2023
89a8cc5
added rudimentary documentation of ENSSTD output into GEOSldas_HIST.rc
gmao-rreichle Mar 22, 2023
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
68 changes: 46 additions & 22 deletions src/Applications/LDAS_App/GEOSldas_HIST.rc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -407,29 +407,43 @@ COLLECTIONS:

# For lndfcstana, *.frequency and *.ref_time must be consistent with the LDAS.rc resource
# parameters LANDASSIM_DT and LANDASSIM_T0.
# Output of the ensemble std-dev (ENSSTD) requires simultaneous output of the ensemble mean. If the
# ensemble mean is not written out for a given field, that field's ENSSTD output will be MAPL_UNDEF.

inst3_1d_lndfcstana_Nt.descr: 'Tile-space,3-Hourly,Instantaneous,Single-Level,Assimilation,Ensemble-Average Land Forecast and Analysis Diagnostics',
inst3_1d_lndfcstana_Nt.descr: 'Tile-space,3-Hourly,Instantaneous,Single-Level,Assimilation,Ensemble Land Forecast and Analysis Diagnostics',
inst3_1d_lndfcstana_Nt.nbits: 12,
inst3_1d_lndfcstana_Nt.template: '%y4%m2%d2_%h2%n2z.bin',
inst3_1d_lndfcstana_Nt.mode: 'instantaneous',
inst3_1d_lndfcstana_Nt.frequency: 030000,
inst3_1d_lndfcstana_Nt.ref_time: 000000,
inst3_1d_lndfcstana_Nt.fields: 'WCSF' , 'ENSAVG' , 'SFMC_FCST' ,
'WCRZ' , 'ENSAVG' , 'RZMC_FCST' ,
'WCPR' , 'ENSAVG' , 'PRMC_FCST' ,
'TPSURF' , 'ENSAVG' , 'TSURF_FCST' ,
'TSOIL1TILE' , 'ENSAVG' , 'TSOIL1_FCST' ,
'WCSF_ANA' , 'LANDASSIM' , 'SFMC_ANA' ,
'WCRZ_ANA' , 'LANDASSIM' , 'RZMC_ANA' ,
'WCPR_ANA' , 'LANDASSIM' , 'PRMC_ANA' ,
'TPSURF_ANA' , 'LANDASSIM' , 'TSURF_ANA' ,
'TSOIL1_ANA' , 'LANDASSIM' , 'TSOIL1_ANA' ,
inst3_1d_lndfcstana_Nt.fields: 'WCSF' , 'ENSAVG' , 'SFMC_FCST' ,
'WCRZ' , 'ENSAVG' , 'RZMC_FCST' ,
'WCPR' , 'ENSAVG' , 'PRMC_FCST' ,
'TPSURF' , 'ENSAVG' , 'TSURF_FCST' ,
'TSOIL1TILE' , 'ENSAVG' , 'TSOIL1_FCST' ,
'WCSF_ENSSTD' , 'ENSAVG' , 'SFMC_FCST_ENSSTD' ,
'WCRZ_ENSSTD' , 'ENSAVG' , 'RZMC_FCST_ENSSTD' ,
'WCPR_ENSSTD' , 'ENSAVG' , 'PRMC_FCST_ENSSTD' ,
'TPSURF_ENSSTD' , 'ENSAVG' , 'TSURF_FCST_ENSSTD' ,
'TSOIL1TILE_ENSSTD' , 'ENSAVG' , 'TSOIL1_FCST_ENSSTD' ,
'WCSF_ANA' , 'LANDASSIM' , 'SFMC_ANA' ,
'WCRZ_ANA' , 'LANDASSIM' , 'RZMC_ANA' ,
'WCPR_ANA' , 'LANDASSIM' , 'PRMC_ANA' ,
'TPSURF_ANA' , 'LANDASSIM' , 'TSURF_ANA' ,
'TSOIL1_ANA' , 'LANDASSIM' , 'TSOIL1_ANA' ,
'WCSF_ANA_ENSSTD' , 'LANDASSIM' , 'SFMC_ANA_ENSSTD' ,
'WCRZ_ANA_ENSSTD' , 'LANDASSIM' , 'RZMC_ANA_ENSSTD' ,
'WCPR_ANA_ENSSTD' , 'LANDASSIM' , 'PRMC_ANA_ENSSTD' ,
'TPSURF_ANA_ENSSTD' , 'LANDASSIM' , 'TSURF_ANA_ENSSTD' ,
'TSOIL1_ANA_ENSSTD' , 'LANDASSIM' , 'TSOIL1_ANA_ENSSTD'
::

# For lndfcstana, *.frequency and *.ref_time must be consistent with the LDAS.rc resource
# parameters LANDASSIM_DT and LANDASSIM_T0.
# Output of the ensemble std-dev (ENSSTD) requires simultaneous output of the ensemble mean. If the
# ensemble mean is not written out for a given field, that field's ENSSTD output will be MAPL_UNDEF.

inst3_2d_lndfcstana_Nx.descr: '2d,3-Hourly,Instantaneous,Single-Level,Assimilation,Ensemble-Average Land Forecast and Analysis Diagnostics',
inst3_2d_lndfcstana_Nx.descr: '2d,3-Hourly,Instantaneous,Single-Level,Assimilation,Ensemble Land Forecast and Analysis Diagnostics',
inst3_2d_lndfcstana_Nx.nbits: 12,
inst3_2d_lndfcstana_Nx.template: '%y4%m2%d2_%h2%n2z.nc4',
inst3_2d_lndfcstana_Nx.archive: '%c/Y%y4',
Expand All @@ -441,16 +455,26 @@ COLLECTIONS:
inst3_2d_lndfcstana_Nx.regrid_name: 'GRIDNAME',
inst3_2d_lndfcstana_Nx.grid_label: PC720x361-DC,
inst3_2d_lndfcstana_Nx.deflate: 2,
inst3_2d_lndfcstana_Nx.fields: 'WCSF' , 'ENSAVG' , 'SFMC_FCST' ,
'WCRZ' , 'ENSAVG' , 'RZMC_FCST' ,
'WCPR' , 'ENSAVG' , 'PRMC_FCST' ,
'TPSURF' , 'ENSAVG' , 'TSURF_FCST' ,
'TSOIL1TILE' , 'ENSAVG' , 'TSOIL1_FCST' ,
'WCSF_ANA' , 'LANDASSIM' , 'SFMC_ANA' ,
'WCRZ_ANA' , 'LANDASSIM' , 'RZMC_ANA' ,
'WCPR_ANA' , 'LANDASSIM' , 'PRMC_ANA' ,
'TPSURF_ANA' , 'LANDASSIM' , 'TSURF_ANA' ,
'TSOIL1_ANA' , 'LANDASSIM' , 'TSOIL1_ANA' ,
inst3_2d_lndfcstana_Nx.fields: 'WCSF' , 'ENSAVG' , 'SFMC_FCST' ,
'WCRZ' , 'ENSAVG' , 'RZMC_FCST' ,
'WCPR' , 'ENSAVG' , 'PRMC_FCST' ,
'TPSURF' , 'ENSAVG' , 'TSURF_FCST' ,
'TSOIL1TILE' , 'ENSAVG' , 'TSOIL1_FCST' ,
'WCSF_ENSSTD' , 'ENSAVG' , 'SFMC_FCST_ENSSTD' ,
'WCRZ_ENSSTD' , 'ENSAVG' , 'RZMC_FCST_ENSSTD' ,
'WCPR_ENSSTD' , 'ENSAVG' , 'PRMC_FCST_ENSSTD' ,
'TPSURF_ENSSTD' , 'ENSAVG' , 'TSURF_FCST_ENSSTD' ,
'TSOIL1TILE_ENSSTD' , 'ENSAVG' , 'TSOIL1_FCST_ENSSTD' ,
'WCSF_ANA' , 'LANDASSIM' , 'SFMC_ANA' ,
'WCRZ_ANA' , 'LANDASSIM' , 'RZMC_ANA' ,
'WCPR_ANA' , 'LANDASSIM' , 'PRMC_ANA' ,
'TPSURF_ANA' , 'LANDASSIM' , 'TSURF_ANA' ,
'TSOIL1_ANA' , 'LANDASSIM' , 'TSOIL1_ANA' ,
'WCSF_ANA_ENSSTD' , 'LANDASSIM' , 'SFMC_ANA_ENSSTD' ,
'WCRZ_ANA_ENSSTD' , 'LANDASSIM' , 'RZMC_ANA_ENSSTD' ,
'WCPR_ANA_ENSSTD' , 'LANDASSIM' , 'PRMC_ANA_ENSSTD' ,
'TPSURF_ANA_ENSSTD' , 'LANDASSIM' , 'TSURF_ANA_ENSSTD' ,
'TSOIL1_ANA_ENSSTD' , 'LANDASSIM' , 'TSOIL1_ANA_ENSSTD'
::

# ========================== EOF ==============================================================
33 changes: 23 additions & 10 deletions src/Applications/LDAS_App/tile_bin2nc4.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -399,16 +399,29 @@ FUNCTION getAttribute (SHORT_NAME, LNAME, UNT) result (str_atr)

! land assimilation forecast and analysis for Catchment model diagnostics

case ('SFMC_FCST'); LONG_NAME = 'soil_moisture_surface_forecast'; UNITS = 'm3 m-3'
case ('RZMC_FCST'); LONG_NAME = 'soil_moisture_rootzone_forecast'; UNITS = 'm3 m-3'
case ('PRMC_FCST'); LONG_NAME = 'soil_moisture_profile_forecast'; UNITS = 'm3 m-3'
case ('TSURF_FCST'); LONG_NAME = 'ave_catchment_temp_incl_snw_forecast'; UNITS = 'K'
case ('TSOIL1_FCST'); LONG_NAME = 'soil_temperatures_layer_1_forecast'; UNITS = 'K'
case ('SFMC_ANA'); LONG_NAME = 'soil_moisture_surface_analysis'; UNITS = 'm3 m-3'
case ('RZMC_ANA'); LONG_NAME = 'soil_moisture_rootzone_analysis'; UNITS = 'm3 m-3'
case ('PRMC_ANA'); LONG_NAME = 'soil_moisture_profile_analysis'; UNITS = 'm3 m-3'
case ('TSURF_ANA'); LONG_NAME = 'ave_catchment_temp_incl_snw_analysis'; UNITS = 'K'
case ('TSOIL1_ANA'); LONG_NAME = 'soil_temperatures_layer_1_analysis'; UNITS = 'K'
case ('SFMC_FCST'); LONG_NAME = 'soil_moisture_surface_forecast'; UNITS = 'm3 m-3'
case ('RZMC_FCST'); LONG_NAME = 'soil_moisture_rootzone_forecast'; UNITS = 'm3 m-3'
case ('PRMC_FCST'); LONG_NAME = 'soil_moisture_profile_forecast'; UNITS = 'm3 m-3'
case ('TSURF_FCST'); LONG_NAME = 'ave_catchment_temp_incl_snw_forecast'; UNITS = 'K'
case ('TSOIL1_FCST'); LONG_NAME = 'soil_temperatures_layer_1_forecast'; UNITS = 'K'

case ('SFMC_FCST_ENSSTD'); LONG_NAME = 'soil_moisture_surface_forecast_ensstd'; UNITS = 'm3 m-3'
case ('RZMC_FCST_ENSSTD'); LONG_NAME = 'soil_moisture_rootzone_forecast_ensstd'; UNITS = 'm3 m-3'
case ('PRMC_FCST_ENSSTD'); LONG_NAME = 'soil_moisture_profile_forecast_ensstd'; UNITS = 'm3 m-3'
case ('TSURF_FCST_ENSSTD'); LONG_NAME = 'ave_catchment_temp_incl_snw_forecast_ensstd'; UNITS = 'K'
case ('TSOIL1_FCST_ENSSTD'); LONG_NAME = 'soil_temperatures_layer_1_forecast_ensstd'; UNITS = 'K'

case ('SFMC_ANA'); LONG_NAME = 'soil_moisture_surface_analysis'; UNITS = 'm3 m-3'
case ('RZMC_ANA'); LONG_NAME = 'soil_moisture_rootzone_analysis'; UNITS = 'm3 m-3'
case ('PRMC_ANA'); LONG_NAME = 'soil_moisture_profile_analysis'; UNITS = 'm3 m-3'
case ('TSURF_ANA'); LONG_NAME = 'ave_catchment_temp_incl_snw_analysis'; UNITS = 'K'
case ('TSOIL1_ANA'); LONG_NAME = 'soil_temperatures_layer_1_analysis'; UNITS = 'K'

case ('SFMC_ANA_ENSSTD'); LONG_NAME = 'soil_moisture_surface_analysis_ensstd'; UNITS = 'm3 m-3'
case ('RZMC_ANA_ENSSTD'); LONG_NAME = 'soil_moisture_rootzone_analysis_ensstd'; UNITS = 'm3 m-3'
case ('PRMC_ANA_ENSSTD'); LONG_NAME = 'soil_moisture_profile_analysis_ensstd'; UNITS = 'm3 m-3'
case ('TSURF_ANA_ENSSTD'); LONG_NAME = 'ave_catchment_temp_incl_snw_analysis_ensstd'; UNITS = 'K'
case ('TSOIL1_ANA_ENSSTD'); LONG_NAME = 'soil_temperatures_layer_1_analysis_ensstd'; UNITS = 'K'

! other land assimilation fields

Expand Down
Loading