Check the snow analysis results in the reanalysis application

[jiaruidong2017]

From Joao Souza's email as below:

The results for the 2023 stream are being uploaded to:
https://noaa-reanalyses-pds.s3.amazonaws.com/index.html#analyses/scout_runs/3dvar_coupledreanl_scoutrun_v1.1/

Please, only look at the results until 202301040000 (as of 02/12 10 am Mountain time) since I am overwriting old results from a previous test.

[jiaruidong2017]

After reviewing the results, it appears that no snow observations are being assimilated. The default configurations are set up for GFSv17, and some modifications are needed for the reanalysis application.

0: QC ims_snow totalSnowDepth: 0 passed out of 0 observations.
0: QC ims_snow totalSnowDepth: 0 assigned the background_check diagnostic flag.
0: QC ims_snow totalSnowDepth: 0 assigned the invalid_snowdepth diagnostic flag.
0: QC ims_snow totalSnowDepth: 0 assigned the land_check diagnostic flag.
0: QC ims_snow totalSnowDepth: 0 assigned the landice_check diagnostic flag.
0: QC ims_snow totalSnowDepth: 0 assigned the seaice_check diagnostic flag.
0: QC ims_snow totalSnowDepth: 0 assigned the spatial_thinning diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 passed out of 0 observations.
0: QC sfcsno totalSnowDepth: 0 assigned the background_check diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 assigned the elevation_bkgdiff diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 assigned the invalid_elevation diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 assigned the invalid_snowdepth diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 assigned the land_check diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 assigned the landice_check diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 assigned the missing_elevation diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 assigned the missing_snowdepth diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 assigned the rejectlist diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 assigned the seaice_check diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 assigned the temporal_thinning diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 passed out of 0 observations.
0: QC madis_snow totalSnowDepth: 0 assigned the background_check diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 assigned the elevation_bkgdiff diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 assigned the invalid_elevation diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 assigned the invalid_snowdepth diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 assigned the land_check diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 assigned the landice_check diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 assigned the missing_elevation diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 assigned the missing_snowdepth diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 assigned the rejectlist diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 assigned the seaice_check diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 assigned the temporal_thinning diagnostic flag.

[jiaruidong2017]

Their obs configure is still GFSv17 as:

0: QC ims_snow snowCoverFraction: 47323 rejected as processed but not assimilated.
0: QC ims_snow snowCoverFraction: 0 passed out of 47323 observations.
0: QC ims_snow snowCoverFraction: 0 assigned the background_check diagnostic flag.
0: QC ims_snow snowCoverFraction: 1374 assigned the invalid_elevation diagnostic flag.
0: QC ims_snow snowCoverFraction: 0 assigned the invalid_snowdepth diagnostic flag.
0: QC ims_snow snowCoverFraction: 3167 assigned the land_check diagnostic flag.
0: QC ims_snow snowCoverFraction: 604 assigned the landice_check diagnostic flag.
0: QC ims_snow snowCoverFraction: 2040 assigned the seaice_check diagnostic flag.
0: QC ims_snow snowCoverFraction: 0 assigned the spatial_thinning diagnostic flag.
0: QC ims_snow totalSnowDepth: 390 missing values.
0: QC ims_snow totalSnowDepth: 3716 out of bounds.
0: QC ims_snow totalSnowDepth: 6239 out of domain of use.
0: QC ims_snow totalSnowDepth: 604 black-listed.
0: QC ims_snow totalSnowDepth: 285 removed by thinning.
0: QC ims_snow totalSnowDepth: 10248 rejected by first-guess check.
0: QC ims_snow totalSnowDepth: 25841 passed out of 47323 observations.
0: QC ims_snow totalSnowDepth: 10248 assigned the background_check diagnostic flag.
0: QC ims_snow totalSnowDepth: 1373 assigned the invalid_elevation diagnostic flag.
0: QC ims_snow totalSnowDepth: 3716 assigned the invalid_snowdepth diagnostic flag.
0: QC ims_snow totalSnowDepth: 3035 assigned the land_check diagnostic flag.
0: QC ims_snow totalSnowDepth: 604 assigned the landice_check diagnostic flag.
0: QC ims_snow totalSnowDepth: 1831 assigned the seaice_check diagnostic flag.
0: QC ims_snow totalSnowDepth: 285 assigned the spatial_thinning diagnostic flag.
0: QC sfcsno totalSnowDepth: 149207 missing values.
0: QC sfcsno totalSnowDepth: 25 out of bounds.
0: QC sfcsno totalSnowDepth: 260 out of domain of use.
0: QC sfcsno totalSnowDepth: 21 black-listed.
0: QC sfcsno totalSnowDepth: 5383 removed by thinning.
0: QC sfcsno totalSnowDepth: 60 rejected by first-guess check.
0: QC sfcsno totalSnowDepth: 31 rejected by difference check.
0: QC sfcsno totalSnowDepth: 1099 passed out of 156086 observations.
0: QC sfcsno totalSnowDepth: 60 assigned the background_check diagnostic flag.
0: QC sfcsno totalSnowDepth: 31 assigned the elevation_bkgdiff diagnostic flag.
0: QC sfcsno totalSnowDepth: 34 assigned the invalid_elevation diagnostic flag.
0: QC sfcsno totalSnowDepth: 25 assigned the invalid_snowdepth diagnostic flag.
0: QC sfcsno totalSnowDepth: 79 assigned the land_check diagnostic flag.
0: QC sfcsno totalSnowDepth: 1 assigned the landice_check diagnostic flag.
0: QC sfcsno totalSnowDepth: 51 assigned the missing_elevation diagnostic flag.
0: QC sfcsno totalSnowDepth: 0 assigned the missing_snowdepth diagnostic flag.
0: QC sfcsno totalSnowDepth: 20 assigned the rejectlist diagnostic flag.
0: QC sfcsno totalSnowDepth: 96 assigned the seaice_check diagnostic flag.
0: QC sfcsno totalSnowDepth: 5383 assigned the temporal_thinning diagnostic flag.
0: QC madis_snow totalSnowDepth: 0 passed out of 0 observations.

I suggest they just keep ims and ghcn for their reanalysis runs.

[jiaruidong2017]

The current GW supports both ASCII and NetCDF formats for IMS snow cover data. The filename itself is not important, as the code automatically detects the data format. For convenience, you may name your NetCDF files with the “.nc” extension and place them in your obs directory.