#
# This script helps "automate" the production of plots as well as simplify
# the process of getting uniform labels.
#
# To access the NCAR escape sequences, use \$ in place of $ and use single 
# quotes.  Ex:
#
#    set line2 = 'Surface Temperature  (\$S1PGU1\$Z\$PRU\$ C)'

#-------------------------------------
#--- Set fairly common parameters. ---
#-------------------------------------

set pebgn = "2003 227"

set line1 = "AOSN-II  (Run 50)"
set line2 = "Synop(27.5,70) Mean(110.0,280)"
set line3 = "August GDEM Climatology; Stab"
set line4 = "unused"

set day0 = "12867.0"
set dayN = "12867.0"
set dday = "1.0"

set veclat = "43.0"
set veclon = "-137.0"

set gridres = "3"
set latlim = "  25.09792519    44.39816666"
set lonlim = "-137.86181641  -115.27670288"

set pfile = "pi_cc_ini.nc"
set sfile = "pi_cc_ini.nc"

set palette = "/home/projects/ASCOT/HOPS/Plot/Palette/gs2.pal.white"
set pltdt = "/home/projects2/ASCOT02/PE/plot_pe_bigarrow.dat"
set pltcst = "yes"
set coast = "/home/projects2/AOSN2/usW_i.dat"
set pltdtl  = "yes"
set details = "/home/projects2/AOSN2/OAG/aosn2_big9km2_sec.dat"
set gfile = "/home/projects2/AOSN2/Grids/Big9km2/grids_AOSN2_big_noishcfsm.nc"

#-----------------------------------
#--- Standard multi-field plots. ---
#-----------------------------------

set nfld = "4"
set fld = "5,6,4,10"
set nlvl = "9"
set lvl = "1 -10 -50 -100 -200 -500 -1000 -3000 20"
set pmin = "0.0"
set pmax = "0.0"
set vscl = "50"
set iref = "0"
set idover = "10"
set nppg = "4"
set gxtnd = "ccnt"

/home/projects/PressTest/Bin_Sun5/pe_ccnt_job "${pebgn}" "${line1}" "${line2}" \
                                     "${line3}" "${line4}" ${nfld} \
            "${fld}" ${nlvl} "${lvl}" ${day0} ${dayN} ${dday} ${pmin} ${pmax} \
            ${vscl} ${veclat} ${veclon} ${iref} ${idover} ${gridres} ${nppg} \
            "${latlim}" "${lonlim}" "${palette}" "${pltdt}" "${pfile}" \
            "${sfile}" "${pltcst}" "${coast}" "${pltdtl}" "${details}" ${gxtnd}

#------------------------
#--- MBAY EW Section. ---
#------------------------

set nfld = "4"
set fld = "5,6,91,92"
set lon1 = "-133.0"
set lat1 = "32.0"
set lon2 = "-122.0"
set lat2 = "37.0"
set sztop = "0.0"
set szbot = "5500.0"
set kms = "150"
set nppg = "4"
set gxtnd = "csec_MBew"

/home/projects/PressTest/Bin_Sun5/pe_csec_job "${pebgn}" "${line1}" \
                                     "${line2}" "${line3}" "${line4}" ${nfld} \
            "${fld}" ${day0} ${dayN} ${dday} ${lon1} ${lat1} ${lon2} ${lat2} \
           ${sztop} ${szbot} ${kms} ${pmin} ${pmax} ${iref} ${idover} ${nppg} \
            "${palette}" "${pltdt}" "${pfile}" "${sfile}" "${gfile}" ${gxtnd}

#----------------------------
#--- Northern EW Section. ---
#----------------------------

set lon1 = "-136.0"
set lat1 = "35.5"
set lon2 = "-124.0"
set lat2 = "41.5"
set szbot = "5500.0"
set gxtnd = "csec_New"

/home/projects/PressTest/Bin_Sun5/pe_csec_job "${pebgn}" "${line1}" \
                                     "${line2}" "${line3}" "${line4}" ${nfld} \
            "${fld}" ${day0} ${dayN} ${dday} ${lon1} ${lat1} ${lon2} ${lat2} \
           ${sztop} ${szbot} ${kms} ${pmin} ${pmax} ${iref} ${idover} ${nppg} \
            "${palette}" "${pltdt}" "${pfile}" "${sfile}" "${gfile}" ${gxtnd}

#----------------------------
#--- Southern EW Section. ---
#----------------------------

set lon1 = "-128.0"
set lat1 = "26.0"
set lon2 = "-117.0"
set lat2 = "32.0"
set szbot = "5500.0"
set gxtnd = "csec_Sew"

/home/projects/PressTest/Bin_Sun5/pe_csec_job "${pebgn}" "${line1}" \
                                     "${line2}" "${line3}" "${line4}" ${nfld} \
            "${fld}" ${day0} ${dayN} ${dday} ${lon1} ${lat1} ${lon2} ${lat2} \
           ${sztop} ${szbot} ${kms} ${pmin} ${pmax} ${iref} ${idover} ${nppg} \
            "${palette}" "${pltdt}" "${pfile}" "${sfile}" "${gfile}" ${gxtnd}