add sum of glacier elevation classes

postproc/e3sm1p0/lnd.sh

@@ -46,6 +46,9 @@ ncwa -O -a lat,lon -w gris_area -B "gris_mask == 1" $outpath/lnd/climos/${case}_
 
 elif [ "$version" == "e3sm" ]; then
 
+# sum% glacier coverage over all elevations to create total % coverage
+ncwa -N -a glc_nec -v PCT_GLC_MEC $outpath/lnd/climos/${case}_ANN_climo.nc $outpath/lnd/climos/${case}_ANN_climo_pct_glc.nc
+
 climos_path=$outpath/lnd/climos
 
 echo "average land monthly climos over ice sheet region"