#!/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, 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 = " "
set line2 = " "
set line3 = " "
set line4 = "unused"

set day0 = "14299.000"
set dayN = "14299.000"
set dday = "1.0"

set veclat = "35.0"
set veclon = "110.5"

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

set sfile = "/projects/philex/PE_initial/2008/GridTest/arch0/pi_ini.nc"

set palette = "/home/phaley/PLUSNet07/HOPS/Plot/Palette/gs2.pal.white_all"
set pltdt = "/home/phaley/PLUSNet07/HOPS/Plot/PlotData/plot_pe.dat"
set pltcst = "yes"
set coast = "/home/phaley/PhilEx/Data/Coasts/archipelago_h.dat"
set pltdtl  = "yes"
set details = "/home/phaley/PhilEx/PE/mindoro.dat"
set gfile = "/projects/philex/Grids/Rsph_CSS/Arch8/grids_Archipelago9.nc"

#-------------------------
#--- 002m Temperature. ---
#-------------------------

set line4 = "2m Temperature Difference"
set nfld = "1"
set fld = "5"
set nlvl = "1"
set lvl = "-2"
set pmin = "0.0"
set pmax = "0.0"
set vscl = "50"
set iref = "3"
set idover = "10"
set nppg = "1"
set gxtnd = "ccntdT002m"

#-------------------------
#--- 100m Temperature. ---
#-------------------------

set line4 = "100m Temperature Error (36 levels)"
set pfile = "/projects/philex/PE_initial/2008/GridTest/arch4/pi_ini.nc"
set lvl = "-100"
set pmin = "-0.33"
set pmax = "0.23"
set gxtnd = "ccnt36dT100m"

/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}

#-------------------------
#--- 100m Temperature. ---
#-------------------------

set line4 = "100m Temperature Error (44 levels)"
set pfile = "/projects/philex/PE_initial/2008/GridTest/arch1/pi_ini.nc"
set lvl = "-100"
set pmin = "-0.33"
set pmax = "0.23"
set gxtnd = "ccnt44dT100m"

/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}

#-------------------------
#--- 100m Temperature. ---
#-------------------------

set line4 = "100m Temperature Error (70 levels)"
set pfile = "/projects/philex/PE_initial/2008/GridTest/arch2/pi_ini.nc"
set lvl = "-100"
set pmin = "-0.33"
set pmax = "0.23"
set gxtnd = "ccnt70dT100m"

/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}

#-------------------------
#--- 350m Temperature. ---
#-------------------------

set line4 = "350m Temperature Error (36 levels)"
set pfile = "/projects/philex/PE_initial/2008/GridTest/arch4/pi_ini.nc"
set lvl = "-350"
set pmin = "-0.06"
set pmax = "0.23"
set gxtnd = "ccnt36ldT350m"

/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}

#-------------------------
#--- 350m Temperature. ---
#-------------------------

set line4 = "350m Temperature Error (70 levels - 300m over deep)"
set pfile = "/projects/philex/PE_initial/2008/GridTest/arch2/pi_ini.nc"
set lvl = "-350"
set pmin = "-0.06"
set pmax = "0.23"
set gxtnd = "ccnt70l300dT350m"

/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}

#-------------------------
#--- 350m Temperature. ---
#-------------------------

set line4 = "350m Temperature Error (70 levels - 400m over deep)"
set pfile = "/projects/philex/PE_initial/2008/GridTest/arch3/pi_ini.nc"
set lvl = "-350"
set pmin = "-0.06"
set pmax = "0.23"
set gxtnd = "ccnt70l400dT350m"

/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}