#!/bin/sh -evx

# load modules
# . /opt/modules/default/init/sh
module load netcdf nco cdo



# merge variable files for all timesteps

# if 3D variable name=sst,PS

cdo select,name=U,PS N20TRAERCN_f19_g16_TAUOCNNEW0_MASKICE/atm/hist/N20TRAERCN_f19_g16_TAUOCNNEW0_MASKICE.cam2.h0.* N20TRAERCN_f19_g16_TAUOCNNEW0_MASKICE/atm/hist/temp.nc
cdo ml2pl,92500,85000,70000,50000,30000,20000 -chname,PS,sp N20TRAERCN_f19_g16_TAUOCNNEW0_MASKICE/atm/hist/temp.nc N20TRAERCN_f19_g16_TAUOCNNEW0_MASKICE/atm/hist/N20TRAERCN_f19_g16_TAUOCNNEW0_MASKICE.u.cam2.h0.1900-2000.nc
rm N20TRAERCN_f19_g16_TAUOCNNEW0_MASKICE/atm/hist/temp.nc