vsdb summary 1

1. Compile vsdb package 

$ cd /data/users/dxu/vsdb_pkg/vsdb_v17

$ ./build.sh

2. Configure two shell scripts: 
#            setup_envs.sh 
#            vsdbjob_submit.sh
$  cd /data/users/dxu/vsdb_pkg/vsdb_v17
$  vi setup_envs.sh
$  vi vsdbjob_submit.sh

3. Run VSDB
$  cd /data/users/dxu/vsdb_pkg/vsdb_v17
$  ./vsdbjob_submit.sh

4. Input, output, log and run dir

steps	run dir
1	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy26072/stats
2	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy47090/acrms47090
3	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy19091/mkup_precip
4	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy22490/plot_pcp
5	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy32596/fit
6	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy20273/2dmaps
steps	log file
1	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy26072/vstep1.out
2	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy47090/vstep2.out
3	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy19091/mkup_rain_stat.out
4	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy22490/plot_pcp.out
5	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy32596/fit2obs.out
6	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy20273/2dmaps[1-4].out
steps	in dir
1	/data/users/dxu/vsdb_workspace/data/input/fcst_data
2	/data/users/dxu/vsdb_workspace/data/output/vsdb_data
3	/data/users/dxu/vsdb_workspace/data/input/fcst_data
4	/data/users/dxu/vsdb_workspace/data/intermediate/dxu/archive
5	/data/users/dxu/vsdb_workspace/data/input/f2o
6	/data/users/dxu/vsdb_workspace/data/input/fcst_data
	/data/users/dxu/vsdb_workspace/data/input/plot2d/obdata
steps	out dir
1	/data/users/dxu/vsdb_workspace/data/output/vsdb_data
2	/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy44333/acrms44333/G2/anom/HGT               score text files /data/users/dxu/vsdb_workspace/data/stmp/dxu/web/allmodel               PNG
3	/data/users/dxu/vsdb_workspace/data/intermediate/dxu/archive
4	/data/users/dxu/vsdb_workspace/data/stmp/dxu/web/rain                       PNG
5	/data/users/dxu/vsdb_workspace/data/stmp/dxu/web/fits                        GIF
6	/data/users/dxu/vsdb_workspace/data/stmp/dxu/web/2D/d[1-4]             GIF

 5. script that contains sub_cardinal
export SUBJOB=$vsdbhome/bin/sub_cardinal
File 1 is :  ./map_util/sfcfcst_1cyc.sh
File 2 is :  ./map_util/allcenters_rmsmap.sh
File 3 is :  ./map_util/allcenters_1cyc.sh
File 4 is :  ./setup_envs.sh
File 5 is :  ./vsdbjob.sh
File 6 is :  ./precip/plot_pcp.sh
File 7 is :  ./precip/precip_score_vsdb.sh
File 8 is :  ./fit2obs/plotall.sh
File 9 is :  ./fit2obs/fit2obs.sh
File 10 is :  ./plot2d/maps2d_new.sh
File 11 is :  ./vsdbjob_submit.sh
File 12 is :  ./grid2obs/grid2obs_plot.sh
File 13 is :  ./grid2obs/grid2obs_driver.sh
File 14 is :  ./grid2obs/scripts/get_opsgfs_data.sh
File 15 is :  ./grid2obs/scripts/get_paragfs_data.sh
File 16 is :  ./grid2obs/scripts/g2o_sfcmap.sh
File 17 is :  ./grid2obs/scripts/grid2obssfc.fits.sh
File 18 is :  ./grid2obs/scripts/g2o_airmap.sh
File 19 is :  ./grid2obs/grid2obs_opsdaily.sh
File 20 is :  ./grid2obs/grid2obs.sh
File 21 is :  ./verify_exp_step2.sh

6. all options used in SUBJOB
1) Flag used in SUBJOB
 $SUBJOB -a $ACCOUNT -q $CUE2FTP -g $GROUP -p 1/1/S -t 1:00:00 -r 128/1 -j ftpcard -o ftpcard$$.out ${rundir}/ftpcard$$.sh

-e : EVN variable list
$SUBJOB -e $listvar -a $task  -q $cue -g $GROUP -p 1/1/S -r 512/1 -t 3:00:00 -o  $SUBJOB -a $ACCOUNT -q $CUE2FTP -g $GROUP -p 1/1/S -r 256/1 -w +${waitfits} -t 1:00:00 -j ftpfits -o $mapdir/ftpfits.out ${mapdir}/ftpfits.sh

