#!/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