mkdir -p tmp
Prefix=LH_SH_OLR_TP_cf_
for iy in {2001..2020}; do
    for idate in `cat hdate.txt`; do
        fname=${Prefix}${iy}-${idate}.nc
        #if [ -e ${fname} ] ; then
        #   echo ${fnae} " exist"
        #else
        #   echo ${fname} " does not exists"
        #fi
        # File ${fname} tat ca cac bien deu chua tong tich luy tu t=1 den t=47. De tinh daily values thi lay ngay sau tru ngay truoc: 
        cdo -f nc -b F32 sub -seltimestep,2/47 ${fname} -seltimestep,1/46 ${fname} tmp.nc
        # Cac bien sshf, slhf, ttr co Unit la J/m^2 cua daily (da tinh o tren). Can doi sang don vi W/m^2:
        cdo expr,'sshf=-sshf/86400; slhf=-slhf/86400; ttr=-ttr/86400; tp=tp' tmp.nc tmp1.nc
        # Doi lai Unit cho phu hop:
        ncatted -a units,sshf,m,c,"W m-2" -a units,slhf,m,c,"W m-2" -a units,ttr,m,c,"W m-2" tmp1.nc  
        # Luu vao Thu muc tmp voi cung ten ban dau:
        mv tmp1.nc ${fname}
    done
done


#cdo -f nc -b F32 sub -seltimestep,2/47 $fname -seltimestep,1/46 $fname tmp/tp_daily.nc # Lay hieu gia tri lien tiep
#cdo expr,'sshf=sshf/86400' output.nc final_output.nc   # Doi sang w/m2
#ncatted -a units,sshf,m,c,"W m-2" final_output.nc      # Doi don vi do