#!/bin/ksh
# @ job_type = parallel
# @ output = /ptmp/wx23ja/out_170_era
# @ error = /ptmp/wx23ja/err170_era
## @ output = /ptmp/wx23ja/out_62_er
## @ error = /ptmp/wx23ja/err62_er
# @ total_tasks = 42
# @ node = 14
## @ total_tasks = 1
## @ node = 1
## @ total_tasks = 28
## @ node = 14
## @ total_tasks = 64
## @ node = 32
# @ preferences = Feature == "dev"
## @ requirements = Feature == "dev"
# @ class = dev
## @ node = 10
## @ class = 1
# @ wall_clock_limit = 01:50:00
## @ wall_clock_limit = 03:30:00
# @ notification = never
# @ network.MPI=switch,shared,us
# @ queue
set -x
#
#    Initial condition directory
#
export MP_PULSE=0
export MP_TIMEOUT=600
#datic=/gloptmp/datprx
#datic=/ptmp/wx23sm/datprx
#datic=/nfsuser/g01/wx23ss/cold/colddata
#datic=/com/fnl/prod/fnl.20010517
 datic=/com/avn/prod/avn.20030101
#datic=/gpfsuser/g01/wx23sm/praaug
#datic=/gpfsuser/g01/wx23sm/para_exp/clw
#datic=/nfsuser/g01/wx23sm/data/para_exp/clw
#datic=/ptmp/wx23sm/clwd
#
#    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=datpro
#export expt=clw_r1
#export expt=clw_tt62b
export RUNDIR=/ptmp/wx23sm/$expt
export COMOUT=$RUNDIR
export SCRIPTS=/nfsuser/g01/wx23sm/gsm/scripts
#export FCSMEAN=/gpfsuser/g01/wx23sm/data
export FCSMEAN=/nfsuser/g01/wx23sm/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/wx23sm/exec/fcst6228_clw_og1.x
#export fname=/nfsuser/g01/wx23sm/exec/fcst6228_clw_ng_3d.x
#export fname=/nfsuser/g01/wx23sm/exec/fcst6228_clw_ogem3d.x
#export fname=/nfsuser/g01/wx23sm/exec/fcst6228_clw_ogem3d.x
#export fname=/nfstmp/wx23sm/clw3d/fcst6228_clw_ogem3d.x
#export fname=/nfsuser/g01/wx23sm/exec/fcst6228_clw_og_test1.x
#export fname=/nfsuser/g01/wx23sm/exec/fcst6228_clw_og_tstyh2.x
#export fname=/nfsuser/g01/wx23sm/exec_01/fcst6228_clw_ng0_3d.x
#export fname=/nfsuser/g01/wx23sm/exec_01/fcst6228_clw_ng5_3d.x
#export fname=/nfsuser/g01/wx23sm/exec_01/fcst${R}_clw_ng5_3dt2.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/wx23sm/exec_01/fcst6228_clw_ng5_3dto.x
#export fname=/nfsuser/g01/wx23sm/exec/fcst6228_clw_og_tstyh.x
#export fname=/nfsuser/g01/wx23sm/exec/fcst6242_clw_og_tstyh.x
#export fname=/nfsuser/g01/wx23sm/exec/fcst6228_clw_ng_tstyh2.x
#export fname=/nfsuser/g01/wx23sm/exec/fcst6228_og.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=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=2001052100
 export sdate_avg=2000020100
 export edate_avg=2000030100
#
#
 export ndays=3
#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=24
 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 SLMASK=/nwprod/fix/global_slmask.t${wave}.grb
#export CO2CON=/nfsuser/g01/wx20mi/newres/global_co2con.l64.f77
#

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/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"
# export CHGRESVARS="ntrac=2"
# export SIGLEVEL=/nfsuser/g01/wx20mi/newres/global_siglevel.l64.txt
#
# /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.186:0
# export DISPLAY=140.90.98.139: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
#
# 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=
  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
  sub -e $listvar -q dev -p 4/4 -t 1:00:00 -o $RUNDIR/post.$CDATE.out $SCRIPTS/global_post_mpi
 fi