Update GDAS INIT utility to optionally use fracoro data

util/gdas_init/run_v17.chgres.sh

@@ -0,0 +1,121 @@
+#!/bin/bash
+
+copy_data()
+{
+
+mkdir -p $SAVEDIR
+cp gfs_ctrl.nc $SAVEDIR
+
+for tile in 'tile1' 'tile2' 'tile3' 'tile4' 'tile5' 'tile6'
+do
+  cp out.atm.${tile}.nc  ${SAVEDIR}/gfs_data.${tile}.nc
+  cp out.sfc.${tile}.nc  ${SAVEDIR}/sfc_data.${tile}.nc 
+done
+}
+
+#---------------------------------------------------------------------------
+# Run chgres using v17 netcdf history data as input.  These history
+# files are part of the OPS v16 gfs/gdas/enkf tarballs, and the
+# v16 retro parallel gfs tarballs.  To run using the v16 retro
+# gdas tarballs (which contain warm restart files), the 
+# run_v16retro.chgres.sh is used.
+#---------------------------------------------------------------------------
+
+set -x
+
+MEMBER=$1
+
+FIX_FV3=$UFS_DIR/fix
+FIX_ORO=${FIX_FV3}/orog
+FIX_AM=${FIX_FV3}/am
+
+WORKDIR=${WORKDIR:-$OUTDIR/work.${MEMBER}}
+
+if [ ${MEMBER} == 'gdas' ] || [ ${MEMBER} == 'gfs' ] ; then
+  CTAR=${CRES_HIRES}
+#---------------------------------------------------------------------------
+# Some gfs tarballs from the v16 retro parallels dont have 'atmos'
+# in their path.  Account for this.
+#---------------------------------------------------------------------------
+  INPUT_DATA_DIR="${EXTRACT_DIR}/${MEMBER}.${yy}${mm}${dd}/${hh}/atmos"
+  if [ ! -d ${INPUT_DATA_DIR} ]; then
+    INPUT_DATA_DIR="${EXTRACT_DIR}/${MEMBER}.${yy}${mm}${dd}/${hh}"
+  fi
+  ATMFILE="${MEMBER}.t${hh}z.atmanl.nc"
+  SFCFILE="${MEMBER}.t${hh}z.sfcanl.nc"
+else  
+  date10=`$NDATE -6 $yy$mm$dd$hh`
+  yy_d=$(echo $date10 | cut -c1-4)
+  mm_d=$(echo $date10 | cut -c5-6)
+  dd_d=$(echo $date10 | cut -c7-8)
+  hh_d=$(echo $date10 | cut -c9-10)
+  CTAR=${CRES_ENKF}
+  INPUT_DATA_DIR="${EXTRACT_DIR}/enkfgdas.${yy_d}${mm_d}${dd_d}/${hh_d}/atmos/mem${MEMBER}"
+  ATMFILE="gdas.t${hh_d}z.atmf006.nc"
+  SFCFILE="gdas.t${hh_d}z.sfcf006.nc"
+fi
+
+if  [ ${CTAR} == 'C48' ] ; then
+   OCNRES='400'
+elif [ ${CTAR} == 'C96' ] ; then
+   OCNRES='100'
+elif [ ${CTAR} == 'C192' ] ; then
+   OCNRES='050'
+elif [ ${CTAR} == 'C384' ] || [ ${CTAR} == 'C768' ] || [ ${CTAR} == 'C1152' ]; then
+   OCNRES='025'
+fi
+
+rm -fr $WORKDIR
+mkdir -p $WORKDIR
+cd $WORKDIR
+
+cat << EOF > fort.41
+
+&config
+ fix_dir_target_grid="${FIX_ORO}/${CTAR}.mx${OCNRES}_frac/fix_sfc"
+ mosaic_file_target_grid="${FIX_ORO}/${CTAR}.mx${OCNRES}_frac/${CTAR}_mosaic.nc"
+ orog_dir_target_grid="${FIX_ORO}/${CTAR}.mx${OCNRES}_frac"
+ orog_files_target_grid="oro_${CTAR}.mx${OCNRES}.tile1.nc","oro_${CTAR}.mx${OCNRES}.tile2.nc","oro_${CTAR}.mx${OCNRES}.tile3.nc","oro_${CTAR}.mx${OCNRES}.tile4.nc","oro_${CTAR}.mx${OCNRES}.tile5.nc","oro_${CTAR}.mx${OCNRES}.tile6.nc"
+ data_dir_input_grid="${INPUT_DATA_DIR}"
+ atm_files_input_grid="${ATMFILE}"
+ sfc_files_input_grid="${SFCFILE}"
+ vcoord_file_target_grid="${FIX_AM}/global_hyblev.l${LEVS}.txt"
+ cycle_mon=$mm
+ cycle_day=$dd
+ cycle_hour=$hh
+ convert_atm=.true.
+ convert_sfc=.true.
+ convert_nst=.true.
+ input_type="gaussian_netcdf"
+ tracers="sphum","liq_wat","o3mr","ice_wat","rainwat","snowwat","graupel"
+ tracers_input="spfh","clwmr","o3mr","icmr","rwmr","snmr","grle"
+/
+EOF
+
+$APRUN $UFS_DIR/exec/chgres_cube
+rc=$?
+
+if [ $rc != 0 ]; then
+  exit $rc
+fi
+
+if [ ${MEMBER} == 'gdas' ] || [ ${MEMBER} == 'gfs' ]; then
+  SAVEDIR=$OUTDIR/${MEMBER}.${yy}${mm}${dd}/${hh}/atmos/INPUT
+  copy_data
+  touch $SAVEDIR/../${MEMBER}.t${hh}z.loginc.txt
+  if [ ${MEMBER} == 'gdas' ]; then
+    cp ${INPUT_DATA_DIR}/*abias* $SAVEDIR/..
+    cp ${INPUT_DATA_DIR}/*radstat $SAVEDIR/..
+  fi
+else  
+  SAVEDIR=$OUTDIR/enkfgdas.${yy}${mm}${dd}/${hh}/atmos/mem${MEMBER}/INPUT
+  copy_data
+  touch $SAVEDIR/../enkfgdas.t${hh}z.loginc.txt
+fi
+
+rm -fr $WORKDIR
+
+set +x
+echo CHGRES COMPLETED FOR MEMBER $MEMBER
+
+exit 0