#!/bin/sh
#=======================================================================
# Radiative forcing plotmix                           --- rf_plotmix ---
# $Id: rf_plotmix_jobdef 618 2010-02-11 01:57:41Z acrnrls $
#=======================================================================
#     keyword :: rf_plotmix
# description :: generate radaitive forcing plots from diagnostic files
#******************************

 set -a
 . betapath2

#  * ........................... Parmsub Parameters ............................

 # These variables are set when the job string is created
 run_start_year=NotSet
 run_start_month=NotSet
 run_stop_year=NotSet
 run_stop_month=NotSet
 current_year=NotSet
 current_month=NotSet
 previous_year=NotSet
 previous_month=NotSet
 next_year=NotSet
 next_month=NotSet

 oldiag="diag4";
 lopgm="lopgmhr";

 jobname=rf_plotmix;
 runid="job000"; uxxx='uxxx'; pool_uxxx=$uxxx; rf_plotmix_uxxx=$pool_uxxx;
 crawork="${runid}_job"; username="acrnxxx"; user="XXX";
 RUNID=`echo "${runid}"|tr '[a-z]' '[A-Z]'`;
 nqsprfx="x${runid}_"; nqsext='';

 # run and days are used as part of several plot titles
 run=$RUNID            # memory99=1
 rf_plotmix_days=NotSet   # memory99=1
 days="DJF"            # memory99=1
 if [ x"$rf_plotmix_days" != "xNotSet" ]; then
   eval days\=$rf_plotmix_days   # memory99=1
 fi
 lc_days=`echo $days|tr '[A-Z]' '[a-z]'`
 case $lc_days in
   jan) lc_days=m01 ;;
   feb) lc_days=m02 ;;
   mar) lc_days=m03 ;;
   apr) lc_days=m04 ;;
   may) lc_days=m05 ;;
   jun) lc_days=m06 ;;
   jul) lc_days=m07 ;;
   aug) lc_days=m08 ;;
   sep) lc_days=m09 ;;
   oct) lc_days=m10 ;;
   nov) lc_days=m11 ;;
   dec) lc_days=m12 ;;
 esac

 year=yyy; mon=mm   # memory99=1

 # yrange is used in the definition of input (flabel) and
 # output (plpfn) file names as well as in the title of most plots.
 yrange=$year   # memory99=1

 # Set yrange to a range of years rather than a single year when
 #   1) current_year and previous_year are defined and not equal or
 #   2) current_year and next_year are defined and not equal
 if [ `expr "$current_year" : '[0-9][0-9]*$'` -gt 0 ]; then
   if [ `expr "$previous_year" : '[0-9][0-9]*$'` -gt 0 ]; then
     if [ $current_year -ne $previous_year ]; then
       # special treatment for djf
       if [ x"$lc_days" = "xdjf" ]; then
         # set if current and previous years are more than 1 year apart
         ym1=`echo $current_year|awk '{y=$1-1;printf "%3.3d",y}' -`
         if [ $previous_year -ne $ym1 ]; then
           yrange=${previous_year}_$current_year   # memory99=1
         fi
       else
         yrange=${previous_year}_$current_year     # memory99=1
       fi
     fi
   fi
   if [ `expr "$next_year" : '[0-9][0-9]*$'` -gt 0 ]; then
     if [ $current_year -ne $next_year ]; then
       # special treatment for djf
       if [ x"$lc_days" = "xdjf" ]; then
         # set if current and next years are more than 1 year apart
         ym1=`echo $next_year|awk '{y=$1-1;printf "%3.3d",y}' -`
         if [ $current_year -ne $ym1 ]; then
           yrange=${current_year}_$next_year       # memory99=1
         fi
       else
         yrange=${current_year}_$next_year         # memory99=1
       fi
     fi
   fi
 fi

 rf_plotmix_prefix="${rf_plotmix_uxxx}_${runid}"                     # memory99=1
 rf_plotmix_prefix_=${rf_plotmix_prefix}_                            # memory99=1
 rf_plotmix_flabel_prefix=${rf_plotmix_prefix_}                      # memory99=1
 rf_plotmix_flabel=${rf_plotmix_flabel_prefix}${yrange}_${lc_days}_  # memory99=1

 # flabel will determine the names of the gp and xp files used
 # ie: files used are ${flabel}gp and ${flabel}xp
 flabel=$rf_plotmix_flabel

 # plpfn is used as part of the resulting gmeta file name
 # The gmeta file name will be of the form uxxx_${plpfn}_${jobname}...
 plpfn=${run}_${lc_days}_${yrange};

 # plunit is used to determine the disposition of the gmeta file
 plunit="VIC"

 # pldir is the full pathname into which the gmeta file will be copied
 pldir=''

 # resol is used as part of some file names defined here
 resol="128_64";
 case $resol in
    96_48) mask=gcm3_landmask_96x48v5  ;;
   128_64) mask=gcm3_landmask_128x64v5 ;;
   192_96) mask=gcm3_landmask_192x96v5 ;;
        *) echo "Unknown resolution $resol"
           exit 1 ;;
 esac

 # obsday is used as part of some file names here as well as in xfind
 # below, where it forms part of the superlabel that xfind looks for
 obsday=$days
 lc_obsday=`echo $obsday|tr '[A-Z]' '[a-z]'`
 case $lc_obsday in
   jan) lc_obsday=m01 ;;
   feb) lc_obsday=m02 ;;
   mar) lc_obsday=m03 ;;
   apr) lc_obsday=m04 ;;
   may) lc_obsday=m05 ;;
   jun) lc_obsday=m06 ;;
   jul) lc_obsday=m07 ;;
   aug) lc_obsday=m08 ;;
   sep) lc_obsday=m09 ;;
   oct) lc_obsday=m10 ;;
   nov) lc_obsday=m11 ;;
   dec) lc_obsday=m12 ;;
 esac

 obsfile="pd_era40_${resol}_198101_200012_${lc_obsday}_";
 obspcp_xa=pd_xie_arkin_1979_1998_pcpn_${resol}_${lc_obsday};
 obspcp_gpcp=pd_gpcp_month_1979_2005_pcpn_${resol}_${lc_obsday};
 obsisccp=pd_isccp_d2_1983_1997_cloud_data_${resol}_${lc_obsday};
 obsuwisc=pd_uwisc_1988_2007_lwp_${resol}_${lc_obsday};
 obsuwisc_msk=pd_uwisc_1988_2007_lwp_tmask_${resol}_${lc_obsday};
 obsssmi=pd_ssmi_1987_2007_${resol}_${lc_obsday};
 obsssmi_msk=pd_ssmi_1987_2007_tmask_${resol}_${lc_obsday};
 obsceres=pd_ceres_ebaf_toa_200003_200510_grid_${resol}_${lc_obsday};
 obsceres_msk=pd_ceres_ebaf_toa_200003_200510_grid_tmask_${resol}_${lc_obsday};
 ceres_cld=pd_ceres_2d_srbavg2_month_2000_2003_data_${resol}_${lc_obsday};
 warren_cld=pd_warren_cld_${resol};

 t1="        1"; t2="999999999"; t3="   1"; s3="   1";
 r1="        1"; r2="999999999"; r3="   1";
 g1="        1"; g2="999999999"; g3="   1";
 a1="        1"; a2="999999999"; a3="   1";
 lmx="  995";  d="B";
 lay="    2";   coord=" ET15"; topsig="-1.00"; moist=" QHYB"; plid="       50.";
 m01="  012";
 plv="   17"; p01="   10"; p02="   20"; p03="   30"; p04="   50"; p05="   70";
            p06="  100"; p07="  150"; p08="  200"; p09="  250"; p10="  300";
            p11="  400"; p12="  500"; p13="  600"; p14="  700"; p15="  850";
            p16="  925"; p17=" 1000"; p18="     "; p19="     "; p20="     ";
            p21="     "; p22="     "; p23="     "; p24="     "; p25="     ";
            p26="     "; p27="     "; p28="     "; p29="     "; p30="     ";
            p31="     "; p32="     "; p33="     "; p34="     "; p35="     ";
            p36="     "; p37="     "; p38="     "; p39="     "; p40="     ";
            p41="     "; p42="     "; p43="     "; p44="     "; p45="     ";
            p46="     "; p47="     "; p48="     "; p49="     "; p50="     ";

 kax="    1"; kin="    1"; lxp="    0"; map="    0"; b="  +";

 dtime="100"; gptime="600"; stime="600";

 memory1="100mb"; memory2="100mb"; memory3="120mb";

  case $resol in
     96_48)
      lrt="   47"; lmt="   47"; typ="    2"
      lon="   96"; lat="   48"; npg="    2"
      ncx="    2"; ncy="    2"
      delt="   1200.0";
      ;;
    128_64)
      lrt="   63"; lmt="   63"; typ="    2"
      lon="  128"; lat="   64"; npg="    2"
      ncx="    2"; ncy="    2"
      delt="    900.0"
      ;;
    192_96)
      lrt="   95"; lmt="   95"; typ="    2"
      lon="  192"; lat="   96"; npg="    2"
      ncx="    2"; ncy="    2"
      delt="    600.0"
      ;;
         *) echo "Unknown resolution $resol"
            exit 1 ;;
 esac

 # ptitle is a common prefix for all plot titles
 yrange=`echo $yrange|sed 's/_/-/g'`
 ptitle="$days $yrange  $run "

#  * ............................ Condef Parameters ............................

 nocldo=on
 auto=off
 join=3
 gcm2plus=on
 gcmtsav=on
 datatype=specsig
 wxstats=on
 blckcld=on
 noprint=on
 nextjob=on
 debug=off
 obsdat=on
 xtradif=on
 colourplots=on
 shade=off
 stat2nd=off

#  * ............................. Deck Definition .............................

 . rf_plotmix.dk

#end_of_job