#!/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 oabgn = "2003 227"

set line1 = " "
set line2 = " "
set line3 = " "
set line4 = " "

set day0 = "12092.000"
set dayN = "12092.000"
set dday = "1.0"

set veclat = "19.0"
set veclon = "113.25"

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

set pfile = "oag_1250m.nc"
set sfile = "oag_1250m.nc"

set palette = "/home/phaley/PLUSNet07/HOPS/Plot/Palette/gs2.pal.white_all"
set pltdt = "/home/phaley/PLUSNet07/HOPS/Plot/PlotData/plot_qg.dat"
set cstflg = "yes"
#set cstfile = "/home/phaley/PhilEx/Data/Coasts/e_asia_i.dat"
set cstfile = "/home/phaley/PhilEx/Data/Coasts/archipelago_h.dat"
#set cstfile = "/home/phaley/PhilEx/Data/Coasts/e_asia_h.dat"
#set cstfile = "/home/phaley/PhilEx/Data/Coasts/archipelago_f.dat"
set detflg  = "no"
set details = "/data/projects2/AOSN2/OAG/aosn2_sml1km2_sec.dat"
set gfile = "/home/phaley/PhilEx/Grids/Rsph_CSS/Archipelago/grids_Archipelago9_sa.nc"

#-----------------------------------
#--- Temp Surface. ---
#-----------------------------------

set line4 = "Surface Temperature, WOA-05"
set pfile = "/projects/philex/OAG/2008/May04/arch0/oag_1250m.nc"
set day0 = "12092.000"
set dayN = "12092.000"

set nfld = "1"
set fld = "11"
set nlvl = "1"
set lvl = "1"
set pmin = "28.9"
set pmax = "30.0"
set vscl = "25"
set iref = "0"
set idover = "10"
set nppg = "1"
set gxtnd = "ccntWoaT0000"

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

./gmetconvjob ${gxtnd}

#-----------------------------------
#--- Temp Surface. ---
#-----------------------------------

set line4 = "Surface Temperature, WOA-05 + False Pts"
set pfile = "/projects/philex/OAG/2008/May04/arch1/oag_1250m.nc"
set day0 = "12092.000"
set dayN = "12092.000"

set gxtnd = "ccntWoaFpT0000"

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

./gmetconvjob ${gxtnd}

#-----------------------------------
#--- Temp Surface. ---
#-----------------------------------

set line4 = "Surface Temperature, Exploratory"
set pfile = "/projects/philex/OAG/2008/May04/arch2/oag_1250m.nc"
set day0 = "14299.000"
set dayN = "14299.000"

set gxtnd = "ccntExplT0000"

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

./gmetconvjob ${gxtnd}

#-----------------------------------
#--- 1000m Temp. ---
#-----------------------------------

set line4 = "1000m Temperature, WOA-05"
set pfile = "/projects/philex/OAG/2008/May04/arch0/oag_1250m.nc"
set day0 = "12092.000"
set dayN = "12092.000"

set lvl = "33"
set pmin = "4.1"
set pmax = "11.1"
set gxtnd = "ccntWoaT1000"

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

./gmetconvjob ${gxtnd}

#-----------------------------------
#--- Temp Surface. ---
#-----------------------------------

set line4 = "1000m Temperature, WOA-05 + False Pts"
set pfile = "/projects/philex/OAG/2008/May04/arch1/oag_1250m.nc"
set day0 = "12092.000"
set dayN = "12092.000"

set gxtnd = "ccntWoaFpT1000"

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

./gmetconvjob ${gxtnd}

#-----------------------------------
#--- Temp Surface. ---
#-----------------------------------

set line4 = "1000m Temperature, Exploratory"
set pfile = "/projects/philex/OAG/2008/May04/arch2/oag_1250m.nc"
set day0 = "14299.000"
set dayN = "14299.000"
set lvl = "33"

set gxtnd = "ccntExplT1000"

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

./gmetconvjob ${gxtnd}