#! /usr/bin/csh -f
#
# Make multiple runs of tgcmproc_f90.
# Copy this file up to IBM and submit with "llsubmit tgcmproc_multi.job"
# (Do NOT submit with the hao "submit" command.)
#
# @ step_name        = tgcmproc_wickwar
# @ class            = dv_reg
# @ wall_clock_limit = 06:00:00
# @ job_type         = serial
# @ notification     = complete
# @ output           = $(step_name).out
# @ error            = $(step_name).out
# @ queue
#
# Build input file and execute for each iteration
# of the while loop:
#
##@ day0 = 1
##@ day1 = 30
##@ day0 = 151
##@ day1 = 180
@ day0 = 301
@ day1 = 330
while ($day1 <= 400)
  set daybeg = $day0
  set dayend = $day1
  if ($day0 < 100) set daybeg = 0$day0
  if ($day1 < 100) set dayend = 0$day1
  if ($day0 < 10) set daybeg = 00$day0
  if ($day1 < 10) set dayend = 00$day1
  if ($day1 == 390) set dayend = 365
  echo daybeg = $daybeg dayend = $dayend
  ed << EOF
a
&proc_input                                                         
;
; Make ascii hourly TN profiles for Vince Wickwar.
; tgcm24r year 2002 run (see /hao/tgcm/tgcm24yr02/sy2n)
;
 tmpdir = '/hao/tgcm/tgcm24yr02/sy2n'
 histvols = '/ROBLE/tgcm24r/sy2n$daybeg.nc','to',
            '/ROBLE/tgcm24r/sy2n$dayend.nc','by','1'
 mtimes = $day0,0,0,'to',$day1,23,0,'by',60
 maxdiskvols = 35
 iprint_fldminmax=0	; print 3d min,max of fields every hist
 modelhts = 0
 cfields = 'TN'
 ipltxyloc=1		; flag for vertical profile xy plots
 xylocs = 41.74,-111.81
 xyloc_locname = 'BEAR LAKE'
 xyloc_htscale = 45.,90.,1.
 senddat = 
   'vishnu.hao:/e/foster/tgcmproc/wickwar/tgcm24r_sy2n_d$daybeg-$dayend.dat'
/
.
w tgcmproc.inp
q
EOF
#
# Execute:
#
#/fs/projects/tgcm/bin.AIX/tgcmproc_f90 < tgcmproc.inp >&! d$daybeg-$dayend.out
./tgcmproc_f90 < tgcmproc.inp >&! d$daybeg-$dayend.out
mv tgcmproc.dat tgcm24r_sy2n_d$daybeg-$dayend.dat
rm -f ROBLE*nc
@ day0 = $day0 + 30
@ day1 = $day1 + 30
end # while loop
exit