#!/bin/tcsh
#========================================================================================
# name: ncjja.tcsh
#
# category: TC-shell script
#
# description
#   This script constructs time series of JJA mean (different from sequential JJA time
#   series) for specified variables of netCDF data. The resulting files may be useful
#   when testing the statistical significance for the climatological difference of two
#   experiments.
#
# input: netCDF monthly time series (e.g., 1979.01-2016.12)
#
# output: netCDF JJA time series
#
# call: NCO commands (ncrcat) and CDO commands (settaxis, mulc, add, divc)
#
# usage
#   ./ncjja.tcsh: to execute
#
# note: works on denali@EES
#
# history
#   2017.11.22: originally written
#
# author
#   m.yoshimori (myoshimo AT ees.hokudai.ac.jp)
#========================================================================================

set diri = '/mnt/logan/data/ERAI/y1979-y2016'  # input file path
set diro = '/mnt/logan/myoshimo/tmp'         # output file path
set fin = 'ts_an_sfc_y1979-y2016_mon.nc'       # input file name
set syear = '1979'                             # start year
set eyear = '1993'                             # end year

cd $diro

foreach var ("t2m" "stl1")

# unpacking netCDF
#if (-f tmp0) /bin/rm -f tmp0
#ncpdq -v ${var} -U ${diri}/${fin} tmp0

# extract JJA
ncrcat -F -d time,6,180,12 ${diri}/${fin} ts_${var}_y${syear}-y${eyear}_jun.nc
ncrcat -F -d time,7,180,12 ${diri}/${fin} ts_${var}_y${syear}-y${eyear}_jul.nc
ncrcat -F -d time,8,180,12 ${diri}/${fin} ts_${var}_y${syear}-y${eyear}_aug.nc

# JJA mean
cdo settaxis,1850-07-01,00:00,1year ts_${var}_y${syear}-y${eyear}_jun.nc tmp1.nc
cdo settaxis,1850-07-01,00:00,1year ts_${var}_y${syear}-y${eyear}_jul.nc tmp2.nc
cdo settaxis,1850-07-01,00:00,1year ts_${var}_y${syear}-y${eyear}_aug.nc tmp3.nc
cdo -b 64 mulc,30. tmp1.nc tmp4.nc  # multiply 30 for Jun.
cdo -b 64 mulc,31. tmp2.nc tmp5.nc  # multiply 31 for Jul.
cdo -b 64 mulc,31. tmp3.nc tmp6.nc  # multiply 31 for Aug.
cdo -b 64 add tmp4.nc tmp5.nc tmp7.nc
cdo -b 64 add tmp6.nc tmp7.nc tmp8.nc
cdo -b 64 divc,92. tmp8.nc tmp9.nc  # devide by 92 for JJA average
cdo settaxis,1850-07-01,00:00,1year tmp9.nc ${diro}/ts_${var}_y${syear}-y${eyear}_jja.nc

# delete intermediate files
/bin/rm -f ts_${var}_y${syear}-y${eyear}_jun.nc ts_${var}_y${syear}-y${eyear}_jul.nc ts_${var}_y${syear}-y${eyear}_aug.nc tmp?.nc

end

exit