#!/bin/ksh
# @ job_type = parallel
# @ output = /ptmp/wx23ja/fcst62/out
## @ output = /ptmp/wx23ja/fcst62.out
# @ error = /ptmp/wx23ja/fcst62/err
## @ total_tasks = 42
## @ node = 21
# @ total_tasks = 28
# @ node = 14
# @ preferences = Feature == "dev"
## @ requirements = Feature == "dev"
# @ class = dev
## @ node = 10
## @ class = 1
## @ wall_clock_limit = 00:10:00
# @ wall_clock_limit = 03:30:00
# @ notification = never
# @ network.MPI=switch,shared,us
# @ queue
set -x
#
#    Forecast model resolution
#
export wave=62
export im=192
export jm=94
export lm=64
#export lm=42
export iop=144
export jop=73
#
#   Control for post and averaging
#
 export post=YES
#export post=NO
export ipgb=1
export iflx=1
export index=1
export avg=1
export iclw=1
#
#    Directory in which the forecast is run
#
export expt=vec_mbja
export RUNDIR=/ptmp/wx23ja/$expt
export COMOUT=$RUNDIR
export SCRIPTS=/nfsuser/g01/wx23ja/gsm2003/scripts
#export FCSMEAN=/gpfsuser/g01/wx23sm/data
export FCSMEAN=/nfsuser/g01/wx23ja/data
#export fname=/ptmp/wx23sm/exec/f6228_$expt.x
#export fname=/ptmp/wx23sm/exec/f6228_clw0.x
#export fname=/nfsuser/g01/wx23sm/exec/fcst6228_clw_og3.x
 export fname=/nfsuser/g01/wx23ja/exec/gf$wave$lm${expt}.x
#export fname=/nfsuser/g01/wx20mi/2jif/global_fcst.fd/global_fcst6228rg
#export fname=/ptmp/wx23sm/exec/f6242_clw2.x
#export fname=/nfstmp/wx23sm/clw0/f6228_$expt.x
mkdir -p $RUNDIR
cd $RUNDIR
#
cp /nfsuser/g01/wx23sm/gsm/sfccyc/aeros_data/aeropac3a* $RUNDIR
#
#   Forecast begining date, and length and averageing dates
#
 export sdate=2003022700
 export CDATE=$sdate
 dt=`echo $CDATE | cut -c1-8`
echo " dt=$dt "
#
#    and Initial condition directory
#
export MP_PULSE=0
export MP_TIMEOUT=600
#datic=/gloptmp/datsum
#datic=/gpfsuser/g01/wx23ss/cold/colddata
#datic=/com/fnl/prod/fnl.19990826
#datic=/gpfsuser/g01/wx23sm/praaug
#datic=/gpfsuser/g01/wx23sm/para_exp/clw
#  datic=/nfsuser/g01/wx23sm/data/para_exp/clw
datic=/com/fnl/prod/fnl.$dt
export ndays=2
#  export ndays=39
#export ndays=48
#export nhourb=1000
 export nhourb=0
#export sdate_avg=1998070100
#export edate_avg=1998073112
#
##  export sdate_avg=2000020100
##  export edate_avg=2000020100
 export sdate_avg=2003022700
 export edate_avg=2003022800
##  export edate_avg=2000030100
#
nhours=`expr $ndays \* 24`
#nhours=168
#nhours=12
#
export fmax=$nhours
#export fout=1
#export fzer=3
 export fout=12
 export fzer=12
#export fout=6
#export fzer=6
#export fcyc=0
 export fcyc=24
#
#export OROGRAPHY=/gpfsuser/g01/wx23sm/topo/global_orography.t${wave}.grb
#export SLMASK=/gpfsuser/g01/wx23sm/topo/global_slmask.t${wave}.grb
#
export OROGRAPHY=/gloptmp/fix/global_orography.t${wave}rg1.grb
export SLMASK=/gloptmp/fix/global_slmask.t${wave}rg1.grb
#

export listvar=CDATE,COMOUT,FCSMEAN,SCRIPTS,fexec,fname,fmax,fout,listvar,fzer,wave,im,jm,lm,iop,jop,OROGRAPHY,SLMASK,fcyc,RUNDIR,sdate_avg,edate_avg,expt,post,index,avg,ipgb,iflx,iclw
#
export JCAP=$wave
export LEVS=$lm
export LONB=$im
export LATB=$jm
export VERBOSE=YES
export FIXDIR='/gloptmp/fix'
#export ANADIR='/gpfsuser/g01/wx23sm/sfcanl'
export ANADIR='/nfsuser/g01/wx23sm/data/sfcanl'
cp /nfsuser/g01/wx20mi/scripts/sighdr .
#
if [ $nhourb -eq 0 ] ; then
  dtg=`echo $CDATE | cut -c3-10`