2) Original meaning of flags.
      where the options are:
      -a account        account (default:none)
      -e envars         copy comm-separated environment variables
      -g group          group name
      -j jobname        specify jobname (default: executable basename)
      -n                write command file to stdout rather than submitting it
      -o output         specify output file (default: jobname.out)
      -p procs[/nodes]
                        number of MPI tasks and number of nodes
      -q queue          queue name
      -r nodetype       node type (harp or neha)
      -v                verbose mode
      -t timew          wall time limit in [[hh:]mm:]ss format (default: 900)
      -w when           when to run, in yyyymmddhh[mm], +hh[mm], thh[mm], or
                        Thh[mm] (full, incremental, today or tomorrow) format
                        (default: now)

3) sub_badger is a wrapper to translate above options into options that scheduler on badger can recognize.
on badger: 
 qsub -V : pass all the EVN variables.

4) sub_cardinal is a wrapper to translate above options into options that scheduler on cardinal can recognize.
on cardinal:
sbatch --export=<environment variables | ALL | NONE>
              Identify which environment variables are propagated to  the  batch  job.
              Multiple environment variable names should be comma separated.  Environ-
              ment variable names may be specified to propagate the current  value  of
              those  variables  (e.g.  "--export=EDITOR")  or  specific values for the
              variables  may  be  exported  (e.g..  "--export=EDITOR=/bin/vi").   This
              option particularly important for jobs that are submitted on one cluster
              and execute on a different  cluster  (e.g.  with  different  paths).  By
              default  all  environment  variables  are propagated. If the argument is
              NONE or specific environment variable  names,  then  the  --get-user-env
              option  will implicitly be set to load other environment variables based
              upon the user’s configuration on the cluster which executes the job.

7. score card related 
File 1 is :  ./run_scorecard.sh
     file3=${scoredir}/score_${stat}_conflimit_${namedaily}_${mdnamec2}_day${dd}
File 2 is :  ./map_util/allcenters_rmsmap.sh
          '${vsdbhome}/map_util/grads/fprintf.gs 'tv.n' score_${var}_conflimit_'
  cp score_${var}_conflimit*.txt $scoredir
File 3 is :  ./map_util/allcenters_1cyc.sh
          '${vsdbhome}/map_util/grads/fprintf.gs 'tv.n' score_cor_conflimit_'${n.


 1). Create score card text files:
              ./vsdbjob_submit.sh  
      -->   ./verify_exp_step2.sh 
      -->   ./map_util/allcenters_rmsmap.sh and ./map_util/allcenters_1cyc.sh       

A) anomaly correlation on single pressure layer
    $SUBJOB -e $listvar -a $ACCOUNT  -q $CUE2RUN -g $GROUP -p $penode -t $cputime -j HGTanom${narea} -o $rundir/HGT_anom.out \
           ${sorcdir}/allcenters_1cyc.sh $edate $ndays $fdays
    if [ $? -ne 0 ]; then ${sorcdir}/allcenters_1cyc.sh $edate $ndays $fdays  ; fi
    sleep 3
B) rms and bias    $SUBJOB -e $listvar -a $ACCOUNT  -q $CUE2RUN -g $GROUP -p $penode -t $cputime -j HGTpres${narea} -o $rundir/HGT_pres.out \
           ${sorcdir}/allcenters_rmsmap.sh $edate $ndays $fdays
    if [ $? -ne 0 ]; then ${sorcdir}/allcenters_rmsmap.sh $edate $ndays $fdays  ; fi
    sleep 3

2) Location of score card files?
 Details of ./verify_exp_step2.sh 
  step 1:   move vsdb status files to     /data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy9457/acrms9457/vsdb_data
  step 2 :  3 variable types:    anom, pres, sfc
  step 3: 5 regions:     G2NHX ,G2SHX, G2TRO, G2,G2PNA
     Eg: 
     /data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy9457/acrms9457/G2NHX
     /data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy9457/acrms9457/G2SHX
     /data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy9457/acrms9457/G2TRO
     /data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy9457/acrms9457/G2
     /data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy9457/acrms9457/G2PNA

  step 4 : Each region has two types:  pres & anom 

                    /data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy9457/acrms9457/G2NHX

├── anom

│   ├── HGT

│   ├── HGTWV0-3

│   ├── HGTWV10-20

│   ├── HGTWV4-9

