##################################################################### echo "-----------------------------------------------------" echo " exglobal_post_gfs.sh.sms" echo " Feb 06 - Wobus - new script to convert gfs to pgrba" echo " Jul 11 - Wobus - rename to exgefs_prdgen_gfs.sh.sms" echo "-----------------------------------------------------" ##################################################################### echo `date` $0 begin set -x cd $DATA msg="HAS BEGUN on `hostname`" postmsg "$jlogfile" "$msg" export MP_LABELIO=YES #poe hostname echo cyc=$cyc echo cycle=$cycle export WGRIB=${WGRIB:-$EXECUTIL/wgrib} export GRBIDX=${GRBIDX:-$EXECUTIL/grbindex} export COPYGB=${COPYGB:-$EXECUTIL/copygb} export CNVGRIB=${CNVGRIB:-$EXECUTIL/cnvgrib} export GRB2IDX=${GRB2IDX:-$EXECUTIL/grb2index} export WGRIB2=${WGRIB2:-$EXECUTIL/wgrib2} export ENSADD=${ENSADD:-$USHGLOBAL/global_ensadd.sh} echo WGRIB=$WGRIB echo COPYGB=$COPYGB echo CNVGRIB=$CNVGRIB echo GRBINDEX=$GRBINDEX echo GRBIDX=$GRBIDX echo GRB2INDEX=$GRB2INDEX echo ENSADD=$ENSADD # Turn off the generation of the grib2 index files for this job: export makegrb2i=no # set variables for ensemble PDS header (( e1 = 1 )) (( e2 = 1 )) ############################################################ # Define Variables: # ----------------- # SHOUR is the starting forecast hour. normally 0 except for restarts. # FHOUR is the ending forecast hour. # FHINC is the increment hour for each forecast steps. # FH is the current forecast hour. # SLEEP_TIME is the number of seconds to sleep before exiting with error. # SLEEP_INT is the number of seconds to sleep between restrt file checks. # restart_file is the name of the file to key off of to kick off post. ############################################################ ############################################################ # Post Analysis Files before starting the Forecast Post ############################################################ if test -f $COMIN/${RUN}.t${cyc}z.master.grbanl -a ${SHOUR} -eq 0 then # RLW 20110725 more complete cleanup of temporary files rm -f master_grbfile rm -f tmpfile rm -f pgbafile pgbaifile rm -f pgb2afile pgb2aifile parmlist=$PARMgefs/gefs_pgrba_f00.parm ln -s $COMIN/${RUN}.t${cyc}z.master.grbanl master_grbfile $WGRIB -s master_grbfile |grep -F -f $parmlist |$WGRIB -i -grib -o tmpfile master_grbfile $COPYGB -g3 -i0 -x tmpfile pgbafile $GRBINDEX pgbafile pgbaifile $ENSADD $e1 $e2 pgbafile pgbaifile epgbafile if [[ "$makepgrb2a" = "yes" ]]; then $CNVGRIB -g12 -p40 epgbafile pgb2afile fi if [[ "$makegrb2i" = "yes" ]]; then if [[ "$makepgrb2a" = "yes" ]]; then $GRB2INDEX pgb2afile pgb2aifile fi fi # # control the inclusion of perturbation identifiers # in the GRIB1 ensemble PDS header extension # if [[ "$addgrb1id" = "yes" ]]; then mv epgbafile pgbafile if [[ "$makegrb1i" = "yes" ]]; then $GRBINDEX pgbafile pgbaifile fi fi if test "$SENDCOM" = 'YES' then # # Save Pressure GRIB/Index files # mv pgbafile $COMOUT/$cyc/pgrba/ge${RUN}.${cycle}.pgrbaanl if [[ "$makegrb1i" = "yes" ]]; then mv pgbaifile $COMOUT/$cyc/pgrba/ge${RUN}.${cycle}.pgrbaianl fi if [[ "$makepgrb2a" = "yes" ]]; then mv pgb2afile $COMOUT/$cyc/pgrb2a/ge${RUN}.${cycle}.pgrb2aanl fi if [[ "$makegrb2i" = "yes" ]]; then if [[ "$makepgrb2a" = "yes" ]]; then mv pgb2aifile $COMOUT/$cyc/pgrb2a/ge${RUN}.${cycle}.pgrb2aianl fi fi #################################### CHECK DBN ALERTS ########### #################################### CHECK DBN ALERTS ########### #################################### CHECK DBN ALERTS ########### if test "$SENDDBN" = 'YES' then if test "$NET" = 'gens' then MEMBER=GFS $DBNROOT/bin/dbn_alert MODEL ENS_PGBA_$MEMBER $job $COMOUT/$cyc/pgrba/ge${RUN}.${cycle}.pgrbaanl if [[ "$makegrb1i" = "yes" ]]; then $DBNROOT/bin/dbn_alert MODEL ENS_PGBAI_$MEMBER $job $COMOUT/$cyc/pgrba/ge${RUN}.${cycle}.pgrbaianl fi if [[ "$makepgrb2a" = "yes" ]]; then $DBNROOT/bin/dbn_alert MODEL ENS_PGB2A_$MEMBER $job $COMOUT/$cyc/pgrb2a/ge${RUN}.${cycle}.pgrb2aanl fi if [[ "$makegrb2i" = "yes" ]]; then if [[ "$makepgrb2a" = "yes" ]]; then $DBNROOT/bin/dbn_alert MODEL ENS_PGB2AI_$MEMBER $job $COMOUT/$cyc/pgrb2a/ge${RUN}.${cycle}.pgrb2aianl fi fi fi fi fi fi SLEEP_LOOP_MAX=`expr $SLEEP_TIME / $SLEEP_INT` (( fhr=SHOUR)) if (( fhr == 0 )); then fhr=00 fi export fhr if test $fhr -lt 10 -a $fhr -gt 0 then export fhr="0$fhr" fi ############################################################ # Loop Through the Post Forecast Files ############################################################ while test $fhr -le $FHOUR do ############################### # Start Looping for the # existence of the restart files ############################### #set +x set -x export pgm="postcheck" ic=1 while [ $ic -le $SLEEP_LOOP_MAX ] do if [[ -s $restart_file_a$fhr ]] && [[ -s $restart_file_b$fhr ]]; then found=yes break else ic=`expr $ic + 1` sleep $SLEEP_INT fi if (( fhr > 180 )) && (( fhr % 12 == 6 )); then if (( ic > 2 )) && (( ic < SLEEP_LOOP_MAX )); then echo fhr=$fhr not expected found=no break fi fi ############################### # If we reach this point assume # fcst job never reached restart # period and error exit ############################### if [ $ic -eq $SLEEP_LOOP_MAX ] then export err=9 err_chk fi done #set -x if [[ "$found" = yes ]]; then # RLW 20110725 more complete cleanup of temporary files rm -f master_grbfile rm -f tmpfile rm -f pgbafile pgbaifile rm -f pgb2afile pgb2aifile msg="Starting prdgen for fhr=$fhr" postmsg "$jlogfile" "$msg" if [ $fhr -eq 0 ] then parmlist=$PARMgefs/gefs_pgrba_f00.parm else parmlist=$PARMgefs/gefs_pgrba_fhh.parm fi ln -s $COMIN/${RUN}.t${cyc}z.master.grbf${fhr} master_grbfile $WGRIB -s master_grbfile |grep -F -f $parmlist |$WGRIB -i -grib -o tmpfile master_grbfile if [[ x$fhoroglist != x ]]; then for fhorog in $fhoroglist do if (( fhr == fhorog )); then $WGRIB -s master_grbfile |grep 'HGT:sfc' |$WGRIB -i -grib -append -o tmpfile master_grbfile fi done fi $COPYGB -g3 -i0 -x tmpfile pgbafile $GRBINDEX pgbafile pgbaifile $ENSADD $e1 $e2 pgbafile pgbaifile epgbafile if [[ "$makepgrb2a" = "yes" ]]; then $CNVGRIB -g12 -p40 epgbafile pgb2afile fi if [[ "$makegrb2i" = "yes" ]]; then if [[ "$makepgrb2a" = "yes" ]]; then $GRB2INDEX pgb2afile pgb2aifile fi fi if [[ "$addgrb1id" = "yes" ]]; then mv epgbafile pgbafile if [[ "$makegrb1i" = "yes" ]]; then $GRBINDEX pgbafile pgbaifile fi fi # RLW 20110725 more complete cleanup of temporary files rm -f tmpfile.2 rm -f pgbafile.2 pgbaifile.2 rm -f pgb2afile.2 pgb2aifile.2 if test "$DO_LOW_RES" = 'YES' -a `expr $fhr % 6 ` -eq 0 then $WGRIB -s master_grbfile |grep -F -f $parmlist |$WGRIB -i -grib -o tmpfile.2 master_grbfile if [[ x$fhoroglist != x ]]; then for fhorog in $fhoroglist do if (( fhr == fhorog )); then $WGRIB -s master_grbfile |grep 'HGT:sfc' |$WGRIB -i -grib -append -o tmpfile.2 master_grbfile fi done fi $COPYGB -g2 -i0 -x tmpfile.2 pgbafile.2 $GRBINDEX pgbafile.2 pgbaifile.2 $ENSADD $e1 $e2 pgbafile.2 pgbaifile.2 epgbafile.2 if [[ "$makepgrb2a" = "yes" ]]; then $CNVGRIB -g12 -p40 epgbafile.2 pgb2afile.2 fi if [[ "$makegrb2i" = "yes" ]]; then if [[ "$makepgrb2a" = "yes" ]]; then $GRB2INDEX pgb2afile.2 pgb2aifile.2 fi fi if [[ "$addgrb1id" = "yes" ]]; then mv epgbafile.2 pgbafile.2 if [[ "$makegrb1i" = "yes" ]]; then $GRBINDEX pgbafile.2 pgbaifile.2 fi fi fi $USHGLOBAL/gefs_transfer_gfs.sh if test $SENDCOM = "YES" then if test $fhr -lt 100 then pad="0" else pad="" fi echo "$PDY$cyc$pad$fhr" > $COMOUT/$cyc/misc/ge${RUN}.t${cyc}z.control fi fi # if not found, come here to increment export fhr=`expr $fhr + $FHINC` if test $fhr -lt 10 then export fhr="0$fhr" fi done cat $pgmout ######################################################## echo `date` $0 end msg='ENDED NORMALLY.' postmsg "$jlogfile" "$msg" ################## END OF SCRIPT #######################