From 21fe38162c800f8d0c9632accef19846763850a3 Mon Sep 17 00:00:00 2001
From: rhaesung <rhaesung.kim@noaa.gov>
Date: Thu, 21 Sep 2023 14:09:48 +0000
Subject: [PATCH 1/3] update do_landDA.sh

---
 do_landDA.sh | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/do_landDA.sh b/do_landDA.sh
index 4c1fd03..2a5a97e 100755
--- a/do_landDA.sh
+++ b/do_landDA.sh
@@ -97,12 +97,15 @@ do
 cp ${RSTRDIR}/${FILEDATE}.sfc_data.tile${tile}.nc  ${RSTRDIR}/${FILEDATE}.sfc_data_back.tile${tile}.nc
 done
 fi 
+
 #stage restarts for applying JEDI update (files will get directly updated)
 for tile in 1 2 3 4 5 6 
 do
   ln -fs ${RSTRDIR}/${FILEDATE}.sfc_data.tile${tile}.nc ${JEDIWORKDIR}/${FILEDATE}.sfc_data.tile${tile}.nc
 done
+
 cres_file=${JEDIWORKDIR}/${FILEDATE}.coupler.res
+
 if [[ -e  ${RSTRDIR}/${FILEDATE}.coupler.res ]]; then 
     ln -sf ${RSTRDIR}/${FILEDATE}.coupler.res $cres_file
 else #  if not present, need to create coupler.res for JEDI 

From 31fff859b0b1d2ce8128781b8abffa8d27e916c5 Mon Sep 17 00:00:00 2001
From: rhaesung <rhaesung.kim@noaa.gov>
Date: Tue, 14 Nov 2023 22:17:11 +0000
Subject: [PATCH 2/3] add the fake observations

---
 do_landDA.sh | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/do_landDA.sh b/do_landDA.sh
index 2a5a97e..0d397b4 100755
--- a/do_landDA.sh
+++ b/do_landDA.sh
@@ -137,7 +137,7 @@ do
   # GHCN are time-stamped at 18. If assimilating at 00, need to use previous day's obs, so that
   # obs are within DA window.
   elif [ ${OBS_TYPES[$ii]} == "GHCN" ]; then 
-     obsfile=$OBSDIR/snow_depth/GHCN/data_proc/v3/${YYYY}/ghcn_snwd_ioda_${YYYP}${MP}${DP}.nc
+     obsfile=$OBSDIR/snow_depth/GHCN/data_proc/v3/${YYYY}/fake_ghcn_snwd_ioda_${YYYP}${MP}${DP}.nc
   elif [ ${OBS_TYPES[$ii]} == "SYNTH" ]; then 
      obsfile=$OBSDIR/synthetic_noahmp/IODA.synthetic_gswp_obs.${YYYY}${MM}${DD}${HH}.nc
   else

From a81c8fd5ec4a6c571f76ce3426f2b0777a6d2349 Mon Sep 17 00:00:00 2001
From: rhaesung <rhaesung.kim@noaa.gov>
Date: Sat, 18 Nov 2023 02:28:12 +0000
Subject: [PATCH 3/3] update the GHCN path

---
 do_landDA.sh | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/do_landDA.sh b/do_landDA.sh
index 0d397b4..7fefae9 100755
--- a/do_landDA.sh
+++ b/do_landDA.sh
@@ -136,8 +136,10 @@ do
      obsfile=$OBSDIR/snow_depth/GTS/data_proc/${YYYY}${MM}/adpsfc_snow_${YYYY}${MM}${DD}${HH}.nc4
   # GHCN are time-stamped at 18. If assimilating at 00, need to use previous day's obs, so that
   # obs are within DA window.
-  elif [ ${OBS_TYPES[$ii]} == "GHCN" ]; then 
-     obsfile=$OBSDIR/snow_depth/GHCN/data_proc/v3/${YYYY}/fake_ghcn_snwd_ioda_${YYYP}${MP}${DP}.nc
+  elif [ $atmos_forc == "era5" ] && [ ${OBS_TYPES[$ii]} == "GHCN" ]; then
+     obsfile=$OBSDIR/snow_depth/GHCN/data_proc/v3/${YYYY}/ghcn_snwd_ioda_${YYYP}${MP}${DP}.nc
+  elif [ $atmos_forc == "gswp3" ] && [ ${OBS_TYPES[$ii]} == "GHCN" ]; then
+     obsfile=$OBSDIR/snow_depth/GHCN/data_proc/v3/${YYYY}/fake_ghcn_snwd_ioda_${YYYP}${MP}${DP}.nc  
   elif [ ${OBS_TYPES[$ii]} == "SYNTH" ]; then 
      obsfile=$OBSDIR/synthetic_noahmp/IODA.synthetic_gswp_obs.${YYYY}${MM}${DD}${HH}.nc
   else