#
# Make a sample tgcmproc.inp template file.
#
set outfile = tgcmproc_example.inp
echo ' '
if (-e $outfile) then
  echo -n "Overwrite $outfile? (y/n): "
  set ans = $<
  if ($ans == 'y') then
    rm $outfile
    echo "Overwriting $outfile"
  else
    echo -n "Enter file name to create: "
    set outfile = $<
    echo "Will create file $outfile"
  endif
endif
touch $outfile
#
set maps      = 'y'
set lonslices = 'y'
set latslices = 'y'
set xyloc     = 'y'
set utvert    = 'y'
set utlat     = 'y'
set xyut      = 'y'
set tgcmdata = '$TGCMDATA'
set mydata = '$MYDATA'
set logname_uc = `echo $LOGNAME | tr "[a-z]" "[A-Z]"`
#
echo ";"                                           >> $outfile
echo "; This sample input file for tgcmproc_f90 was generated by " >> $outfile
echo "; $0 at `date`"                              >> $outfile
echo ";"                                           >> $outfile
echo "&proc_input"                                 >> $outfile
echo ";"                                           >> $outfile
echo "; History file(s) to be read:"               >> $outfile
echo "histvols = '$tgcmdata/TGCM.tiegcm.p001-2002-080.nc'" >> $outfile
echo ";"                                           >> $outfile
echo "; Model time(s) (ddd,hh,mm):"                >> $outfile
echo "mtimes = 80,0,0"                             >> $outfile
;
echo ";"                                           >> $outfile
echo "; Examples of history file lists (env vars may be imbedded in paths):" >> $outfile
echo ";histvols = '$TGCMDATA/TGCM.tiegcm.p001-2002-080.nc'"            >> $outfile
echo ";histvols = '/mydir/dir1/file_001.nc','/mydir/dir1/file_002.nc'" >> $outfile
echo ";histvols = '$mydata/file_001.nc','$mydata/file_002.nc'"         >> $outfile
echo ";histvols = '/$logname_uc/dir/file001.nc','to','/$logname_uc/dir/file008.nc','by','1'" >> $outfile
;
echo ";"                                           >> $outfile
echo "; Examples of multiple model-time requests:" >> $outfile
echo ";mtimes = 80,0,0, 80,6,0, 80,12,0    ; day 80, hours 0, 6, and 12" >> $outfile 
echo ";mtimes = 80,0,0,'to',81,0,0,'by',60 ; day 80 to 81 every 60 minutes" >> $outfile
echo "maxdiskvols = 365"                           >> $outfile
echo "cfields = 'TN','UN','VN','O2','O1','N2','Z'" >> $outfile
echo ";"                                           >> $outfile
echo "; modelhts = 0 -> calculate geopotential height (non-constant gravity)" >> $outfile
echo "; modelhts = 1 -> use geopotential height on the history" >> $outfile
echo ";" >> $outfile
echo "modelhts = 0"                                >> $outfile
echo "icolor = 1"                                  >> $outfile
echo ";"                                           >> $outfile
echo "; iden=0 -> leave species as on history (most are mass mix ratios)" >> $outfile
echo "; iden=1 -> convert species to number densities (cm3)"              >> $outfile
echo "; iden=2 -> convert species to number density mixing ratios"        >> $outfile
echo "; iden=3 -> convert species to mass density (gm/cm3)"               >> $outfile
echo ";"                                           >> $outfile
echo "iden = 1"                                    >> $outfile
echo "iplot = 1"                                   >> $outfile
echo "outplt = 'cgm','ps'"                         >> $outfile
#
if ($maps == 'y') then
  echo ";"                                         >> $outfile
  echo "; Make maps (global and/or polar):"        >> $outfile
  echo ";"                                         >> $outfile
  echo "ipltmaps = 1"                              >> $outfile
  echo "  map_global = 1"                          >> $outfile
  echo "  map_polar = 0"                           >> $outfile
  echo "  fmap_polar_perimlat = -40.,40."          >> $outfile
  echo "  fmap_zpht = -7., -4., 0., 2."            >> $outfile
