#!/bin/bash

#SBATCH --job-name=do_tasks
#SBATCH --nodes=1
#SBATCH --ntasks-per-node=20
#SBATCH --time=240:00:00
#SBATCH --output=task_out
#SBATCH --error=task_err

. configure_regcm; pwd=`pwd`

idir=$DATA_DIR/CFS/00

for d in `ls -d $idir/????????`; do d=`basename $d`

    f1=`ls $idir/$d/f*2|head -1`
    f2=`ls $idir/$d/f*2|tail -1`

    S_DATE=`echo ${f1##*/}|cut -d. -f1|sed s/flxf//g`
    E_DATE=`echo ${f2##*/}|cut -d. -f1|sed s/flxf//g`

#
# just convert to nc
#
    if [[ $S_DATE == 2025121100 && ! -d Dataset/CFS/$S_DATE && ! -d wdir/$S_DATE && ! -d wdir2/$S_DATE && ! -d done/$S_DATE ]]; then
        echo ''
        echo ' !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! '
        echo ' ! converting *.grb2 to netcdf         ! '
        echo ' ! S_DATE='$S_DATE'                    ! '
        echo ' !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! '
        echo ''
        #new_grid="-set_grib_type jpeg -new_grid_winds earth -new_grid `grid_defn.pl $pwd/Dataset/flxf.grb2`"
        CFSDIR=$idir
        emem=01               # The ensemble member (one in 01 02 03 04)
        do_pack=1             # Keep this to lower disk usage
        cmdline="-nc_table $pwd/Dataset/prc.table -append -nc4 -netcdf"

        today=$S_DATE
        START_TIME=${today}
        PLVDEST=$pwd/Dataset/CFS/${START_TIME}/${emem}/PLEV; mkdir -p $PLVDEST
        SSTDEST=$pwd/Dataset/CFS/${START_TIME}/${emem}/SST ; mkdir -p $SSTDEST

        cd $CFSDIR/${today:0:8}
        for flxf in `ls flxf*g*2`; do
            wgrib2 $flxf -match 'TMP:surface' $cmdline ${SSTDEST}/sst.${START_TIME}.nc
        done
        for pgbf in `ls pgbf*g*2`; do
            wgrib2 $pgbf -match ':TMP'  -match 'mb:' $cmdline ${PLVDEST}/air.${START_TIME}.nc
            wgrib2 $pgbf -match ':RH'   -match 'mb:' $cmdline ${PLVDEST}/rhum.${START_TIME}.nc
            wgrib2 $pgbf -match ':UGRD' -match 'mb:' $cmdline ${PLVDEST}/uwnd.${START_TIME}.nc
            wgrib2 $pgbf -match ':VGRD' -match 'mb:' $cmdline ${PLVDEST}/vwnd.${START_TIME}.nc
            wgrib2 $pgbf -match ':HGT'  -match 'mb:' $cmdline ${PLVDEST}/hgt.${START_TIME}.nc
        done
        cd $pwd
    fi
    
done

exit 0