│   ├── PMSL

│   ├── T

│   ├── U

│   ├── V

│   └── WIND

└── pres

    ├── HGT

    ├── O3

    ├── T

    ├── U

    ├── V

    └── WIND

A#) anomaly correlation on single pressure layer:
vtype=anom  + map_util/allcenters_1cyc.sh
-------                                      ------------
 Level                                      Parameters
-------                                       ------------

"P1000 P700 P500 P250"    "HGT HGT_WV1/0-3 HGT_WV1/4-9 HGT_WV1/10-20"

"P850 P500 P250"               "WIND U V T"
"MSL"                                 "PMSL"

 log files:
/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy9457/acrms9457/G2NHX

PMSL_anom.out
UVT_anom.out        # combine U, V and T
HGT_anom.out        # combine all HGT including layered HGT.

B#) rms and bias :

vtype=pres   + map_util/allcenters_rmsmap.sh
-------                       ------------
 Level                       Parameters
-------                       ------------if maptop = "10" 
"P1000 P925 P850
 P700 P500 P400
 P300 P250 P200          "HGT WIND U V T"
 P150 P100 P50
 P20 P10"                
if maptop = "50" 
"P1000 P925 P850
 P700 P500 P400          "HGT WIND U V T"   
 P300 P250 P200
 P150 P100 P50"
if maptop = "100"
"P1000 P925 P850
 P700 P500 P400       "HGT WIND U V T"
 P300 P250 P200
 P150 P100"

"P100 P70 P50 P30 P20 P10" : "O3"

log files: 
 /data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy9457/acrms9457/G2NHX

O3_pres.out
HGT_pres.out 
WIND_pres.out
U_pres.out
V_pres.out
T_pres.out

step 5 :  Process "sfc" separately.

vtype=sfc  + map_util/sfcfcst_1cyc.shreglist="G2 G2/NHX G2/SHX G2/TRO G2/N60 G2/S60 G2/NPO G2/SPO G2/NAO G2/SAO G2/CAM G2/NSA"

-------                   ------------
 Level                     Parameters
-------                   ------------

"SL1L2"              "CAPE CWAT PWAT HGTTRP"

                          "TMPTRP HPBL PSFC PSL"   
                          "RH2m SPFH2m T2m TOZNE TG"
                          "U10m V10m WEASD TSOILT WSOILT"   

/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy9457/acrms9457/sfc

└── sfc
    ├── CAPE
    ├── CWAT
    ├── HGTTRP
    ├── HPBL
    ├── PSFC
    ├── PSL
    ├── PWAT
    ├── RH2m
    ├── SPFH2m
    ├── T2m
    ├── TG
    ├── TMPTRP
    ├── TOZNE
    ├── TSOILT
    ├── U10m
    ├── V10m
    ├── WEASD
    └── WSOILT

3) Final html output files
eg: output location: /scratch2/portfolios/NESDIS/h-sandy/noscrub/Deyong.Xu/vsdb_workspace/data/stmp/Deyong.Xu/nwpvrfy29853/acrms29853/score
scorecard.html
scorecard.css
mainindex.html
legend.html




8. How VSDB status files are used? 
verify_exp_step2.sh
-->  map_util/allcenters_1cyc.sh    &     map_util/allcenters_rmsmap.sh
--> map_util/gen_scal.sh
      map_util/gen_scal_pres.sh
      map_util/gen_wind.sh
      map_util/gen_wind_pres.sh
      map_util/gen_sfc.sh

Detail of "./gen_scal.sh"
1) convert vsdb file into txt file
while [ $cdate -le $edate ]; do
  fhour=00; vhr=$cyc
  while [ $fhour -le $vlength ]; do
    datadir=${vsdb_data}/${vtype}/${vhr}Z/${model}
    vsdbname=${datadir}/${model}_${cdate}.vsdb               ### vsdb files
    string=" $mdl $fhour ${cdate}${vhr} $mdl $reg SAL1L2 $vnam $lev "
    mycheck=$( grep "$string" $vsdbname )
    if [ $? -ne 0 ]; then
      echo "missing" >>$outname.txt
    else
      grep "$string" $vsdbname |cat >>$${outname}.txt           ### grep data out and save it into text file
    fi
    fhour=` expr $fhour + $fhout `
    if [ $fhour -lt 10 ]; then fhour=0$fhour ; fi
    vhr=` expr $vhr + $fhout `
    if [ $vhr -ge 24 ]; then vhr=`expr $vhr - 24 `; fi
    if [ $vhr -lt 10 ]; then vhr=0$vhr ; fi
  done

