4. NEXRAD Experiment - delippi/PhD-globalRadarOSSE GitHub Wiki

THE NEXRAD EXPERIMENT

You will need to run the cycled DA workflow.
You will also need to use the initial conditions generated for the NODA experiment and create initial conditions for the C384 ensemble.

1. Ensure that reflectivity is in the diag table.

See "add reflectivity to the diag table" in the setup FV3GFS wiki

!set up the directories (only do this the very first time)
cd /gpfs/hps3/emc/meso/save/Donald.E.Lippi/PhD-globalRadarOSSE
ksh run_NEXRAD_288hrrestart.ksh

!download ICs
mkdir -p /gpfs/hps2/ptmp/Donald.E.Lippi/fv3gfs_dl2rw/2018092300
cd /gpfs/hps2/ptmp/Donald.E.Lippi/fv3gfs_dl2rw/2018092300
htar -xvf /NCEPDEV/emc-meso/5year/Donald.E.Lippi/rw_FV3GFS/FV3ICS/NODA-2018091100-2018092300_ICS288.tar
cd /gpfs/hps3/emc/meso/save/Donald.E.Lippi/PhD-globalRadarOSSE.Surge/get_ics
ksh efcspull.ksh
ksh get_ics.ksh.nexrad

cd /gpfs/hps2/ptmp/Donald.E.Lippi/fv3gfs_dl2rw/2018092300/NEXRAD-2018092300-2018100700
mv gdas.2/NODA-2018091100-2018092300/gfs.20180911 .
rm -rf gdas.2/
cd /gpfs/hps2/ptmp/Donald.E.Lippi/fv3gfs_dl2rw/2018092300/NEXRAD-2018092300-2018100700/gfs.20180911/00
rm -f gfs.t00z.*
mv RERUN_RESTART/ RESTART

cd /gpfs/hps2/ptmp/Donald.E.Lippi/fv3gfs_dl2rw/2018092300/NEXRAD-2018092300-2018100700/
mv enkfgdas.20180911/ enkfgdas.20180922
cd enkfgdas.20180922
mv 00 18
mv gfs.20180911 gdas.20180922
cd gdas.20180922
mv 00 18

!modify config.base
cd /gpfs/hps3/emc/meso/noscrub/Donald.E.Lippi/fv3gfs_dl2rw/2018092300/NEXRAD-2018092300-2018100700

sed -i "s/HPSS_PROJECT=emc-global/HPSS_PROJECT=emc-meso/g" config.base
sed -i "s/export FHCYC=1/export FHCYC=00/g" config.base
sed -i "s/export FHCYC=24/export FHCYC=00/g" config.base
sed -i "s/export FHMAX_GFS_00=240/export FHMAX_GFS_00=168/g" config.base
sed -i "s/export FHMAX_GFS_06=180/export FHMAX_GFS_06=168/g" config.base
sed -i "s/export FHMAX_GFS_12=180/export FHMAX_GFS_12=168/g" config.base
sed -i "s/export FHMAX_GFS_18=180/export FHMAX_GFS_18=168/g" config.base
sed -i "s/export FHOUT_GFS=3/export FHOUT_GFS=1/g" config.base
sed -i "s/export DONST=\"YES\"/export DONST=\"NO\"/g" config.base
sed -i "s/export nst_anl=.true./export nst_anl=.false./g" config.base

!modify config.fcst
sed -i "s/export hord_mt_nh_nonmono=6/export hord_mt_nh_nonmono=5/g" config.fcst
sed -i "s/export hord_xx_nh_nonmono=6/export hord_xx_nh_nonmono=5/g" config.fcst

!modify config.post
sed -i "s/export GOESF=YES/export GOESF=NO/g" config.post
sed -i "s/export FLXF=YES/export FLXF=NO/g" config.post
sed -i "s/export PGB1F=YES/export PGB1F=NO/g" config.post
sed -i "s/export downset=2/export downset=1/g" config.post

!download radial wind obs
cd /gpfs/hps3/emc/meso/noscrub/Donald.E.Lippi/fv3gfs_dl2rw/2018092300/NEXRAD-2018092300-2018100700
htar -xvf /NCEPDEV/emc-meso/5year/Donald.E.Lippi/rw_FV3GFS/drwsim_obs/nexrad_1p00_13km.tar

!edit /gpfs/hps3/emc/meso/save/Donald.E.Lippi/global-workflow-20190306/scripts/exglobal_fcst_nemsfv3gfs.sh
  # Link sfcanl_data restart files from $memdir
- for file in $memdir/RESTART/${PDY}.${cyc}0000.*.nc; do
+ for file in $gmemdir/RESTART/${PDY}.${cyc}0000.*.nc; do
    file2=$(echo $(basename $file))
    file2=$(echo $file2 | cut -d. -f3-) # remove the date from file
    fsufanl=$(echo $file2 | cut -d. -f1)
-   if [ $fsufanl = "sfcanl_data" ]; then
+   if [ $fsufanl = "sfc_data" ]; then
-     file2=$(echo $file2 | sed -e "s/sfcanl_data/sfc_data/g")
+     #file2=$(echo $file2 | sed -e "s/sfcanl_data/sfc_data/g")
      $NLN $file $DATA/INPUT/$file2
    fi
  done

