mkdir -p tmp1
Prefix=TX_TM_cf_
for iy in {2003..2022}; do
    for idate in `cat hdate_2023.txt`; do
        fname=${Prefix}${iy}-${idate}.nc
        #if [ -e ${fname} ] ; then
        #   echo ${fname} " exist"
        #else
        #   echo ${fname} " does not exists"
        #fi
        # Daily maximum temp: 
        cdo timselmax,4 -selvar,mx2t6 $fname tx.nc 
        # Daily minimum temp: 
        cdo timselmin,4 -selvar,mn2t6 $fname tm.nc 
        # Ghep Tx va Tm vao cung mot file:
        cdo merge tx.nc tm.nc txm.nc
        # Convert do K sang do C:
        cdo expr,'mx2t6=mx2t6-273.15; mn2t6=mn2t6-273.15' txm.nc txm1.nc
        # Doi Unit K->C
        ncatted -a units,mx2t6,m,c,"C" -a units,mn2t6,m,c,"C" txm1.nc  
        # Luu vao Thu muc tmp voi cung ten ban dau:
        mv txm1.nc tmp1/${fname}
        rm -f tx.nc tm.nc txm.nc
    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