#!/bin/ksh

edatei=20210216
edatef=20210216

dir_in=/home/ioper/data/external
dir_out=/pesq/share/oceanwave/dist/Barbara/gfs_0p25

edate=$edatei
while [ $edate -le $edatef ]; do
echo ${edate}
mkdir -p ${dir_out}/${edate:0:4}/${edate:4:2}

cdo -f nc import_binary ${dir_in}/${edate}00/dataout/GFS_0p25/WW_${edate}00.ctl ${dir_out}/${edate:0:4}/${edate:4:2}/temp.nc
cdo select,name=u10,v10 -seltimestep,1/9 ${dir_out}/${edate:0:4}/${edate:4:2}/temp.nc ${dir_out}/${edate:0:4}/${edate:4:2}/gfs_0p25_${edate}00_wind_hourly.nc

rm -f ${dir_out}/${edate:0:4}/${edate:4:2}/temp.nc

edate=$(date -u "+%Y%m%d" -d "${edate:0:8} 00:00 24 hour")
done

#  cdo select,name=uwnd,vwnd,hs WW_2021062500.nc WW_2021062500_hs_wind.nc\n
#  cdo -f nc import_binary binary/WW_2021062500.ctl netcdf/WW_2021062500.nc


#cdo -f nc copy /home/ioper/data/external/2021021600/dataout/GFS_0p25/gfs.t00z.pgrb2.0p25.f000.2021021600.grib2 file.nc
#cdo select,name=\10u,\10v temp.nc temp2.nc
#cdo chname,\10u,u10,\10v,v10 file2.nc file3.nc#