2) convert txt file into binary file via convert.f that is generated on the fly. 
  ( within gen* script such as gen_scal.sh script )
       open(9,file="modelname.txt",form="formatted",status="old")
       open(10,file="${outname}.txt",form="formatted",status="old")        ## text file generated above
       open(11,file="tmp.txt",form="formatted",status="new")
       open(20,file="${outname}.bin",form="unformatted",status="new")      ## output binary file

$FC $FFLAG -o convert.x convert.f
./convert.x
meantxt=${vnam1}_${lev}_${reg1}_${yyyymm}
mv fort.13 meancor_${meantxt}.txt
mv fort.14 meanrms_${meantxt}.txt
mv fort.15 meanbias_${meantxt}.txt

3) create grads control file
cat >${outname}.ctl <<EOF1
dset ^${outname}.bin                                 ## binary file generated above
undef -99.9
options big_endian sequential
title scores
xdef   $nmdcyc linear 1  1


4) output ctl files 
HGT_P1000_G2NHX_2014020120140228.ctl
HGT_P700_G2NHX_2014020120140228.ctl
HGT_P500_G2NHX_2014020120140228.ctl
HGT_P250_G2NHX_2014020120140228.ctl


Details of "allcenters_1cyc.sh"

0)# -- search data for all models; write out binary data, create grads control file
/gen_wind.sh
/gen_scal.sh

1) # ----- PLOT TYPE 1:  time series of anomaly correlations ----
cat >acz_${outname1}.gs <<EOF1
'reinit'; 'set font 1'
'open ${outname}.ctl'
mdc.1=${mdnamec[0]}
* Create verification scorecard text files  ( score card files )
  '${vsdbhome}/map_util/grads/fprintf.gs 'sc.i' score_cor_'${namedaily}'_'mdc.i'_day'%day'.txt %-7.6f'          ### GOLD 
...
$GRADSBIN/grads -bcp "run acz_${outname1}.gs" 
2) # ----- PLOT TYPE 2:  Die-off plot for mean correlation over $ndays days----
ndaysp1=`expr $ndays + 1`
fdaysp1=`expr $fdays + 1`
cat >cordieoff_${outname1}.gs <<EOF1
'reinit'; 'set font 1'
'open ${outname}.ctl'
** Create verification scorecard text files  ( conflimit files  between models )
'${vsdbhome}/map_util/grads/fprintf.gs 'tv.n' score_cor_conflimit_'${namedaily}'_'mdc.i'_day'n-1'.txt %-7.6f'       ## GOLD:  generate score_cor_conflimit* files needed to run score card.
...
$GRADSBIN/grads -bcp "run cordieoff_${outname1}.gs"


3)  # ----- PLOT TYPE 3:  difference of AC, other models minus first model  ----
cat >cordiff_${outname1}.gs <<EOF1
'reinit'; 'set font 1'
'run $sorcdir/grads/white.gs'
'open ${outname}.ctl'
 ...
$GRADSBIN/grads -bcp "run cordiff_${outname1}.gs"


4) # ----- PLOT TYPE 4:  frequency distribution of anomaly correlations ----
ndayfq=$ndays
if [ $ndayfq -gt 20 ]; then ndayfq=20; fi
nday05=`expr $ndayfq \/ 2 `
cat >freq_${outname1}.gs <<EOF1
'reinit'; 'set font 1'
'open ${outname}.ctl'
...
$GRADSBIN/grads -bcp "run freq_${outname1}.gs"

5) output grads script (.gs) files
freq_HGT_P700_G2NHX_2014020120140228.gs
freq_HGT_P250_G2NHX_2014020120140228.gs
freq_HGT_P500_G2NHX_2014020120140228.gs
freq_HGT_P1000_G2NHX_2014020120140228.gs

cordiff_HGT_P700_G2NHX_2014020120140228.gs
cordiff_HGT_P250_G2NHX_2014020120140228.gs
cordiff_HGT_P500_G2NHX_2014020120140228.gs
cordiff_HGT_P1000_G2NHX_2014020120140228.gs

acz_HGT_P700_G2NHX_2014020120140228.gs
acz_HGT_P250_G2NHX_2014020120140228.gs
acz_HGT_P500_G2NHX_2014020120140228.gs
acz_HGT_P1000_G2NHX_2014020120140228.gs

