#!/bin/sh -evx

FCMIP='vmro3_input4MIPs_ozone_CMIP_UReading-CCMI-1-0_gn_185001-185012.nc';
FCAM6='tracer_cnst_CAM6chem_piControl_cycle_3D_monthly_c170901.nc';

# change to ozone directory 
IDIR='/cluster/shared/noresm/inputdata/atm/cam/ozone';
cd $IDIR

# create output file  
FOUT=`basename $FCMIP .nc`_CAM4.nc 
cp -f $FCMIP $FOUT 

# change name of level coordinate
ncrename -d plev,lev -d bnds,nbnd -v plev,lev -v vmro3,O3 $FOUT
ncpdq -O -a -lev $FOUT $FOUT
ncatted -a units,O3,o,c,'mol/mol' $FOUT


# add time information 
ncks -C -A -v time,time_bnds,date,datesec -o $FOUT $FCAM6