#!/bin/ksh 
set -ex

#CFSR grib2 data:
dir0=/cpc/cfsr/archive/month_h/pgb

dir=/cpc/cfsr_ipvcdb/ipv/CFSv2/ypan/SFS/CFSR
dir2=$dir/data

wgrib2=/usr/local/bin/wgrib2

i=0
#ymd=19790101
ymd=19910101

while [[ ${ymd} -le 20201201 ]]; do
#while [[ ${ymd} -le 19790101 ]]; do

yr=`echo $ymd | cut -c1-4`
mon=`echo $ymd | cut -c5-6`
ym=${yr}${mon}

echo "ym=", ${ym}

dfa=${dir0}/pgb.gdas.${ym}.grb2

cd ${dir2}

i=$((i + 1))
  if [[ $i -ge 1 && $i -le 9 ]] ; then
    N=000$i
  fi

  if [[ $i -ge 10 && $i -le 99 ]] ; then
    N=00$i
  fi

  if [[ $i -ge 100 && $i -le 999 ]] ; then
    N=0$i
  fi

  if [[ $i -ge 1000 ]] ; then
    N=$i
  fi
echo "N=", $N

$wgrib2 $dfa -match "SOILW:0-0.1 m" | $wgrib2 -i $dfa -bin sm.$N

cd $dir

gfortran x_read_sm_1deg.f
./a.out $N

rm -f ${dir2}/sm.$N

cat ${dir2}/sm.$N.gr >> ${dir2}/sm.199101-202012.gr

rm -f ${dir2}/sm.$N.gr

cd ${dir2}

i=$((i + 1))
  if [[ $i -ge 1 && $i -le 9 ]] ; then
    N=000$i
  fi

  if [[ $i -ge 10 && $i -le 99 ]] ; then
    N=00$i
  fi

  if [[ $i -ge 100 && $i -le 999 ]] ; then
    N=0$i
  fi

  if [[ $i -ge 1000 ]] ; then
    N=$i
  fi
echo "N=", $N

$wgrib2 $dfa -match "SOILW:0.1-0.4 m" | $wgrib2 -i $dfa -bin sm.$N

cd $dir

gfortran x_read_sm_1deg.f
./a.out $N

rm -f ${dir2}/sm.$N

cat ${dir2}/sm.$N.gr >> ${dir2}/sm.199101-202012.gr

rm -f ${dir2}/sm.$N.gr

cd ${dir2}

i=$((i + 1))
  if [[ $i -ge 1 && $i -le 9 ]] ; then
    N=000$i
  fi

  if [[ $i -ge 10 && $i -le 99 ]] ; then
    N=00$i
  fi

  if [[ $i -ge 100 && $i -le 999 ]] ; then
    N=0$i
  fi

  if [[ $i -ge 1000 ]] ; then
    N=$i
  fi
echo "N=", $N

$wgrib2 $dfa -match "SOILW:0.4-1 m" | $wgrib2 -i $dfa -bin sm.$N

cd $dir

gfortran x_read_sm_1deg.f
./a.out $N

rm -f ${dir2}/sm.$N

cat ${dir2}/sm.$N.gr >> ${dir2}/sm.199101-202012.gr

rm -f ${dir2}/sm.$N.gr

cd ${dir2}

i=$((i + 1))
  if [[ $i -ge 1 && $i -le 9 ]] ; then
    N=000$i
  fi

  if [[ $i -ge 10 && $i -le 99 ]] ; then
    N=00$i
  fi

  if [[ $i -ge 100 && $i -le 999 ]] ; then
    N=0$i
  fi

  if [[ $i -ge 1000 ]] ; then
    N=$i
  fi
echo "N=", $N

$wgrib2 $dfa -match "SOILW:1-2 m" | $wgrib2 -i $dfa -bin sm.$N

cd $dir

gfortran x_read_sm_1deg.f
./a.out $N

rm -f ${dir2}/sm.$N

cat ${dir2}/sm.$N.gr >> ${dir2}/sm.199101-202012.gr

rm -f ${dir2}/sm.$N.gr


ymd=$(date -d "${ymd} +1 month" +%Y%m01)

done

exit