#!/bin/bash

# Load nco
module load nco

YEAR=${1}
echo $YEAR

# Get time from the original high res file
ncks -v time global_30s_${YEAR}.nc time${YEAR}.nc

# Make a back-up copy just in case
cp LAI.01x01.${YEAR}.nc LAI.01x01.${YEAR}.nc.sav

# Remove time coordinate from regridded file as we will replace it
ncks -C -x -v time LAI.01x01.${YEAR}.nc LAI.01x01.${YEAR}.nc.notime

# Add original file time coordinate to the regridded file
ncks -A -v time time${YEAR}.nc LAI.01x01.${YEAR}.nc.notime 

# Rename regridded  file
mv LAI.01x01.${YEAR}.nc.notime LAI.01x01.${YEAR}.nc

# change time units from days since 1/1/${YEAR} to minutes since 1/1/${YEAR}
ncap2 -s 'time=time*1440' LAI.01x01.${YEAR}.nc tmp1.nc

# fix time units
ncatted -a units,time,m,c,"minutes since ${YEAR}-01-01 00:00:00 GMT" tmp1.nc tmp2.nc
ncatted -a long_name,time,m,c,"time" tmp2.nc tmp3.nc
ncatted -a begin_date,time,a,c,"${YEAR}0101" tmp3.nc tmp4.nc
ncatted -a begin_time,time,a,c,"000000" tmp4.nc tmp5.nc

# convert LAI from ubyte to float
ncap2 -s 'lai=float(lai)' tmp5.nc tmp6.nc

# edit LAI attributes
ncatted -a units,lai,o,c,'cm2 cm-2' tmp6.nc tmp7.nc
ncatted -a long_name,lai,o,c,'MODIS leaf area index' tmp7.nc tmp8.nc
ncatted -a gamap_category,lai,o,c,'MODIS' tmp8.nc tmp9.nc

# chunk and compress
nccopy -c lon/3600,lat/1800,time/1 -d1 tmp9.nc MODIS.LAI.vBNU.generic.01x01.${YEAR}.nc

# change variable name
ncrename -v lai,MODIS MODIS.LAI.vBNU.generic.01x01.${YEAR}.nc

# remove temporary files
rm -f tmp?.nc