#
export SIGINP=$datic/gdas2.t00z.sanl
export SFCINP=$datic/gdas2.t00z.sfcanl
#
# export SIGINP=$datic/siganl.$dtg
# export SFCINP=$datic/sfcanl.$dtg
#
# export SIGINP=$datic/siganl.$CDATE
# export SFCINP=$datic/sfcanl.$CDATE
#
  export SIGOUT=$RUNDIR/siganl.$CDATE
  export SFCOUT=$RUNDIR/sfcanl.$CDATE
#
  export CHGRESVARS="ntrac=3,ncldt=1"
# export CHGRESVARS="ntrac=2"
#
# /nfsuser/g01/wx20mi/scripts/global_chgres.sh
cp $SIGINP $SIGOUT
cp $SFCINP $SFCOUT
#
  export SIGI=$RUNDIR/siganl.$CDATE
  export SFCI=$RUNDIR/sfcanl.$CDATE
  export FHINI=0$nhourb
fi
if [ $nhourb -lt $nhours ] ; then
  if [ $nhourb -gt 0 ] ; then
    cp sigr1 sigr1a
    cp sigr2 sigr2a
    cp sfcr  sfcra
    export SIGI=$RUNDIR/sigr1a
    export SIGI2=$RUNDIR/sigr2a
    export SFCI=$RUNDIR/sfcra
    export FHINI=$nhourb
  fi
#
#  Boundary analyses files for cycle
#
# export sstgrb=$ANADIR/sst9899
# export snogrb=$ANADIR/snow1999
  export snogrb='        '
# export engicegrb=$ANADIR/ice9899
  export sstgrb=$ANADIR/sst980607
  export engicegrb=$ANADIR/ice980607
#
  export FCSTEXEC=$fname
# export FCSTEXEC="pedb  $fname"
# export FCSTEXEC="totalview poe -a  $fname"
  export DISPLAY=140.90.192.158:0
# export SIGHDR=$RUNDIR/sighdr
  export SIGO='sigf$(/nwprod/util/exec/ndate $FH $CDATE)'
  export SFCO='sfcf$(/nwprod/util/exec/ndate $FH $CDATE)'
  export FLXO='flxf$(/nwprod/util/exec/ndate $FH $CDATE)'
  export LOGO='logf$(/nwprod/util/exec/ndate $FH $CDATE)'
  export FIXGLOBAL=/gloptmp/fix
#
# cp /nfsuser/g01/wx23ja/terrain/s99terrain/data/mtn.T${wave}.ieee .
# export MTNVAR=mtn.T${wave}.ieee
# export MTNVAR=/gpfsuser/g01/wx23sm/topo/global_mtnvar.t${wave}.f77
#   export MTNVAR=/gloptmp/fix/global_mtnvar.t${wave}rg1.f77
  export MTNVAR=/nfsuser/g01/wx23ja/data/terr01/ml/mtn.T${wave}.ieee
#
# sst, ice and snow analyses file for cycle
  export FNTSFA=/nfsuser/g01/wx23sm/data/sfcanl/sstgrb.2000011900.2000031100
  export FNACNA=/nfsuser/g01/wx23sm/data/sfcanl/icegrb.2000011900.2000031100
# export FNTSFA=/gpfsuser/g01/wx23sm/sfcanl/sstgrb.2000011900.2000031100
# export FNACNA=/gpfsuser/g01/wx23sm/sfcanl/icegrb.2000011900.2000031100
# export FNSNOA=/gpfsuser/g01/wx23sm/sfcanl/snogrb.2000011900.2000031100
#
  export EXECGLOBAL=$RUNDIR
  export FHOUT=$fout
  export FHZER=$fzer
  export FHCYC=$fcyc
  export FHLWR=3
  export FHSWR=1
  export FHMAX=$nhours
# export FHMAX=$((45*24))
  export SIGR1=$RUNDIR/sigr1
  export SIGR2=$RUNDIR/sigr2
  export SFCR=$RUNDIR/sfcr
# export FHRES=$((6))
  export FHRES=$((1*24))
  export DELTIM=1200.
# export DELTIM=600.
  export VERBOSE=YES
# export FCSTVARS='ras=.true.'
# echo $FCSTVARS
#
  export FTSFS=0.0
  export FAISS=0.0
# export FSNOL=0.0
  export CYCLVARS='FALBL=0.0,FALBS=0.0,'
#
#  /nfsuser/g01/wx23sm/gsm/scripts/global_forecast.sh_cyc
#  /nfsuser/g01/wx20mi/scripts/global_forecast.sh3
#  prof -s -m mon.out.? mon.out.?? > tem
fi
if [ $sdate_avg -lt $edate_avg ] ; then
  #llsubmit /nfsuser/g01/wx23sm/gsm/scripts/post62_mpi3
 sub -e $listvar -q 1 -p 4/1 -t 1:00:00 -o $RUNDIR/post.$CDATE.out $SCRIPTS/global_post_mpi
fi