#!/bin/csh -f
# +-======-+ 
#  Copyright (c) 2003-2007 United States Government as represented by 
#  the Admistrator of the National Aeronautics and Space Administration.  
#  All Rights Reserved.
#  
#  THIS OPEN  SOURCE  AGREEMENT  ("AGREEMENT") DEFINES  THE  RIGHTS  OF USE,
#  REPRODUCTION,  DISTRIBUTION,  MODIFICATION AND REDISTRIBUTION OF CERTAIN 
#  COMPUTER SOFTWARE ORIGINALLY RELEASED BY THE UNITED STATES GOVERNMENT AS 
#  REPRESENTED BY THE GOVERNMENT AGENCY LISTED BELOW ("GOVERNMENT AGENCY").  
#  THE UNITED STATES GOVERNMENT, AS REPRESENTED BY GOVERNMENT AGENCY, IS AN 
#  INTENDED  THIRD-PARTY  BENEFICIARY  OF  ALL  SUBSEQUENT DISTRIBUTIONS OR 
#  REDISTRIBUTIONS  OF THE  SUBJECT  SOFTWARE.  ANYONE WHO USES, REPRODUCES, 
#  DISTRIBUTES, MODIFIES  OR REDISTRIBUTES THE SUBJECT SOFTWARE, AS DEFINED 
#  HEREIN, OR ANY PART THEREOF,  IS,  BY THAT ACTION, ACCEPTING IN FULL THE 
#  RESPONSIBILITIES AND OBLIGATIONS CONTAINED IN THIS AGREEMENT.
#  
#  Government Agency: National Aeronautics and Space Administration
#  Government Agency Original Software Designation: GSC-15354-1
#  Government Agency Original Software Title:  GEOS-5 GCM Modeling Software
#  User Registration Requested.  Please Visit http://opensource.gsfc.nasa.gov
#  Government Agency Point of Contact for Original Software:  
#  			Dale Hithon, SRA Assistant, (301) 286-2691
#  
# +-======-+ 

#######################################################################
#             Experiment Specific Environment Variables
#######################################################################

setenv  ARCH   `uname`
setenv  SITE    @SITE

setenv  EXPID   @EXPID
setenv  EXPDIR  @EXPDIR
setenv  HOMDIR  @HOMDIR
setenv GEOSUTIL @GEOSUTIL

#######################################################################
#                    Set Parameters for Forecasts
#
#    FCST_TYPE: STRAIGHT
#               Forecast using archived 21z Restarts including
#               AGCM_IMPORT (if available)
#
#    FCST_TYPE: REPLAY
#               Forecast using archived 21z Restarts, REPLAY to 
#               archived ANA.ETA to create internal AGCM_IMPORT
#
#    FCST_TYPE: CYCLED_REPLAY
#               Forecast using a single 21z Restart to initialze and cycle
#               subsequent restarts, REPLAY to archived ANA.ETA to create 
#               internal AGCM_IMPORT
#
#######################################################################

set    BDATE = $1
set    EDATE = $2
set FSEGMENT = $3
set     QSUB = $4

if( .$QSUB == . ) set QSUB = FALSE

set       NX = @NX
set       NY = @NY
set  AGCM_IM = @AGCM_IM
set  AGCM_JM = @AGCM_JM
set  AGCM_LM = @AGCM_LM

set HEARTBEAT = @DT
set       ANA = GMAO
set     NFMAX = 1

set    FCST_TYPE = STRAIGHT
set    ANA_EXPID = x0018
set ANA_LOCATION = /discover/nobackup/projects/gmao/g6dev/ltakacs/x0018

#######################################################################
#                   Create Forecast SubScripts
#######################################################################

     @ NPES   = $NX * $NY
     @ count  = 0
set     date  = $BDATE
while( $date <= $EDATE )

   @ NF = 1
   set     test  = $date
   while( $test <= $EDATE & $NF <= $NFMAX )
   set nymd = `$GEOSUTIL/post/tick $test 000000 86400`
   set test = $nymd[1]
   @ NF = $NF + 1
   end
   set nymde = `$GEOSUTIL/post/tick $test 000000 -86400`
   set nymde =  $nymde[1]

   @ NF = $NF - 1
   @ count = $count + 1
   echo $date $NF
   /bin/rm -f gcm_forecast.j$count
   /bin/rm -f sedfile
   /bin/rm -f gcm_forecast.tmp
   /bin/cp -f $EXPDIR/forecasts/gcm_forecast.tmpl $EXPDIR/forecasts/gcm_forecast.tmp
cat > sedfile << EOF
s?&NPES?$NPES?g
s?&DATE?$date?g
s?&AGCM_IM?$AGCM_IM?g
s?&AGCM_JM?$AGCM_JM?g
s?&NX?$NX?g
s?&NY?$NY?g
s?&HEARTBEAT?$HEARTBEAT?g
s?&NFMAX?$NF?g
s?&FCST_TYPE?$FCST_TYPE?g
s?&ANA_EXPID?$ANA_EXPID?g
s?&ANA_LOCATION?$ANA_LOCATION?g
s?&FSEGMENT?$FSEGMENT?g
s?FCSTOUT?$EXPDIR/forecasts/gcm_${FCST_TYPE}_forecast.o${date}-${nymde}?g
s?&ANA?$ANA?g
EOF
sed -f sedfile $EXPDIR/forecasts/gcm_forecast.tmp > $EXPDIR/forecasts/gcm_${FCST_TYPE}_forecast.j${date}-${nymde}
if( $QSUB != 'FALSE' )                        qsub  $EXPDIR/forecasts/gcm_${FCST_TYPE}_forecast.j${date}-${nymde}

set   date = $test
end

/bin/rm -f sedfile
/bin/rm -f $EXPDIR/forecasts/gcm_forecast.tmp