atm_plot_bias.sh

source config.diag
module load netcdf
mkdir -p tmp
cp -pr grads-scripts tmp/.
cd tmp
lenexp=${#EXPLIST[@]}
(( nplots = lenexp + 2 ))
echo $lenexp

for CDATE in $CDATELIST; do
  for i in "${!VARLIST[@]}"; do
    VAR=${VARLIST[$i]}
    VARANA=${VARLIST_ERA[$i]}
    VAR=$(echo "$VAR" | tr  '[:lower:]' '[:upper:]')

    echo "Model variable: $VAR"
    echo "ERA variable:   $VARANA"
echo $VAR ' ' plot_${VAR}_${CDATE}.gs
cat >plot_${VAR}_${CDATE}.gs <<EOF
EXPLIST="$EXPLIST"
CDATE=$CDATE
VAR=$VAR
VARA='$VARANA'
VARA_name="$VARANA"
Nmonth=$Nmonth
cmap_field="$cmap_field"
cmap_diff="$cmap_diff"
t=1
lenexp=$lenexp
nplots=$nplots
while(t<=Nmonth)
tt=t-1
iexp=1
while (iexp<=2)
exp.iexp=subwrd(EXPLIST,iexp)
exp=subwrd(EXPLIST,iexp)
'sdfopen  $DATAOUT/'VAR'/'exp.iexp'.'CDATE'.'VAR'.1p0.monthly.nc'
'set t 't
'q time'
say result
ctime=subwrd(result,3)
mm=substr(ctime,6,3)
yyyy=substr(ctime,9,42)
say mm
say yyyy
'q file'
say result
if ( VAR !='V10M' & VAR != 'V200' )
res=sublin(result,7)
else
res=sublin(result,8)
endif
say res
varname=subwrd(res,1)
say varname
'define fcst'iexp'='varname'.'iexp
iexp=iexp+1
endwhile
if ( VARA !='none')
'sdfopen $ANADATADIR/$ERAFILENAME'
'set dfile 'iexp
'set time 00Z01'mm''yyyy
'define ana='VARA_name
endif
say result
'$GRADSDIR/subplot.gs 'nplots'  1'
'set gxout shaded'
'set grid off'
'set grads off'
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
'set stat on'
'd fcst1'
say result
range=sublin(result,9)
cmin=subwrd(range,5)
cmax=subwrd(range,6)
cint=subwrd(range,7)
cint=cint/2
'd fcst1-fcst2'
say result
range=sublin(result,9)
cmind=subwrd(range,5)
cmaxd=subwrd(range,6)
cintd=subwrd(range,7)
cintd=cintd/2
if (math_abs(cmind)>math_abs(cmaxd))
cmaxd=math_abs(cmind)
else
cmind=-1*cmaxd
endif

'set stat off'
*'$GRADSDIR/colormaps.gs -map 'cmap_field'  'flipped_field' -levels 'cmin' 'cmax' 'cint
'$GRADSDIR/color.gs -kind 'cmap_field' 'cmin' 'cmax' 'cint
'set grid off'
'set grads off'
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
'd fcst1'
'd aave(fcst1,g)'
say result
ress=subwrd(result,4)
res=substr(ress,1,10)
'$GRADSDIR/cbarm.gs 1 0 1'
'draw title 'exp.1' 'VAR'\ mean='res' 'mm' 'yyyy

'$GRADSDIR/subplot.gs 'nplots'  2'
'set grid off'
'set grads off'
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
if ( VARA !='none')
'$GRADSDIR/color.gs -kind 'cmap_diff' 'cmind' 'cmaxd' 'cintd
'd fcst2-ana'
'd aave(fcst2-ana,g)'
say result
ress=subwrd(result,4)
res=substr(ress,1,10)
'$GRADSDIR/cbarm.gs 1 0 1'
'draw title 'exp.2' 'VAR' bias\ mean='res' 'mm' 'yyyy
else
'$GRADSDIR/color.gs -kind 'cmap_field' 'cmin' 'cmax' 'cint
'd fcst2'
'd aave(fcst2,g)'
say result
ress=subwrd(result,4)
res=substr(ress,1,10)
'$GRADSDIR/cbarm.gs 1 0 1'
'draw title 'exp.2' 'VAR'\ mean='res' 'mm' 'yyyy
endif

'$GRADSDIR/subplot.gs 'nplots'  3'
'set grid off'
'set grads off'
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
if ( VARA !='none')
'$GRADSDIR/color.gs -kind 'cmap_diff' 'cmind' 'cmaxd' 'cintd
'd fcst1-ana'
'd aave(fcst1-ana,g)'
say result
ress=subwrd(result,4)
res=substr(ress,1,10)
'draw title 'exp.1' 'VAR' bias \ mean='res' 'mm' 'yyyy
endif
'$GRADSDIR/cbarm.gs 1 0 1'

'$GRADSDIR/subplot.gs 'nplots'  4'
'set grid off'
'set grads off'
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
'$GRADSDIR/color.gs -kind 'cmap_diff' 'cmind' 'cmaxd' 'cintd
'set grid off'
'set grads off'
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
'd fcst2-fcst1'
'd aave(fcst2-fcst1,g)'
say result
ress=subwrd(result,4)
res=substr(ress,1,10)

'$GRADSDIR/cbarm.gs 1 0 1'
'draw title 'exp.2'-'exp.1' 'VAR' DIFF mean='res' \ 'mm' 'yyyy



'printim 'exp'.bias_'VAR'_'CDATE'_leadmonth'tt'.png x1000 y1000  white'
'close 1'
t=t+1
*pull dummy
*'c'
reinit
endwhile
'quit'
EOF
grads -blc "run  plot_${VAR}_${CDATE}.gs"
done
done