#!/bin/csh -f
#
# 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 dirstr = $cwd:t

set pebgn = "2003 231.0"

set line1 = "PhilEx (${dirstr})"
set line2 = "FsPs + Top(8)+Assi(Expl+IC)+Tides; BC(7373)"
set line3 = "WOA-05+PbOp+TpSm3sm+45Vlvl; DenExt(1900m)"
set line4 = "unused"

set day0 = "5.0"
set dayN = "20.0"
set dday = "5.0"

set veclat = "19.0"
set veclon = "115.2"

set gridres = "2"
set latlim = "  1.58213854  20.82702827"
set lonlim = "112.81455231 132.85191345"

set pfile = "/projects/philex/PE/2008/Aug23/arch11/pe_out.nc"
set sfile = "/projects/philex/PE/2008/Aug23/arch11/pe_out.nc"

set palette = "/home/phaley/PLUSNet07/HOPS/Plot/Palette/gs2.pal.white_all"
set pltdt = "/projects/philex/HOPS/Plot/Data/plot_pe_25aug2007.dat"
set pltcst = "yes"
set coast = "/home/phaley/PhilEx/Data/Coasts/archipelago_h.dat"
set pltdtl  = "no"
set details = "/projects/philex/PE/mindoro_pac.dat"
set gfile = "/projects/philex/Grids/Rsph_CSS/Arch8/grids_Archipelago9.nc"
set pltdft = "no"
set drifters = "driftout.dat"

#----------------------------
#--- Surface Temperature. ---
#----------------------------

set oline1 = "${line1}"
set oline2 = "${line2}"
set oline3 = "${line3}"
set oline4 = "${line4}"

set line1 = " "
set line3 = " "
set line2 = " "
set line4 = "20m Temperature"
set nfld = "1"
set fld = "5"
set nlvl = "1"
set lvl = "-20"
set pmin = "27.5"
set pmax = "29.5"
set vscl = "30"
set iref = "1"
set idover = "10"
set nppg = "1"
set gxtnd = "ccntT020_${dirstr}"
#set pltdtl  = "yes"

/home/phaley/PLUSNet07/Bin/pe_ccnt_job_at_nl "${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}

../gmetconvjob ${gxtnd}

#----------------------------
#--- Upper 200m velocity. ---
#----------------------------

set vpalette = "/home/phaley/PLUSNet07/HOPS/Plot/Palette/gs2v.pal.white_all"
set opalette = "${palette}"
set palette = "${vpalette}"
set oiref = "${iref}"

set line4 = "Upper 500m velocity"
set fld = "10"
set lvl = "-500"
set pmin = "0.0"
set pmax = "35.0"
set iref = "0"
set gxtnd = "ccntVAV500m_${dirstr}"

#/home/phaley/PLUSNet07/Bin/pe_ccnt_asapvp_va "${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}

#../gmetconvjob ${gxtnd}

set palette = "${opalette}"
set iref = "${oiref}"

#----------------------------
#--- Surface velocity. ---
#----------------------------

set palette = "${vpalette}"

set line4 = "Surface velocity"
set fld = "10"
set lvl = "1"
set pmin = "0.0"
set pmax = "35.0"
set iref = "0"
set gxtnd = "ccntV000m_${dirstr}"

#/home/phaley/PLUSNet07/Bin/pe_ccnt_job_at_nl "${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}

#../gmetconvjob ${gxtnd}

set palette = "${opalette}"
set iref = "${oiref}"

#----------------------------
#--- 100m velocity. ---
#----------------------------

set palette = "${vpalette}"

set line4 = "100m velocity"
set fld = "10"
set lvl = "-100"
set pmin = "0.0"
set pmax = "35.0"
set iref = "0"
set gxtnd = "ccntV100m_${dirstr}"

#/home/phaley/PLUSNet07/Bin/pe_ccnt_job_at_nl "${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}

set palette = "${opalette}"
set iref = "${oiref}"

#----------------------------
#--- Barotropic velocity. ---
#----------------------------

set palette = "${vpalette}"

set line4 = "Barotropic velocity"
set fld = "72"
set lvl = "1"
set pmin = "0.0"
set pmax = "35.0"
set iref = "0"
set gxtnd = "ccntVbar_${dirstr}"

#/home/phaley/PLUSNet07/Bin/pe_ccnt_job_at_nl "${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}

set palette = "${opalette}"
set iref = "${oiref}"

#---------------------------------------------
#--- Northern Temperature Mindoro Section. ---
#---------------------------------------------

set line4 = "Temperature - N Mindoro"