endif
if ($lonslices == 'y') then
  echo ";"                                         >> $outfile
  echo "; Make longitude slices (lat vs zp/ht at selected lons):" >> $outfile
  echo ";"                                         >> $outfile
  echo "ipltlon = 1"                               >> $outfile
  echo "  flons = 0., 180.,'zm'"                   >> $outfile
  echo "  fslts = 0.,6.,12.,18."                   >> $outfile
  echo "  flon_zprange = -99.,99."                 >> $outfile
  echo "; flon_htscale = 75.,150.,5."              >> $outfile
endif
if ($latslices == 'y') then
  echo ";"                                         >> $outfile
  echo "; Make latitude slices (lon vs zp/ht at selected lats):" >> $outfile
  echo ";"                                         >> $outfile
  echo "ipltlat = 1"                               >> $outfile
  echo "  flats = -40.,-60.,40.,60."               >> $outfile
  echo "  flat_zprange = -99.,99."                 >> $outfile
  echo "; flat_htscale = 75.,150.,5."              >> $outfile
endif
if ($xyloc == 'y') then
  echo ";"                                         >> $outfile
  echo "; Make xy location line plots (vertical profiles):" >> $outfile
  echo ";"                                         >> $outfile
  echo "ipltxyloc = 1"                             >> $outfile
  echo ";"                                         >> $outfile
  echo "; xylocs = lat,lon locations"              >> $outfile
  echo "; Note 'zm' means plot zonal means"        >> $outfile
  echo "; Note 'zm','zm' means plot global means"  >> $outfile
  echo ";"                                         >> $outfile
  echo "  xylocs = 40.,'zm',60.,'slt12.',-60.,180.,'zm','zm'"   >> $outfile
  echo "  xyloc_locname = 'loc1','loc2','Global Mean'" >> $outfile
  echo "  xyloc_zprange = -999.,999."                >> $outfile
  echo "; xyloc_htscale = 50.,500.,10."              >> $outfile
endif
if ($utvert == 'y') then
  echo ";"                                         >> $outfile
  echo "; Make ut vs zp/ht contours at selected locations:" >> $outfile
  echo ";"                                         >> $outfile
  echo "ipltutvert = 0"                            >> $outfile
  echo "  utvert_locs = 2.0,-157., 20.8,-156.5, 41.9,-111.4, 69.6,19.2" >> $outfile
  echo "  utvert_locname = 'CHRISTMAS','MAUI','BEAR LAKE','EISCAT'"     >> $outfile
  echo "  utvert_zprange = -999.,999."             >> $outfile
  echo "; utvert_htscale = 50.,100.,1."            >> $outfile
endif
if ($utlat == 'y') then
  echo ";"                                         >> $outfile
  echo "; Make ut vs latitude contours at selected zp/ht and longitudes" >> $outfile
  echo ";"                                         >> $outfile
  echo "ipltutlat = 0"                             >> $outfile
  echo "  utlat_zphtlon = -7.,0., -7.,180., 2.,0., 2.,180." >> $outfile
endif
if ($xyut == 'y') then
  echo ";"                                         >> $outfile
  echo "; Make ut vs field line plots at selected locations:" >> $outfile
  echo ";"                                         >> $outfile
  echo "ipltxyut = 0"                             >> $outfile
  echo "  xyut_locs = 40.,-90., 40.,'zm', 40.,'slt12.', 'zm','zm'" >> $outfile
  echo "; xyut_locname = 'LOC 1','LOC 2','LOC 3','LOC 4'" >> $outfile
  echo "  xyut_zpht = -4.,120.,'hmf2'"            >> $outfile
endif
echo ";"                                           >> $outfile
echo "; Example destinations for scp of cgm and ps output plot files:" >> $outfile
echo "; (If sendcgm and sendps are not provided, plots will be output" \
                                                   >> $outfile
echo ";   to local files tgcmproc.cgm and tgcmproc.ps)" >> $outfile
echo ";"                                           >> $outfile
echo ";sendcgm = 'machine.hao:/e/$LOGNAME/tgcmproc/tgcmproc.cgm'" >> $outfile
echo ";sendps  = 'machine.hao:/e/$LOGNAME/tgcmproc/tgcmproc.ps'"  >> $outfile
echo "/"                                           >> $outfile
#
echo ' '
echo "Created file $outfile."
#echo "See hao:/home/tgcm/tgcmproc/master.inp for explanations and additional options."
echo "Please edit $outfile for your own needs before executing." 
echo ' '
exit
