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

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

set oabgn = "2003 227"

set line1 = "AOSN-II  (Run 20)"
set line2 = "Synop(27.5,70) Mean(110.0,280)"
set line3 = "AugGDEM Infl(100,556,1000) Shap(2,1)"
set line4 = "unused"

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

set veclat = "42.5"
set veclon = "-136.5"

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

set pfile = "oa_gdem_aug_2000m.nc"
set sfile = "oa_gdem_aug_2000m.nc"

set palette = "/home/projects/ASCOT/HOPS/Plot/Palette/gs2.pal.white"
set pltdt = "/home/projects/ASCOT/HOPS/Plot/Data/plot_qg.dat"
set cstflg = "yes"
set cstfile = "/home/projects2/AOSN2/usW_i.dat"
set detflg  = "yes"
set details = "/home/projects2/AOSN2/OAG/aosn2_big9km2_sec.dat"
set gfile = "/home/projects2/AOSN2/Grids/Big9km2/grids_AOSN2_big_nois.nc"

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

set nfld = "4"
set fld = "11,12,306,5"
set nlvl = "11"
set lvl = "1,2,3,5,7,10,14,19,26,29,31"
set pmin = "0.0"
set pmax = "0.0"
set vscl = "25"
set iref = "0"
set idover = "10"
set nppg = "4"
set gxtnd = "ccnt"

/home/projects/PressTest/Bin_Sun5/oa_ccnt_job "${oabgn}" "${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}" "${cstflg}" "${cstfile}" "${detflg}" "${details}" ${gxtnd}

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

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 gxtnd = "csec_MBew"

/home/projects/PressTest/Bin_Sun5/oa_csec_job "${oabgn}" "${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/oa_csec_job "${oabgn}" "${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/oa_csec_job "${oabgn}" "${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}