cordieoff_HGT_P700_G2NHX_2014020120140228.gs
cordieoff_HGT_P250_G2NHX_2014020120140228.gs
cordieoff_HGT_P500_G2NHX_2014020120140228.gs
cordieoff_HGT_P1000_G2NHX_2014020120140228.gs


6) Sample score files needed to run scorecard.

[dxu@s4-cardinal HGT]$ pwd/data/users/dxu/vsdb_workspace/data/stmp/dxu/nwpvrfy_deyong/acrms_deyong/G2NHX/anom/HGT
[dxu@s4-cardinal HGT]$ tree -f |grep score_cor  |grep 250
├── ./score_cor_conflimit_HGT_P250_G2NHX_ECM_day1.txt  # for 2nd model only
├── ./score_cor_conflimit_HGT_P250_G2NHX_ECM_day2.txt
├── ./score_cor_conflimit_HGT_P250_G2NHX_ECM_day3.txt
├── ./score_cor_conflimit_HGT_P250_G2NHX_ECM_day4.txt
├── ./score_cor_conflimit_HGT_P250_G2NHX_ECM_day5.txt
├── ./score_cor_HGT_P250_G2NHX_ECM_day1.txt
├── ./score_cor_HGT_P250_G2NHX_ECM_day3.txt
├── ./score_cor_HGT_P250_G2NHX_ECM_day5.txt
├── ./score_cor_HGT_P250_G2NHX_GFS_day1.txt    # Reference model
├── ./score_cor_HGT_P250_G2NHX_GFS_day3.txt
├── ./score_cor_HGT_P250_G2NHX_GFS_day5.txt


Details of "allcenters_rmsmap.sh"
varlist="rms bias pcor emd epv rsd msess"

0) # -- search data for all models; write out binary data, create grads control file
gen_wind_pres.sh
gen_scal_pres.sh

# ----- PLOT TYPE 1:  maps of $var as a function of calendar day and pressure for each forecast time ----
cat >${var}p_${outname1}.gs <<EOF1
'reinit'; 'set font 1'
'run $sorcdir/grads/white.gs'
'open ${outname1}.ctl'
...
EOF1
$GRADSBIN/grads -bcp "run ${var}p_${outname1}.gs"


eg: 
rmsp_WIND_G2NHX_2014020120140228.gs
biasp_WIND_G2NHX_2014020120140228.gs
pcorp_WIND_G2NHX_2014020120140228.gs
emdp_WIND_G2NHX_2014020120140228.gs
epvp_WIND_G2NHX_2014020120140228.gs
rsdp_WIND_G2NHX_2014020120140228.gs
msessp_WIND_G2NHX_2014020120140228.gs


# ----- PLOT TYPE 2:  maps of mean ${var} as a function of forecast time and pressure  ----ndaysp1=`expr $ndays + 1`
fdaysp1=`expr $fdays + 1`
cat >${var}pmean_${outname1}.gs <<EOF1
'reinit'; 'set font 1'
'run $sorcdir/grads/white.gs'
'open ${outname1}.ctl'
...
EOF1
$GRADSBIN/grads -bcp "run ${var}pmean_${outname1}.gs"
eg: 
rmspmean_WIND_G2NHX_2014020120140228.gs
biaspmean_WIND_G2NHX_2014020120140228.gs
pcorpmean_WIND_G2NHX_2014020120140228.gs
emdpmean_WIND_G2NHX_2014020120140228.gs
epvpmean_WIND_G2NHX_2014020120140228.gs
rsdpmean_WIND_G2NHX_2014020120140228.gs
msesspmean_WIND_G2NHX_2014020120140228.gs


# ----- PLOT TYPE 3:  time series of ${var} Errors----cat >${var}_${outname1}${lev}.gs <<EOF1
'reinit'; 'set font 1'
'open ${outname1}.ctl'
'${vsdbhome}/map_util/grads/fprintf.gs 'sc.i' score_${var}_${namedaily1}_'mdc.i'_day'%day'.txt %-7.6f'       ## GOLD: generate single score file for each model 
...
EOF1
$GRADSBIN/grads -bcp "run ${var}_${outname1}${lev}.gs"
eg: 
rsddieoff_WIND_G2NHX_2014020120140228P200.gs
rsddieoff_WIND_G2NHX_2014020120140228P10.gs
rsddieoff_WIND_G2NHX_2014020120140228P100.gs


