Add all changes for global_enspqpf

1) 6c32ae9 on Nov 25, 2020 by Dingchen

 In GEFS v12 implementation, the application of pgpf procudure was accidently
 removed from the ex-script of the atmos_enspost job and consquently, the
 follwing files, their grb2 version as well as the corresponding index
 files (20 files in total) are missing from the the atmos/ensstat directory
 gefs.${PDY}/${cyc}/atmos/ensstat/ensstat.t00z.pqff
 gefs.${PDY}/${cyc}/atmos/ensstat/ensstat.t00z.pqif
 gefs.${PDY}/${cyc}/atmos/ensstat/ensstat.t00z.pqpf
 gefs.${PDY}/${cyc}/atmos/ensstat/ensstat.t00z.pqrf
 gefs.${PDY}/${cyc}/atmos/ensstat/ensstat.t00z.pqsf

 Scripts of the atmos_enspost jobs are modified to fix this problem and
 other minor issues:
 In the J-job script, $COMIN is distinguished from $COMOUT
 In the ex-script, the pqpf block is added back for hr job (0-16d, ext_h= )
 In the ush script, nepert is added to the input parameter list

2) c8f2ca9 on Nov 30, 2020 by Xianwu Xue

 Delete unused statements in sorc/build_global_enspqpf.sh

3) 1451baf on Nov 30, 2020

 In the scripts/exgefs_atmos_enspost.sh script, improved indentation
 and variable name boundary
 in the pqpf code file, updated information about ensemble size

 On branch hotfix/pqpf_ops
	modified:   jobs/JGEFS_ATMOS_ENSPOST
	modified:   scripts/exgefs_atmos_enspost.sh
	modified:   sorc/build_global_enspqpf.sh
	modified:   ush/global_enspqpf.sh

Refs: #2

jobs/JGEFS_ATMOS_ENSPOST

@@ -72,7 +72,7 @@ export PARMgefs=${PARMgefs:-$HOMEgefs/parm}
 ##############################################
 # Define COM directories
 ##############################################
-export COMIN=${COMOUT:-${COMROOT}/${NET}/${envir}/${RUN}.${PDY}/$cyc}
+export COMIN=${COMIN:-${COMROOT}/${NET}/${envir}/${RUN}.${PDY}/$cyc}
 export COMOUT=${COMOUT:-${COMROOT}/${NET}/${envir}/${RUN}.${PDY}/$cyc}
 
 #############################

scripts/exgefs_atmos_enspost.sh

@@ -173,7 +173,7 @@ while [[ $fh -le $FHOUR ]]; do
 			previncr=no
 			for mem in $memberlist; do
 				(( nmem = nmem + 1 ))
-				testfile=$COMIN/$COMPONENT/pgrb2$lr/ge${mem}.${cycle}.pgrb2$FXT\f$fh$EXT.idx
+				testfile=$COMIN/$COMPONENT/pgrb2$lr/ge${mem}.${cycle}.pgrb2${FXT}f$fh$EXT.idx
 				if [ -f $testfile ]; then
 					echo testfile=$testfile found
 					(( nfiles = nfiles + 1 ))
@@ -428,6 +428,14 @@ if [ "$SENDCOM" = "YES" ]; then
 	done # for file in $postvarlist
 fi # [ "$SENDCOM" = "YES" ]
 
+#
+# create pqpf 24h and probabilistic precip forecast files
+## write out with ensemble extended messages
+#
+
+## part (1a): calculate PQPF (named ensppf ) for each 24 hours period,
+#
+
 if [[ "$cycle" == "t00z" ]] && [[ -z $ext_h ]]; then
 	$ENSPPF $COMOUT/$COMPONENT/ensstat/enspost_grb2${ext_h}.$cycle.prcp ensppf${ext_h}.$PDY$cyc.grib2 $npert
 	$WGRIB2 ensppf${ext_h}.$PDY$cyc.grib2 -s >ensppf${ext_h}.$PDY$cyc.grib2.idx 
@@ -460,7 +468,24 @@ if [[ "$cycle" == "t00z" ]] && [[ -z $ext_h ]]; then
 fi # test "$cycle" = "t00z"
 
 # part (1b): probabilistic forecasts ( PQPF, PQRF, PQFF, PQSF and PQIF )
-export CDATE=$PDY$cyc; 
+if [[ -z $ext_h ]]; then
+
+	export CDATE=$PDY$cyc; 
+	$ENSPQPF
+
+	for file in pqpf pqrf pqff pqsf pqif; do
+		$CNVGRIB -g21 $DATA/$file $DATA/${file}_grb1
+		$WGRIB2 $DATA/$file -s >$DATA/${file}.idx
+		$GRBINDEX $DATA/${file}_grb1 $DATA/${file}i_grb1
+		if [ $SENDCOM = "YES" ];        then
+			mv $DATA/${file} $COMOUT/$COMPONENT/ensstat/ensstat_grb2${ext_h}.$cycle.$file
+			mv $DATA/${file}.idx $COMOUT/$COMPONENT/ensstat/ensstat_grb2${ext_h}.$cycle.${file}.idx
+			mv $DATA/${file}_grb1 $COMOUT/$COMPONENT/ensstat/ensstat${ext_h}.$cycle.$file
+			mv $DATA/${file}i_grb1 $COMOUT/$COMPONENT/ensstat/ensstat${ext_h}.$cycle.${file}i
+		fi # [ $SENDCOM = "YES" ]
+	done # for file in pqpf pqrf pqff pqsf pqif
+
+fi
 
 ############################################################################
 ###########  ADD DBN ALERTS FOR PPF AND PQPF FILES IF NEEDED  ##############

sorc/build_global_enspqpf.sh

@@ -32,20 +32,6 @@ export INCSM="-I ${G2_INC4}"
 
 export LIBSM="${G2_LIB4} ${W3NCO_LIB4} ${BACIO_LIB4} ${JASPER_LIB} ${PNG_LIB} ${Z_LIB}"
 
-
-# If you want to get the original size of excutable file, 
-GetOriginal=${GetOriginal:-false}
-
-if $GetOriginal; then
-    if [ $target == wcoss_cray ]; then
-        echo "This is on wcoss_cray"
-        export INCSM="-I ${G2_INC4} -I/opt/cray/iobuf/2.0.5/include"
-        export LIBSM="${G2_LIB4} ${W3NCO_LIB4} ${BACIO_LIB4} ${JASPER_LIB} ${PNG_LIB} ${Z_LIB} -Wl,/opt/cray/iobuf/2.0.5/lib/iobuf.o"
-    elif [ $target == wcoss_dell_p3 ]; then
-        echo "This is on wcoss_dell_p3"
-    fi
-fi
-
 make -f Makefile clobber
 make -f Makefile
 make -f Makefile install

ush/global_enspqpf.sh

@@ -40,12 +40,13 @@ cat <<- EOF >inputpqpf
 		cfrzr='$CFRZR',cfrzro='$CFRZRO'
 		cicep='$CICEP',cicepo='$CICEPO'
 		csnow='$CSNOW',csnowo='$CSNOWO'
+		npert=$npert
 	/
 	EOF
 
 cat inputpqpf
 
-rm $CPGO $CRAINO $CFRZRO $CICEPO $CSNOWO
+#rm $CPGO $CRAINO $CFRZRO $CICEPO $CSNOWO
 
 export pgm=$enspqpf
 source prep_step
@@ -66,4 +67,4 @@ fi
 
 echo "$(date -u) end ${.sh.file}"
 
-exit $err
+exit $err