#!/bin/ksh
# @ job_type = parallel
## @ output = /ptmp/wx23ja/out170
## @ error = /ptmp/wx23ja/err170
# @ output = /ptmp/wx23ja/std28/out62
# @ error = /ptmp/wx23ja/std28/err62
## @ total_tasks = 42
## @ node = 21
## @ total_tasks = 1
## @ node = 1
# @ total_tasks = 28
## @ total_tasks = 32
# @ blocking = unlimited
## @ node = 8
## @ node = 32
## @ requirements = Feature == "dev"
# @ class = dev
## @ class = 1
## @ wall_clock_limit = 01:50:00
## @ wall_clock_limit = 00:10:00
# @ wall_clock_limit = 00:30:00
# @ notification = never
# @ network.MPI=csss,shared,us
# @ queue
set -x
#   Forecast begining date, and length and averageing dates
#
 export sdate=2002102400
 export CDATE=$sdate
 dt=`echo $CDATE | cut -c1-8`
echo " dt=$dt "
#
#    Initial condition directory
#
export MP_PULSE=0
export MP_TIMEOUT=600
#datic=/gloptmp/datprx
#datic=/ptmp/wx23sm/datprx
#datic=/com/fnl/prod/fnl.20010517
# datic=/com/avn/prod/avn.20030309
## datic=/com/fnl/prod/fnl.$dt
# datic=/com/avn/prod/avn.$dt
#datic=/gpfsuser/g01/wx23sm/praaug
#datic=/nfsuser/g01/wx23sm/data/para_exp/clw
datic=/nfsuser/g01/wx23ys/data.6228
#
#    Forecast model resolution
#
 export wave=62
#export wave=170
export lm=28
#export lm=42
# export lm=64
#
if [ $wave -eq 62 ] ; then
  export im=192
  export jm=94
  export iop=144
  export jop=73
fi
#
if [ $wave -eq 170 ] ; then
  export im=512
  export jm=256
  export iop=360
  export jop=181
fi
export R=${wave}${lm}
#
#
#   Control for post and averaging
#
 export post=YES
#export post=NO
export ipgb=1
export iflx=1
export index=1
#export avg=1
# export avg=0
export iclw=1
#
#    Directory in which the forecast is run
#
# export expt=datprx_r8
export expt=std28
#export expt=datpro
#export expt=clw_tt62b
export RUNDIR=/ptmp/wx23ja/$expt
export COMOUT=$RUNDIR
export SCRIPTS=/nfsuser/g01/wx23ja/gsm2003/scripts
export SCRIPTS2=/nfsuser/g01/wx23sm/gsm/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/wx23sm/exec_01/fcst6228_clw_ng5_3d.x
#export fname=/nfsuser/g01/wx23sm/exec_03/fcst${R}_clw_og_ras1.x
# export fname=/nfsuser/g01/wx23sm/exec_03/global_fcst${R}a.x
#export fname=/nfsuser/g01/wx23sm/exec_04/global_fcst${R}
#export fname=/nfsuser/g01/wx23sm/exec_03/global_fcst${R}_rrtm.x
#export fname=/nwprod/exec/global_fcst${R}
#export fname=/nfsuser/g01/wx23sm/exec_01/global_fcst6228
#export fname=/nfsuser/g01/wx23sm/exec_02/global_fcst6228
#export fname=/nfsuser/g01/wx20mi/2jif/global_fcst.fd/global_fcst6228rg
mkdir -p $RUNDIR
cd $RUNDIR
#
# cp /nfsuser/g01/wx23sm/gsm/sfccyc/aeros_data/aeropac3a* $RUNDIR
# cp /nwprod/fix/global_aeropac3a* $RUNDIR
cp /nfsuser/g01/wx23ys/data_fix/aeropac3a* $RUNDIR
#
#   Forecast begining date, and length and averageing dates
#
#export sdate=1997121500
#
#                          Summer 1998 case
#export sdate=1998061500
#export sdate_avg=1998070100
#export edate_avg=1998073112
#
#                          Summer 2000 case
#export sdate=2000061500
#export sdate_avg=2000070100
#export edate_avg=2000073112
#export cycle_date=2000060100.2000080100
#                          Winter 2001 case
#export sdate=2000121500
#export sdate_avg=2001010100
#export edate_avg=2001013112
#export cycle_date=2000121500.2001020100
#
#export sdate=2000012300
#export sdate_avg=2000020100
#export edate_avg=2000030100
#export cycle_date=2000011900.2000031100
#
#                          Just a Random 2001 case
#export sdate=2001042000
 export sdate_avg=2002102400
 export edate_avg=2002102900
