#!/bin/ksh export PATH=/bin:/usr/bin:/usr/bsd:/etc:/usr/ucb:/usr/bin/X11:/usr/lpp/X11/Xamples/bin:$HOME/bin:$HOME/mcidas/bin:~mcidas/bin:. export MCPATH=$HOME/mcidas/data:$HOME/mcidas/help:~mcidas/data:~mcidas/help export SHELL=/bin/ksh cat << 'EOF' | ~mcidas/bin/mcenv adtdir="$HOME/aodt/ADTV7.2.3" scripts_dir="$adtdir/scripts" fcstfiletype[0]="ATCF" fcstfiletype[1]="DISC" fcstfiletype[2]="PACWARN" fcstfiletype[3]="GENERIC" fcstfiletype[4]="ICAO" fcstfiletype[5]="RSMC" fcstfiletype[6]="AXAU" wmoheader=WTNT stormnumber=44 stormname=19L ifcsttype=1 let ttime=24 let ic_flag=1 # first time use on new storm. day=`date +%Y%j` time=`date +%T` minute=`date +%M` let timex=10#`date +%H%M%S` origstormname=$stormname echo "Begin getforecast_met10_91L.sh" ls -la $adtdir/history/$stormname'.ODT' if [ -f $adtdir/history/$stormname'.ODT' ] then echo " Storm is already there - Don't use IC " ic_flag=-1 fi echo "ifcsttype == $ifcsttype" echo "day == $day" echo "time == $timex" echo "ttime == $ttime" # tlo added 3/10/09 start sizefile=`ls -l $adtdir/forecasts/$stormname.FORECAST | awk {'print $5'}` echo "sizefile == $sizefile" if [ $sizefile -lt 34 ] then echo " Storm too weak or forecast non-existant" let valid=-99 else echo "$adtdir/forecasts/$stormname.FORECAST" echo "$scripts_dir/interpll-all $adtdir/forecasts/$stormname.FORECAST $ifcsttype $day $timex $ttime >$scripts_dir/interpll_output" # $scripts_dir/interpll $adtdir/forecasts/$stormname.FORECAST $ifcsttype $day $timex $ttime >$scripts_dir/interpll_output $scripts_dir/interpll-all $adtdir/forecasts/$stormname.FORECAST $ifcsttype $day $timex $ttime >$scripts_dir/interpll_output echo interpll_output OUTPUT cat $scripts_dir/interpll_output read valid lat lon strength < $scripts_dir/interpll_output echo $valid $lat $lon $strength rm -r -f $scripts_dir/interpll_output fi # tlo added 3/10/09 end if [ $valid -lt 0 ] then if [ $valid -eq -1 ] then echo "INVALID FORECAST FILE, IMAGE BEFORE FIRST FORECAST TIME" CSTAT="INVALID_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "INVALID FORECAST FILE, IMAGE BEFORE FIRST FORECAST TIME" echo " " elif [ $valid -eq -2 ] then echo "INVALID FORECAST FILE, IMAGE AFTER LAST FORECAST TIME" CSTAT="INVALID_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "INVALID FORECAST FILE, IMAGE AFTER LAST FORECAST TIME" echo " " elif [ $valid -eq -3 ] then echo "FORECAST FILE MORE THAN $ttime HOURS OLD" CSTAT="24HRS_OLD_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "FORECAST FILE MORE THAN 24 HOURS OLD" echo " " # Not run the mklist command from the back up TN 03/23/09 # echo "The run_mklist.sh script is called !!!!" # $HOME/aodt/ADTV7.2.3/scripts/run_mklist.sh $stormname # elif [ $valid -eq -3 ] # then # # echo "dates are invalid. Either image date is before forecast (month switch) # # or is well beyond forecast availability # echo "Invalid forecast file (JTWC forecast)" # CSTAT="INVALJTWC_FCST" elif [ $valid -eq -4 ] then echo "INTERPOLATION FAILED" CSTAT="INTERP_FAILED" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "INTERPOLATION FAILED" echo " " elif [ $valid -eq -5 ] then echo "FORECAST FILE DOES NOT EXIST... EXITING" CSTAT="NOT_EXIST_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "FORECAST FILE DOES NOT EXIST... EXITING" echo " " # tlo added 3/10/09 start elif [ $valid -eq -99 ] then echo "STORM TOO WEAK... EXITING" CSTAT="INVALID_FCST" # tlo added 3/10/09 end echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "STORM TOO WEAK...... EXITING" echo " " else CSTAT="UNKNOWN" fi curimgtime=999999 # $adtdir/scripts/emailnote.sh $CSTAT $stormname $scripts_dir/adtstats.sh $day $timex $curimgtime $wmoheader $stormnumber $stormname $valid $CSTAT exit else # RUN ODT FOR STORM wind=$strength echo "wind = $wind" ####------- Convert the lat/lon from float to into ----------- #### let ilat=`$HOME/aodt/ADTV7.2.3/scripts/floattoint $lat` let ilon=`$HOME/aodt/ADTV7.2.3/scripts/floattoint $lon` echo lat=$lat $ilat lon=$lon $ilon ####------- Convert the Max wind to Tno ----------- #### let w1or10=1 # added 5/6/09 tlo if [[ $wmoheader = "WTIO" && $stormnumber -eq 30 ]] then ic=$wind elif [ $wmoheader = "AXAU" ] then ic=$wind else ic=`$HOME/aodt/ADTV7.2.3/scripts/windtotno $wind $w1or10` # modified 5/6/09 tlo echo IC=$ic fi ####------- Download image from server ----------- #### ldsname="LOC/AREA" let image=9989 mceval.k IMGDEL $ldsname.$image echo lat=$lat echo lon=$lon ####------- Determine the Satellite ----------- #### #if [[ $ilon -le -3000 && $ilon -gt -10000 ]] #then ## INDIAN OCEAN #satellite_data='IND/INDOEXIR' #bsat='MET-7' #fcsttype='ICAO' #fi fcsttype=${fcstfiletype[${ifcsttype}]} if [ $ilat -gt 0 ] then # NORTHERN HEMISPHERE if [ $ilon -le -10000 ] then # WEST PACIFIC satellite_data='MTS/MTNHEM04I2' bsat='MTSAT' #fcsttype='PACWARN' fi if [[ $ilon -gt -10000 && $ilon -lt 0 ]] then # INDIAN OCEAN satellite_data='IND/INDOEXIR' bsat='MET-7' #fcsttype='PACWARN' fi if [[ $ilon -ge 0 && $ilon -lt 3500 ]] then # EAST ATLANTIC satellite_data='MSG/MSGLOB09I' bsat='MET-9' #fcsttype='DISC' fi if [[ $ilon -ge 3500 && $ilon -lt 10500 ]] then # WEST ATLANTIC satellite_data='GCR/GCNHEM04I4' bsat='GOES-13' #fcsttype='DISC' fi if [ $ilon -ge 10500 ] then # EAST PACIFIC satellite_data='GCR/GCNHEM04I4' bsat='GOES-14' #fcsttype='DISC' fi else # SOUTHERN HEMISPHERE if [ $ilon -le -10000 ] then # WEST PACIFIC ## Switch from using MTGLOB to using MTSHEM. YL 05/08/09 satellite_data='MTS/MTSHEM04I2' bsat='MTSAT' #fcsttype='PACWARN' fi if [[ $ilon -gt -10000 && $ilon -lt 0 ]] then # INDIAN OCEAN satellite_data='IND/INDOEXIR' bsat='MET-7' #fcsttype='ICAO' fi if [[ $ilon -ge 0 && $ilon -lt 6000 ]] then # EAST ATLANTIC satellite_data='MSG/MSGLOB09I' bsat='MET-9' #fcsttype='DISC' fi if [ $ilon -ge 10000 ] then # EAST PACIFIC satellite_data='GCR/GCSHEM04I4' bsat='GOES-14' #fcsttype='DISC' fi fi echo "fcsttype= $fcsttype" mceval.k IMGLIST $satellite_data.-5 ## Take the current image instead of the image minus one. YL 04/30/09 mceval.k IMGCOPY $satellite_data $ldsname.$image SIZE=500 500 LATLON=$lat $lon mceval.k IMGLIST $ldsname.$image ####------- check the image time for the current image. ----------- #### imglist.k $ldsname.$image | grep -i "IMG" > /dev/null if [ $? = 0 ] then curimghh=`imglist.k $ldsname.$image | grep $image | awk '{print $7}' |cut -c1-2` curimgmm=`imglist.k $ldsname.$image | grep $image | awk '{print $7}' |cut -c4-5` curimgss=`imglist.k $ldsname.$image | grep $image | awk '{print $7}' |cut -c7-8` curimgtime=$curimghh$curimgmm$curimgss fi imglist.k $ldsname.$image | grep -i "MET" > /dev/null if [ $? = 0 ] then curimghh=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c1-2` curimgmm=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c4-5` curimgss=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c7-8` curimgtime=$curimghh$curimgmm$curimgss fi imglist.k $ldsname.$image | grep -i "MSG" > /dev/null if [ $? = 0 ] then curimghh=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c1-2` curimgmm=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c4-5` curimgss=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c7-8` curimgtime=$curimghh$curimgmm$curimgss fi imglist.k $ldsname.$image | grep -i "MTS" > /dev/null if [ $? = 0 ] then curimghh=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c1-2` curimgmm=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c4-5` curimgss=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c7-8` curimgtime=$curimghh$curimgmm$curimgss fi imglist.k $ldsname.$image | grep -i "GMS" > /dev/null if [ $? = 0 ] then curimghh=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c1-2` curimgmm=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c4-5` curimgss=`imglist.k $ldsname.$image | grep $image | awk '{print $6}' |cut -c7-8` curimgtime=$curimghh$curimgmm$curimgss fi echo " Current image time == $curimgtime" ####------- check the previous image time from History list file ----------- #### if [ -f $HOME/aodt/ADTV7.2.3/history/$stormname'-list.txt' ] then preimagetime=`tail -1 $HOME/aodt/ADTV7.2.3/history/$stormname'-list.txt' | awk '{print $2}'` echo " previous image time == $preimagetime" ####------- Run ADT on the new image only ----------- #### if [ $preimagetime -eq $curimgtime ] then echo " The current image is already processed -- Exit ADT" CSTAT="OLDIMG" $scripts_dir/adtstats.sh $day $timex $curimgtime $wmoheader $stormnumber $stormname $valid $CSTAT echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " tail -10 $adtdir/history/$stormname-list.txt echo " " exit else echo RUNNING ODT ON $stormname echo $ic_flag $scripts_dir/runaodt.sh $fcsttype $stormname $lat $lon $ilon $ic $ic_flag $bsat $origstormname ####------- Run ADT again if the output record contains [EXTRP] value ----------- #### tail -1 $HOME/aodt/ADTV7.2.3/history/$stormname'-list.txt' | grep -i "EXTRP" if [ $? = 0 ] then echo " [EXTRP] return - Run ADT again " cp $adtdir/forecasts/$stormname.FORECAST.pregood $adtdir/forecasts/$stormname.FORECAST echo "$scripts_dir/interpll-all $adtdir/forecasts/$stormname.FORECAST $ifcsttype $day $timex $ttime" $scripts_dir/interpll-all $adtdir/forecasts/$stormname.FORECAST $ifcsttype $day $timex $ttime >$scripts_dir/interpll_output read valid lat lon strength < $scripts_dir/interpll_output echo $valid $lat $lon $strength rm -r -f $scripts_dir/interpll_output if [ $valid -lt 0 ] then if [ $valid == -1 ] then echo "INVALID FORECAST FILE, IMAGE BEFORE/AFTER FORECAST TIMES" CSTAT="INVALID_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "INVALID FORECAST FILE, IMAGE BEFORE FIRST FORECAST TIME" echo " " elif [ $valid -eq -2 ] then echo "INVALID FORECAST FILE, IMAGE AFTER LAST FORECAST TIME" CSTAT="INVALID_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "INVALID FORECAST FILE, IMAGE AFTER LAST FORECAST TIME" echo " " elif [ $valid == -3 ] then echo "FORECAST FILE MORE THAN 24 HOURS OLD" CSTAT="24HRS_OLD_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "FORECAST FILE MORE THAN 24 HOURS OLD" echo " " # Not run the mklist command from the back up TN 03/23/09 # echo "The run_mklist.sh script is called !!!!" # $HOME/aodt/ADTV7.2.3/scripts/run_mklist.sh $stormname #elif [ $valid == -3 ] #then # # echo "dates are invalid. Either image date is before forecast (month switch) # # or is well beyond forecast availability # echo "Invalid forecast file (JTWC forecast)" # CSTAT="INVALJTWC_FCST" elif [ $valid == -4 ] then echo "INTERPOLATION FAILED" CSTAT="INTERP_FAILED" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "INTERPOLATION FAILED" echo " " elif [ $valid == -5 ] then echo "FORECAST FILE DOES NOT EXIST... EXITING" CSTAT="NOT_EXIST_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "FORECAST FILE DOES NOT EXIST... EXITING" echo " " # tlo added 3/10/09 start elif [ $valid -eq -99 ] then echo "STORM TOO WEAK... EXITING" CSTAT="INVALID_FCST" # tlo added 3/10/09 end echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "STORM TOO WEAK...... EXITING" echo " " else echo "UNKNOWN" CSTAT="UNKNOWN" fi # $adtdir/scripts/emailnote.sh $CSTAT $stormname $scripts_dir/adtstats.sh $day $timex $curimgtime $wmoheader $stormnumber $stormname $valid $CSTAT exit else $scripts_dir/runaodt.sh $fcsttype $stormname $lat $lon $ilon $ic $ic_flag $bsat $origstormname fi fi fi else echo RUNNING ODT ON $stormname echo $ic_flag $scripts_dir/runaodt.sh $fcsttype $stormname $lat $lon $ilon $ic $ic_flag $bsat $origstormname tail -1 $HOME/aodt/ADTV7.2.3/history/$stormname'-list.txt' | grep -i "EXTRP" if [[ $? = 0 ]] then echo " [EXTRP] return - Run ADT again " cp $adtdir/forecasts/$stormname.FORECAST.pregood $adtdir/forecasts/$stormname.FORECAST echo "$scripts_dir/interpll-all $adtdir/forecasts/$stormname.FORECAST $ifcsttype $day $timex $ttime" $scripts_dir/interpll-all $adtdir/forecasts/$stormname.FORECAST $ifcsttype $day $timex $ttime >$scripts_dir/interpll_output read valid lat lon strength < $scripts_dir/interpll_output echo $valid $lat $lon $strength rm -r -f $scripts_dir/interpll_output if [ $valid -lt 0 ] then if [ $valid == -1 ] then echo "INVALID FORECAST FILE, IMAGE BEFORE/AFTER FORECAST TIMES" CSTAT="INVALID_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "INVALID FORECAST FILE, IMAGE BEFORE/AFTER FORECAST TIMES" echo " " elif [ $valid == -2 ] then echo "INVALID FORECAST FILE, IMAGE AFTER LAST FORECAST TIME" CSTAT="INVALID_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "INVALID FORECAST FILE, IMAGE AFTER LAST FORECAST TIME" echo " " elif [ $valid -eq -3 ] then echo "FORECAST FILE MORE THAN $ttime HOURS OLD" CSTAT="24HRS_OLD_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "FORECAST FILE MORE THAN 24 HOURS OLD" echo " " # # run the mklist command TN 03/23/09 # echo "The run_mklist.sh script is called !!!" # $HOME/aodt/ADTV7.2.3/scripts/run_mklist.sh $stormname elif [ $valid == -4 ] then echo "INTERPOLATION FAILED" CSTAT="INTERP_FAILED" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "INTERPOLATION FAILED" echo " " elif [ $valid == -5 ] then echo "FORECAST FILE DOES NOT EXIST... EXITING" CSTAT="NOT_EXIST_FCST" echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "FORECAST FILE DOES NOT EXIST... EXITING" echo " " # tlo added 3/10/09 start elif [ $valid -eq -99 ] then echo "STORM TOO WEAK... EXITING" CSTAT="INVALID_FCST" # tlo added 3/10/09 end echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "STORM TOO WEAK...... EXITING" echo " " else CSTAT="UNKNOWN" fi # $adtdir/scripts/emailnote.sh $CSTAT $stormname $scripts_dir/adtstats.sh $day $timex $curimgtime $wmoheader $stormnumber $stormname $valid $CSTAT exit else echo $fcsttype $stormname $lat $lon $ilon $ic $ic_flag $bsat $origstormname $scripts_dir/runaodt.sh $fcsttype $stormname $lat $lon $ilon $ic $ic_flag $bsat $origstormname fi fi fi grep -i "AODT Error" $HOME/aodt/ADTV7.2.3/scripts/xxxx > /dev/null if [ $? = 0 ] then echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " echo "NAODT command ERROR" echo "Please contact the programmer and open a help desk ticket" echo " " CSTAT="ERR_NAODT" else echo "LIST TEXT FOR ACTIVE STORM NAME = $stormname" echo " " tail -10 $adtdir/history/$stormname-list.txt echo " " CSTAT="GOODRUN" fi fi $scripts_dir/adtstats.sh $day $timex $curimgtime $wmoheader $stormnumber $stormname $valid $CSTAT ####------- Arhive ADT oupputs to the final direcotry ----------- #### cp $adtdir/forecasts/$stormname'.FORECAST' $adtdir/forecasts/$stormname'.FORECAST.pregood' cp $adtdir/forecasts/$stormname'.FORECAST' $adtdir/final/fcst/$stormname'.FORECAST'_$wmoheader$stormnumber'_'$day'_'$timex cp $HOME/aodt/ADTV7.2.3/scripts/xxxx $HOME/aodt/ADTV7.2.3/final/text/xxxx.$stormname'_'$day'_'$timex ####------- Send ADT products to the Web Server ----------- #### #### Stop to send to the Web Server (GP5) from the back up TN 03/23/09 ### Bob Potash activate this send 2012JUNE30 #$HOME/aodt/ADTV7.2.3/scripts/send_adt2webserv.sh $stormname #mv $HOME/mcidas/data/$stormname'.GIF' $HOME/aodt/ADTV7.2.3/final/gif/'0'$stormname'_'$day'_'$timex'.GIF' #cp $HOME/aodt/ADTV7.2.3/history/$stormname'.ODT' $HOME/aodt/ADTV7.2.3/final/history/'0'$stormname'.ODT' #cp $HOME/aodt/ADTV7.2.3/history/$stormname'-list.txt' $HOME/aodt/ADTV7.2.3/final/history/'0'$stormname'-list.txt' ####------- Post ADT products to the Web page ----------- #### #### Stop to Post to the web site from the back up TN 03/23/09 ### $HOME/aodt/ADTV7.2.3/scripts/post_adt2webpage.pl $stormname-list.txt echo "End getforecast_met10_91L.sh" EOF