#
# 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 zmslt     = '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 "cfields = 'TN','UN','VN','O2','O1','N2','N4S','NO','OP','N2D'," >> $outfile
echo "          'TI','TE','NE','O2P','OMEGA','POTEN','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 = 'ps'"                               >> $outfile
echo ";"                                           >> $outfile
echo "; Local ps output file:"                     >> $outfile
echo "flnm_ps = 'tgcmproc.ps'"                     >> $outfile
echo ";"                                           >> $outfile
echo "; Local ascii output file:"                  >> $outfile
echo "flnm_dat = 'tgcmproc.dat'"                   >> $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 = 'zm'"                            >> $outfile
  echo "  fslts = 0.,12."                          >> $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 = -60.,-40.,0.,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 = 'zm','zm'"                      >> $outfile
  echo "  xyloc_locname = 'Global Mean'"           >> $outfile
  echo "  xyloc_zprange = -999.,999."              >> $outfile
  echo "; xyloc_htscale = 50.,500.,10."            >> $outfile
endif
if ($zmslt == 'y') then
  echo ";"                                               >> $outfile
  echo "; Make lat vs zp/ht plots of zm at selected slt" >> $outfile
  echo "; (must have 1 day time series):"                >> $outfile
  echo ";"                                               >> $outfile
  echo "ipltzmslt = 0"                                   >> $outfile
  echo "; zmslt = 0.,6.,12.,18."                         >> $outfile
  echo "zmslt = 0.,12."                                  >> $outfile
  echo "zmslt_zprange = -999.,999."                      >> $outfile
  echo "zmslt_htscale = 90.,400.,5."                     >> $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 = -60.,180., -40.,180., 0.,180., 40.,180., 60.,180." >> $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.,180., -4.,180., 0.,180., 2.,180., 6.,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 = 'zm','zm'" >> $outfile
  echo "; xyut_locname = 'LOC 1','LOC 2','LOC 3','LOC 4'" >> $outfile
  echo "  xyut_zpht = -7.,-4.,0.,2.,6."            >> $outfile
endif
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
