Homec4science
Diffusion solar_potential (master) Browse

solarPOT/Potential_hourly/tilt_aspect_aggregates/postprocess_communes_gt_tilt_aspect.sh
67c48c8e30edmaster

Last Change
Raw File

postprocess_communes_gt_tilt_aspect.sh
View Options

#!/bin/bash
echo STARTING AT $(date)
echo
## ====== Get slurm array ID
WORKDIR=/scratch/walch/workdir_solar_hourly # CHANGE
CODEDIR=~/code/solarPOT/Potential_hourly/tilt_aspect_aggregates # CHANGE
## ====== User input
PATH_TO_SCRIPT_AGGREGATE=$CODEDIR\/extract_aggregate_subsets.py
PATH_TO_SCRIPT_COMMUNE=$CODEDIR\/commune_batches_gt.py
PATH_TO_SCRIPT_COMM_SUM=$CODEDIR\/merge_communes_gt.py
N_BATCHES=2 ## 95
LAST_BATCH=`expr $N_BATCHES - 1`
START_BATCH=0
ROOF_AGGREGATE_FP=$CODEDIR\/roof_aggregates.csv # CHANGE
## ===== Execute script
source activate py3_geo
for YEAR in 2016 # 2011 2012 2013 2014 2015
do
# Save annual data
AREA_FP=/work/hyenergy/output/solar_potential/hourly_potential/$YEAR\_annual_roofs.csv # CHANGE
TMP_PATH=/scratch/walch/scratch/tmp_comm_aggr_$YEAR # CHANGE
mkdir $TMP_PATH
GT_PATH=$WORKDIR\/tmp_$YEAR
for ASPECT in S # SE E NE N NW W SW
do
for TILT in 10 # 15 20 25 30 35 40 45 50
do
python $PATH_TO_SCRIPT_AGGREGATE $ROOF_AGGREGATE_FP $GT_PATH $TMP_PATH $N_BATCHES $ASPECT $TILT
for BATCH in `seq $START_BATCH $LAST_BATCH`
do
GT_FP=$TMP_PATH\/TMP_tilted_irrad_$ASPECT\_$TILT\_$BATCH\.nc
OUTFP=$TMP_PATH\/$YEAR\_commune_hourly_$ASPECT\_$TILT\_$BATCH\.nc
python $PATH_TO_SCRIPT_COMMUNE $GT_FP $AREA_FP $OUTFP
done
OUT_ALL=$TMP_PATH\/$YEAR\_commune_hourly_$ASPECT\_$TILT\.nc
python $PATH_TO_SCRIPT_COMM_SUM $OUT_ALL $N_BATCHES
if test -f $TMP_PATH\/$YEAR\_commune_hourly_$ASPECT\_$TILT\_gt.csv; then
echo $TMP_PATH\/$YEAR\_commune_hourly_$ASPECT\_$TILT\_gt.csv exists
rm $TMP_PATH\/$YEAR\_commune_hourly_$ASPECT\_$TILT\_*.nc
rm $TMP_PATH\/TMP_tilted_irrad_$ASPECT\_$TILT\_*.nc
fi
done
done
done
conda deactivate
echo
echo FINISHED AT $(date)
c4science · Help