#!/bin/bash

GCMs=(
    access1.0_RCP45
    access1.0_RCP85
    access1.3_RCP85
    bcc-csm1.1_RCP85
    canesm2_RCP85
    ccsm4_RCP85
    csiro-mk3.6.0_RCP85
    fgoals-g2_RCP85
    gfdl-cm3_RCP85
    giss-e2-h_RCP85
    miroc5_RCP85
    mri-cgcm3_RCP85
    noresm1-m_RCP85
)

VARs=(
    PREC
    PSFC
    GLW
    SWDOWN
    Q2
    T2
    U10
    V10
)

#GCMs=(access1.0_RCP45)
#VARs=(LWDNB
#     SWDNB)

bdir="/home/disk/rocinante/DATA/temp/WRF/BIA-WRF-forcing/"
scpt="/home/disk/rocinante/DATA/temp/WRF/BIA-WRF-forcing/scripts/regrid.py"
ddir="/home/disk/rocinante/DATA/temp/WRF/var/"



for g in ${GCMs[@]}; do
    for v in ${VARs[@]}; do
	#echo $g - $v
	ifmt="${ddir}/${g}/2d/${v}_year/"	
	
	# Interpolate yearly
	#find ${ifmt}/*.nc | parallel -j10 --eta python $scpt {}
	#find ${ifmt}/*.nc | head -n1 | parallel -j10 --eta python $scpt {}
	#echo "${ifmt}/*.nc"

	
    done
    echo $g
    ifmt=${ddir}/$g/

    # Merge yearly
    #idir=${bdir}/nc/${g}/
    #python merge.py $idir
    
    
    # Create DHSVM forcings with bias-correction
    ifile=${bdir}/merge/${g}.nc
    python nc2ascii.py ${ifile}
    
done





#find "${ifmt}/${g}_${v}_d02_*.nc" | parallel echo {}