#!/bin/ksh -x

########## config.anal ##########
# Analysis specific

echo "BEGIN: config.anal"

# Get task specific resources
. $EXPDIR/config.resources anal

if [ $DONST = "YES" ]; then
    . $EXPDIR/config.nsst
fi

if [[ "$CDUMP" = "gfs" ]] ; then
    export USE_RADSTAT="NO" # This can be only used when bias correction is not-zero.
    export GENDIAG="NO"
    export SETUP='diag_rad=.false.,diag_pcp=.false.,diag_conv=.false.,diag_ozone=.false.,write_diag(3)=.false.,niter(2)=100,'
    export DIAG_TARBALL="NO"
fi

export npe_gsi=$npe_anal

if [[ "$CDUMP" == "gfs" ]] ; then
   export npe_gsi=$npe_anal_gfs
   export nth_anal=$nth_anal_gfs
fi

# Set parameters specific to L127
if [ $LEVS = "128" ]; then
    export GRIDOPTS="nlayers(63)=1,nlayers(64)=1,"
    export SETUP="gpstop=55,nsig_ext=45,$SETUP"
fi

# Set namelist option for LETKF
export lobsdiag_forenkf=".false."  # anal does not need to write out jacobians
                                   # set to .true. in config.eobs and config.eupd

if [ $OUTPUT_FILE = "nemsio" ]; then
    export DO_CALC_INCREMENT="YES"
    export DO_CALC_ANALYSIS="NO"
fi

# Do not process the following datasets
export GSNDBF=${GSNDBF:-/dev/null}
export AMSREBF=${AMSREBF:-/dev/null}
export SSMITBF=${SSMITBF:-/dev/null}
export AMSR2BF=${AMSR2BF:-/dev/null}


# Set default values for info files and observation error
#  NOTE:  Remember to set PRVT in config.prep as OBERROR is set below
export CONVINFO=$FIXgsi/global_convinfo.txt
export OZINFO=$FIXgsi/global_ozinfo.txt
export SATINFO=$FIXgsi/global_satinfo.txt
export OBERROR=$FIXgsi/prepobs_errtable.global


# Use experimental dumps in EMC GFS v16 parallels
export COMPONENT=atmos
if [[ $RUN_ENVIR == "emc" ]]; then
    export ABIBF="/dev/null"
    if [[ "$CDATE" -ge "2019022800" ]] ; then
	export ABIBF="$DMPDIR/${CDUMP}x.${PDY}/${cyc}/${COMPONENT}/${CDUMP}.t${cyc}z.gsrcsr.tm00.bufr_d"
	if [[ "$CDATE" -ge "2019111000" && "$CDATE" -le "2020052612" ]]; then
	    export ABIBF="$DMPDIR/${CDUMP}y.${PDY}/${cyc}/${COMPONENT}/${CDUMP}.t${cyc}z.gsrcsr.tm00.bufr_d"
	fi
    fi

    export AHIBF="/dev/null"
    if [[ "$CDATE" -ge "2019042300" ]]; then
	export AHIBF="$DMPDIR/${CDUMP}x.${PDY}/${cyc}/${COMPONENT}/${CDUMP}.t${cyc}z.ahicsr.tm00.bufr_d"
    fi

    export HDOB=$DMPDIR/${CDUMP}x.${PDY}/${cyc}/${COMPONENT}/${CDUMP}.t${cyc}z.hdob.tm00.bufr_d

#   Use dumps from NCO GFS v16 parallel
    if [[ "$CDATE" -ge "2020103012" ]]; then
	export ABIBF=""
	export AHIBF=""
	export HDOB=""
    fi

#   Set info files and prepobs.errtable.global for GFS v16 retrospective parallels
    if [[ "$CDATE" -ge "2019021900" && "$CDATE" -lt "2019110706" ]]; then
	export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2019021900
	export OBERROR=$FIXgsi/gfsv16_historical/prepobs_errtable.global.2019021900
    fi

#   Place GOES-15 AMVs in monitor, assimilate GOES-17 AMVs, assimilate KOMPSAT-5 gps
    if [[ "$CDATE" -ge "2019110706" && "$CDATE" -lt "2020040718" ]]; then
	export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2019110706
	export OBERROR=$FIXgsi/gfsv16_historical/prepobs_errtable.global.2019110706
    fi