! edit /gpfs/hps3/emc/meso/save/Donald.E.Lippi/global-workflow-20190306/sorc/gsi.fd/scripts/exglobal_analysis_fv3gfs.sh.ecf.NEXRAD
export ndate=/u/Donald.E.Lippi/bin/ndate
valtime_tm00=${PDY}${cyc}
valtime_tm01=`${ndate} -1 $valtime_tm00`
valtime_tm02=`${ndate} -2 $valtime_tm00`
valtime_tm03=`${ndate} -3 $valtime_tm00`
tm00=`echo $valtime_tm00 | cut -c 9-10`
tm01=`echo $valtime_tm01 | cut -c 9-10`
tm02=`echo $valtime_tm02 | cut -c 9-10`
tm03=`echo $valtime_tm03 | cut -c 9-10`

PATH1="/gpfs/hps3/emc/meso/noscrub/Donald.E.Lippi/fv3gfs_dl2rw/2018092300/NEXRAD-2018092300-2018100700"
nexrad_L2RWBUFR_tm00="${PATH1}/nexrad_1p00_13km/simf90/nexrad_${valtime_tm00}_fv3.t${tm00}z_drw.bin"
nexrad_L2RWBUFR_tm01="${PATH1}/nexrad_1p00_13km/simf90/nexrad_${valtime_tm00}_fv3.t${tm01}z_drw.bin"
nexrad_L2RWBUFR_tm02="${PATH1}/nexrad_1p00_13km/simf90/nexrad_${valtime_tm00}_fv3.t${tm02}z_drw.bin"
nexrad_L2RWBUFR_tm03="${PATH1}/nexrad_1p00_13km/simf90/nexrad_${valtime_tm00}_fv3.t${tm03}z_drw.bin"


[[ $DOL2RW = "YES" ]] && $NLN $nexrad_L2RWBUFR_tm00 l2rwbufr_tm00
[[ $DOL2RW = "YES" ]] && $NLN $nexrad_L2RWBUFR_tm01 l2rwbufr_tm01
[[ $DOL2RW = "YES" ]] && $NLN $nexrad_L2RWBUFR_tm02 l2rwbufr_tm02
[[ $DOL2RW = "YES" ]] && $NLN $nexrad_L2RWBUFR_tm03 l2rwbufr_tm03


&GRIDOPTS
  JCAP_B=$JCAP,JCAP=$JCAP_A,NLAT=$NLAT_A,NLON=$NLON_A,nsig=$LEVS,
  regional=.false.,global_l2rw=.true.,nlayers(63)=3,nlayers(64)=6,
  $GRIDOPTS
/


   dfile          dtype       dplat       dsis                dval    dthin dsfcalc
   l2rwbufr_tm00  rw          null        l2rw                0.0     0     0
   l2rwbufr_tm01  rw          null        l2rw                0.0     0     0
   l2rwbufr_tm02  rw          null        l2rw                0.0     0     0
   l2rwbufr_tm03  rw          null        l2rw                0.0     0     0
::


!edit /gpfs/hps3/emc/meso/noscrub/Donald.E.Lippi/fv3gfs_dl2rw/2018092300/NEXRAD-2018092300-2018100700/NEXRAD-2018092300-2018100700.xml
<task name="gdasfcst" cycledefs="gdas" maxtries="&MAXTRIES;">

        <command>&JOBS_DIR;/fcst.sh</command>

        <jobname><cyclestr>&PSLOT;_gdasfcst_@H</cyclestr></jobname>
        <account>&ACCOUNT;</account>
        <queue>&QUEUE_FCST_GDAS;</queue>
        &RESOURCES_FCST_GDAS;
        <walltime>&WALLTIME_FCST_GDAS;</walltime>
        <memory>&MEMORY_FCST_GDAS;</memory>
        <native>&NATIVE_FCST_GDAS;</native>

        <join><cyclestr>&ROTDIR;/logs/@Y@m@d@H/gdasfcst.log</cyclestr></join>

        <envar><name>RUN_ENVIR</name><value>&RUN_ENVIR;</value></envar>
        <envar><name>HOMEgfs</name><value>&HOMEgfs;</value></envar>
        <envar><name>EXPDIR</name><value>&EXPDIR;</value></envar>
        <envar><name>CDATE</name><value><cyclestr>@Y@m@d@H</cyclestr></value></envar>
        <envar><name>CDUMP</name><value>gdas</value></envar>
        <envar><name>PDY</name><value><cyclestr>@Y@m@d</cyclestr></value></envar>
        <envar><name>cyc</name><value><cyclestr>@H</cyclestr></value></envar>

        <dependency>
+           <datadep><cyclestr offset="-6:00:00">&ROTDIR;/gdas.@Y@m@d/@H/RESTART</cyclestr></datadep>
        </dependency>
</task>

!Run experiment
cd /gpfs/hps3/emc/meso/noscrub/Donald.E.Lippi/fv3gfs_dl2rw/2018092300/NEXRAD-2018092300-2018100700
ksh run.ksh
⚠️ **GitHub.com Fallback** ⚠️