1) Updates regarding MODIS obs. reader
  -  Added the description of MODIS obs. reader at the subroutine
     (Table of longitude band is moved up)
  -  Changed the hdf-related function to read (vdata vs. SDS dataset)
     : Previous hdf-related functions are cleaned up
  -  Updated and checkced the QC procedure by comparing against Matlab result

2) Updates on subroutine cat_enkf_increment
  - Removing the cat_diagS at ths subroutine
    :Cat_diagS is also removed when the subroutine is called (e.g., ensupd_enkf_update.F90)
  - Computed the snow diagnostics to apply condition when to add/remove snow or no further action needed

: This commit has been made to update the error status located @ subroutine get_ind_obs()
: Basically, not significant change made on the code compared to the previous commit. All the print statements for debugging purposes are commented out
: Minor change made on the MODIS SCA reader(in clsm_ensupd_read_obs.F90) to correctly read the MODIS observation
: Things need to be done
  1) Check the compressed vector ('Observation') lat, long, and tilenum to see whether it correctly calculates or not
  2) Check the tile_grid properties

** Brief updates compared to commits made on 11/08/2021]
- No more errors related to the subroutine get_ind_obs()
- Checked the details of tilenum, Observation vector, and Obs_pred vector
  (meeting slide at 11/19/2021)
- Another floating error occured when calculating the snow heat content (located within the subroutine clsm_ensupd_upd_routine.F90)
  code line where showed error: TSN = min(met_force(n)%Tair-MAPL_TICE,0.0)
- I've tried to calculate the different between met_force%Tair -MAPL_TICE before entering the loop, and i defined it as deduction (array) and call   nth deduction value in the if loop. Still, it revealed a problem.
- For the review purpose, It would be great help if you can look close into the subroutine clsm_ensupd_upd_routine.F90 (from line 4076 ot 4261)

- fixed typos
- fixed indentation
- added/edited comments
- removed redundant initialization of cat_progn_incr
- removed redundant cat_progn_incr%wesn=0 etc in "do nothing" case of snow analysis
- undid partial (!) deletion of an old, commented-out subroutine

…_fraction_assimilation

land_assim: yes; pertubation:0; no mwRTM file

…ensup_read_obs.F90)

…n_assimilation

…n_assimilation

…n_assimilation

Bug fixes:
- Toure et al 2018 equation 1
  - Added missing parentheses in ADD SNOW case
  - Changed 1 to 1. (real number) in REMOVE SNOW case
- Changed declaration of isnow from real to integer
Cleanup:
- Simplified if block logic for computing layer-based analysis
- Replaced loop index n with kk for clarity and searchability
- Removed unused variables

…n_assimilation

…n_assimilation

…n_assimilation

…n_assimilation

…n_assimilation

…n_assimilation

…n_assimilation

…n_assimilation

…n_assimilation

…n_assimilation

…n_assimilation

…ion_assimilation'; manually resolved conflict in clsm_ensupd_enkf_update.F90

…n_assimilation

…n_assimilation

…n_assimilation

…n_assimilation

…R21C (tile_bin2nc4.F90)

…g back to "develop" version

…reverting back to "develop" version

… diffs with develop (clsm_ensupd_enkf_update.F90)

…mments, white-space changes to minimize diffs with "develop" branch (clsm_ensupd_enkf_update.F90, clsm_ensupd_read_obs.F90, clsm_ensupd_upd_routines.F90)

…0, clsm_ensupd_upd_routines.F90)

…m_ensupd_upd_routines.F90, clsm_ensupd_glob_param.F90, clsm_ensupd_read_obs.F90)

…clsm_ensupd_glob_param.F90)

…ions, etc) (clsm_ensupd_read_obs.F90)

…ensupd_read_obs.F90)

…n_assimilation

Added '/' after path and before YYYY dir so path does not have to be specified without trailing '/'
This has no effect if the path does include the trailing '/'

…le/MODIS_snow_cover_fraction_assimilation; manually resolved conflict in tile_bin2nc4.F90

…op" branch (LDASsa_DEFAULT_inputs_ensupd.nml)

…h "develop" branch (LDASsa_DEFAULT_inputs_ensupd.nml)"

This reverts commit 5844ce9.

…aram

…aram.F90)

…n_assimilation

…now_fcst); cleanup of SCF analysis (clsm_ensupd_upd_routines.F90)

…ment calcs (clsm_ensupd_upd_routines.F90)

…_asnow (clsm_ensupd_upd_routines.F90)

… feature/rreichle/MODIS_snow_cover_fraction_assimilation (manually resolved conflicts in components.yaml, clsm_ensupd_enkf_update.F90)

… feature/rreichle/MODIS_snow_cover_fraction_assimilation (manually resolved conflict in tile_bin2nc4.F90)

…upd_upd_routines.F90)

…tant (rather than by using "swe_ratio") (clsm_ensupd_upd_routines.F90)

…_ensupd_upd_routines.F90)

…d_upd_routines.F90)

…) for computation of snow analysis increments (clsm_ensupd_upd_routines.F90)

…h CPP directive for now)