#   Assimilate 135 (T) & 235 (uv) Canadian AMDAR observations
    if [[ "$CDATE" -ge "2020040718" && "$CDATE" -lt "2020052612" ]]; then
	export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2020040718
	export OBERROR=$FIXgsi/gfsv16_historical/prepobs_errtable.global.2020040718
    fi

#   Assimilate COSMIC-2
    if [[ "$CDATE" -ge "2020052612" && "$CDATE" -lt "2020082412" ]]; then
        export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2020052612
	export OBERROR=$FIXgsi/gfsv16_historical/prepobs_errtable.global.2020040718
    fi

#   Assimilate HDOB
    if [[ "$CDATE" -ge "2020082412" && "$CDATE" -lt "2020091612" ]]; then
        export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2020082412
    fi

#   Assimilate Metop-C GNSSRO
    if [[ "$CDATE" -ge "2020091612" && "$CDATE" -lt "2021031712" ]]; then
        export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2020091612
    fi

#   Assimilate DO-2 GeoOptics
    if [[ "$CDATE" -ge "2021031712" && "$CDATE" -lt "2021091612" ]]; then
        export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2021031712
    fi

#   Assimilate DO-3 Spire
    if [[ "$CDATE" -ge "2021091612" && "$CDATE" -lt "2022031612" ]]; then
        export CONVINFO=$FIXgsi/gfsv16_historical/gfsv16.3/global_convinfo.txt.2021091612
    fi


#   NOTE:
#   As of 2022021812, gfsv16_historical/global_convinfo.txt.2022031612 is
#   identical to ../global_convinfo.txt.  Thus, the logic below is not
#   needed at this time.
#   Assimilate DO-4 (Spire and GeoOptics) and turn off uv 224 VADWND
#   if [[ "$CDATE" -ge "2022031612" && "$CDATE" -lt "YYYYMMDDHH" ]]; then
#     export CONVINFO=$FIXgsi/gfsv16_historical/global_convinfo.txt.2022031612
#   fi


#   Turn off assmilation of OMPS during period of bad data
    if [[ "$CDATE" -ge "2020011600" && "$CDATE" -lt "2020011806" ]]; then
	export OZINFO=$FIXgsi/gfsv16_historical/global_ozinfo.txt.2020011600
    fi


#   Set satinfo for start of GFS v16 parallels
    if [[ "$CDATE" -ge "2019021900" && "$CDATE" -lt "2019110706" ]]; then
	export SATINFO=$FIXgsi/gfsv16_historical/global_satinfo.txt.2019021900
    fi

#   Turn on assimilation of Metop-C AMSUA and MHS
    if [[ "$CDATE" -ge "2019110706" && "$CDATE" -lt "2020022012" ]]; then
	export SATINFO=$FIXgsi/gfsv16_historical/global_satinfo.txt.2019110706
    fi

#   Turn off assimilation of Metop-A MHS
    if [[ "$CDATE" -ge "2020022012" && "$CDATE" -lt "2021052118" ]]; then
	export SATINFO=$FIXgsi/gfsv16_historical/global_satinfo.txt.2020022012
    fi

#   Turn off assimilation of S-NPP CrIS
    if [[ "$CDATE" -ge "2021052118" && "$CDATE" -lt "2021092206" ]]; then
	export SATINFO=$FIXgsi/gfsv16_historical/global_satinfo.txt.2021052118
    fi

#   Turn off assimilation of MetOp-A IASI
    if [[ "$CDATE" -ge "2021092206" && "$CDATE" -lt "2021102612" ]]; then
	export SATINFO=$FIXgsi/gfsv16_historical/gfsv16.3/global_satinfo.txt.2021092206
    fi

#   NOTE:  
#   As of 2021110312, gfsv16_historical/global_satinfo.txt.2021102612 is
#   identical to ../global_satinfo.txt.  Thus, the logic below is not
#   needed at this time
#
#   Turn on assimilation of MetOp-C IASI
#   if [[ "$CDATE" -ge "2021102612" && "$CDATE" -lt "YYYYMMDDHH" ]]; then
#       export SATINFO=$FIXgsi/gfsv16_historical/global_satinfo.txt.2021102612
#   fi

fi

echo "END: config.anal"