#!/bin/sh # # Metafile Script : gefs_meta_lows.sh_comb # # Log : # J. Carr/HPC 3/10/2003 Moved script from prodsp to hdtb1/dtb. # F. Achorn/NCO 05/18/2006 modify script to run for more than one garea. # F. Achorn/OPC 06/03/2008 Changed the gdfile for members 11-20 and cntrl # (previous run) from pXX_06 to pXX_6 # C. Magee/NCO 10/06/2008 Changed to use COMINs and COMIN for input file # locations (to make testing easier). # # Set Up Local Variables # set -x export PS4='gefs_lows_spag:$SECONDS + ' mkdir $DATA/gefs_meta_lows_spag cd $DATA/gefs_meta_lows_spag sh $utilscript/setup.sh cp $FIXgempak/datatype.tbl datatype.tbl mdl=gefs MDL=GEFS #ddate=`echo $PDY | cut -c3-8` #ddatem1=`echo $PDYm1 | cut -c3-8` fcsthrs="012 024 036 048 060 072 084 096 108 120 132 144 156 168 180 192 204" # SET VARIOUS TIME/DATE PARAMETERS case $cyc in 00)cyc_6=18 PDY_6=$PDYm1 ECM_cyc=12 ECM_PDY=$PDYm1 ECM_fcsthrs=12 # CMC_cyc=00 # CMC_fcsthrs=00 ;; 06)cyc_6=00 PDY_6=$PDY ECM_cyc=12 ECM_PDY=$PDYm1 ECM_fcsthrs=18 # CMC_cyc=00 # CMC_fcsthrs=06 ;; 12)cyc_6=06 PDY_6=$PDY ECM_cyc=00 ECM_PDY=$PDY ECM_fcsthrs=12 # CMC_cyc=12 # CMC_fcsthrs=00 fcsthrs="012 024 036 048 060 072 084 096 108 120 132 144 156 168 180 192 204 216 228" ;; 18)cyc_6=12 PDY_6=$PDY ECM_cyc=00 ECM_PDY=$PDY ECM_fcsthrs=18 # CMC_cyc=12 # CMC_fcsthrs=06 ;; esac metaname="gefs_${PDY}_${cyc}_lows" device="nc|$metaname" # GENERATE THE PMSL LOW CENTERS for area in nam us trop sam hi do proj=" " if [ ${area} = "nam" ]; then num=" " name="NAM" garea="nam" proj=" " elif [ ${area} = "us" ]; then num="#" name="US" garea="bwus" proj=" " elif [ ${area} = "trop" ]; then num=" " name="TROP" garea="5;-100;48;-25" proj="MER" elif [ ${area} = "sam" ]; then num=" " name="SAM" garea="sam" proj=" " else num=" " name="HI" garea="-16;173;34;-67" proj="STR/90.0;-155.0;0.0" fi prodname="LOW CNTRS" ######################################################## # for testing only, restrict to one lvl and fcsthr. will comment these out # when implementing. #fcsthrs="024 048" ######################################################## # START PRODUCTION OF METAFILE for fcsthr in ${fcsthrs} do fcsthr_6=`expr ${fcsthr} + 6` fcsthr_ecm=`expr ${fcsthr} + $ECM_fcsthrs` # fcsthr_cmc=`expr ${fcsthr} + $CMC_fcsthrs` if [ ${fcsthr_6} -lt 10 ]; then fcsthr_6="00${fcsthr_6}" elif [ ${fcsthr_6} -lt 100 ]; then fcsthr_6="0${fcsthr_6}" fi if [ ${fcsthr_ecm} -lt 10 ]; then fcsthr_ecm="00${fcsthr_ecm}" elif [ ${fcsthr_ecm} -lt 100 ]; then fcsthr_ecm="0${fcsthr_ecm}" fi # rm -rf gefs_avg gefs_avg_6 gfs gfs_6 ecmwf ukmet cmc nogaps dgex rm -rf gefs_avg gefs_avg_6 gfs gfs_6 ecmwf ukmet dgex grids="p01 p02 p03 p04 p05 p06 p07 p08 p09 p10 p11 p12 p13 p14 p15 p16 p17 p18 p19 p20 c00" for fn in `echo $grids` do rm -rf $fn ${fn}_6 if [ -r $COMIN/ge${fn}_${PDY}${cyc}f${fcsthr} ] then ln -s $COMIN/ge${fn}_${PDY}${cyc}f${fcsthr} $fn fi if [ -r $COMINs/gefs.${PDY_6}/ge${fn}_${PDY_6}${cyc_6}f${fcsthr_6} ] then ln -s $COMINs/gefs.${PDY_6}/ge${fn}_${PDY_6}${cyc_6}f${fcsthr_6} ${fn}_6 fi done if [ -r $COMIN/geavg_${PDY}${cyc}f${fcsthr} ] then ln -s $COMIN/geavg_${PDY}${cyc}f${fcsthr} gefs_avg fi if [ -r $COMINs/gefs.${PDY_6}/geavg_${PDY_6}${cyc_6}f${fcsthr_6} ] then ln -s $COMINs/gefs.${PDY_6}/geavg_${PDY_6}${cyc_6}f${fcsthr_6} gefs_avg_6 fi if [ -r $COMINs/gfs.${PDY}/gfs_${PDY}${cyc}f${fcsthr} ] then ln -s $COMINs/gfs.${PDY}/gfs_${PDY}${cyc}f${fcsthr} gfs fi if [ -r $COMINs/gfs.${PDY_6}/gfs_${PDY_6}${cyc_6}f${fcsthr_6} ] then ln -s $COMINs/gfs.${PDY_6}/gfs_${PDY_6}${cyc_6}f${fcsthr_6} gfs_6 fi if [ -r $COMINs/ecmwf.${ECM_PDY}/ecmwf_glob_${ECM_PDY}${ECM_cyc} ] then ln -s $COMINs/ecmwf.${ECM_PDY}/ecmwf_glob_${ECM_PDY}${ECM_cyc} ecmwf fi if [ -r $COMINs/ukmet.${ECM_PDY}/ukmet_${ECM_PDY}${ECM_cyc}f${fcsthr_ecm} ] then ln -s $COMINs/ukmet.${ECM_PDY}/ukmet_${ECM_PDY}${ECM_cyc}f${fcsthr_ecm} ukmet fi # if [ -r $COMINs/cmc.${PDY}/cmc${PDY}${CMC_cyc}f${fcsthr_cmc} ] # then # ln -s $COMINs/cmc.${PDY}/cmc${PDY}${CMC_cyc}f${fcsthr_cmc} cmc # fi # # if [ -r $COMINs/nogaps.${PDY}/nogaps${PDY}${CMC_cyc}f${fcsthr_cmc} ] # then # ln -s $COMINs/nogaps.${PDY}/nogaps${PDY}${CMC_cyc}f${fcsthr_cmc} nogaps # fi # if [ -r $COMINs/dgex.${PDY}/dgex_${PDY}${cyc}f${fcsthr} ] then ln -s $COMINs/dgex.${PDY}/dgex_${PDY}${cyc}f${fcsthr} dgex fi cat > cmdfilelows << EOF DEVICE = ${device} PANEL = 0 TEXT = s/22/1/1/hw CONTUR = 2 MAP = 1 CLEAR = yes GAREA = ${garea} PROJ = ${proj} LATLON = 1/10/1/2/10;10 GLEVEL = 0 GVCORD = none SKIP = 0 SCALE = 0 GDPFUN = sm5s(pmsl) TYPE = c !CINT = 4/1/8 CINT = 4 LINE = 0 FINT = FLINE = HILO = 25/L${num}/900-1016/5/50/y HLSYM = l/22/3/hw CLRBAR = 0 WIND = 0 REFVEC = GDFILE = gfs HILO = 22/L${num}/900-1016/5/50/y TITLE = 22/-2/~ ? GFS ${cyc}Z|~${name} GDATTIM = F${fcsthr} run MAP = 0 LATLON = 0 CLEAR = no GDFILE = gfs_6 HILO = 21/L${num}/900-1016/5/50/y TITLE = 21/-1/~ ? GFS ${cyc_6}Z|~${name} GDATTIM = F${fcsthr_6} GDFILE = p01_6 ! p02_6 ! p03_6 ! p04_6 ! p05_6 ! p06_6 ! p07_6 ! p08_6 ! p09_6 ! p10_6 HILO = 25/L${num}/900-1016/5/50/y TITLE = 25/+1/~ ? ${cyc_6}Z ENS MEMBERS|~${name} GDATTIM = F${fcsthr_6} run GDFILE = p11_6 ! p12_6 ! p13_6 ! p14_6 ! p15_6 ! p16_6 ! p17_6 ! p18_6 ! p19_6 ! p20_6 ! c00_6 run EOF grids="C00 P01 P02 P03 P04 P05 P06 P07 P08 P09 P10 P11 P12 P13 P14 P15 P16 P17 P18 P19 P20" line_count=2 color_number=9 for grid in ${grids} do gridl=`echo $grid | tr [A-Z] [a-z]` cat >> cmdfilelows << EOF !GDFILE = \$COMIN/ge${gridl}_${PDY}${cyc}f${fcsthr} GDFILE = $gridl HILO = ${color_number}/L${num}/900-1016/5/50/y TITLE = ${color_number}/+${line_count}/~ ? ${cyc}Z ${grid}|~${name} GDATTIM = F${fcsthr} run EOF let line_count=$line_count+1 let color_number=$color_number+1 done cat >> cmdfilelows << EOF HILO = 21/L${num}/900-1016/5/50/y ! 22/L${num}/900-1016/5/50/y GDFILE = ukmet GDATTIM = F${fcsthr_ecm} HILO = 7/L${num}/900-1016/5/50/y TITLE = 7/-5/~ ? UKMET ${ECM_cyc}Z|~${name} run GDFILE = ecmwf GDATTIM = F${fcsthr_ecm} HILO = 6/L${num}/900-1016/5/50/y TITLE = 6/-4/~ ? ECMWF ${ECM_cyc}Z|~${name} run GDFILE = gefs_avg HILO = 1/L${num}/900-1016/5/50/y TITLE = 1/-7/~ ? ENSMEAN ${cyc}Z |~${name} GDATTIM = F${fcsthr} run GDFILE = gefs_avg_6 HILO = 31/L${num}/900-1016/5/50/y TITLE = 31/-6/~ ? ENSMEAN ${cyc_6}Z|~${name} GDATTIM = F${fcsthr_6} run GDFILE = dgex GDATTIM = F${fcsthr} HILO = 13/L${num}/900-1016/5/50/y TITLE = 13/-8/~ ? ${cyc}Z DGEX|~${name} run EOF echo cdmfile gdplot2_nc < cmdfilelows done done if [ $SENDCOM = "YES" ] ; then mv ${metaname} ${COMOUT}/$metaname if [ $SENDDBN = "YES" ] ; then $DBNROOT/bin/dbn_alert MODEL ${DBN_ALERT_TYPE} $job ${COMOUT}/$metaname fi fi gpend rm gemglb.nts last.nts exit