#!/bin/bash

ml CDO/1.9.8-intel-2019b

#6-hourly everything

#rlds mean
#rsds mean 
#pr, mm/s -> mm/6h mean

#ps	instantaneous
#huss	instantaneous
#tas	instantaneous
#uas	instantaneous
#vas	instantaneous

# bash script to compute 6-hourly values from 1-hourly RCM files (input as separate files, one file per year per variable):

output_dir="/nird/projects/6hr/"

cd /nird/projects/

# set a single year manually
#yyyy=2000

# or loop over multiple years:
for yyyy in {1985..2014}; do
    # 6-hourly means
    for var in rlds rsds pr; do
	cdo timselmean,6 ${var}_ARC-11_NorESM2-MM_historical_r1i1p1f1_AWI-ICON-2024-01_v1_1hr_${yyyy}*-${yyyy}*.nc ${output_dir}${var}_6hr_${yyyy}.nc
    done
    # instantaneous 6-hourly
    for var in ps huss uas vas tas; do
	cdo selhour,0,6,12,18 ${var}_ARC-11_NorESM2-MM_historical_r1i1p1f1_AWI-ICON-2024-01_v1_1hr_${yyyy}*-${yyyy}*.nc ${output_dir}${var}_6hr_${yyyy}.nc
    done
done