# ----- PLOT TYPE 4:  mean ${var} error growth curve over $ndays days----ndaysp1=`expr $ndays + 1`
fdaysp1=`expr $fdays + 1`
cat >${var}dieoff_${outname1}${lev}.gs <<EOF1
'reinit'; 'set font 1'
'open ${outname1}.ctl'
'${vsdbhome}/map_util/grads/fprintf.gs 'tv.n' score_${var}_conflimit_'${namedaily1}'_'mdc.i'_day'n-1'.txt %-7.6f'          ## GOLD  *conflimit*
...
EOF1
$GRADSBIN/grads -bcp "run ${var}dieoff_${outname1}${lev}.gs"
eg: 
rsddieoff_WIND_G2NHX_2014020120140228P10.gs
rsddieoff_WIND_G2NHX_2014020120140228P100.gs
msessdieoff_WIND_G2NHX_2014020120140228P850.gs
msessdieoff_WIND_G2NHX_2014020120140228P700.gs


log files: 
/scratch2/portfolios/NESDIS/h-sandy/noscrub/Deyong.Xu/vsdb_workspace/data/stmp/Deyong.Xu/nwpvrfy774/acrms774/G2NHX/pres/WIND
./score_rsd_WIND_P1000_G2NHX_GFS_day3.txt
./score_epv_WIND_P500_G2NHX_GFS_day3.txt
./score_pcor_WIND_P700_G2NHX_ECM_day5.txt
./score_pcor_WIND_P50_G2NHX_ECM_day5.txt


8. Region control / distribution 

Go to here:  /scratch2/portfolios/NESDIS/h-sandy/noscrub/Deyong.Xu/vsdb_pkg/vsdb_v17/exe
on zeus and do following command and find the lines attached below.
cntl_sfc.sh seems to be what you're trying to find.

$ vi  cntl_sfc.sh
    12 ${gd}
       ${gd}/NHX 0  20  360  80
       ${gd}/SHX 0 -80  360 -20
       ${gd}/TRO 0 -20  360  20
       ${gd}/N60 0  60  360  90
       ${gd}/S60 0 -90  360 -60
       ${gd}/NPO
       ${gd}/SPO
       ${gd}/NAO
       ${gd}/SAO
       ${gd}/CAM
       ${gd}/NSA

$ vi cntl_pres.sh
    5  ${gd}
       ${gd}/NHX 0  20  360  80
       ${gd}/SHX 0 -80  360 -20
       ${gd}/TRO 0 -20  360  20
       ${gd}/PNA 180 20 320  75

$ vi cntl_anom.sh
        5  ${gd}
       ${gd}/NHX 0  20  360  80
       ${gd}/SHX 0 -80  360 -20
       ${gd}/TRO 0 -20  360  20
       ${gd}/PNA 180 20 320  75


9. Account related issue: 
On some platforms such as zeus, account setting becomes matter because not all of accounts are allowed to submit jobs on zeus. This account information is part of parameters of job submission.

So for zeus, I set it to "h-sandy"
export ACCOUNT=h-sandy

gsi diagnostic data information

==============================
diagnostic data info:
==============================

module read_diag
    subroutine read_radiag_header
    subroutine read_radiag_data

### 1.1) header fix ####
  type diag_header_fix_list
     character(len=20) :: isis           ! sat and sensor type
     character(len=10) :: id             ! sat type
     character(len=10) :: obstype        ! observation type
     integer(i_kind) :: jiter            ! outer loop counter
     integer(i_kind) :: nchan            ! number of channels in the sensor
     integer(i_kind) :: npred            ! number of updating bias correction predictors
     integer(i_kind) :: idate            ! time (yyyymmddhh)
     integer(i_kind) :: ireal            ! # of real elements in the fix part of a data record
     integer(i_kind) :: ipchan           ! # of elements for each channel except for bias correction terms
     integer(i_kind) :: iextra           ! # of extra elements for each channel
     integer(i_kind) :: jextra           ! # of extra elements
     integer(i_kind) :: idiag            ! first dimension of diag_data_chan
     integer(i_kind) :: angord           ! order of polynomial for adp_anglebc option
     integer(i_kind) :: iversion         ! radiance diagnostic file version number
     integer(i_kind) :: inewpc           ! indicator of newpc4pred (1 on, 0 off)
  end type diag_header_fix_list

