#!/bin/bash

gcms=(
    #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
    #access1.0_RCP45
)

#basin="chehalis"
basin="snoho"
#basin="skkm"
#basin="kndy"
#basin='puyallup'

idir="/home/disk/rocinante/DATA/temp/kcp3/forcings/"
odir="/home/disk/becassine/jswon11/forcing/"
#mkdir -p $idir/forc_nc
#mkdir -p $idir/forc_ascii

for gcm in ${gcms[@]}; do
    echo $gcm

    #echo "Split netcdf"
    #python save_forc_nc.py ${idir}/3.interp/${gcm}_interp.nc
    
    #echo "Create ascii file"
    #find ${idir}/4.forc_nc/${gcm}_interp/*.nc | parallel -j10 --eta python save_forc_ascii.py {}

    #echo "Create delta"
    #python calc_delta.py ${gcm} ${basin}
    
    #echo "Create delta-applied forcing"
    #find ${odir}/6.delta/${basin}/${gcm}/delta*
    #echo "Fix leap year" 
    #find ${odir}/7.final/${basin}/${gcm}/data* | parallel -j10 --eta python fix_leap.py {} ${basin}

    #echo "apply delta adjustment"
    #find ${odir}/6.delta/${basin}/${gcm}/delta* | parallel -j10 --eta python apply_ptmin_delta.py {} ${basin}
    #find ${odir}/6.delta/${basin}/${gcm}/delta_47.93794_-122.18076.csv | parallel -j10 --eta python adjust_ptmin_grid.py {} ${basin}
    find ${odir}/6.delta/${basin}/${gcm}/delta_47.93794_-122.18076.csv | parallel -j10 --eta python apply_prism_delta.py {} ${basin}
    

    
    ##-----------------------------------------
    # Use existing prism deltas
    #delta_dir="/home/disk/tsuga2/jswon11/workdir/2021_09_KingCounty-Phase3/scripts/lelv/delta/prism_delta_snoho_fixed_tz/"
    
    #echo "Create delta-applied forcing"
    #find ${delta_dir}/delta* | parallel -j10 --eta python apply_prism_delta.py {} ${delta_dir}/low_elev_delta.csv $gcm $basin
    
    
    ##--------- Extra sites
    #echo "Split netcdf"
    #python save_forc_nc.py ${idir}/xtra/3.interp/${gcm}_interp.nc
    #python save_forc_nc.py ${idir}/3.interp/${gcm}_interp.nc
    
    #echo "Create ascii file"
    #find ${idir}/xtra/4.forc_nc/${gcm}_interp/*.nc | parallel -j10 --eta python save_forc_ascii.py {}
    #find ${idir}/4.forc_nc/${gcm}_interp/*.nc | parallel -j10 --eta python save_forc_ascii.py {}


    ##------------------------------------------------
    # Final forcings
    #pnnl=/home/disk/rocinante/DATA/temp/WRF/scripts/regrid_c2c/longwave/output/prism_grid_pyl_tlapse_1C_px.90/
    #python wrf_delta.py ${gcm} ${basin} ${pnnl}
    #find ${odir}/6.delta/${basin}/${gcm}/delta* | parallel -j10 --eta python adjust_ptmin_grid.py {} ${basin}

    # Apply wind fix, shortwave fix, longwave fix, fix leap year, add tlapse
    #find ${odir}/7.final/${basin}/${gcm}/data* | parallel -j10 --eta python final_forc.py {}
    
    
    echo '-----------------------------------'
    date
    echo '-----------------------------------'; echo; echo
    
done