#!/bin/sh 

FIN=surfdata_1.9x2.5_simyr1850_c110114.nc
FOUT=`basename $FIN .nc`_CMIP6.nc
FCMIP6=surfdata_1.9x2.5_16pfts_Irrig_CMIP6_simyr1850_c170824.nc

# create output file 
cd /cluster/shared/noresm/inputdata/lnd/clm2/surfdata
cp -f $FIN $FOUT 
chmod +xwr $FOUT

# add old PCT_PFT 
ncks -A -v PCT_NAT_PFT,PCT_CFT,PCT_NATVEG,PCT_CROP ../surfdata_map/$FCMIP6 $FOUT