### 1.2) header chan ####
  type diag_header_chan_list
     real(r_single) :: freq              ! frequency (Hz)
     real(r_single) :: polar             ! polarization
     real(r_single) :: wave              ! wave number (cm^-1)
     real(r_single) :: varch             ! error variance (or SD error?)
     real(r_single) :: tlapmean          ! mean lapse rate
     integer(i_kind):: iuse              ! use flag
     integer(i_kind):: nuchan            ! sensor relative channel number
     integer(i_kind):: iochan            ! satinfo relative channel number
  end type diag_header_chan_list

### 2.1) data name ####
  type diag_data_name_list
     character(len=10),dimension(ireal_radiag) :: fix
     character(len=10),dimension(:),allocatable :: chn
  end type diag_data_name_list


### 2.2) data fix ####
  type diag_data_fix_list
     real(r_single) :: lat               ! latitude (deg)
     real(r_single) :: lon               ! longitude (deg)
     real(r_single) :: zsges             ! guess elevation at obs location (m)
     real(r_single) :: obstime           ! observation time relative to analysis
     real(r_single) :: senscn_pos        ! sensor scan position (integer(i_kind))
     real(r_single) :: satzen_ang        ! satellite zenith angle (deg)
     real(r_single) :: satazm_ang        ! satellite azimuth angle (deg)
     real(r_single) :: solzen_ang        ! solar zenith angle (deg)
     real(r_single) :: solazm_ang        ! solar azimumth angle (deg)
     real(r_single) :: sungln_ang        ! sun glint angle (deg)
     real(r_single) :: water_frac        ! fractional coverage by water
     real(r_single) :: land_frac         ! fractional coverage by land
     real(r_single) :: ice_frac          ! fractional coverage by ice
     real(r_single) :: snow_frac         ! fractional coverage by snow
     real(r_single) :: water_temp        ! surface temperature over water (K)
     real(r_single) :: land_temp         ! surface temperature over land (K)
     real(r_single) :: ice_temp          ! surface temperature over ice (K)
     real(r_single) :: snow_temp         ! surface temperature over snow (K)
     real(r_single) :: soil_temp         ! soil temperature (K)
     real(r_single) :: soil_mois         ! soil moisture
     real(r_single) :: land_type         ! land type (integer(i_kind))
     real(r_single) :: veg_frac          ! vegetation fraction
     real(r_single) :: snow_depth        ! snow depth
     real(r_single) :: sfc_wndspd        ! surface wind speed
     real(r_single) :: qcdiag1           ! ir=cloud fraction, mw=cloud liquid water
     real(r_single) :: qcdiag2           ! ir=cloud top pressure, mw=total column water
     real(r_single) :: tref              ! reference temperature (Tr) in NSST
     real(r_single) :: dtw               ! dt_warm at zob
     real(r_single) :: dtc               ! dt_cool at zob
     real(r_single) :: tz_tr             ! d(Tz)/d(Tr)
  end type diag_data_fix_list

### 2.3) data chan ####
  type diag_data_chan_list
     real(r_single) :: tbobs              ! Tb (obs) (K)
     real(r_single) :: omgbc              ! Tb_(obs) - Tb_(simulated w/ bc)  (K)
     real(r_single) :: omgnbc             ! Tb_(obs) - Tb_(simulated_w/o bc) (K)
     real(r_single) :: errinv             ! inverse error (K**(-1))
     real(r_single) :: qcmark             ! quality control mark
     real(r_single) :: emiss              ! surface emissivity
     real(r_single) :: tlap               ! temperature lapse rate
     real(r_single) :: tb_tz              ! d(Tb)/d(Tz)
     real(r_single) :: bicons             ! constant bias correction term
     real(r_single) :: biang              ! scan angle bias correction term
     real(r_single) :: biclw              ! CLW bias correction term
     real(r_single) :: bilap2             ! square lapse rate bias correction term
     real(r_single) :: bilap              ! lapse rate bias correction term
     real(r_single) :: bicos              ! node*cos(lat) bias correction term
     real(r_single) :: bisin              ! sin(lat) bias correction term
     real(r_single) :: biemis             ! emissivity sensitivity bias correction term
     real(r_single),dimension(:),allocatable :: bifix          ! angle dependent bias
     real(r_single) :: bisst              ! SST bias correction term
  end type diag_data_chan_list

### 2.4) data extra ####
  type diag_data_extra_list
     real(r_single) :: extra              ! extra information
  end type diag_data_extra_list