GetGFSanalysis.csh

#!/bin/csh -f
# Get GFS analysis (0-h forecast) for cold start initial conditions

date

# Setup environment
# =================
source config/model.csh
source config/experiment.csh
source config/builds.csh
set yymmdd = `echo ${CYLC_TASK_CYCLE_POINT} | cut -c 1-8`
set yy = `echo ${CYLC_TASK_CYCLE_POINT} | cut -c 1-4`
set hh = `echo ${CYLC_TASK_CYCLE_POINT} | cut -c 10-11`
set thisCycleDate = ${yymmdd}${hh}
set thisValidDate = ${thisCycleDate}
source ./getCycleVars.csh

# static work directory
set WorkDir = ${InitICWorkDir}/${thisValidDate}
echo "WorkDir = ${WorkDir}"
mkdir -p ${WorkDir}
cd ${WorkDir}

# ================================================================================================

set res = 0p25
set fhour = 000
set linkWPS = link_grib.csh
ln -sfv ${WPSBuildDir}/${linkWPS} .
rm -rf GRIBFILE.*

if ( ${model__AnalysisSource} == "GFSfromRDAOnline" ) then
  echo "Getting GFS analysis from RDA"
  # RDA GFS forecasts directory
  set GFSgribdirRDA = /gpfs/fs1/collections/rda/data/ds084.1
  # link ungribbed GFS
  ./${linkWPS} ${GFSgribdirRDA}/${yy}/${yymmdd}/gfs.${res}.${yymmdd}${hh}.f${fhour}.grib2
else if ( ${model__AnalysisSource} == "GFSfromNCEPFTPOnline" ) then
  echo "Getting GFS analysis from the NCEP FTP"
  # url for GFS data
  set gfs_ftp = https://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.${yymmdd}/${hh}/atmos
  set gfs_file = gfs.t${hh}z.pgrb2.${res}.f${fhour}
  # check if the GFS analysis is available
  if ( ! -e ${gfs_file}) then
    set gfs_ftp_file = ${gfs_ftp}/${gfs_file}
    wget -S --spider $gfs_ftp_file >&! log_check_gfs_f000
    grep "HTTP/1.1 200 OK" log_check_gfs_f000
    # if the file exists then download it
    # otherwise, exit with failure
    if ( $status == 0 ) then
     echo "Downloading $gfs_ftp_file ..."
     wget -r -np -nd $gfs_ftp_file
    else
     echo "$gfs_file not available yet -- exiting"
     exit 1
    endif
  endif
  # link ungribbed GFS
  ./${linkWPS} $gfs_file
endif

date

exit 0