set pmin = "3.77363038"
set pmax = "29.7978706"
set fld = "5"
set lon1 = "120.0525"
set lat1 = "13.0445"
set lon2 = "121.7796"
set lat2 = "11.3671"
set dels = "0.01"
set sztop = "0.0"
set szbot = "500.0"
set itopo = "1"
set kms = "150"
set nppg = "1"
set iref = "0"
set gxtnd = "csec_TNmdor_${dirstr}"

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

#../gmetconvjob ${gxtnd}

#---------------------------------------------
#--- Southern Temperature Mindoro Section. ---
#---------------------------------------------

set line4 = "Temperature - S Mindoro"

set lon1 = "121.7796"
set lat1 = "11.3671"
set lon2 = "121.6433"
set lat2 = "9.2871"
set dels = "0.01"
set szbot = "500.0"
set gxtnd = "csec_TSmdor_${dirstr}"

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

#../gmetconvjob ${gxtnd}

#--------------------------------------------------
#--- Temperature Sulawesi-Sulu. ---
#--------------------------------------------------

set line4 = "Temperature - Sulawesi-Sulu"

set lon1 = "119.4915"
set lat1 = "5.0927"
set lon2 = "119.5925"
set lat2 = "6.7859"
set pmin = "5.2"
set pmax = "29.2"
set dels = "0.01"
set szbot = "500.0"
set itopo = "1"
set gxtnd = "csec_Tsulsul_${dirstr}"

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

#../gmetconvjob ${gxtnd}

#--------------------------------------------------
#--- Temperature Along Bohol. ---
#--------------------------------------------------

set line4 = "Temperature - Along Bohol"

set lon1 = "123.9563"
set lat1 = "9.0210"
set lon2 = "125.1079"
set lat2 = "9.6080"
set pmin = "5.2"
set pmax = "29.2"
set dels = "0.01"
set szbot = "500.0"
set itopo = "1"
set gxtnd = "csec_Tbohol_${dirstr}"

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

#../gmetconvjob ${gxtnd}

#-------------------------
#--- Surface Salinity. ---
#-------------------------

set line4 = "Surface Salinity"
set nfld = "1"
set fld = "6"
set nlvl = "1"
set lvl = "1"
set pmin = "33.30"
set pmax = "34.50"
set vscl = "40"
set iref = "1"
set nppg = "1"
set gxtnd = "ccntS000_${dirstr}"
#set pltdtl  = "no"

#/home/phaley/PLUSNet07/Bin/pe_ccnt_job_at_nl "${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}

#../gmetconvjob ${gxtnd}

#------------------------------------------
#--- Northern Salinity Mindoro Section. ---
#------------------------------------------

set line4 = "Salinity - N Mindoro"

set pmin = "33.3"
set pmax = "34.5"
set lon1 = "120.0525"
set lat1 = "13.0445"
set lon2 = "121.7796"
set lat2 = "11.3671"
set dels = "0.01"
set szbot = "500.0"
set itopo = "1"
set iref = "0"
set gxtnd = "csec_SNmdor_${dirstr}"

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

#../gmetconvjob ${gxtnd}

#------------------------------------------
#--- Southern Salinity Mindoro Section. ---
#------------------------------------------

set line4 = "Salinity - S Mindoro"

set lon1 = "121.7796"
set lat1 = "11.3671"
set lon2 = "121.6433"
set lat2 = "9.2871"
set dels = "0.01"
set szbot = "500.0"
set gxtnd = "csec_SSmdor_${dirstr}"

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

#../gmetconvjob ${gxtnd}

#--------------------------------------------------
#--- Salinity Sulawesi-Sulu. ---
#--------------------------------------------------

set line4 = "Salinity - Sulawesi-Sulu"

set lon1 = "119.4915"
set lat1 = "5.0927"
set lon2 = "119.5925"
set lat2 = "6.7859"
set pmin = "33.7"
set pmax = "34.5"
set dels = "0.01"
set szbot = "500.0"
set itopo = "1"
set gxtnd = "csec_Ssulsul_${dirstr}"

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

#../gmetconvjob ${gxtnd}

#--------------------------------------------------
#--- Salinity Along Bohol. ---
#--------------------------------------------------

set line4 = "Salinity - Along Bohol"

set lon1 = "123.9563"
set lat1 = "9.0210"
set lon2 = "125.1079"
set lat2 = "9.6080"
set pmin = "33.9"
set pmax = "34.5"
set dels = "0.01"
set szbot = "500.0"
set itopo = "1"
set gxtnd = "csec_Sbohol_${dirstr}"

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

#../gmetconvjob ${gxtnd}