#!/bin/bash
source ~/envs/py3env3.bash
## extract NorCPM DCPP hindcast data for Yihan Tseng from CWA

varlist='SST,TS,Z200,Z500,U200,U850,V200,V850'
num=01
freq=h0 ## monthly
freq=h1 ## daily

basepath=/nird/datalake/NS9039K/projects/CMIP6_DCPP/cases/norcpm-cmip6_hindcast
caseprefix=norcpm-cmip6_hindcast_YYYYMMDD
nmem=10

MM=10
DD=15
for y in {1960..2021} ; do
    case=$(echo $caseprefix|sed -e"s/YYYY/$y/" -e"s/MM/$MM/" -e"s/DD/$DD/")
    for m in $(seq 1 $nmem); do
        mem=$(printf "%2.2d" $m)
        ofn=${case}_mem${mem}.cam2.${freq}.extract${num}.nc
        casepath=$basepath/$case/${case}_mem$mem
        printf "making $ofn\r"
        cdo -O -s -mergetime -apply,-selname,$varlist [ $casepath/atm/hist/${case}_mem${mem}.cam2.${freq}.*.nc ] $ofn 2>/dev/null
    done
done
echo DONE