Add soil temp QC to IMS pre-processing (NOAA-EMC#1988)

# Description

This PR adds soil temperature QC to IMS pre-processing by removing the
snow obs if the IMS derived snow depth > 0 and model soil temp in layer
1 is > 273.155.

# Issues
Resolves NOAA-EMC#1986 

# Automated CI tests to run in Global Workflow
<!-- Which Global Workflow CI tests are required to adequately test this
PR? -->
- [ ] atm_jjob <!-- JEDI atm single cycle DA !-->
- [ ] C96C48_ufs_hybatmDA <!-- JEDI atm cycled DA !-->
- [ ] C96C48_hybatmsnowDA <!-- JEDI snow cycled DA !-->
- [ ] C96_gcafs_cycled <!-- JEDI aerosol cycled DA !-->
- [ ] C48mx500_3DVarAOWCDA <!-- JEDI low-res marine 3DVar cycled DA !-->
- [ ] C48mx500_hybAOWCDA <!-- JEDI marine hybrid envar cycled DA !-->
- [ ] C96C48_hybatmDA <!-- GSI atm cycled DA !-->

Author: jiaruidong2017

utils/land/gdas_fv3jedi_calc_scf_to_ioda.cc

@@ -654,7 +654,8 @@ void gdasapp::CalcSCFtoIODA::IMSscf::updateIMSsd(fv3jedi::State &state,
   // convert the state to an atlas fieldset
   atlas::FieldSet xBfs;
   state.toFieldSet(xBfs);
-  // Get the snow cover fraction and depth fields from the state
+  // Get the necessary fields from the state
+  auto bkg_stc = atlas::array::make_view<double, 2>(xBfs["stc"]);
   auto bkg_scf = atlas::array::make_view<double, 2>(xBfs["surface_snow_area_fraction"]);
   auto bkg_snd = atlas::array::make_view<double, 2>(xBfs["totalSnowDepth"]);
   const auto bkg_idx =
@@ -675,6 +676,11 @@ void gdasapp::CalcSCFtoIODA::IMSscf::updateIMSsd(fv3jedi::State &state,
          // set the IMS snow depth to a fixed value to QC in JEDI
         this->sndIMS[fv3_j][fv3_i] = -10.0f;
       }
+      if ((bkg_stc(jnode, 0) > 273.155) && (this->sndIMS[fv3_j][fv3_i] > 0)) {
+         // if model soil is too warm to add snow, and the observed snow depth is >0,
+         // set the IMS snow depth to a fixed value to QC in JEDI
+        this->sndIMS[fv3_j][fv3_i] = -20.0f;
+      }
       if (this->sndIMS[fv3_j][fv3_i] > sndIMS_max) {
         // if the IMS snow depth is greater than the maximum, set to nodata
         this->sndIMS[fv3_j][fv3_i] = nodata_float;