#
#
 export ndays=5
#export ndays=39
#export ndays=48
#export nhourb=2000
#export nhourb=72
#export nhourb=48
 export nhourb=0
 export CDATE=$sdate
#
#
nhours=`expr $ndays \* 24`
#nhours=168
#nhours=12
#
export fmax=$nhours
 export fout=3
 export fzer=3
#export fout=12
#export fzer=12
#export fout=6
#export fzer=6
#export fcyc=0
export fcyc=0
# fcyc=0 cycle not called - inital values used
# export fcyc=24
#
# set fdfi=0 for no digital filter  
export fdfi=0
# export fdfi=3
#export fdfi=12
#
#export OROGRAPHY=/gloptmp/fix/global_orography.t${wave}rg1.grb
#export SLMASK=/gloptmp/fix/global_slmask.t${wave}rg1.grb
## export OROGRAPHY=/nwprod/fix/global_orography.t${wave}.grb
export OROGRAPHY=/nfsuser/g01/wx23ys/t6228/fix/global_orography.t62rg.grb
export SLMASK=/nwprod/fix/global_slmask.t${wave}.grb
#export CO2CON=/nfsuser/g01/wx20mi/newres/global_co2con.l64.f77
# above for l64
#export CO2CON=/gpfstmp/wx23sm/clm1/co2l28.sp.new
## export CO2CON=/nwprod/fix/global_co2con.l${lm}.f77
export CO2CON=/nfsuser/g01/wx23ys/t6228/fix/co2l28
#

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='/nwprod/fix'
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/gdas1.t00z.sanl
#  export SFCINP=$datic/gdas1.t00z.sfcanl
#
# ------------------------------------
## export SIGINP=$datic/gdas1.t00z.sanl
## export SFCINP=$datic/gdas1.t00z.sfcanl
# ------------------------------------
export SIGINP=$datic/siganl
export SFCINP=$datic/sfcanl
echo $SIGINP
ls -l $SIGINP
#
# export SIGINP=$datic/siganl.$dtg.new
# export SFCINP=$datic/sfcanl.$dtg.new
#
# 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 -  number of cloud species: set to 1
##   export CHGRESVARS="ntrac=3,ncldt=1"
# export CHGRESVARS="ntrac=2"
#  
# export SIGLEVEL=/nfsuser/g01/wx20mi/newres/global_siglevel.l64.txt
# CHGRES chgres points to ...l64 for 64 levels
##  /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
#
  export FCSTEXEC=$fname
# export FCSTEXEC="pedb  $fname"
# export FCSTEXEC="xldb  $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 D3DO='d3df$(/nwprod/util/exec/ndate $FH $CDATE)'
  export FIXGLOBAL=/nwprod/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=/nwprod/fix/global_mtnvar.t${wave}.f77
#  export MTNVAR=/emc1/wx23ja/terr01/ml/mtnvar14_${wave}
export MTNVAR=/nfsuser/g01/wx23ja/data/terr01/mtnvar14_${wave}
ls -l $MTNVAR
#
#  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

# sst, ice and snow analyses file for cycle
#  export FNTSFA=/nfsuser/g01/wx23sm/data/sfcanl/sstgrb.$cycle_date
#  export FNACNA=/nfsuser/g01/wx23sm/data/sfcanl/icegrb.$cycle_date
# export FNTSFA=/nfsuser/g01/wx23sm/data/sfcanl/sstgrb.2000060100.2000080100
# export FNACNA=/nfsuser/g01/wx23sm/data/sfcanl/icegrb.2000121500.2001020100
#
   export FNTSFA=
   export FNACNA=
#
  export EXECGLOBAL=$RUNDIR
  export FHOUT=$fout
  export FHZER=$fzer
  export FHCYC=$fcyc
  export FHDFI=$fdfi
  export FHLWR=3
  export FHSWR=1
  export FHMAX=$fmax
# 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=450.
# export FCSTVARS='emc=.true.'
  export FCSTVARS='LDIAG3D=.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/wx23sm/gsm/scripts/global_forecast.sh_3d
#  /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
echo " sub -e $listvar -q dev -p 4/4 -t 1:00:00 -o $RUNDIR/post.$CDATE.out $SCRIPTS2/global_post_mpi "
  sub -e $listvar -q dev -p 4/4 -t 1:00:00 -o $RUNDIR/post.$CDATE.out $SCRIPTS2/global_post_mpi
#  fi